高分三号NSC模式SAR图像舰船目标检测初探

刘泽宇 柳 彬* 郭炜炜 张增辉 张 波 周月恒 马 高 郁文贤

①(上海交通大学智能探测与识别上海市重点实验室 上海 200240)

②(中国空间技术研究院西安分院 西安 710100)

高分三号NSC模式SAR图像舰船目标检测初探

刘泽宇①柳 彬*①郭炜炜①张增辉①张 波②周月恒②马 高②郁文贤①

①(上海交通大学智能探测与识别上海市重点实验室 上海 200240)

②(中国空间技术研究院西安分院 西安 710100)

高分三号卫星是我国首颗分辨率达到1m的C波段多极化合成孔径雷达(Synthetic Aperture Radar,SAR)卫星，拥有多种成像模式。该文针对高分三号NSC模式SAR图像提出一种海上舰船目标检测方法，其核心为基于贝叶斯框架的像素分类以实现目标筛选，并根据数据特点设计有效的图像降质条件下的性能提升方法。该文提出的检测算法与多种恒虚警率(Constant False Alarm Rate, CFAR)检测算法进行对比实验分析，实验结果证明了该文所提方法的有效性与性能优势。

高分三号卫星；合成孔径雷达；舰船目标检测；像素分类

1 引言

高分三号卫星(GF-3 satellite)是我国最新自主研发的多成像模式、多极化C波段合成孔径雷达(Synthetic Aperture Radar, SAR)卫星，其设计寿命长达8年。高分三号的SAR成像模式覆盖聚束、条带、扫描、超宽幅、扩展入射角等共12种。目前，针对高分三号，发表的研究包括系统技术[1]、几何定位[2]、海浪定量遥感[3]、海表面流场测量[4]、海洋内波测量[5]等。

在高分三号的多种成像模式中，窄幅扫描(Narrow Scan, NSC)模式成像幅宽可以达到300 km以上，其标准产品的距离向空间分辨率30～60 m，方位向分辨率50～60 m[1]。这种模式的产品结合中等分辨率与大幅宽的特点，适合海洋监测和舰船监视等任务。结合高分三号NSC模式数据特点，本文选取舰船目标检测问题作为主要关注点。

经过多年的发展，SAR图像舰船检测方法层出不穷，其中最为广泛研究和应用的是恒虚警率(Constant False Alarm Rate, CFAR)检测算法。CFAR检测算法利用的是SAR图像数据中背景像素(相对于舰船目标而言)的统计特征，对于图像中背景统计特征描述的准确程度是影响CFAR检测性能的关键。具体地，影响背景建模精度的因素主要有：(1)描述背景统计特性所用的概率分布模型；(2)用于背景建模的像素选择方法。以上两个因素也正是SAR图像舰船检测CFAR方法近年来的主要研究方向。为了适应不同海况的统计特征，多种统计分布模型，如G0分布、K分布[6]、Weibull分布、伽马分布[7]、广义伽马分布(Generalized Gamma Distribution, GΓD)[8]、有限混合分布模型，常被用于对SAR图像中背景部分统计建模。其中，基于K分布[9]和GΓD的CFAR方法[10,11]被证明可适用于多种海况的舰船目标检测。为了去除目标像素对背景统计模型参数估计的影响，基于局部窗口内背景像素进行统计建模的CFAR检测器得以发展，其中的典型应用是各类自适应滑窗CFAR检测器，如艾加秋等[12]提出的改进双参数CFAR检测器，种劲松等[13]将局部窗口引入K-CFAR并证明对于海况复杂的图像检测性能有所提升。

SAR图像舰船检测方法还包括将计算机视觉领域其他方法迁移到SAR图像中进行应用。例如，形状信息先验方法[14]，视觉显著度衡量方法[15]，利用上下文信息[16]进行目标检测等。通常情况下，这类特征应用于高分辨率SAR图像目标检测，并且为适应SAR图像特点需要事先进行特征改进。除此之外，基于机器学习，特别是深度学习SAR图像舰船检测方法是目前研究前沿方向，已有研究结合卷积神经网络对SAR图像舰船检测开展研究[17]，但目标样本不足仍是SAR图像目标认知中学习方法亟待克服的障碍。

