基于FPGA的P-H法星上解算卫星相对姿态

周国清, 黄景金, 舒磊

周国清, 黄景金, 舒磊. 基于FPGA的P-H法星上解算卫星相对姿态[J]. 武汉大学学报 ( 信息科学版), 2018, 43(12): 1838-1846. DOI: 10.13203/j.whugis20180248
引用本文: 周国清, 黄景金, 舒磊. 基于FPGA的P-H法星上解算卫星相对姿态[J]. 武汉大学学报 ( 信息科学版), 2018, 43(12): 1838-1846. DOI: 10.13203/j.whugis20180248
ZHOU Guoqing, HUANG Jingjin, SHU Lei. An FPGA-Based P-H Method On-Board Solution for Satellite Relative Attitude[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12): 1838-1846. DOI: 10.13203/j.whugis20180248
Citation: ZHOU Guoqing, HUANG Jingjin, SHU Lei. An FPGA-Based P-H Method On-Board Solution for Satellite Relative Attitude[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12): 1838-1846. DOI: 10.13203/j.whugis20180248

基于FPGA的P-H法星上解算卫星相对姿态

基金项目: 

国家自然科学基金 41431179

国家重点研发计划 2016YFB0502501

广西创新驱动发展专项(科技重大专项) GuikeAA18118038

详细信息
    作者简介:

    周国清, 博士, 教授, 主要从事智能卫星观测系统、激光雷达研究。gzhou@glut.edu.cn

    通讯作者:

    黄景金, 博士生。jingjin_huang@tju.edu.cn

  • 中图分类号: P23

An FPGA-Based P-H Method On-Board Solution for Satellite Relative Attitude

Funds: 

The National Natural Science Foundation of China 41431179

the National Key Research and Development Program of China 2016YFB0502501

Guangxi Innovative Development Grand GuikeAA18118038

