Influence of MGEX Differential Code Bias Products on BDS-3 Pseudorange Single Point Positioning
-
摘要:
为研究不同差分码偏差(differential code bias,DCB)产品对BDS-3伪距单点定位的影响,推导了BDS-3单频和双频无电离层组合伪距单点定位卫星端DCB改正模型,基于多模全球导航卫星系统(global navigation satellite system,GNSS)实验跟踪网发布的中国科学院(Chinese Academy of Sciences, CAS)和德国宇航中心(Deutsches Zentrum für Luftund Raumfahrt,DLR)两种DCB产品,并结合全球不同纬度均匀分布的国际GNSS服务测站连续7 天的BDS-3实测数据对B1I、B1C和B2a单频与B1C/B2a、B1I/B3I两种常用无电离层组合伪距单点定位精度进行对比分析。实验结果表明,在B1I、B1C、B2a 3种频率和B1C/B2a、B1I/B3I两种组合中,CAS产品改正后的定位精度整体上优于DLR产品,具体来说,定位精度的平均值在E、N、U 3个方向上分别提高了3 cm、1 cm和6 cm。
Abstract:ObjectivesIn order to study the influence of different differential code bias (DCB) products on BDS-3 pseudorange single point positioning, the BDS-3 single-frequency and dual-frequency ionospheric-free combined pseudorange single point positioning satellite-side DCB correction model was derived.
MethodsBased on two DCB products of Chinese Academy of Sciences (CAS) and Deutsches Zentrum für Luftund Raumfahrt(DLR) released by the multi-global navigation satellite system(GNSS) experiment, combined with the BDS-3 measured data from the International GNSS Service stations evenly distributed in different latitudes around the world for 7 days of B1I, B1C and B2a single frequency and B1C/B2a, B1I/B3I two commonly used ionospheric-free combined pseudorange single point positioning accuracy for comparison and analysis.
ResultsExperimental results show that in the three frequencies of B1I, B1C, and B2a and the two combinations of B1C/B2a and B1I/B3I, the positioning accuracy of CAS products after correction is better than that of DLR products as a whole. Specifically, the average of positioning accuracy has increased 3 cm, 1 cm, and 6 cm in the directions of E, N, and U,respectively.
ConclusionsIt can be concluded that for BDS-3 pseudorange single point positioning, the positioning accuracy of CAS products after correction is generally better than that of DLR products.
-
舰船目标作为河流与海洋的主要运输载体,对其进行定位和监测在军事侦察、资源生产管理、环境保护等领域都具有重要意义[1-5]。舰船目标检测作为海洋监测的一项重要技术,是进行舰船目标跟踪、识别的前提,在军事及民用领域都有着广袤的应用前景。随着航天技术的飞速发展,利用遥感进行舰船目标检测具有大范围任意区域观测、高时效性、数据获取高性价比等优点,是进行舰船目标感知的绝佳手段。
舰船目标实时、准确的态势感知需要结合海、陆、空大量传感数据进行探测与分析。舰船自动识别系统(automatic identification system,AIS)是最常见的舰船跟踪系统[5],它能够将舰船的位置与属性信息长时间持续地记录传输,是监测合作区域运动舰船目标的重要工具。此外,基于图像信息进行舰船目标解译,例如合成孔径雷达(synthetic aperture radar,SAR)[6]、热红外传感器[7]、光学传感器[8-9]等数据,是进行大范围舰船目标探测的重要手段。不同传感器类型都有其自身的优点和缺点,能够在不同时间、天气条件、监测场景、覆盖范围、更新频率、检测粒度等方面发挥各自的优势。SAR图像具有不受天气、云雾干扰,图像中船身金属反射特性区别于周围水面的特性,已成为海上舰船目标探测的主要手段。热红外图像通过成像对象本身的发射值来区分目标,在夜间的舰船目标监测中发挥着一定的作用。光学床图像中基于序列图像、视频图像以及基于成像延迟的单幅凝视卫星图像进行运动舰船目标检测已在动态舰船目标感知中得到了成功应用;基于高分辨率可见光的图像数据对近岸复杂区域的舰船目标进行定位识别也在舰船目标分析识别中发挥着重要作用。随着大量光学传感器的广泛使用,高空间分辨率(以下统称高分)的可见光图像数据量以几何级数增加。由于其具有获取经济便捷、细节清晰、成像与自然图像相似等优势,既能开展大规模数据下的验证评估,又能够借鉴成熟的计算机视觉目标检测方法论,基于高分可见光影像进行舰船目标检测已经成为当前主要的热点领域,相关研究呈现爆发式的增长,见图 1。
本文聚焦高分图像舰船目标检测技术,对近十年国内外有影响力文献(中文文献包括中国知网上核心期刊、硕博士论文,共718篇;英文文献包括科学引文索引(science citation index,SCI)论文、有影响力的会议论文,共234篇)的研究主题进行了统计分析。从研究主题来看,随着获取到图像的数量与质量提升,舰船目标从特殊场景目标变为现在的作为典型遥感目标检测的一部分开展研究;从早期的聚焦于海面影像的舰船定位与提取(包括云、海面、海浪等主题)逐渐转向复杂环境下的舰船检测(包括复杂场景、多尺度、类内变化多样等主题)研究;检测方法从人工设计、层层筛选定位检测(如疑似区域、部件、虚警等主题)发展为当前主流的基于卷积神经网络、一体化定位检测(包括深层特征、包围框、锚点等主题)。近年来,已有学者开展了有关可见光遥感影像舰船目标检测综述分析,普遍认为高分可见光遥感图像可以为精确的舰船感知提供高价值的特征及定位识别信息,但同时表明当前仅利用高分可见光图像进行舰船监测仍以检测为主,进行大范围、长时间、高精度、细粒度的舰船属性分析非常困难,需要结合多传感器数据进一步探究。现有方法模型聚焦于特定数据、特定场景下的目标检测研究,在实际不同传感器、不同复杂环境中开展鲁棒的舰船目标检测依然面临严重挑战。
为了比较全面地回顾近十多年高分可见光图像舰船目标检测领域的研究进展,本文从舰船检测的方法体系与发展历程概述、代表性方法归纳、检测精度评估等方面对高分可见光图像舰船检测的进展进行回顾和总结,并对高分可见光舰船检测未来研究趋势进行了展望。
1 高分可见光遥感图像舰船检测的方法体系
由于高分可见光遥感图像空间分辨率不同、覆盖的场景变化多样,许多学者为此开发了多种多样的舰船目标检测方法。本文从检测粒度、建模方法两方面对这些方法进行总结,如图 2所示。
随着图像空间分辨率的提升,舰船目标所呈现的结构纹理信息更加丰富,目标解析层次从区域级识别(图像块分类识别是否是舰船目标)、像素级提取(按一定规则获取目标像素)到对象级目标定位(获取实例级对象的位置),在舰船目标自动化解译中得到有效应用。从检测模型尤其是舰船区域/目标特征建模来分,高分可见光图像舰船目标检测可以分为基于人工特征建模的方法与基于深度特征学习的方法两类。前者利用专家知识规则对舰船目标进行特征建模表达(大量的低层、中层特征),实现舰船目标的有效判别。后者依托于人工智能技术,利用深度神经网络自动学习目标不同层次的特征进行舰船目标的检测定位,目前在遥感图像目标检测中得到越来越广泛的应用。由于舰船目标分布的特殊性和广泛性,从检测场景来分,舰船目标检测可以分为面向特殊场景的舰船目标检测与一般场景的舰船目标检测两类。
2 高分可见光遥感图像舰船检测发展历程
高分可见光图像舰船目标检测的任务是从图像中定位识别出舰船目标实例,需要解决目标“是什么”“在哪里”两大问题,遵循目标疑似区域定位、目标精细定位检测由粗到细的检测流程。在数据预处理完成后,首先进行大范围的图像舰船目标疑似区域定位,快速获取有限候选区域;然后通过对舰船目标特征精细建模在候选区域内进行进一步分析确认,完成对舰船目标的定位和判别;最后经过后处理步骤进一步剔除冗余目标,并对检测结果进行精度评价。
遥感图像目标检测方法的演化史就是对地观测技术、信息技术的发展史,遥感图像空间分辨率的提升、计算机视觉的发展、并行计算的突破,都对高分可见光遥感图像目标检测技术的发展产生了深入影响。本文对2010—2020年高分可见光图像舰船检测方法从年代、图像分辨率、检测场景、处理流程等维度进行了分析归纳,如图 3所示,分为以下3个典型的发展阶段:(1)基于中高分可见光图像舰船目标判别提取时期。2010年左右5 m分辨率的图像是用于可见光图像舰船检测的主要对象,基于底层特征的对象分析是舰船目标检测的主要思路,利用舰船目标区别于水面的独特亮度、纹理与形状特性进行统计、学习与判别、分类,从而实现中高分辨率、较为简单场景的舰船目标检测。(2)基于高分可见光图像舰船目标特征建模检测时期。2013年左右,2 m分辨率的可见光图像获取难度降低,图像中大量中小型的舰船目标信息变得丰富、容易识别,视觉注意机制、中层特征编码等机理研究成为视觉领域的热点主题,前者在大幅面遥感影像的舰船目标定位,尤其是疑似区域定位方面得到了成功应用;后者通过对底层特征的有效高层建模弥补了语义鸿沟,广泛地使用在复杂背景、类内变化多样的舰船目标建模中,从而实现高分图像中较为复杂区域的舰船检测。(3)基于深度学习的高分可见光图像舰船目标检测时期。随着深度学习在计算机视觉领域取得了革命性的成功[10],以及2016年左右大规模高分可见光图像舰船目标检测数据的涌现,舰船目标检测的整体研究思路迅速向基于深度学习的方法演变。舰船检测流程也由分步进行向当前的一体化定位检测的方向发展,卷积神经网络(convolutional neural networks,CNN)被广泛应用,其中基于分类的检测框架[11]与基于回归的检测框架[12-14]引领着目标检测体系发展潮流。
目前,基于高分可见光遥感图像的舰船检测在速度、精度、场景适应性上都得到了巨大的提升,但距离实际复杂场景图像中高效的舰船目标检测仍然有一定差距。高分可见光图像舰船目标检测存在以下难点:(1)舰船目标区域定位困难。在大范围遥感图像上,舰船目标通常只占很小一部分,冗余、复杂的背景环境(如大量无关场景、破碎岛屿、凸堤、陆地目标等)干扰给舰船目标区域有效定位带来了巨大挑战。(2)舰船目标表征难。舰船目标周围复杂的背景(云、浪、岛屿、海杂波、尾杂波、黑白极性、沿岸泥沙、光照、阴影等干扰)导致舰船目标被遮挡、弱化,舰船目标本身分布的多样性(多尺度、多类别、多形态、密集排布等干扰)导致目标特征对齐困难,这些都使得舰船目标的特征类间区分性、类内统一性建模非常困难。(3)舰船检测方法普适性不高。现有的高分可见光图像舰船检测标注数据集中,目标在尺度、类型、场景、视角分布上相对有限,舰船目标类型呈长尾特性分布,当前检测模型对难分目标,以及分布外的不同传感器、不同拍摄环境图像目标检测适应性差,难以满足大范围的自动化检测应用需求。
3 高分可见光遥感图像舰船检测方法归纳分析
3.1 基于人工特征建模的舰船检测方法
此类舰船目标检测方法逐步去除干扰因素,对于大幅面的高分可见光图像,一般需要进行舰船疑似候选区域提取从而获取目标区域,并在此基础上进行舰船目标检测定位得到精确的目标属性与位置信息。也有一些研究基于图像切片直接进行舰船目标检测定位。
3.1.1 舰船疑似候选区域提取
基于先验信息的海陆分离针对大范围水陆图像的舰船检测在有影像位置信息和相应的地理先验信息时,利用精确的地理信息辅助进行水陆分离是获取舰船目标候选区域的有效方法。其中先验信息包括港口地理信息、海岸线信息等[15],并通过精确的模板匹配或图像配准后分离出水域图像。此类方法的优点是分割效果准确稳定,不受场景、光照等因素影响。但同时缺点也很明显,当图像位置信息缺失或不准确、地理先验信息缺乏时,配准结果不再稳定,此过程就无法保证顺利进行。
基于分割的疑似区域提取对于灰度分布较为一致的水面区域,利用分割算法能够有效地实现大面积云雾去除、水陆区域分离、水面与舰船疑似区域分离,因此被广泛应用于舰船疑似区域提取。具体方法有基于全局阈值分割方法[16]、基于能量函数分割方法[17]、基于局部多阈值分割方法[18-21]等。此类方法能够有效地从简单场景图像中提取疑似舰船区域,但对于海杂波、大量碎云、海浪、小岛干扰较多的海面区域以及复杂的近岸区域,往往不能有效地进行分割,导致结果虚警和漏检都很高。
基于特征分析的疑似区域提取对于较为复杂的海面和近岸场景需要对目标进行有效特征建模来区分背景与疑似舰船区域。提取的特征大多具有一定的局部或全局不变性,如灰度[18]、颜色直方图[19]、纹理[20-22]、频域[23]、几何形状特征[24]等。此外基于生物视觉启发的视觉显著特征是一种特别的对比差异性特征,能够快速将少量舰船疑似区域从大范围的海面背景中区分出来。具体包括生物视觉注意模型[24-25,39]、频域显著模型[26]以及多种显著性融合特征[27]。特征判别方式包括统计判别方法[22,28,41-42]、决策方法[29]、异常分析方法[30]、特征学习分类方法[26-27,31-32]等。基于特征统计判别(尤其是基于显著特征判别)、异常分析的方法大多以像素为单位分析,易受不同环境因素变化的影响,同时对目标类内变化的适应能力较弱,适合分辨率较低、船体稀疏分布的目标候选区域提取。利用特征投票、决策或学习分类的方法能够一定程度上以对象为单位进行特征建模,但其提取结果受特征设计及数据影响大,很难产生较为通用的大范围复杂场景舰船候选区域提取方法。
3.1.2 舰船目标检测定位
在经过大量背景、虚警去除获取到有限的疑似舰船候选区域后,需要进一步对其进行精细分析,输出更为精细的舰船目标定位结果。此阶段面临的主要问题来自于与舰船区域相似的背景干扰(类间差异小),以及舰船内部多尺度、多类型分布偏差(类内差异大)。主要的方法包括基于统计判别的方法、基于投票判别的方法、基于特征学习的方法等。
基于统计判别的方法通过对舰船与非舰船目标的特征进行统计,基于聚类或相似性度量的方法进行舰船目标判别[20-22,33-34]。常用的特征有灰度、纹理、几何结构(长宽比、圆度、角点、边缘、形状轮廓等)等特征。此类方法基于对单个舰船的目标区域块进行分析,因此仅适用于海面场景、稀疏分布的舰船目标,对于复杂的舰船场景无法进行有效的统计判别。
基于投票判别的方法对于高分辨率的复杂近岸场景图像,基于部件表示的模型将舰船目标检测转化为利用一组对象组件、形状等特征进行相似性测量。在舰船检测中大多利用广义霍夫变换 [35-37]进行特征度量,再通过后续的投票获取最终的定位检测结果。此类方法对于平静水面场景的靠岸舰船能取得良好的检测提取结果,但并不能适用于各种复杂背景干扰如阴影、光照、海浪、云雾等,以及舰船目标多尺度、多姿态、多方向分布复杂的情况。
基于特征学习的方法为了适用于大范围、大规模的舰船目标定位检测,基于大量数据学习分类的方法被广泛使用。一些具有不变性质的特征[22,25,38-39]表示,如尺度不变特征(scale invariant feature transform,SIFT)、Gabor特征、梯度直方图(histogram of oriented gradient,HOG)特征、视觉显著特征等被广泛用于特征提取。为了弥补底层视觉特征与高层语义之间的鸿沟,基于中层特征编码建模的目标特征表征模型被广泛地使用,包括词袋模型(bag of words,BOW)[23,40]、稀疏编码[32,41-44]等。这类基于特征编码的方法能够较好地适用于复杂背景下的舰船目标判别,但大多是基于图像块的建模,因此不能进行目标的有效定位。后续基于部件模型(deformable part models,DPM)的方法[45-47]能够有效对目标进行各个部件的组合建模,从而实现目标的精确定位与检测。但这种方法需要将目标进行统一尺度的对齐,人工参与较多,对于旋转、尺度、类型差异大的舰船目标检测效果差。
3.2 基于深度特征学习的舰船目标检测方法
尽管基于深度CNN模型与传统的舰船目标检测模型在特征设计上有巨大差异,但整个检测流程后者对前者的影响深远。虽然基于深度学习的目标检测流程逐渐形成了一体化的目标检测体系,但传统方法尤其是后期特征编码学习的方法中的滑动窗口策略、多层次信息整合等思想依然在基于深度学习的舰船目标检测中起着重要作用。当前基于深度学习的高分可见光图像舰船目标方法依托于计算机视觉中目标检测领域的方法体系,在此基础上针对可见光遥感图像本身的特殊性、舰船目标的独特性进行进一步深化改进,推动着高分可见光遥感图像的舰船检测向着更加自动化、高精度、实用性强的方向发展。增强目标特征表征的有效性、多尺度定位的适应性、多姿态位置回归的稳定性、不同场景范围的鲁棒性等都是有效提升舰船目标检测性能的关键,现有的检测方法往往针对其中一点或多个点综合开展研究,本文将这4个方面的具体方法分别进行总结分析。
3.2.1 基于深度特征增强的目标特征表征方法
高分可见光遥感图像可以清晰地呈现一定范围内的目标及其背景,在带给目标更丰富细节特征的同时引入了大量的冗余信息。为了能将目标从大量复杂的负样本中分离出来,同时尽可能满足不同类型、尺度、姿态的目标表征,就需要提取更为有效的目标特征来提升目标判别的准确性。
除了加深网络层数外,增大卷积核的感受野(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],以此作为进一步约束检测框架产生更高精度的目标定位检测结果。
此类方法能充分应用目标及其上下文信息进行丰富的特征建模,但同时也不可避免引入新的问题,如新引入的超参数的敏感性与鲁棒性,新增模块与检测定位网络如何联合优化等。
3.2.2 基于多尺度特征融合的候选区域提取方法
为适应不同尺度的目标定位与识别,便于并行运算,充分利用CNN网络特征层的多尺度特性代替多尺度图像多次进行不同尺度的统一区域定位与检测是一种更为快速有效的方法。特征金字塔网络(feature pyramid network,FPN)[67]是一种典型的联合自上而下的多尺度特征增强网络,通过使用更深的层探测较大的物体,使用较浅层在较大的尺度上检测小目标被广泛用于目标检测研究中。Li等[68]构造的HSF-Net将不同尺度的特征映射到同一尺度空间进行采样回归检测,高层特征能够有效检测大尺度目标的同时,底层尺度空间特征经过与上层特征的融合提升了对于小目标的感知能力。文献[69-71]基于FPN的结构或改进不同尺度特征层的连接方式如稠密连接、跳跃连接等进行多尺度特征融合的舰船检测。
与深度特征增强的目标特征表征方法尤其是上下文信息融合的方法不同,FPN系列的多尺度特征融合模型更强调在解码阶段将不同层的特征经过采样后与原先层进行融合,能够有效地将图像中不同尺度的目标进行统一学习,提升了舰船目标检测的速度;同时通过合理分配相应尺度的目标回归检测,有效地提升了目标检测的性能。如何构建强健的多尺度特征融合结构是此类方法的重点与难点,后续很多研究在此基础上探索更好的连接方式设计[72]或进行网络架构搜索(NAS-FPN[73])学习,以追求更高效、更高精度的目标检测性能。
3.2.3 基于目标精确位置回归的检测方法
遥感目标不同于一般自然图像与人脸目标检测数据集中的竖直排布,其目标位置分布是任意方向的。舰船目标的长宽比范围比一般的目标更为多样,且目标与目标之间、目标与背景之间干扰都非常严重,影响着舰船目标的准确定位与判别。使用传统的一般水平框回归会导致目标正样本区域存在大量背景干扰使得回归位置不准确,而在后处理过程中也可能因为水平框重叠度过高导致部分检测框被删除而导致漏检。基于此,遥感领域的学者开展了大量的基于更准确目标位置回归的研究[62,74-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]来灵活学习目标的位置。对于常用的位置回归损失与真实目标预测回归框不完全对应导致回归偏差的问题,提出了基于改进目标检测分类与回归的损失函数[79,82-83]、多次逐步精细回归[80,84-85]以及改进的非最大值抑制(non-maximum suppression,NMS)[59,86]等方法进行更为有效的特征学习、模型训练以及结果优化,从而达到更高的检测精度。
此类方法通过深入分析影响舰船目标定位精度的因素,利用增强目标区域对齐的方法进行检测性能优化,有效地提升了目标定位与检测的上限。如何保证回归过程的稳定性以及提升定位的泛化能力是此类方法面临的主要问题。
3.2.4 舰船目标检测模型自适应检测方法
尽管多尺度特征金字塔的方法能够一定程度上解决目标的尺度问题,但由于遥感图像舰船目标尺度差异极大,当前将图像裁切缩放后统一尺寸地输入检测网路不可避免地导致部分过大或极小的目标无法得到有效的学习。为了缓解这种目标尺度分布规模不平衡问题,如何合理地进行多尺度候选区域获取,使得不同尺寸的目标都能够尽可能在一定的尺度范围内被检测到是此类方法的重点研究目标。Li等[87]以及You等[88]在统计小尺寸舰船目标分布的基础上进行基本图像块的确定,并基于显著性评估的方法对图像中的大尺寸船图像块缩小后进行统一目标尺度分布下的定位与检测。Yu等[89]在获取较为准确的舰船候选区后,对每个区域基于中心裁剪的方法获取多尺度图像,并以此输入卷积网络中进一步进行目标判别。Wang等[90]采用不同的块尺度对大规模舰船数据集中的目标进行检测,计算不同块尺度下不同目标检测结果的平均精度,并统计最优值获取对应的图像尺度进一步用于自适应的舰船检测以提升整体检测精度。此外,无论是模型训练还是测试阶段,随机裁剪、多尺度图像训练/测试都能极为有效地提升检测方法对不同尺度分布的目标适应性[91]。
针对实际应用特别是可见光卫星图像上的舰船检测范围相当大,真实舰船分布的区域往往只占很小一部分,直接进行目标检测由于大量的冗余背景使得正负样本分布极度不平衡,严重影响检测的准确性和时效性。因此如何自适应提取大范围图像上的舰船目标疑似区域尤为关键。此类型方法目前的研究思路都是利用舰船目标本身的场景关联特性,利用先验地理信息[88]进行港口场景区域的定位实现疑似图像区域的提取,在此基础上进一步执行舰船目标检测。
此类方法可以显式地适应不同大小、场景的舰船目标以提升检测精度,但需要有专业知识辅助学习,而对于不同场景图像之间尺度参数的适应性及迁移性研究需要进一步探索。
3.3 高分可见光遥感图像舰船目标检测精度评估
对于舰船目标的检测结果评估是高分辨率可见光遥感图像目标检测的重要一环。目前对于遥感图像舰船检测性能的定量评估主要基于对象级统计,检测的结果可能会出现正确检测、错误检测和漏检3种情况,主要用以下4个指标来评价检测的结果:查全率(Recall)、查准率(Precision)、平均精度(average precision,AP)、平均均值精度(mean average precision,MAP)。查全率R指预测值中被正确识别到的正样本数占所有正样本数量的比例;查准率P表示被正确识别的正样本在所有预测到为正样本中的占比;AP表示在不同的测试召回率下的平均精度,对应着不同的IOU(intersection over union)阈值;MAP首先对每一个类别都计算出AP,然后再计算平均值,是对检测到的目标平均精度的一个综合度量。计算如下:
式中,
(true positives)、 (false positives)、 (false negatives)分别表示检测到的目标与真值比对后产生的真阳性、假阳性以及假阴性目标个数;m表示所有样本的类别个数。预测框(Bp)位置是否预测正确(实际标注真值框(Bt))通过IOU判定: 式中,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 从表 1中可以看出,近几年高分可见光遥感图像舰船检测相关数据集发展迅速且备受瞩目,亚米级的图像被广泛应用于检测任务中,并且有由少量类型专用数据集向综合多类数据集发展的趋势,目标标注粒度也从水平框的标注向更为精确的旋转框、实例掩码标注发展。HRSC2016和Kaggle-ship是专门针对舰船目标的数据集,前者是目前公开的唯一可以进一步识别舰船类型的高分辨率图像数据集,而后者是舰船目标掩码标注的数据集。
4 结语
随着对地观测技术和信息技术的蓬勃发展,研究人员在高分可见光遥感图像舰船检测技术方面做了大量的研究工作,2010—2020年基于深度学习的高分可见光图像舰船检测无论在检测规模、检测性能还是检测自动化程度上都得到了巨大的提升。但目前的检测技术面对数据源的不确定性、算法的适应性、稳定性与可靠性等难题依然存在挑战,距离自适应、高精度、实时舰船目标检测及智能解译依然存在很大的差距,有待进一步的研究。本文对基于高分可见光遥感图像舰船目标检测的展望如下:
1)舰船目标检测性能进一步增强。目前利用高分可见光遥感图像的舰船目标检测在大规模数据集上的检测精度有待进一步提高,实际应用于大范围、多云雾干扰的近岸及海上图像舰船检测时,存在大量虚检与漏检,需要更加有效的检测算法。对于背景与目标分布极度不平衡问题、不同尺度目标尤其是小目标(二十个像素内)、极限长宽比目标、难分目标等的定位判别问题需要进一步探索。
2)舰船目标检测算法适应性提升。目前基于深度学习的舰船检测框架种类繁多、设计复杂,包含的超参数非常多且严重依赖于经验,对于检测性能影响很大。发展自动的网络结构设计、参数初始化与调整策略优化等技术是未来自适应舰船目标检测的重要方向。对于不同的数据集、不同的应用任务,目前方法需要经过部分的新数据标注及模型迁移后才能进行目标检测,过程成本开销大且效率低下,因此需要开发基于弱监督甚至无监督的图像标注与目标检测,提高舰船目标检测对不同数据源、任务的适应能力。开展受限资源下如移动平台、在轨卫星数据的舰船模型轻量化、实时检测技术研究,为实际舰船目标动态监测应用提供支撑。
3)舰船目标检测识别能力深化。实际舰船检测应用中对于舰船目标的定位及进一步的类别识别有着迫切需求,当前舰船目标类型识别主要还是依赖人工判读解译舰船类别属性。高分可见光遥感图像仅仅提供了有限视角下(下视与侧下视)的目标形态,且图像中的目标往往存在模糊、扭曲、遮挡的情况,仅依赖高分可见光图像进行舰船目标识别相当困难。建立权威、完善的大规模舰船目标分类数据集是实现自动化舰船目标识别的基础。同时需要融合多源遥感信息(如AIS、凝视卫星图像、SAR图像、光学图像等)实现对舰船目标多维度特征互补关联,提升任意时间、场景下的智能化舰船目标解译能力。
http://ch.whu.edu.cn/cn/article/doi/10.13203/j.whugis20210454
-
表 1 CAS产品和DLR产品改正后单频伪距单点定位精度RMS均值/m
Table 1 RMS Mean Values of Single-Frequency Pseudorange Single Point Positioning Accuracy After the Correction of CAS Products and DLR Products /m
测站 CAS产品改正后 DLR产品改正后 B1I B1C B2a B1I B1C B2a E N U E N U E N U E N U E N U E N U CHPG 0.79 0.98 2.20 0.69 0.97 1.85 1.25 1.74 4.78 0.80 0.98 2.23 0.72 0.94 1.90 1.27 1.71 4.80 FAA1 0.81 2.30 2.90 0.78 2.25 2.68 1.23 4.19 5.87 0.82 2.30 2.89 0.81 2.23 2.78 1.25 4.16 5.95 KOUG 1.04 0.98 1.92 0.99 0.93 1.78 1.75 1.62 5.16 1.04 0.98 1.92 1.01 0.94 1.86 1.75 1.62 5.20 MAS1 0.72 2.38 2.53 0.67 2.34 2.43 1.04 4.23 5.38 0.72 2.39 2.53 0.70 2.35 2.50 1.05 4.25 5.45 NKLG 0.79 1.06 3.35 0.73 1.03 3.21 1.19 1.82 6.40 0.79 1.06 3.34 0.77 1.05 3.22 1.21 1.84 6.41 THTG 0.79 2.29 3.21 0.73 2.24 2.97 1.23 4.18 6.53 0.80 2.28 3.21 0.77 2.22 3.04 1.25 4.15 6.60 FFMJ 0.51 0.84 1.95 0.45 0.79 1.86 0.44 1.38 1.30 0.52 0.84 1.96 0.50 0.80 1.88 0.48 1.37 1.43 BRST 0.55 0.89 1.75 0.48 0.85 1.52 0.51 1.44 2.36 0.56 0.90 1.78 0.53 0.86 1.67 0.56 1.44 2.56 ENAO 0.52 0.8 3.17 0.45 0.72 2.78 0.46 1.35 2.10 0.53 0.80 3.19 0.50 0.75 3.01 0.51 1.38 2.48 WTZZ 0.89 1.29 2.54 0.72 1.15 2.18 0.69 1.63 1.84 0.91 1.31 2.56 0.84 1.25 2.36 0.84 1.68 2.16 RIO2 0.48 0.73 2.25 0.48 0.67 2.26 0.56 1.24 1.58 0.48 0.74 2.27 0.54 0.70 2.36 0.62 1.26 1.67 TIT2 0.52 0.85 1.97 0.45 0.79 1.91 0.45 1.37 1.29 0.53 0.86 1.98 0.50 0.80 1.93 0.49 1.37 1.42 WUH2 0.63 1.06 2.49 0.59 1.00 2.34 0.77 1.98 3.21 0.64 1.07 2.49 0.64 1.04 2.41 0.79 2.02 3.24 KIRU 0.65 0.62 2.27 0.61 0.64 2.26 0.57 0.62 1.44 0.64 0.62 2.27 0.62 0.68 2.34 0.60 0.64 1.57 MET3 0.51 0.69 2.25 0.45 0.58 2.16 0.44 0.88 1.35 0.52 0.70 2.27 0.48 0.61 2.18 0.47 0.87 1.45 METG 0.48 0.58 1.18 0.45 0.52 1.07 0.43 0.86 1.84 0.49 0.58 1.17 0.48 0.55 1.14 0.46 0.87 1.93 NYA2 0.47 0.45 2.96 0.45 0.43 2.92 0.63 0.51 2.02 0.47 0.45 2.95 0.47 0.47 2.98 0.64 0.54 2.08 SOD3 0.57 0.67 2.46 0.54 0.58 2.31 0.57 0.62 1.55 0.57 0.68 2.47 0.56 0.62 2.37 0.59 0.64 1.67 YEL2 0.44 0.64 2.00 0.41 0.63 1.98 0.52 0.71 1.42 0.44 0.64 2.00 0.44 0.65 2.01 0.52 0.73 1.52 表 2 CAS产品和DLR产品改正后双频无电离层组合伪距单点定位精度RMS均值/m
Table 2 RMS Mean Values of Single Point Positioning Accuracy of Dual-Frequency Ionospheric-Free Combined Pseudorange After the Correction of CAS Products and DLR Products /m
测站 CAS产品改正后 DLR产品改正后 B1C/B2a B1I/B3I B1C/B2a B1I/B3I E N U E N U E N U E N U CHPG 0.63 0.62 2.06 1.19 1.17 4.69 0.66 0.63 2.13 1.28 1.21 4.78 FAA1 0.85 0.81 3.06 1.19 1.11 4.10 0.86 0.80 3.11 1.19 1.10 4.09 KOUG 0.60 0.48 1.78 1.01 0.82 2.55 0.64 0.47 1.79 1.00 0.8 2.52 MAS1 0.65 0.66 1.78 0.96 0.95 2.58 0.70 0.68 1.82 0.96 0.94 2.55 NKLG 0.60 0.46 4.03 0.88 0.77 4.10 0.64 0.46 4.05 0.88 0.78 4.11 THTG 0.65 0.60 3.19 1.07 1.01 4.28 0.67 0.60 3.19 1.06 1.00 4.26 FFMJ 0.51 0.61 1.51 0.84 0.98 2.13 0.58 0.70 1.66 0.87 1.02 2.20 BRST 0.62 0.83 2.91 1.08 1.40 3.99 0.66 0.88 2.97 1.13 1.45 4.17 ENAO 0.48 0.61 2.09 1.06 1.12 3.79 0.54 0.63 2.25 1.07 1.12 3.89 WTZZ 0.98 1.26 2.51 1.58 1.98 3.77 1.07 1.38 2.77 1.62 2.01 3.83 RIO2 0.47 0.63 1.29 0.51 0.66 1.41 0.55 0.66 1.34 0.50 0.63 1.36 TIT2 0.51 0.63 1.40 0.83 0.97 1.96 0.58 0.72 1.53 0.85 1.00 2.02 WUH2 0.69 0.76 1.82 1.41 1.18 3.23 0.73 0.75 1.85 1.43 1.19 3.23 KIRU 0.65 0.75 1.83 0.86 1.08 2.57 0.66 0.78 1.88 0.87 1.08 2.58 MET3 0.53 0.67 1.26 0.97 1.35 2.58 0.54 0.69 1.30 0.99 1.37 2.65 METG 0.46 0.58 2.36 0.82 0.84 2.49 0.49 0.61 2.39 0.83 0.83 2.45 NYA2 0.44 0.47 1.70 0.51 0.71 2.35 0.50 0.54 1.77 0.53 0.72 2.33 SOD3 0.47 0.62 1.42 0.83 1.21 2.92 0.49 0.66 1.49 0.85 1.23 3.00 YEL2 0.48 0.50 1.17 0.57 0.75 1.53 0.53 0.51 1.21 0.57 0.74 1.53 -
[1] Dach R, Lutz S, Walser P, et al. Bernese GNSS Software Version 5.2[M]. Berne: Astronomical Institute of University of Berne, 2015.
[2] Håkansson M,Jensen A B O,Horemuz M,et al. Review of Code and Phase Biases in Multi-GNSS Positioning[J]. GPS Solutions, 2017, 21(3) : 849-860.
[3] Wilson B D , Yinger C H , Feess W A , et al. New and Improved: The Broadcast Interfrequency Biases[J]. GPS World, 1999, 10(9):56-66.
[4] 王宁波,袁运斌,张宝成, 等.GPS民用广播星历中ISC参数精度分析及其对导航定位的影响[J].测绘学报,2016,45(8):919-928. Wang Ningbo,Yuan Yunbin,Zhang Baocheng, et al.Accuracy Analysis of ISC Parameters in GPS Civil Broadcast Ephemeris and Its Influence on Navigation and Positioning[J].Acta Geodaetica et Cartographica Sinica,2016,45(8):919-928.
[5] 曾添,隋立芬,鲍亚东, 等.BDS卫星端差分码偏差对定位的影响及改正模型研究[J].大地测量与地球动力学,2017,37(1):53-57. Zeng Tian,Sui Lifen,Bao Yadong, et al.Study on the Influence of BDS Satellite Differential Code Deviation on Positioning and the Correction Model[J].Journal of Geodesy and Geodynamics,2017,37(1):53- 57.
[6] 张强,赵齐乐,章红平, 等.利用北斗观测实验网解算北斗卫星差分码偏差[J].武汉大学学报(信息科学版),2016,41(12):1649-1655. Zhang Qiang,Zhao Qile,Zhang Hongping, et al.Using BeiDou Observation and Experiment Network to Solve the Difference of BeiDou Satellite Differential Code Deviation[J]. Geomatics and Information Science of Wuhan University,2016,41(12):1649-1655.
[7] Li X , Xie W , Huang J , et al. Estimation and Analysis of Differential Code Biases for BDS3/BDS2 Using iGMAS and MGEX Observations[J]. Journal of Geodesy, 2019, 93(3):419-435.
[8] Montenbruck O , Hauschild A , Steigenberger P . Differential Code Bias Estimation Using Multi‐GNSS Observations and Global Ionosphere Maps[J]. Navigation, 2015, 61(3):191-201.
[9] Wang N , Yuan Y , Li Z , et al. Determination of Differential Code Biases with Multi-GNSS Observations[J]. Journal of Geodesy, 2016, 90(3):209-228.
[10] Li Z, Yuan Y, Hui L, et al. Two-step Method for the Determination of the Differential Code Biases of COMPASS Satellites[J]. Journal of Geodesy, 2012, 86(11):1059-1076.
[11] Lou Y , Zheng F , Gu S , et al. Multi-GNSS Precise Point Positioning with Raw Single-Frequency and Dual-Frequency Measurement Models[J]. GPS Solutions, 2016, 20(4):849-862.
[12] Guo F , Zhang X , Wang J . Timing Group Delay and Differential Code Bias Corrections for BeiDou Positioning[J]. Journal of Geodesy, 2015, 89(5):427-445.
[13] 舒宝,刘晖,张明, 等.北斗系统硬件延迟解算及精度分析[J].武汉大学学报(信息科学版),2016,41(2):279-284. Shu Bao,Liu Hui,Zhang Ming, et al.BeiDou System Hardware Delay Calculation and Accuracy Analysis[J]. Geomatics and Information Science of Wuhan University,2016,41(2):279-284.
[14] 张辉, 郝金明, 刘伟平, 等. 估计接收机差分码偏差的GPS/BDS非组合精密单点定位模型[J]. 武汉大学学报 ( 信息科学版), 2019, 44(4): 495-500, 592. Zhang Hui, Hao Jinming, Liu Weiping, et al. GPS/BDS Precise Point Positioning Model with Receiver DCB Parameters for Raw Observations[J]. Geomatics and Information Science of Wuhan University, 2019, 44(4): 495-500, 592.
[15] 梅登奎,闻德保.MGEX北斗差分码偏差两种精确处理方法对比分析[J].空间科学学报,2019,39(5):662-669. Mei Dengkui,Wen Debao.Comparison Analysis of Two Accurate Processing Methods for MGEX BeiDou Differential Code Deviation[J].Journal of Space Science,2019,39(5):662-669.
[16] 谷世铭,党亚民,王虎, 等.北斗差分码偏差改正对单点定位的影响[J].测绘科学,2020,45(10):10-15. Gu Shiming,Dang Yamin,Wang Hu, et al.The Influence of BeiDou Differential Code Deviation Correction on Single Point Positioning[J].Science of Surveying and Mapping,2020,45(10):10-15.
[17] 王利军,焦文海,贾小林, 等.BDS-3精密单点定位性能比较分析[J].大地测量与地球动力学,2021,41(4):357-361. Wang Lijun,Jiao Wenhai,Jia Xiaolin, et al.Comparative Analysis of BDS-3 Precision Single Point Positioning Performance[J].Journal of Geodesy and Geodynamics,2021,41(4):357-361.
[18] 邓远帆,郭斐,张小红, 等.北斗三号卫星多频多通道差分码偏差估计与分析[J].测绘学报,2021,50(4):448-456. Deng Yuanfan,Guo Fei,Zhang Xiaohong, et al.BeiDou-3 Satellite Multi-frequency and Multi-channel Differential Code Bias Estimation and Analysis[J]. Acta Geodaetica et Cartographica Sinica,2021,50(4):448-456.
[19] Zhang B C, Teunissen P J G, Yuan Y B,et al. A Modified Carrier-to-Code Leveling Method for Retrieving Ionospheric Observables and Detecting Short-Term Temporal Variability of Receiver Differential Code Biases[J]. Journal of Geodesy,2019,93(1):19-28.
[20] Zhang B C, Teunissen P J G, Yuan Y B. On the Short-Term Temporal Variations of GNSS Receiver Differential Phase Biases[J]. Journal of Geodesy,2017,91(5):563-572.
[21] 袁海军,章浙涛,何秀凤,等.北斗三号卫星差分码偏差稳定性分析及其对单点定位的影响[J].武汉大学学报(信息科学版),2023, 48(3): 425-432. Yuan Haijun, Zhang Zhetao, He Xiufeng, et al. Stability Analysis of Differential Code Deviation of BeidDou-3 Satellite and Its Influence on Single Point Positioning [J]. Geomatics and Information Science of Wuhan University, 2023, 48(3): 425-432.
[22] 徐磊,常国宾,高井祥,等.附加闭合差约束的BDS频间偏差估计模型[J].武汉大学学报(信息科学版),2021,46(4):520-529. Xu Lei, Chang Guobin, Gao Jingxiang, et al.BDS Inter-Frequency Deviation Estimation Model with Closure Constraint [J]. Geomatics and Information Science of Wuhan University, 2021,46 (4): 520-529.
[23] 林高宇,王磊,何飞扬,等.Swarm低轨卫星星座的GPS接收机差分码偏差估计[J/OL].武汉大学学报(信息科学版),2023, 48(1): 119-126. Lin Gaoyu, Wang Lei, He Feiyang,et al . Differential Code Deviation Estimation of GPS Receiver for Swarm LEO Satellite Constellation [J]. Geomatics and Information Science of Wuhan University, 2023, 48(1): 119-126.
[24] 付宁波,陈艳梅,黄霞飞.基于GPS跟踪站的精密单点定位精度分析[J].测绘地理信息,2015,40(4):29-31. Fu Ningbo, Chen Yanmei, Huang Xiafei. Precision Analysis of Precise Point Positioning Based on GPS Tracking Station[J]. Journal of Geomatics, 2015,40(4):29-31.
-
期刊类型引用(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)