各种SAR图像舰船检测方法对不同传感器获取的SAR图像适应能力有所不同，本文主要针对高分三号NSC模式的SAR图像数据进行舰船检测方法的探索。NSC模式L1级产品可从多方面进行图像质量提升，本文针对两种影响舰船检测性能的因素分别设计图像质量提升算法，并将图像质量提升步骤融入检测流程。本文所提出的舰船检测方法也是一种基于SAR图像统计特征的检测方法。与CFAR方法相比，不同点在于本文方法是基于像素分类方法对舰船目标进行筛选，该方法既考虑背景的复杂性，在模型构建中也考虑到目标干扰，提升了稳定性。在真实数据实验中，本文将设计的检测算法与G0-CFAR, K-CFAR, GΓD-CFAR等常用检测方法进行性能比较分析，结果证明本文所提方法对NSC模式的SAR图像舰船检测比常用CFAR方法有更好的适应能力。

2 方法

2.1 方法概要介绍

根据SAR图像的分辨率和图像质量特点，本文提出一种适应NSC模式SAR图像的海面舰船目标检测流程。整个流程如图1所示，包括以下4个主要环节。(1)图像质量提升：本文选取高分三号NSC模式图像中的条带噪声和旁瓣影响因素分别设计图像质量提升方法，提升海面舰船检测流程稳定性和数据适应能力。(2)海陆分割：海陆分割是舰船目标检测的基本前提。复杂海况场景中，依据幅值或强度进行海陆分割在真实数据测试中性能不稳定，可遵循纹理信息进行海陆分割。本文采用局部邻域幅值方差作为衡量纹理强弱的指标，选择混合高斯模型对纹理信息进行建模分类提取海面感兴趣区域(Region Of Interest, ROI)。(3)基于像素分类的目标筛选：本文结合SAR图像统计分布特征和空间邻域相关分析对像素进行分类，并组合目标类像素来筛选得到疑似舰船目标。(4)目标鉴别：在获得海面ROI区域的候选舰船目标后，本文依据目标几何参数(如舰船目标的长、宽等)和散射特性(如幅值，强度等)对虚警进行剔除。对于受旁瓣影响的海域候选目标，我们采图像质量提升步骤中去除旁瓣干扰的目标核心部分来提取几何参数和散射特性，并通过对比实验证明图像质量提升算法对检测性能提高有积极作用。

本文提出的检测方法主要特点为：

(1) 针对高分三号NSC模式L1A级图像方位向上存在规律分布的条带噪声，本文提出基于1维傅里叶变换的频谱分析方法来抑制条带噪声对后续统计建模过程的影响。条带噪声对应频谱中频分量的峰值，本文采用大津法[18]和霍夫变换检测直线[19]的方法自动提取频谱中条带噪声对应的中频分量，对检测到的中频分量采用带阻滤波。最终将经过逆变换恢复成抑制条带噪声后的SAR图像，图2为局部海域的SAR幅值图像，图3即为图2抑制条带噪声后的结果。

图1 高分三号NSC模式L1A级产品海面舰船检测全流程Fig. 1 Flow chart of maritime ship detection on GF3 NSC imaging mode, L1A product

(2) 针对SAR图像舰船目标提取几何参数和散射特性易受旁瓣干扰的问题，本文提出一种“十字”旁瓣定位和提取受旁瓣干扰舰船目标核心区域的质量提升算法，并辅助目标鉴别过程减少漏检。常见方法是基于目标二值分割结果进行形态学滤波或局部插值减弱旁瓣干扰，本文提出一种基于辐射信息排除旁瓣干扰的方法。本文方法排除干扰的同时不改变原本的图像数据，数据不会偏离原有概率分布，对后续统计建模分析有利。SAR图像中影响舰船目标的“十字”旁瓣可以看作是垂直相交的线目标组合。Tupin等[20]设计了一种适合于SAR图像的线检测器，在此基础上我们改进线检测器(如图4(a)所示)，无需形状先验，可直接定位“十字”旁瓣；接着按照辐射信息相似的原则进行区域增长；从而，可提取受“十字”旁瓣干扰舰船目标核心部分的连通域，具体效果如图4(b)， 图4(c)， 图4(d)所示。通过目标核心部分的连通域能够估计得到更准确的几何参数，从而排除旁瓣对舰船目标的干扰，降低漏检和虚警。

(3) 本文采取的目标筛选手段是基于有限混合模型和像素间邻域相关分析的像素分类方法，适合复杂海况场景中舰船目标检测。基于像素分类进行目标筛选的方法是本文所提方法的主要创新部分，下节中将作详细说明。

图2 局部海域的SAR幅值图像Fig. 2 SAR amplitude image of marine area