More Information
    Author Bio:

    ZHOU Guoqing, PhD, professor, specializes in the theories and methods of intelligent satellite observation systems and LiDAR technology. E-mail:gzhou@glut.edu.cn

    Corresponding author:

    HUANG Jingjin, PhD candidate. E-mail: jingjin_huang@tju.edu.cn

  • 摘要: 针对目前星上遥感图像实时处理只能实现低级别算法的情况,提出了基于现场可编程门阵列(field-programmable gate array,FPGA)的P-H法星上相对姿态实时解算模型。该模型不仅避免了传统基于欧拉角的复杂三角函数计算与初值估算,还降低了迭代次数。试验选用FPGA(V7 xc7vx1140t)作为实时解算的硬件平台。在FPGA实现中,采用64位的浮点数据结构和串行/并行相结合策略;并采用LU(Lower-Upper)分解-分块算法实现矩阵求逆。试验结果表明,该模型的迭代次数比基于欧拉角的少了13次。该模型在FPGA和计算机的实现结果相差仅为5.0×10-14,加速度比为10。另外,该模型可广泛适用于实时性要求高的图像处理领域。
    Abstract: Aimed at the situation that the low-level algorithms were implemented in satellite real time processing system for remote sensing image, this paper proposes an FPGA (field programmable gate array) -based P-H method for satellite relative attitude on-board solution. The proposed algorithm not only avoids computations of trigonometric function and estimation of initial value, but also reduces the number of iterations when comparing with the Eulerian angle-based algorithm. The Xilinx FPGA (V7 xc7vx1140tflg1930-1) is selected as the hardware platform for the real-time processing. In FPGA implementation:①We adopt a 64-bit floating point data structure and a strategy of combination of serial and parallel processing; ②a lower-upper (LU) decomposition-block algorithm is adopted for matrix inversion. The experimental results indicate that the number of iterations of the proposed algorithm is 13 less than the Eulerian angle-based algorithm. The difference of FPGA and PC implementation is about 5.0×10-14 and the speedup is about 10, which meets the requirements of precision and speed for on-board image real time processing. The proposed algorithm can be suitable for high real-time image processing field.
  • 舰船目标作为河流与海洋的主要运输载体,对其进行定位和监测在军事侦察、资源生产管理、环境保护等领域都具有重要意义[1-5]。舰船目标检测作为海洋监测的一项重要技术,是进行舰船目标跟踪、识别的前提,在军事及民用领域都有着广袤的应用前景。随着航天技术的飞速发展,利用遥感进行舰船目标检测具有大范围任意区域观测、高时效性、数据获取高性价比等优点,是进行舰船目标感知的绝佳手段。

    舰船目标实时、准确的态势感知需要结合海、陆、空大量传感数据进行探测与分析。舰船自动识别系统(automatic identification system,AIS)是最常见的舰船跟踪系统[5],它能够将舰船的位置与属性信息长时间持续地记录传输,是监测合作区域运动舰船目标的重要工具。此外,基于图像信息进行舰船目标解译,例如合成孔径雷达(synthetic aperture radar,SAR)[6]、热红外传感器[7]、光学传感器[8-9]等数据,是进行大范围舰船目标探测的重要手段。不同传感器类型都有其自身的优点和缺点,能够在不同时间、天气条件、监测场景、覆盖范围、更新频率、检测粒度等方面发挥各自的优势。SAR图像具有不受天气、云雾干扰,图像中船身金属反射特性区别于周围水面的特性,已成为海上舰船目标探测的主要手段。热红外图像通过成像对象本身的发射值来区分目标,在夜间的舰船目标监测中发挥着一定的作用。光学床图像中基于序列图像、视频图像以及基于成像延迟的单幅凝视卫星图像进行运动舰船目标检测已在动态舰船目标感知中得到了成功应用;基于高分辨率可见光的图像数据对近岸复杂区域的舰船目标进行定位识别也在舰船目标分析识别中发挥着重要作用。随着大量光学传感器的广泛使用,高空间分辨率(以下统称高分)的可见光图像数据量以几何级数增加。由于其具有获取经济便捷、细节清晰、成像与自然图像相似等优势,既能开展大规模数据下的验证评估,又能够借鉴成熟的计算机视觉目标检测方法论,基于高分可见光影像进行舰船目标检测已经成为当前主要的热点领域,相关研究呈现爆发式的增长,见图 1

    图  1  2011—2020年高分可见光遥感图像舰船目标检测相关主题文献发表统计
    Figure  1.  Statistics of the Number of Publications on the Subject of Ship Detection Based on High-Resolution Optical Remote Sensing Images from 2011 to 2020

    本文聚焦高分图像舰船目标检测技术,对近十年国内外有影响力文献(中文文献包括中国知网上核心期刊、硕博士论文,共718篇;英文文献包括科学引文索引(science citation index,SCI)论文、有影响力的会议论文,共234篇)的研究主题进行了统计分析。从研究主题来看,随着获取到图像的数量与质量提升,舰船目标从特殊场景目标变为现在的作为典型遥感目标检测的一部分开展研究;从早期的聚焦于海面影像的舰船定位与提取(包括云、海面、海浪等主题)逐渐转向复杂环境下的舰船检测(包括复杂场景、多尺度、类内变化多样等主题)研究;检测方法从人工设计、层层筛选定位检测(如疑似区域、部件、虚警等主题)发展为当前主流的基于卷积神经网络、一体化定位检测(包括深层特征、包围框、锚点等主题)。近年来,已有学者开展了有关可见光遥感影像舰船目标检测综述分析,普遍认为高分可见光遥感图像可以为精确的舰船感知提供高价值的特征及定位识别信息,但同时表明当前仅利用高分可见光图像进行舰船监测仍以检测为主,进行大范围、长时间、高精度、细粒度的舰船属性分析非常困难,需要结合多传感器数据进一步探究。现有方法模型聚焦于特定数据、特定场景下的目标检测研究,在实际不同传感器、不同复杂环境中开展鲁棒的舰船目标检测依然面临严重挑战。

    为了比较全面地回顾近十多年高分可见光图像舰船目标检测领域的研究进展,本文从舰船检测的方法体系与发展历程概述、代表性方法归纳、检测精度评估等方面对高分可见光图像舰船检测的进展进行回顾和总结,并对高分可见光舰船检测未来研究趋势进行了展望。

    由于高分可见光遥感图像空间分辨率不同、覆盖的场景变化多样,许多学者为此开发了多种多样的舰船目标检测方法。本文从检测粒度、建模方法两方面对这些方法进行总结,如图 2所示。

    图  2  高分可见光遥感图像舰船目标检测方法体系
    Figure  2.  Methodology System of Ship Detection in High-Resolution Optical Remote Sensing Images

    随着图像空间分辨率的提升,舰船目标所呈现的结构纹理信息更加丰富,目标解析层次从区域级识别(图像块分类识别是否是舰船目标)、像素级提取(按一定规则获取目标像素)到对象级目标定位(获取实例级对象的位置),在舰船目标自动化解译中得到有效应用。从检测模型尤其是舰船区域/目标特征建模来分,高分可见光图像舰船目标检测可以分为基于人工特征建模的方法与基于深度特征学习的方法两类。前者利用专家知识规则对舰船目标进行特征建模表达(大量的低层、中层特征),实现舰船目标的有效判别。后者依托于人工智能技术,利用深度神经网络自动学习目标不同层次的特征进行舰船目标的检测定位,目前在遥感图像目标检测中得到越来越广泛的应用。由于舰船目标分布的特殊性和广泛性,从检测场景来分,舰船目标检测可以分为面向特殊场景的舰船目标检测与一般场景的舰船目标检测两类。

    高分可见光图像舰船目标检测的任务是从图像中定位识别出舰船目标实例,需要解决目标“是什么”“在哪里”两大问题,遵循目标疑似区域定位、目标精细定位检测由粗到细的检测流程。在数据预处理完成后,首先进行大范围的图像舰船目标疑似区域定位,快速获取有限候选区域;然后通过对舰船目标特征精细建模在候选区域内进行进一步分析确认,完成对舰船目标的定位和判别;最后经过后处理步骤进一步剔除冗余目标,并对检测结果进行精度评价。

    遥感图像目标检测方法的演化史就是对地观测技术、信息技术的发展史,遥感图像空间分辨率的提升、计算机视觉的发展、并行计算的突破,都对高分可见光遥感图像目标检测技术的发展产生了深入影响。本文对2010—2020年高分可见光图像舰船检测方法从年代、图像分辨率、检测场景、处理流程等维度进行了分析归纳,如图 3所示,分为以下3个典型的发展阶段:(1)基于中高分可见光图像舰船目标判别提取时期。2010年左右5 m分辨率的图像是用于可见光图像舰船检测的主要对象,基于底层特征的对象分析是舰船目标检测的主要思路,利用舰船目标区别于水面的独特亮度、纹理与形状特性进行统计、学习与判别、分类,从而实现中高分辨率、较为简单场景的舰船目标检测。(2)基于高分可见光图像舰船目标特征建模检测时期。2013年左右,2 m分辨率的可见光图像获取难度降低,图像中大量中小型的舰船目标信息变得丰富、容易识别,视觉注意机制、中层特征编码等机理研究成为视觉领域的热点主题,前者在大幅面遥感影像的舰船目标定位,尤其是疑似区域定位方面得到了成功应用;后者通过对底层特征的有效高层建模弥补了语义鸿沟,广泛地使用在复杂背景、类内变化多样的舰船目标建模中,从而实现高分图像中较为复杂区域的舰船检测。(3)基于深度学习的高分可见光图像舰船目标检测时期。随着深度学习在计算机视觉领域取得了革命性的成功[10],以及2016年左右大规模高分可见光图像舰船目标检测数据的涌现,舰船目标检测的整体研究思路迅速向基于深度学习的方法演变。舰船检测流程也由分步进行向当前的一体化定位检测的方向发展,卷积神经网络(convolutional neural networks,CNN)被广泛应用,其中基于分类的检测框架[11]与基于回归的检测框架[12-14]引领着目标检测体系发展潮流。

    图  3  高分辨率可见光遥感图像舰船目标检测方法发展历程
    Figure  3.  Development History of Ship Detection Methods in High-Resolution Optical Remote Sensing Images

    目前,基于高分可见光遥感图像的舰船检测在速度、精度、场景适应性上都得到了巨大的提升,但距离实际复杂场景图像中高效的舰船目标检测仍然有一定差距。高分可见光图像舰船目标检测存在以下难点:(1)舰船目标区域定位困难。在大范围遥感图像上,舰船目标通常只占很小一部分,冗余、复杂的背景环境(如大量无关场景、破碎岛屿、凸堤、陆地目标等)干扰给舰船目标区域有效定位带来了巨大挑战。(2)舰船目标表征难。舰船目标周围复杂的背景(云、浪、岛屿、海杂波、尾杂波、黑白极性、沿岸泥沙、光照、阴影等干扰)导致舰船目标被遮挡、弱化,舰船目标本身分布的多样性(多尺度、多类别、多形态、密集排布等干扰)导致目标特征对齐困难,这些都使得舰船目标的特征类间区分性、类内统一性建模非常困难。(3)舰船检测方法普适性不高。现有的高分可见光图像舰船检测标注数据集中,目标在尺度、类型、场景、视角分布上相对有限,舰船目标类型呈长尾特性分布,当前检测模型对难分目标,以及分布外的不同传感器、不同拍摄环境图像目标检测适应性差,难以满足大范围的自动化检测应用需求。

    此类舰船目标检测方法逐步去除干扰因素,对于大幅面的高分可见光图像,一般需要进行舰船疑似候选区域提取从而获取目标区域,并在此基础上进行舰船目标检测定位得到精确的目标属性与位置信息。也有一些研究基于图像切片直接进行舰船目标检测定位。

    基于先验信息的海陆分离针对大范围水陆图像的舰船检测在有影像位置信息和相应的地理先验信息时,利用精确的地理信息辅助进行水陆分离是获取舰船目标候选区域的有效方法。其中先验信息包括港口地理信息、海岸线信息等[15],并通过精确的模板匹配或图像配准后分离出水域图像。此类方法的优点是分割效果准确稳定,不受场景、光照等因素影响。但同时缺点也很明显,当图像位置信息缺失或不准确、地理先验信息缺乏时,配准结果不再稳定,此过程就无法保证顺利进行。

    基于分割的疑似区域提取对于灰度分布较为一致的水面区域,利用分割算法能够有效地实现大面积云雾去除、水陆区域分离、水面与舰船疑似区域分离,因此被广泛应用于舰船疑似区域提取。具体方法有基于全局阈值分割方法[16]、基于能量函数分割方法[17]、基于局部多阈值分割方法[18-21]等。此类方法能够有效地从简单场景图像中提取疑似舰船区域,但对于海杂波、大量碎云、海浪、小岛干扰较多的海面区域以及复杂的近岸区域,往往不能有效地进行分割,导致结果虚警和漏检都很高。

    基于特征分析的疑似区域提取对于较为复杂的海面和近岸场景需要对目标进行有效特征建模来区分背景与疑似舰船区域。提取的特征大多具有一定的局部或全局不变性,如灰度[18]、颜色直方图[19]、纹理[20-22]、频域[23]、几何形状特征[24]等。此外基于生物视觉启发的视觉显著特征是一种特别的对比差异性特征,能够快速将少量舰船疑似区域从大范围的海面背景中区分出来。具体包括生物视觉注意模型[24-2539]、频域显著模型[26]以及多种显著性融合特征[27]。特征判别方式包括统计判别方法[222841-42]、决策方法[29]、异常分析方法[30]、特征学习分类方法[26-2731-32]等。基于特征统计判别(尤其是基于显著特征判别)、异常分析的方法大多以像素为单位分析,易受不同环境因素变化的影响,同时对目标类内变化的适应能力较弱,适合分辨率较低、船体稀疏分布的目标候选区域提取。利用特征投票、决策或学习分类的方法能够一定程度上以对象为单位进行特征建模,但其提取结果受特征设计及数据影响大,很难产生较为通用的大范围复杂场景舰船候选区域提取方法。

    在经过大量背景、虚警去除获取到有限的疑似舰船候选区域后,需要进一步对其进行精细分析,输出更为精细的舰船目标定位结果。此阶段面临的主要问题来自于与舰船区域相似的背景干扰(类间差异小),以及舰船内部多尺度、多类型分布偏差(类内差异大)。主要的方法包括基于统计判别的方法、基于投票判别的方法、基于特征学习的方法等。

    基于统计判别的方法通过对舰船与非舰船目标的特征进行统计,基于聚类或相似性度量的方法进行舰船目标判别[20-2233-34]。常用的特征有灰度、纹理、几何结构(长宽比、圆度、角点、边缘、形状轮廓等)等特征。此类方法基于对单个舰船的目标区域块进行分析,因此仅适用于海面场景、稀疏分布的舰船目标,对于复杂的舰船场景无法进行有效的统计判别。

    基于投票判别的方法对于高分辨率的复杂近岸场景图像,基于部件表示的模型将舰船目标检测转化为利用一组对象组件、形状等特征进行相似性测量。在舰船检测中大多利用广义霍夫变换 [35-37]进行特征度量,再通过后续的投票获取最终的定位检测结果。此类方法对于平静水面场景的靠岸舰船能取得良好的检测提取结果,但并不能适用于各种复杂背景干扰如阴影、光照、海浪、云雾等,以及舰船目标多尺度、多姿态、多方向分布复杂的情况。

    基于特征学习的方法为了适用于大范围、大规模的舰船目标定位检测,基于大量数据学习分类的方法被广泛使用。一些具有不变性质的特征[222538-39]表示,如尺度不变特征(scale invariant feature transform,SIFT)、Gabor特征、梯度直方图(histogram of oriented gradient,HOG)特征、视觉显著特征等被广泛用于特征提取。为了弥补底层视觉特征与高层语义之间的鸿沟,基于中层特征编码建模的目标特征表征模型被广泛地使用,包括词袋模型(bag of words,BOW)[2340]、稀疏编码[3241-44]等。这类基于特征编码的方法能够较好地适用于复杂背景下的舰船目标判别,但大多是基于图像块的建模,因此不能进行目标的有效定位。后续基于部件模型(deformable part models,DPM)的方法[45-47]能够有效对目标进行各个部件的组合建模,从而实现目标的精确定位与检测。但这种方法需要将目标进行统一尺度的对齐,人工参与较多,对于旋转、尺度、类型差异大的舰船目标检测效果差。

    尽管基于深度CNN模型与传统的舰船目标检测模型在特征设计上有巨大差异,但整个检测流程后者对前者的影响深远。虽然基于深度学习的目标检测流程逐渐形成了一体化的目标检测体系,但传统方法尤其是后期特征编码学习的方法中的滑动窗口策略、多层次信息整合等思想依然在基于深度学习的舰船目标检测中起着重要作用。当前基于深度学习的高分可见光图像舰船目标方法依托于计算机视觉中目标检测领域的方法体系,在此基础上针对可见光遥感图像本身的特殊性、舰船目标的独特性进行进一步深化改进,推动着高分可见光遥感图像的舰船检测向着更加自动化、高精度、实用性强的方向发展。增强目标特征表征的有效性、多尺度定位的适应性、多姿态位置回归的稳定性、不同场景范围的鲁棒性等都是有效提升舰船目标检测性能的关键,现有的检测方法往往针对其中一点或多个点综合开展研究,本文将这4个方面的具体方法分别进行总结分析。

    高分可见光遥感图像可以清晰地呈现一定范围内的目标及其背景,在带给目标更丰富细节特征的同时引入了大量的冗余信息。为了能将目标从大量复杂的负样本中分离出来,同时尽可能满足不同类型、尺度、姿态的目标表征,就需要提取更为有效的目标特征来提升目标判别的准确性。

    除了加深网络层数外,增大卷积核的感受野(receptive field,RF)能够对更大范围的图像块(包括目标及其上下文)进行特征描述,从而有效提升CNN网络的特征提取能力。早期的一些方法直接利用大卷积核(如Alexnet[48])进行特征抽取,增强网络表征能力的同时给计算资源带来了极大挑战,因此后续发展的多个小卷积核代替大卷积核(Inception [49]系列)、空洞卷积(dilated CNN[50])、金字塔池化(spatial pyramid pooling,SPP[51])等方法被广泛用于视觉任务中。在舰船目标检测方法中,Zheng等[52]开发了多感受野卷积分组表征聚合模块,用于捕获超尺度特征的尺度不变性,提升模型对于舰船目标尺度变化的适应性。此外,各种上下文关联模块[53-54]被用来增强网络的舰船检测性能尤其是小目标的检测性能。

    另一种能够辅助网络有效忽略无关信息,关注重点信息,从而加强网络对于目标及其相关上下文的模型为注意力机制。不同于单纯的特征设计,注意力机制通过显式的监督,将学习到的权重信息施加在特征上进一步强化对于目标表征更有用的特征。典型的视觉任务有关的注意力机制有通道注意力[55]、空间注意力[56]、混合域注意力[57]、全局注意力[58]等。在遥感目标舰船目标检测中,一些学者基于不同层次特征[59-60]、不同空间尺度的注意力[61-63]机制来提升目标的表征能力,取得了不错的结果。考虑到空间上的和全局的注意力需要的计算量很大,舰船检测中目前更多的是使用分割的结果来定位关键点或区域[64-66],以此作为进一步约束检测框架产生更高精度的目标定位检测结果。

    此类方法能充分应用目标及其上下文信息进行丰富的特征建模,但同时也不可避免引入新的问题,如新引入的超参数的敏感性与鲁棒性,新增模块与检测定位网络如何联合优化等。

    为适应不同尺度的目标定位与识别,便于并行运算,充分利用CNN网络特征层的多尺度特性代替多尺度图像多次进行不同尺度的统一区域定位与检测是一种更为快速有效的方法。特征金字塔网络(feature pyramid network,FPN)[67]是一种典型的联合自上而下的多尺度特征增强网络,通过使用更深的层探测较大的物体,使用较浅层在较大的尺度上检测小目标被广泛用于目标检测研究中。Li等[68]构造的HSF-Net将不同尺度的特征映射到同一尺度空间进行采样回归检测,高层特征能够有效检测大尺度目标的同时,底层尺度空间特征经过与上层特征的融合提升了对于小目标的感知能力。文献[69-71]基于FPN的结构或改进不同尺度特征层的连接方式如稠密连接、跳跃连接等进行多尺度特征融合的舰船检测。

    与深度特征增强的目标特征表征方法尤其是上下文信息融合的方法不同,FPN系列的多尺度特征融合模型更强调在解码阶段将不同层的特征经过采样后与原先层进行融合,能够有效地将图像中不同尺度的目标进行统一学习,提升了舰船目标检测的速度;同时通过合理分配相应尺度的目标回归检测,有效地提升了目标检测的性能。如何构建强健的多尺度特征融合结构是此类方法的重点与难点,后续很多研究在此基础上探索更好的连接方式设计[72]或进行网络架构搜索(NAS-FPN[73])学习,以追求更高效、更高精度的目标检测性能。

    遥感目标不同于一般自然图像与人脸目标检测数据集中的竖直排布,其目标位置分布是任意方向的。舰船目标的长宽比范围比一般的目标更为多样,且目标与目标之间、目标与背景之间干扰都非常严重,影响着舰船目标的准确定位与判别。使用传统的一般水平框回归会导致目标正样本区域存在大量背景干扰使得回归位置不准确,而在后处理过程中也可能因为水平框重叠度过高导致部分检测框被删除而导致漏检。基于此,遥感领域的学者开展了大量的基于更准确目标位置回归的研究[6274-83],引领了独属于遥感领域目标检测的研究热潮。其中最直接的解决方法就是通过将水平的锚点(anchor)框进行旋转增强[74],以此更好地贴合旋转框定位检测任务。此类方法能够改善密集排布的舰船目标检测情况,但由于不可能穷举所有的可能角度,导致不能产生足够的与正样本比较贴合的正例,而且对于参数相当敏感,不容易学习,检测速度缓慢。后续有学者陆续提出了旋转池化区域推荐层[75-76],用于进行旋转特征学习定位。Zhang等 [77]直接利用旋转anchor进行两阶段的回归检测;Yang等[78]研究R3Det、SCRDet[62]等模型在生成初始更容易回归的区域后进一步细化对齐到精确的目标位置;Ding等 [79]通过设计一个可学习的RROI Learner模块,进行更高精度的锚点框生成。针对旋转框角度不连续的问题,Zhu等 [80]基于一种新的表示方法进行目标角度的可连续学习。

    除了水平框不能很好贴合目标的问题外,基于锚点设计(anchor-based)的策略对于舰船这种具有极大长宽比变化的目标也不能很好地贴合(用传统的1∶1、1∶2、2∶1预设锚点产生的正样本很少)。一些学者基于数据集中的舰船尺寸统计[81],并以此选择合适的尺度分配进行定位检测;也有学者基于任意锚点(anchor-free)定位的思路,基于可学习的锚点预测模型[82-83]来灵活学习目标的位置。对于常用的位置回归损失与真实目标预测回归框不完全对应导致回归偏差的问题,提出了基于改进目标检测分类与回归的损失函数[7982-83]、多次逐步精细回归[8084-85]以及改进的非最大值抑制(non-maximum suppression,NMS)[5986]等方法进行更为有效的特征学习、模型训练以及结果优化,从而达到更高的检测精度。

    此类方法通过深入分析影响舰船目标定位精度的因素,利用增强目标区域对齐的方法进行检测性能优化,有效地提升了目标定位与检测的上限。如何保证回归过程的稳定性以及提升定位的泛化能力是此类方法面临的主要问题。

    尽管多尺度特征金字塔的方法能够一定程度上解决目标的尺度问题,但由于遥感图像舰船目标尺度差异极大,当前将图像裁切缩放后统一尺寸地输入检测网路不可避免地导致部分过大或极小的目标无法得到有效的学习。为了缓解这种目标尺度分布规模不平衡问题,如何合理地进行多尺度候选区域获取,使得不同尺寸的目标都能够尽可能在一定的尺度范围内被检测到是此类方法的重点研究目标。Li等[87]以及You等[88]在统计小尺寸舰船目标分布的基础上进行基本图像块的确定,并基于显著性评估的方法对图像中的大尺寸船图像块缩小后进行统一目标尺度分布下的定位与检测。Yu等[89]在获取较为准确的舰船候选区后,对每个区域基于中心裁剪的方法获取多尺度图像,并以此输入卷积网络中进一步进行目标判别。Wang等[90]采用不同的块尺度对大规模舰船数据集中的目标进行检测,计算不同块尺度下不同目标检测结果的平均精度,并统计最优值获取对应的图像尺度进一步用于自适应的舰船检测以提升整体检测精度。此外,无论是模型训练还是测试阶段,随机裁剪、多尺度图像训练/测试都能极为有效地提升检测方法对不同尺度分布的目标适应性[91]

    针对实际应用特别是可见光卫星图像上的舰船检测范围相当大,真实舰船分布的区域往往只占很小一部分,直接进行目标检测由于大量的冗余背景使得正负样本分布极度不平衡,严重影响检测的准确性和时效性。因此如何自适应提取大范围图像上的舰船目标疑似区域尤为关键。此类型方法目前的研究思路都是利用舰船目标本身的场景关联特性,利用先验地理信息[88]进行港口场景区域的定位实现疑似图像区域的提取,在此基础上进一步执行舰船目标检测。

    此类方法可以显式地适应不同大小、场景的舰船目标以提升检测精度,但需要有专业知识辅助学习,而对于不同场景图像之间尺度参数的适应性及迁移性研究需要进一步探索。

    对于舰船目标的检测结果评估是高分辨率可见光遥感图像目标检测的重要一环。目前对于遥感图像舰船检测性能的定量评估主要基于对象级统计,检测的结果可能会出现正确检测、错误检测和漏检3种情况,主要用以下4个指标来评价检测的结果:查全率(Recall)、查准率(Precision)、平均精度(average precision,AP)、平均均值精度(mean average precision,MAP)。查全率R指预测值中被正确识别到的正样本数占所有正样本数量的比例;查准率P表示被正确识别的正样本在所有预测到为正样本中的占比;AP表示在不同的测试召回率下的平均精度,对应着不同的IOU(intersection over union)阈值;MAP首先对每一个类别都计算出AP,然后再计算平均值,是对检测到的目标平均精度的一个综合度量。计算如下:

    R=TPTP+FN×100%
    P=TPTP+FP×100%
    AP=01PRdR
    MAP=mAPm

    式中,TP(true positives)、FP(false positives)、FN(false negatives)分别表示检测到的目标与真值比对后产生的真阳性、假阳性以及假阴性目标个数;m表示所有样本的类别个数。预测框(Bp)位置是否预测正确(实际标注真值框(Bt))通过IOU判定:

    val=true,IOU>Tfalse,
    IOU=BpBtBpBt

    式中,T一般通过具体指定来约束检测框与真值的贴合程度要求。

    对于检测结果的精度评估离不开数据集的构建。随着第一个大规模遥感公开目标场景数据集UCmerged[92]的提出,后续的大规模、细粒度标注数据集不断涌现。本文对位置标注级的可见光图像舰船检测数据集进行介绍,具体包括NWPU VHR-10[93]、HRSC2016[94]、Kaggle-ship[95]、HRRSD[96]、DOTA[97]、iSAID[98]、DIOR[99]、xView[100]等(见表 1)。

    表  1  高分辨率可见光遥感图像舰船目标检测有关公开数据集统计
    Table  1.  Public Detection Datasets Related to Ship Targets in High-Resolution Optical Remote Sensing
    数据集 目标个数(训练+验证) 图像大小/像素 分辨率/m 标注方法 被引用次数
    NWPU VHR-10[93] 约300 450~1 500 0.5~2 水平框 261
    HRSC2016 [94] 约3 000 300~1 500 0.4~2 水平框/旋转框 39
    Kaggle-ship[95] 约80 000 768 约2 掩码 43
    HRRSD[96] 约2 000 150~15 000 0.15~1.2 水平框 20
    DOTA [97] 约45 000 400~15 000 0.1~5 水平框/旋转框 315
    iSAID[98] 约80 000 400~15 000 0.1~5 掩码 11
    DIOR[99] 约60 000 800 0.5~5 水平框 43
    xView[100] 约10 000 2 500~5 000 0.3 水平框 10
    下载: 导出CSV 
    | 显示表格

    表 1中可以看出,近几年高分可见光遥感图像舰船检测相关数据集发展迅速且备受瞩目,亚米级的图像被广泛应用于检测任务中,并且有由少量类型专用数据集向综合多类数据集发展的趋势,目标标注粒度也从水平框的标注向更为精确的旋转框、实例掩码标注发展。HRSC2016和Kaggle-ship是专门针对舰船目标的数据集,前者是目前公开的唯一可以进一步识别舰船类型的高分辨率图像数据集,而后者是舰船目标掩码标注的数据集。

    随着对地观测技术和信息技术的蓬勃发展,研究人员在高分可见光遥感图像舰船检测技术方面做了大量的研究工作,2010—2020年基于深度学习的高分可见光图像舰船检测无论在检测规模、检测性能还是检测自动化程度上都得到了巨大的提升。但目前的检测技术面对数据源的不确定性、算法的适应性、稳定性与可靠性等难题依然存在挑战,距离自适应、高精度、实时舰船目标检测及智能解译依然存在很大的差距,有待进一步的研究。本文对基于高分可见光遥感图像舰船目标检测的展望如下:

    1)舰船目标检测性能进一步增强。目前利用高分可见光遥感图像的舰船目标检测在大规模数据集上的检测精度有待进一步提高,实际应用于大范围、多云雾干扰的近岸及海上图像舰船检测时,存在大量虚检与漏检,需要更加有效的检测算法。对于背景与目标分布极度不平衡问题、不同尺度目标尤其是小目标(二十个像素内)、极限长宽比目标、难分目标等的定位判别问题需要进一步探索。

    2)舰船目标检测算法适应性提升。目前基于深度学习的舰船检测框架种类繁多、设计复杂,包含的超参数非常多且严重依赖于经验,对于检测性能影响很大。发展自动的网络结构设计、参数初始化与调整策略优化等技术是未来自适应舰船目标检测的重要方向。对于不同的数据集、不同的应用任务,目前方法需要经过部分的新数据标注及模型迁移后才能进行目标检测,过程成本开销大且效率低下,因此需要开发基于弱监督甚至无监督的图像标注与目标检测,提高舰船目标检测对不同数据源、任务的适应能力。开展受限资源下如移动平台、在轨卫星数据的舰船模型轻量化、实时检测技术研究,为实际舰船目标动态监测应用提供支撑。

    3)舰船目标检测识别能力深化。实际舰船检测应用中对于舰船目标的定位及进一步的类别识别有着迫切需求,当前舰船目标类型识别主要还是依赖人工判读解译舰船类别属性。高分可见光遥感图像仅仅提供了有限视角下(下视与侧下视)的目标形态,且图像中的目标往往存在模糊、扭曲、遮挡的情况,仅依赖高分可见光图像进行舰船目标识别相当困难。建立权威、完善的大规模舰船目标分类数据集是实现自动化舰船目标识别的基础。同时需要融合多源遥感信息(如AIS、凝视卫星图像、SAR图像、光学图像等)实现对舰船目标多维度特征互补关联,提升任意时间、场景下的智能化舰船目标解译能力。

  • 图  1   FPGA硬件架构

    Figure  1.   FPGA Architecture

    图  2   Rpqr子模块

    Figure  2.   R and p, q, r Sub-modules

    图  3   AL子模块

    Figure  3.   A and L Sub-modules

    图  4   ATAATL子模块

    Figure  4.   ATA and ATL Sub-modules

    图  5   LU分解-分块算法流程图

    Figure  5.   Flowchart of LU Decomposition Block Algorithm

    图  6   L11-1U11-1子模块

    Figure  6.   L11-1 and U11-1 Sub-modules

    图  7   L21U12子模块

    Figure  7.   L21 and U12 Sub-modules

    图  8   L22-1U22-1子模块

    Figure  8.   L22-1 and U22-1 Sub-modules

    图  9   M12-1N21-1子模块

    Figure  9.   M12-1 and N21-1 Sub-modules

    图  10   U11-1L11-1M12-1N21-1子模块

    Figure  10.   U11-1L11-1 and M12-1N21-1 Sub-modules

    图  11   B-1子模块

    Figure  11.   B-1 Sub-module

    图  12   改正数模块

    Figure  12.   Correction Module

    图  13   迭代过程的仿真波形

    Figure  13.   Simulation Waveform of Iteration Processing

    图  14   LU分解-分块算法的仿真波形

    Figure  14.   Simulation Waveform of LU Decomposition Block Algorithm

    表  1   两幅图像点对(f=0.1 m)

    Table  1   Point Correspondences of Two Images (f=0.1 m)

    点号 左图像 右图像
    行号 列号 行号 列号
    1 -37.104 -20.849 -43.682 -21.725
    2 -34.718 -21.669 -41.243 -22.652
    3 -32.359 -22.520 -38.839 -23.610
    4 -29.959 -23.372 -36.396 -24.567
    5 -27.631 -24.187 -34.030 -25.483
    6 -36.121 -18.675 -42.539 -19.529
    7 -33.601 -19.385 -39.958 -20.350
    8 -31.343 -20.141 -37.657 -21.206
    9 -28.945 -20.954 -35.218 -22.125
    下载: 导出CSV

    表  2   结果对比

    Table  2   Results Comparison

    参数 PC1 (欧拉角) PC2 (P-H法) FPGA (P-H法) |PC1-PC2| |PC2-FPGA|
    φ(picth)
    航向倾角
    0.043 373 3
    6 510 974
    0.043 373 3
    6 316 752
    0.043 373 3
    6 316 747
    1.9×10-9 5×10-14
    κ(yaw)
    旁向倾角
    0.044 063 6
    4 319 247
    0.044 063 6
    5 064 706
    0.044 063 6
    5 064 708
    7.5×10-9 2×10-14
    ω(roll)
    像片旋角
    0.021 523 4
    0 090 194
    0.021 547 6
    7 515 801
    0.021 547 6
    7 515 795
    2.4×10-5 5×10-14
    迭代次数 20 7 7 13 0
    运行时间/ms 2.269
    (3.40 GHz)
    2.873
    (3.40 GHz)
    0.308
    (100 MHz)
    / /
    下载: 导出CSV

    表  3   FPGA硬件资源消耗

    Table  3   Utilization of FPGA Resources

    资源 消耗情况 占总资源百分比/%
    FF 656 464 46.1
    LUT 382 344 53.7
    Memory LUT 45 596 16.1
    DSP48 2 712 80.7
    下载: 导出CSV
  • [1]

    Zhou Guoqing, Baysal O, Kaye J. Concept Design of Future Intelligent Earth Observing Satellites[J]. International Journal of Remote Sensing, 2004, 25(14):2667-2685 doi: 10.1080/0143116031000101558

    [2] 李德仁, 沈欣.论智能化对地观测系统[J].测绘科学, 2005, 30(4):9-11 doi: 10.3771/j.issn.1009-2307.2005.04.001

    Li Deren, Shen Xin. On Intelligent Earth Observation Systems[J]. Science of Surveying and Mapping, 2005, 30(4):9-11 doi: 10.3771/j.issn.1009-2307.2005.04.001

    [3] 张兵.智能遥感卫星系统[J].遥感学报, 2011, 15(3):415-431 http://d.old.wanfangdata.com.cn/Periodical/wjsjxx201101013

    Zhang Bing. Intelligent Remote Sensing Satellite System[J]. Journal of Remote Sensing, 2011, 15(3):415-431 http://d.old.wanfangdata.com.cn/Periodical/wjsjxx201101013

    [4] 李德仁, 王密, 沈欣, 等.从对地观测卫星到对地观测脑[J].武汉大学学报·信息科学版, 2017, 42(2):143-149 http://ch.whu.edu.cn/CN/abstract/abstract5653.shtml

    Li Deren, Wang Mi, Shen Xin, et al. From Earth Observation Satellite to Earth Observation Brain[J]. Geomatics and Information Science of Wuhan University, 2017, 42(2):143-149 http://ch.whu.edu.cn/CN/abstract/abstract5653.shtml

    [5] 李德仁, 龚健雅, 邵振峰.从数字地球到智慧地球[J].武汉大学学报·信息科学版, 2010, 35(2):127-132 http://ch.whu.edu.cn/CN/abstract/abstract836.shtml

    Li Deren, Gong Jianya, Shao Zhenfeng. From Digi-tal Earth to Smart Earth[J]. Geomatics and Information Science of Wuhan University, 2010, 35(2):127-132 http://ch.whu.edu.cn/CN/abstract/abstract836.shtml

    [6] 杜列波, 肖学敏, 鲁琴, 等.基于FPGA+多DSP的JPEG2000星载遥感图像压缩实现方案[J].测试技术学报, 2008, 22(6):478-482 doi: 10.3969/j.issn.1671-7449.2008.06.002

    Du Liebo, Xiao Xuemin, Lu Qin, et al. An Implementation Scheme for JPEG2000 Satellite-Borne Remote Sensing Image Compression Based on FPGA + Multi-DSPs[J]. Journal of Test and Measurement Technology, 2008, 22(6):478-482 doi: 10.3969/j.issn.1671-7449.2008.06.002

    [7] 龚健雅, 钟燕飞.光学遥感影像智能化处理研究进展[J].遥感学报, 2016, 20(5):733-747 http://d.old.wanfangdata.com.cn/Periodical/ygxb201605006

    Gong Jianya, Zhong Yanfei. Survey of Intelligent Optical Remote Sensing Image Processing[J]. Journal of Remote Sensing, 2016, 20(5):733-747 http://d.old.wanfangdata.com.cn/Periodical/ygxb201605006

    [8] 高立宁, 龙腾.基于FPGA的星上实时信息处理系统[C].全国信号处理与应用学术会议, 南京, 2008

    Gao Lining, Long Teng. On-Board Spaceborne Real-Time Digital Signal Processing System[C]. Compu-ter Engineering and Applications, Nanjing, China, 2008

    [9] 尤政, 戴汩. "航天清华一号"微小卫星及其图像处理[J].遥感学报, 2001, 5(3):177-182 http://d.old.wanfangdata.com.cn/Periodical/ygxb200103004

    You Zheng, Dai Mi. "Hangtian Tsinghua_1" Micro-Satellite and Its Image Processing[J]. Journal of Remote Sensing, 2001, 5(3):177-182 http://d.old.wanfangdata.com.cn/Periodical/ygxb200103004

    [10]

    Williams J A, Dawood A S, Visser S J. FPGA-Based Cloud Detection for Real-Time Onboard Remote Sensing[C]. IEEE International Conference on Field-Programmable Technology, Hong Kong, China, 2002

    [11] 高昆, 刘迎辉, 倪国强, 等.光学遥感图像星上实时处理技术的研究[J].航天返回与遥感, 2008, 29(1):50-54 doi: 10.3969/j.issn.1009-8518.2008.01.010

    Gao Kun, Liu Yinghui, Ni Guoqiang, et al. Study on On-Board Real-Time Image Processing Technology of Optical Remote Sensing[J]. Spacecraft Recovery & Remote Sensing, 2008, 29(1):50-54 doi: 10.3969/j.issn.1009-8518.2008.01.010

    [12]

    Huang Jingjin, Zhou Guoqing. On-Board Detection and Matching of Feature Points[J]. Remote Sen-sing, 2017, 9:601 doi: 10.3390/rs9060601

    [13]

    Huang Jingjin, Zhou Guoqing, Zhang Dianjun, et al. An FPGA-Based Implementation of Corner Detection and Matching with Outlier Rejection[J]. International Journal of Remote Sensing, 2018, doi: 10.1080101431161.2018.1500728

    [14]

    Zhou Guoqing, Jiang Linjun, Huang Jingjin, et al. FPGA-Based On-Board Geometric Calibration for Linear CCD Array Sensors[J]. Sensors, 2018, 18(6):1794 doi: 10.3390/s18061794

    [15] 闫利, 聂倩, 赵展.利用四元数描述线阵CCD影像的空间后方交会[J].武汉大学学报·信息科学版, 2010, 35(2):201-204 http://ch.whu.edu.cn/CN/abstract/abstract846.shtml

    Yan Li, Nie Qian, Zhao Zhan. Space Resection of Line Scanner CCD Image Based on the Description of Quaternions[J]. Geomatics and Information Science of Wuhan University, 2010, 35(2):201-204 http://ch.whu.edu.cn/CN/abstract/abstract846.shtml

    [16] 周拥军, 邓才华.利用HGA和单位四元数的相对定向解法[J].武汉大学学报·信息科学版, 2011, 36(6):670-673 http://ch.whu.edu.cn/CN/abstract/abstract561.shtml

    Zhou Yongjun, Deng Caihua. A New Method for Relative Orientation with Hybrid Genetic Algorithm and Unit Quaternion[J]. Geomatics and Information Science of Wuhan University, 2011, 36(6):670-673 http://ch.whu.edu.cn/CN/abstract/abstract561.shtml

    [17] 贲进, 童晓冲, 闾海庆.基于Pope-Hinsken算法的空间后方交会[J].测绘科学技术学报, 2011, 28(1):37-41 doi: 10.3969/j.issn.1673-6338.2011.01.009

    Ben Jin, Tong Xiaochong, Lv Haiqing. Space Resection Based on the Pope-Hinsken Algorithm[J]. Journal of Geomatics Science and Technology, 2011, 28(1):37-41 doi: 10.3969/j.issn.1673-6338.2011.01.009

    [18] 江刚武, 姜挺, 王勇, 等.基于单位四元数的无初值依赖空间后方交会[J].测绘学报, 2007, 36(2):169-175 doi: 10.3321/j.issn:1001-1595.2007.02.010

    Jiang Gangwu, Jiang Ting, Wang Yong, et al. Space Resection Independent of Initial Value Based on Unit Quaternions[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(2):169-175 doi: 10.3321/j.issn:1001-1595.2007.02.010

    [19] 王之卓.摄影测量原理[M].武汉:武汉大学出版社, 2007

    Wang Zhizhuo. Photogrammetry Principle[M]. Wuhan:Wuhan University Press, 2007

    [20] 邬贵明. FPGA矩阵计算并行算法与结构[D].长沙: 国防科学技术大学, 2011 http: //cdmd.cnki.com.cn/Article/CDMD-90002-1011303506.htm

    Wu Guiming. Parallel Algorithms and Architectures for Matrix Computations on FPGA[D]. Changsha: National University of Defense Technology, 2011 http: //cdmd.cnki.com.cn/Article/CDMD-90002-1011303506.htm

  • 期刊类型引用(6)

    1. 杨凯钧,雷帆,曹里,唐伟,张哲. 联合GNSS和InSAR地灾区纯净形变场解算研究. 地理空间信息. 2024(06): 113-117 . 百度学术
    2. 赵奇. 基于GNSS和InSAR技术的矿区建筑物形变监测. 测绘通报. 2024(11): 126-132+166 . 百度学术
    3. 陈明锴,许光煜,王乐洋. InSAR同震地表三维形变反演:一种顾及形变梯度的联合解算方法. 武汉大学学报(信息科学版). 2023(08): 1349-1358 . 百度学术
    4. 柯福阳,胡祥祥,明璐璐,刘学武,尹继鑫,刘宇航. 面向地表形变高精度监测的GNSS-InSAR融合方法. 遥感技术与应用. 2023(05): 1028-1041 . 百度学术
    5. 奎达敏,喜文飞,杨正荣,赵子龙. 基于D-InSAR技术的不同裁剪尺度下沉降监测精度对比分析. 测绘与空间地理信息. 2023(11): 39-41 . 百度学术
    6. 何秀凤,高壮,肖儒雅,罗海滨,贾东振,章浙涛. InSAR与北斗/GNSS综合方法监测地表形变研究现状与展望. 测绘学报. 2022(07): 1338-1355 . 百度学术

    其他类型引用(17)

图(14)  /  表(3)
计量
  • 文章访问数:  2136
  • HTML全文浏览量:  195
  • PDF下载量:  234
  • 被引次数: 23
出版历程
  • 收稿日期:  2018-08-01
  • 发布日期:  2018-12-04

目录

/

返回文章
返回