图3 抑制条带噪声后的幅值图像Fig. 3 Amplitude image after suppressing strip noise

图4 改进滤波器和提取受旁瓣影响的舰船目标核心部分示意图Fig. 4 Modified filter and extraction result of ship target core region after region growing based on center region of sidelobes

2.2 基于像素分类的目标筛选方法

CFAR检测方法对目标所在背景进行统计建模，并通过找到偏离背景像素值概率分布模型的像素作为目标，这种方法容易受到多种干扰导致性能降低。在舰船目标近岸、海况复杂、目标密集等情况下背景建模容易受到干扰。考虑到上述各种复杂情况，我们采用有限混合模型对舰船目标所在的复杂场景进行统计建模，并结合空间相关性对于场景像素进行分类，从而实现目标筛选，这一方法的优势在于：(1)将目标作为有限混合模型中的一类，无需预先选择背景像素样本。(2)受相干斑的影响，仅考虑各像素统计特征的分类方法性能会受到限制，本文所提方法将各像素强度统计特征与像素间的空域相关特性相结合，在贝叶斯框架下，在分类算法中融合后向散射与空域相关两类信息进行疑似目标筛选，能更好地适应受相干斑影响的情况。该算法主要包括两个关键部分：一是对SAR图像强度数据进行统计建模，按照有限混合伽马模型进行参数估计，二是按照最大后验概率(MaximumA Posteriori, MAP)准则考虑空间相关对像素进行分类。

对于有限混合模型，定义y为观测数据集，那么混合模型的概率密度函数可以定义为：

其中K为混合模型中单一模型的个数，pi为有限混合模型中第i个混合分量的权重因子，权重因子应满足且为第i个混合分量的单一概率密度函数。表示有限混合模型中第i个混合分量成分中所包含的待估计的参数组合。

在满足完全发展的相干斑[21]噪声条件下，SAR图像(视数为L)的强度值服从伽马分布，其统计模型为：

基于有限混合伽马模型的参数估计得到的是统计意义上与数据相似的图像描述方式，统计建模仅考虑每个像素的强度信息，但像素类别不仅由本身的后向散射强度决定，还与其空间邻域像素相关。令S代表整幅图像网格，IS代表强度图，Y={y，y∈IS}代表强度图像中各像素值，KS代表类别标号图像，X={x，x∈KS}代表各像素的类别标号。本文所提的像素分类过程可视为对强度图像像素类别标号进行估计的问题。根据MAP准则和贝叶斯公式，x的取值应使后验概率：

达到最大值，亦即：

其中K为类别数，根据混合伽马建模的结果，提示x的取值可为为在给像素给定类别编号时强度数据的条件概率密度，为类别标号出现的先验概率。基于像素空间相关性的考虑，像素最终的类别标号需要由像素强度观测值y和该像素邻域标号决定，因此式(4)可改写为：

只考虑像素标号与其邻域标号相关，随机场X符合关于邻域η的MRF。MRF与吉布斯随机场的等效性在文献[23]中经过详细证明，式(5)中对应的先验概率分布可描述为：

式(7)全局最优解需要极大的计算量，本文采用迭代条件模型法(Iterated Conditional Model,ICM)求解MAP准则下局部最优解。ICM算法符合MRF先验和MAP准则，在ICM迭代过程中，K类伽马分布概率密度函数重新估计，更新μi(i=1，2，···，K)，再利用更新的概率密度函数和类别标号共同确定分类的后验概率，在满足像素标号两次迭代之间改变小于一定比例或迭代次数达到固定次数即认为已经收敛，此时迭代终止，像素分类完成。ICM算法分类结果受初值影响较大，因此我们采用混合伽马模型估计结果作为像素类别标号的初始值。本文参考了文献[24]中对于ICM方法与MRF随机场的详细说明。

考虑到SAR图像中至少包括平静海域、起伏海域、平坦陆地和人造目标4类地物，因此将混合模型预设为4类混合(即K=4)。基于像素分类的目标筛选算法中结合了散射信息和空域相关信息，最终提取强度最高的一类像素作为目标筛选的结果，并与海域ROI部分取交集获得海域疑似目标，最终按照预先设定的筛选规则对疑似目标逐个鉴别获得舰船检测结果。

3 实验

3.1 实验数据情况

NSC模式的L1A级产品是单视复数据(Single Look Complex, SLC)图像，方位向像元尺寸20 m左右，地距向像元尺寸5 m左右。本文选取4幅成像区域内海域占比高、舰船目标丰富的图像作为实验图像。本文获得的4幅实验图像中舰船目标真值总计为111个。图像中的真值采用人工标注方法，标注时参考成像区域光学卫星影像排除SAR图像中陆地、港口、岛屿、礁石等与舰船目标易混淆的问题。

3.2 实验结果与分析

受海浪、风速等多重因素影响，SAR图像中海域部分并非完全均匀。为了适应复杂海况的舰船检测，除本文引文部分提到的K分布模型和广义伽马模型外，基于G0分布的CFAR方法也被证实可应用于SAR图像舰船检测[25]。因此，本文选取基于G0分布，K分布，GΓD的3种CFAR检测器(G0-CFAR, K-CFAR, GΓD-CFAR)作为基线方法与本文所提方法进行比较。在整个检测流程中，性能比较的核心在于目标筛选步骤的比较。因此，我们将CFAR目标检测方法与本文所提的融合混合伽马建模和求空域相关的像素分类目标检测方法做性能比较。为了公平起见，检测流程的其他步骤操作和参数一致。除此之外，实验过程中我们对CFAR检测器的设定如下：(1)CFAR检测中包括对海域背景进行建模，性能比较过程中需基于检测流程中获得的海陆分割结果作为先验排除陆地对建模的干扰。(2)除陆地区域外，我们还需要排除海上目标像素对建模的干扰。本文采用阈值筛选方法剔除目标像素点，阈值经过多幅数据调优。(3)各种CFAR检测器设定的恒虚警率相同，虚警率设定为Pfa=10–3，在此情况下，各CFAR检测器均能取得较好性能。

为了比较检测性能，本文采用“Precision-Recall”指标组合来评价检测方法的性能。其中“Precision”又称查准率，计算公式为“正确检测数/检测数”，代表检测结果中是真正舰船目标的比例。“Recall”又称查全率，计算公式为“正确检测数/真值数”，代表检测得到的真正舰船目标能涵盖全部真值的比例。Precision越高，代表检测结果中真实舰船目标的比例越高，检测方法对目标和虚警的辨别更灵敏。Recall越高，代表图像中真实舰船被检测到的比例越高，检测方法在图像中搜寻真实舰船的能力越强。

表1中列出了高分三号NSC模式SAR图像实验数据使用4种检测方法的指标对比，从表1中经过图像质量提升后4种检测方法指标可以发现，本文所提方法的虚警和漏检数目在4种方法中最少，“Precision-Recall”指标也是最优。通过仔细比较不同方法结果的虚警和漏检总结发现：CFAR检测结果中比本文检测方法结果多出的虚警主要是礁石和海浪，而起伏海域区域内的舰船目标则是CFAR检测结果中比本文检测方法结果多出的漏检。本文采用的基于有限混合伽马建模和空间相关信息融合的像素分类方法能够将海域中复杂的地物背景组成情况和目标进行更精细的区分，因此可排除礁石和平静海面上起伏等易混淆的虚警，也可以更好区分起伏海域背景和舰船目标。

本文设计的实验中除CFAR检测方法与本文基于像素分类的检测方法性能比较外，还对本文第2节中所提出的图像质量提升方法对检测性能的影响进行比较实验。表1中包括4种检测方法加入图像质量提升步骤前后的性能对比。性能对比的共同点是“Precision-Recall”指标均有提升。不同点是对于3种CFAR方法，主要提升是漏检数目减少，与之相关的是排除“十字”旁瓣对舰船目标形状参数估计的影响，从而减少因形状参数估计超出范围剔除目标造成的漏检。而对于本文提出的基于混合伽马建模和考虑空域相关性的检测方法，虚警和漏检数目在图像质量提升后均有明显下降。这是因为本文所提方法中混合伽马建模对图像中的条带噪声敏感，抑制条带噪声后的混合建模结果才能够对应预设的4类混合，对各类像素无噪强度真值估计准确程度得以改善，可有效减少虚警，融合排除“十字”旁瓣的干扰可进一步减少漏检。

表1 高分三号NSC模式SAR图像舰船检测指标统计Tab. 1 GF-3 satellite NSC mode SAR image maritime ship detection results

为说明图像质量提升前后本文所提检测方法的性能差距，我们选取了一幅实验图像中典型海域检测结果进行展示。其中的舰船目标真值是40个，图5 为未经过图像质量提升的检测结果，图6 为经过抑制条带噪声和排除“十字”旁瓣干扰后的检测结果。图5和图6中正确检测的舰船目标以绿色矩形框标出，漏检的舰船目标以红色圆圈标出，虚警以黄色矩形框标出。经过图像质量提升后，受“十字”旁瓣干扰造成的漏检减少，仍然存在的漏检是因为靠岸舰船与陆地紧凑相连，在海陆分割的过程中未能将舰船目标与陆地分开。本文主要关注复杂海况条件下海面舰船的检测，针对这点将在结论部分展开阐述。经过图像质量提升后，海面的虚警目标减少，这是因为抑制条带噪声后，混合伽马建模和考虑空域相关的像素分类结果更接近于真实的地物分类情况。

本文所提检测方法的复杂度主要体现在混合伽马建模参数估计和考虑空域相关的迭代优化过程，与本文流程中这一步骤相对应的是CFAR检测方法中的参数估计，建模和阈值求取步骤。我们以4幅尺寸大小接近1500×1500的实验图像来测试本算法的平均时间效率。在CPU为Intel(R) Core-i5 3.30 GHz，内存为2 GB的硬件环境下，采用Windows 10操作系统下MATLAB 2016平台实现算法。G0-CFAR, K-CFAR, GΓD-CFAR的参数估计，建模和阈值求取的平均每幅图像耗时接近1 s，本文所提方法中混合建模和考虑空间相关性的步骤采用EM算法和ICM算法，需要对全图中逐像素点进行迭代运算。其中，混合建模过程平均每幅图像耗时接近3 s，考虑空间相关的迭代优化过程平均每幅图像耗时接近14 s。本文所提方法中混合建模与CFAR过程的耗时接近，考虑空间相关的迭代优化过程是主要耗时步骤，可在未来采用更为高效的优化方案。通过以上实验结果分析可知，本文的算法的检测效果在结合图像质量提升操作后明显优于G0-CFAR, K-CFAR, GΓD-CFAR等公认稳定的检测算法，能够适应复杂海况的舰船目标检测，真实数据实验结果证实了本文所提算法适合高分三号NSC数据。

4 结论

本文提出了一种适用于高分三号NSC成像模式的海面舰船检测方法，整个检测流程的核心包括：(1)根据数据特点设计适应图像降质条件下的性能提升算法，包括抑制条带噪声和排除“十字”旁瓣干扰。算法遵循原始数据损失最小原则，实验结果证明性能提升算法的有效性；(2)提出基于贝叶斯框架像素分类的目标筛选方法。与目前公认稳定的G0-CFAR, K-CFAR, GΓD-CFAR检测器进行性能比较，结果证明本文所提方法对高分三号NSC模式数据适用性更强。

图5 未经过图像质量提升的检测结果Fig. 5 Detection result without image quality enhancement

图6 经过图像质量提升的检测结果Fig. 6 Detection result after image quality enhancement

传统的CFAR检测器和本文所提的融合混合伽马建模和考虑空域相关性的像素分类方法均属于利用统计特征进行目标检测的算法，不同之处在于CFAR需要预先设定规则将用于建模的背景像素筛选出来，需要首先对背景像素有一个先验，而本文方法直接将舰船目标类像素作为有限混合模型中的一类，无需预先筛选背景像素再进行统计建模。在对SAR图像进行建模和描述时，考虑到各种干扰因素的同时存在，本文方法中采用的混合模型精确程度超过单一模型，能够更好的反映真实复杂的海况，结合考虑空域相关特性，本文所提方法可以从SAR图像上提取更为完整的目标连通域，为目标特征提取打下基础。但本文所提方法也存在着以下几点不足：(1)本文所提方法中混合伽马模型预设4类组合，但复杂海况下，背景和目标像素很难用4类全部代表，平静海域中，4类又会存在冗余。因此，对于混合建模需要继续研究自动估计或搜寻合适的类别数目。(2)尽管本文所提的检测方法在复杂海况下舰船检测性能优于几种CFAR检测器，但精确的像素分类结果是以迭代过程的耗时为代价。从实验结果耗时情况分析可知，混合建模与CFAR过程的耗时接近，考虑空间相关时采用的ICM算法占用最多时间。在后续的改进过程中，我们考虑新的寻找最优解的快速算法来取代ICM算法。(3)本文所提流程中海陆分割也是舰船目标检测的关键步骤，在本文流程框架下，海陆分割结果决定了近岸和靠岸舰船是否可以正常检测。尽管本文按照目前实验数据特点关注的是海面舰船检测，但靠岸舰船检测问题仍需要继续利用更多数据来验证本文方法的适用性，并在此基础上提出改进的流程和方法，以期推广到高分三号更多的成像模式数据中。

[1]张庆君. 高分三号卫星总体设计与关键技术[J]. 测绘学报,2017, 46(3): 269–277.Zhang Qingjun. System design and key technologies of the GF-3 satellite[J].Acta Geodaetica et Cartographica Sinica,2017, 46(3): 269–277. DOI: 10.11947/j.AGCS.2017.2017 0049.

[2]丁赤飚, 刘佳音, 雷斌, 等. 高分三号SAR卫星系统级几何定位精度初探[J]. 雷达学报, 2017, 6(1): 11–16.Ding Chibiao, Liu Jiayin, Lei Bin,et al.. Preliminary exploration of systematic geolocation accuracy of GF-3 SAR satellite system[J].Journal of Radars, 2017, 6(1): 11–16.

[3]杨劲松, 任林, 王隽. 高分三号卫星对海浪的首次定量遥感[J].海洋与湖沼, 2017, 48(2): 207–209.Yang Jingsong, Ren Lin, and Wang Juan. The first quantitative remote sensing of ocean surface waves by Chinese GF-3 SAR satellite[J].Oceanologia et Limnologia Sinica, 2017, 48(2): 207–209.

[4]王文煜, 谢春华, 袁新哲, 等. 高分三号卫星ATI模式海表面流场测量性能分析[J]. 航天器工程, 2017, 26(1): 132–139.Wang Wenyu, Xie Chunhua, Yuan Xinzhe,et al..Performance analysis of sea surface current measurements by GF-3 satellite ATI mode[J].Spacecraft Engineering,2017, 26(1): 132–139.

[5]杨劲松, 王隽, 任林. 高分三号卫星对海洋内波的首次定量遥感[J]. 海洋学报, 2017, 39(1): 148.Yang Jingsong, Wang Juan, and Ren Lin. The first quantitative remote sensing of ocean internal wave By Chinese GF-3 SAR satellite[J].Haiyang Xuebao, 2017,39(1): 148.

[6]Li Hengchao, Hong Wen, Wu Yirong,et al.. On the empirical-statistical modeling of SAR images with generalized gamma distribution[J].IEEE Journal of Selected Topics in Signal Processing, 2011, 5(3): 386–397. DOI:10.1109/JSTSP.2011.2138675.

[7]Ni Weiping, Yan Weidong, Wu Junzheng,et al.. Statistical analysis and modeling of TerraSAR-X images for CFAR based target detection[C]. Proceedings of 2013 IEEE International Geoscience and Remote Sensing Symposium,Melbourne, VIC, Australia, 2013: 1983–1986.

[8]高贵, 周蝶飞, 蒋咏梅, 等. SAR图像目标检测研究综述[J]. 信号处理, 2008, 24(6): 971–981.Gao Gui, Zhou Diefei, Jiang Yongmei,et al.. Study on target detection in SAR image: A survey[J].Signal Processing, 2008, 24(6): 971–981.

[9]艾加秋, 齐向阳. 一种基于局部K-分布的新的SAR图像舰船检测算法[J]. 中国科学院研究生院学报, 2010, 27(1): 36–42.Ai Jiaqiu and Qi Xiangyang. A new ship detection algorithm based on local K-distribution in SAR images[J].Journal of the Graduate School of the Chinese Academy of Sciences, 2010, 27(1): 36–42.

[10]Qin Xianxiang, Zhou Shilin, Zou Huanxin,et al.. A CFAR detection algorithm for generalized gamma distributed background in high-resolution SAR images[J].IEEE Geoscience and Remote Sensing Letters, 2013, 10(4):806–810. DOI: 10.1109/LGRS.2012.2224317.

[11]Gao Gui, Ouyang Kewei, Luo Yongbo,et al.. Scheme of parameter estimation for generalized gamma distribution and its application to ship detection in SAR images[J].IEEE Transactions on Geoscience and Remote Sensing,2017, 55(3): 1812–1832. DOI: 10.1109/TGRS.2016.2634862.

[12]艾加秋, 齐向阳, 禹卫东. 改进的SAR图像双参数CFAR舰船检测算法[J]. 电子与信息学报, 2009, 31(12): 2881–2885.Ai Jiaqiu, Qi Xiangyang, and Yu Weidong. Improved two parameter CFAR ship detection algorithm in SAR images[J].Journal of Electronicsamp;Information Technology,2009, 31(12): 2881–2885.

[13]种劲松, 朱敏慧. SAR图像局部窗口K-分布目标检测算法[J].电子与信息学报, 2003, 25(9): 1276–1280.Chong Jinsong and Zhu Minhui. Target detection algorithm of SAR image based on local window K-distribution[J].Journal of Electronics and Information Technology, 2003,25(9): 1276–1280.

[14]Akylmaz E, Demirkesen C, Nar F,et al.. Interactive ship segmentation in SAR images[C]. Proceedings of the 2013 21st Signal Processing and Communications Applications Conference (SIU), Haspolat, 2013: 1–4.

[15]Wang Shigang, Wang Min, Yang Shuyuan,et al.. New hierarchical saliency filtering for fast ship detection in highresolution SAR images[J].IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(1): 351–362. DOI:10.1109/TGRS.2016.2606481.

[16]Zhai Liang, Li Yu, and Su Yi. Inshore ship detection via saliency and context information in high-resolution SAR images[J].IEEE Geoscience and Remote Sensing Letters,2016, 13(12): 1870–1874. DOI: 10.1109/LGRS.2016.2616187.

[17]Zhao J, Guo W, Liu B,et al.. Ship detection for sentinel-1 SAR imagery based on a coupled convolutional neural network[J].IEEE Geoscience and Remote Sensing Letters,submitted for reviewing.

[18]Otsu N. A threshold selection method from gray-level histograms[J].IEEE Transactions on Systems,Man,and Cybernetics, 1979, 9(1): 62–66. DOI: 10.1109/TSMC.1979.4310076.

[19]Duda R O and Hart P E. Use of the Hough transformation to detect lines and curves in pictures[J].Communications of the ACM, 1972, 15(1): 11–15. DOI: 10.1145/361237.361242.

[20]Tupin F, Maitre H, Mangin J F,et al.. Detection of linear features in SAR images: Application to road network extraction[J].IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(2): 434–453. DOI: 10.1109/36.662728.

[21]Goodman J W. Some fundamental properties of speckle[J].Journal of the Optical Society of America, 1976, 66(11):1145–1150. DOI: 10.1364/JOSA.66.001145.

[22]Dempster A P, Laird N M, and Rubin D B. Maximum likelihood from incomplete data via the EM algorithm[J].Journal of the Royal Statistical Society. Series B(Methodological), 1977, 39(1): 1–38.

[23]Wu Yonghui, Ji Kefeng, Yu Wenxian,et al.. Region-based classification of polarimetric SAR images using wishart MRF[J].IEEE Geoscience and Remote Sensing Letters,2008, 5(4): 668–672. DOI: 10.1109/LGRS.2008.2002263.

[24]吴永辉. 极化SAR图像分类技术研究[D]. [博士论文], 国防科学技术大学, 2007: 117–121.Wu Yonghui. Classification of polarimetric SAR images[D].[Ph.D. dissertation], National University of Defense Technology, 2007: 117–121.

[25]曹峰, 邢相薇, 计科峰. 一种SAR图像舰船目标快速检测方法[J]. 雷达科学与技术, 2012, 10(4): 380–386.Cao Feng, Xing Xiangwei, and Ji Kefeng. A fast method for ship detection in SAR imagery[J].Radar Science and Technology, 2012, 10(4): 380–386.

刘泽宇(1993–)，男，黑龙江哈尔滨人，于2014年获上海交通大学学士学位，2014年9月至今在上海交通大学电子信息与电气工程学院攻读博士研究生。研究方向为雷达图像解译、数据挖掘。

E-mail: ribosomal@sjtu.edu.cn

柳 彬(1985–)，男，湖南衡阳人，博士，助理研究员，分别于2007年、2009年和2015年获上海交通大学信息工程、信号与信息处理和信号与信息处理学士、硕士和博士学位。2012年10月至2013年4月在法国巴黎高科电信学院访问研究。2015年12月，任上海交通大学电信学院信息技术与电气工程研究院助理研究员。主要从事雷达图像的分割分类、目标检测识别、多时相分析等方面的研究。

E-mail: bliu.rsti@sjtu.edu.cn

郭炜炜(1983–)，男，江苏南通人，博士，分别于2005年、2007年和2011年获国防科技大学信息与通信工程专业学士、硕士和博士学位。2014年至今，在上海交通大学电子信息与电气工程学院做博士后。主要从事图像理解、模式识别与机器学习等方面的研究。

E-mail: gwnudt@163.com

张增辉(1980–)，男，山东金乡人，博士，副研究员，分别于2001年、2003年和2008年获国防科技大学应用数学、计算数学和信息与通信工程专业学士、硕士和博士学位。2008年6月，任国防科大理学院数学与系统科学系讲师；2014年2月，任上海交通大学电子信息与电气工程学院副研究员。主要从事新体制雷达系统、雷达信号处理、压缩感知理论等方面的研究。

E-mail: zenghui.zhang@sjtu.edu.cn

张 波(1981–)，男，陕西西安人，硕士，工程师，分别于2004年、2009年获西安交通大学信息工程、系统工程学士和硕士学位。2009年8月任中国空间技术研究院西安分院助理工程师；2011年任中国空间技术研究院西安分院工程师。主要从事遥感卫星数据传输系统设计、遥感卫星地面接收处理系统设计、卫星应用技术等方面研究。

E-mail: seipopzb@163.com

周月恒(1990–)，女，陕西西安人，硕士，工程师，于2012年获武汉大学地理信息系统学士学位，2015年获纽约州立大学布法罗分校地理硕士学位。2015年至今在中国空间技术研究院西安分院工作。主要从事遥感数据处理、目标识别等方面研究。

E-mail: zyhzyh360@126.com

马 高(1989–)，男，山西长治人，硕士，工程师，2012年获陕西师范大学地图学与地理信息系统硕士学位。2015年至今在中国空间技术研究院西安分院工作。主要从事遥感数据处理、信息智能提取等方面研究。

E-mail: magao2002@163.com

郁文贤(1964–)，男，上海松江人，博士，教授，博士生导师，上海交通大学讲席教授。中国第2代卫星导航系统重大专项测试评估与试验验证专家组专家，高分辨率对地观测系统重大专项专家委员会地面系统组专家，“十二五”总装备部卫星应用技术专业组顾问，装发部上海市“北斗导航与位置服务”共建重点实验室主任，上海交通大学学术委员会委员，雷达信号处理国防科技重点实验室学术委员会委员，“十一五”国家863计划信息获取与处理技术主题第一、第二届专家组组长，“十一五”总装备部雷达探测技术专业组专家，主要研究方向为先进探测技术和多维信号与信息处理，研究内容包括新型成像系统、微波图像处理和解译、信息融合、目标识别等。

E-mail: wxyu@sjtu.edu.cn

Key Program of National Natural Science Foundation of China—High Resolution SAR Database and Data Quality Evaluation (61331015)

Ship Detection in GF-3 NSC Mode SAR Images

Liu Zeyu①Liu Bin①Guo Weiwei①Zhang Zenghui①Zhang Bo②Zhou Yueheng②Ma Gao②Yu Wenxian①
①(Shanghai Key Laboratory of Intelligent Sensing and Recognition,Shanghai Jiao Tong University,Shanghai200240,China)
②(China Academy of Space Technology-Xi’an,Xi’an710100,China)

GF-3, the first C-band full-polarimetric Synthetic Aperture Radar (SAR) satellite with a space resolution up to 1 m, has multiple strip and scan imaging modes. In this paper, we propose a maritime ship detection algorithm that detects ship targets via pixel classification in a Bayesian framework and employ effective enhancement methods to improve detection performance based on the data characteristics. We compare and analyze the results of detection experiments using the proposed algorithm with those of several Constant False Alarm Rate (CFAR) algorithms. The experimental results verify the effectiveness of the proposed algorithm.

GF-3 satellite; Synthetic Aperture Radar (SAR); Ship detection; Pixel classification

TP75

A

2095-283X(2017)05-0473-10

10.12000/JR17059

刘泽宇, 柳彬, 郭炜炜, 等. 高分三号NSC模式SAR图像舰船目标检测初探[J]. 雷达学报, 2017, 6(5):473–482.

10.12000/JR17059.

Reference format:Liu Zeyu, Liu Bin, Guo Weiwei,et al.. Ship detection in GF-3 NSC mode SAR images[J].Journal of Radars, 2017, 6(5): 473–482. DOI: 10.12000/JR17059.

2017-06-15；改回日期：2017-07-31；网络出版：2017-08-18

*通信作者： 柳彬 bliu.rsti@sjtu.edu.cn

国家自然科学基金重点项目—高分辨率SAR测试库及数据质量评估(61331015)