SAR图像舰船目标去旁瓣处理

0 引言

合成孔径雷达(SAR)具有分辨率高、探测范围广、全天时、全天候的优点，在舰船目标监视、海洋资源保护中起到重要作用。利用SAR获取海上舰船目标的电磁散射特征、结构特征、位置及运动状态，实现舰船目标分类识别，是各海洋大国的研究热点。

舰船目标一般由金属部件构成，当舰船某一结构的后向散射特别强烈或图像聚焦效果较差时，相干叠加形成明显的十字旁瓣效应[1]。在基于SAR图像的舰船检测识别系统中，舰船的长度、宽度、姿态角、面积等参数被作为分类识别的重要依据，而旁瓣的存在严重影响了上述参数的精度，进而影响舰船分类识别准确性。

有不少文献中都提出了去旁瓣的方法，如邢相薇[1]利用Radon变换的特点将旁瓣和舰船目标主轴进行区分，该方法认为当主轴和旁瓣具有不同的方向时，峰值点在Radon变换域中处于不同的位置，可在Radon变换域中对其进行区分；当旁瓣方向与主轴方向一致，或者旁瓣污染严重、存在多个旁瓣时，该算法将无法确定哪些峰值点是旁瓣。蒋少峰等[2]利用沿方位向和距离向积分寻找旁瓣，认为旁瓣宽度一般为1～2个像素，采用局部均值滤波消除其影响，但是，由于舰船内部结构在SAR图像上引起的强弱变化，在方位向和距离向上的积分受舰船本身回波影响，并不能真实反映出旁瓣位置，而且在处理的图像中，有的旁瓣宽度达十几个像素，局部均值滤波区域如何确定是另一个要考虑的问题。甄勇等[3]利用舰船方位角，确定主轴方向，计算目标的外接矩形，去除旁瓣影响，该方法同样对聚焦效果好的图像有效果；当旁瓣对目标整体结构影响较大，主轴无法确定时，该方法将同样不能去除旁瓣影响。

本文基于未经量化的SAR模图像，结合旁瓣由强散射点产生并且只会在方位向和距离向出现的特点，提出一种舰船目标去旁瓣方法，该方法可以去除旁瓣影响，同时最大限度保留船只弱散射区域。

1 旁瓣特征分析

旁瓣能量相比主瓣虽然大大减弱，但仍然高于其他弱散射部位回波信号，如图1所示。图1(a)为舰船SAR图像，图1(b)为旁瓣能量分布示意图。因此，基于阈值删除旁瓣将导致船只内部弱散射区域信息丢失。旁瓣是由舰船强散射结构产生，船只自身的弱散射区域不会引起旁瓣效应，且其形状一般是沿距离向和方位向垂直的“十”字形。通过统计发现，旁瓣能量一般低于舰船目标整体区域能量的均值。

  

图1 旁瓣对舰船成像区域影响

2 旁瓣删除流程

本文提出的旁瓣删除算法流程如图2所示，考虑到量化后的SAR图像丢失过多回波信息，使用未经量化的模图像。

  

图2 旁瓣删除流程图

具体处理步骤如下：

步骤1：对模图像，利用CFAR实现目标区域检测[4-5]，检测结果为二值化图像(定义为B)，如图3所示，其中目标区域置1，其余区域置0。

T=mean(x)+t×var(x)

  

图3 SAR图像CFAR检测结果

步骤2：对二值图像B，用形态学滤波去除孤立小区域。利用邻域搜索确定目标所在区域，考虑到目标可能由数个孤立区域组成，邻域搜索范围应有一定扩大，比如在10个像素范围内搜索目标点。对目标所在区域计算最小外接矩形M1，确定M1的4个角点坐标[6]，二值图像B在M1区域内部分定义为D，如图4中方框所示。

  

图4 区域最小外接矩形示意图

农业废弃物资源化利用是省农业现代化的重要考核目标，也是农业生态文明建设的重要内容。各地要认真总结前几年的先进经验和做法，围绕中央、省市政府下达的农业废弃物资源化利用工作目标和任务，继续把农业废弃物资源化利用列入农业农村重点考核目标作。同时结合本地实际，制定明确年度工作任务和实施方案、督查方案和考核办法，做到督查能发现问题和不足，考核能拉开分数和等次，进一步加强这项工作的督查考核力度，保证全市农业废弃物资源化利用的有效推进。

式中，x为目标像素点模值，mean(x)为D中非零区域能量均值，var(x)为D中非零区域能量方差，t为阈值调节系数，T为检测门限。M1内的模数据，大于阈值T的像素点为强散射点，置1，其余置0，强散射区域二值化结果定义为S，如图5所示。

步骤3：读入M1内的模数据，统计目标能量均值mean(x)和方差var(x)。

图7(b)中墩柱均有耗能钢筋的参与，其中墩柱③，耗能钢筋部分裸露在外，墩柱的耗能能力以及抗侧力较墩柱①有明显的提升.滞回曲线饱满，耗能能力优异，滞回曲线中的各个滞回环界限分明，相邻滞回环的后一等级覆盖上一等级的滞回环.当加载位移达到100 mm时，局部节段墩柱间局部混凝土出现应力集中，应力迅速增加，墩柱的抗侧力稍有提升.墩柱④有现浇混凝土包裹耗能钢筋，耗能钢筋与现浇混凝土共同作用，使墩柱的抗侧力大幅提升，耗能能力明显提高.

(1)突出环境照明。所谓的突出环境照明，又可以称之为基础照明，是指整个展示现场的空间照明。整个环境明亮，展示空间和展品的细节都能清晰的展现。展区的面貌和气氛可以由灯光改变，不同的灯光营造不同的展区效果，柔和灯光让人舒适，延缓视觉疲劳；强烈灯光突出展品，让参展这印象深刻。突出环境照明的基础照明可用直接光源，也可以用间接光源，需根据展示的需要和类型而定。

(1)

步骤4：设定阈值，用于确定目标强散射区域。

 

(2)

  

图5 目标强散射区域

2)统计矩阵D在L区域内，第i列非0像素点个数CDi，则有数组CD=[CD1,CD2,…,CDw]。

实现社会治理的有效、有序进行，基层警务工作有着不可替代的作用，充分、合法履行职能是其发挥作用的前提。在现有法律、法规基础上，推动社会治理重心向基层下移，对基层警务工作进行职能优化，不断满足人民日益增长的美好生活需要，不断促进社会公平正义，形成有效的社会治理、良好的社会秩序。

革命爆发之后，德伐日太太由昔日的被压迫者转为新的压迫者，她可以随意的将贵族送上断头台，为了报仇她要“斩草除根”。“复仇女神”和“雅克”们是她的拥护者，替她实现她的复仇计划。她逐渐变为没有人性，惨无人道的统治者，像昔日的统治阶级一样残酷无情。新的压迫者代替了旧的压迫者，而新的压迫者也必将灭亡，德伐日太太在成为新的统治者时，就注定了灭亡的结局。历史的进步无人能阻挡，新统治阶级的出现无法阻挡，德伐日太太的最终失败无法改变。这又是其悲剧性的一方面。

  

图6 强散射区域最小外接矩形

步骤6：假设目标区域宽度为w，高度为h，方位向为图像行向，距离向为图像列向，定义区域L=M1-M2，如图7所示。在L，M1和M2内，进行如下信息统计。

  

图7 区域示意图

1)L区域内：统计区域L在第i列像素点个数，定义为CLi，i=1,…,w，则有数组CL=[CL1,CL2,…,CLw]。

3)统计矩阵S在M2区域内，第i列非0像素点个数CSi，则有数组CS=[CS1,CS2,…,CSw]。

步骤5：对S，重新计算最小外接矩形M2，如图6中方框所示。

在当前的互联网时代，积极推进教育信息化在护理专业中的应用和推广，可以有效推进护理专业的内涵式发展，提高高职院校护理专业教学的信息化以及护理专业毕业生理论能力、信息操作能力的提升，更好地满足我国当前医院临床护理岗位的建设和发展需求。

4)由于旁瓣由强散射点产生，故只有CSi>0时，列向存在旁瓣。相比舰船自身的弱散射区域，旁瓣在距离向空间上有更长延伸，利用弱散射点在距离向所占比例作为旁瓣判断依据，记为数组PC，PC=[Pc1,Pc2,…,Pcw]，当CSi>0时，时，PCi=0，i=1,…,w。

5)对数组PC中PCi>0的值，统计均值，记为Avgc，如果PCi>Avgc，则认为第i列存在旁瓣，在L区域内，将矩阵D第i列的非0值置0。

不同时间点10 μmol/L GSK126对细胞增殖的影响见图1b。结果显示在3个不同时间点，正常对照组细胞活力明显高于激素性股骨头坏死组（P＜0.001），用10 μmol/L GSK126干预后，细胞活力较未干预前明显提高（P＜0.01），但仍低于对照组。表明10 μmol/L GSK126能明显促进激素性股骨头坏死MSCs的增殖。

6)同样方法，统计方位向(行向)L区域像素点个数RLj，D在L区域内像素点个数RDj，S在M2区域内像素点个数RSj，当RSj>0时，时，PRj=0，j=1,…,h，对数组PR中PRj>0的值，统计均值，记为AvgR，PRj>AvgR，则认为第j行存在旁瓣，将矩阵D在L区域中，第j行的非0值置0。经过此步处理后，目标旁瓣删除结果如图8所示。

  

图8 旁瓣删除示意图

步骤7：对旁瓣删除后的区域，用形态学滤波去除孤立小区域，剩余区域重新计算外接矩形M3，从模图像中取出M3区域内的数据，用于进一步目标识别处理，如图9所示。

  

图9 旁瓣删除处理后的目标成像区域

3 实验及结果分析

为验证本算法的适用性，使用20幅机载SAR图像进行测试，图像分辨率1～3 m，每幅图像中包含40～60艘船只，本文的算法均较好删除了目标旁瓣，更精确地获取目标真实成像区域，为舰船目标进一步分类识别奠定扎实基础。图10为旁瓣删除前后获取的目标成像区域比较。使用同样的算法，对经过旁瓣删除处理和未经旁瓣删除处理的图像进行目标几何参数提取，结果如图11所示。图11(a)是对未经旁瓣删除处理的图像进行目标几何参数提取的结果，图11(b)是经过旁瓣删除处理后舰船几何参数提取结果。可以看出，经过本文算法处理获取的几何参数更接近于目标真实尺寸，在旁瓣删除基础上进行目标检测识别，准确率将有大幅度提高。

  

图10 旁瓣删除前后目标成像区域比较

  

图11 旁瓣删除对目标几何特征提取的改善

4 结束语

本文分析SAR图像舰船旁瓣成像机理，基于旁瓣能量分布特征，提出了一种SAR图像舰船目标去旁瓣算法，基于未经量化的SAR模图像，统计成像区域能量分布，以及旁瓣与强散射点的位置关系，进而确定旁瓣区域。对比实验证明，该算

法可以在剔除旁瓣对目标影响的同时，尽可能保留目标弱散射区域，精确分割出目标成像区域，为准确提取目标几何参数、提高识别精度提供有力支撑。

参考文献：

[1] 邢相薇. HRWS SAR图像舰船目标监视关键技术研究[D]. 长沙: 国防科学技术大学, 2014.

[2] 蒋少峰, 王超, 吴樊, 等. 基于结构特征分析的COSMO-SkyMed图像商用船舶分类算法[J]. 遥感技术与应用, 2014, 29(4):607-614.

[3] 甄勇, 刘伟, 陈建宏, 等. 高分辨率SAR图像舰船目标几何结构特征提取[J].信号处理, 2016, 32(4):424-429.

[4] CERUTI-MAORI D, SIKANETA I, GIERULL C. Detection and Imaging of Moving Objects with High-Resolution Wide-Swath SAR Systems[C]∥ 10th European Conference on Synthetic Aperture Radar, Berlin: VDE, 2014:977-980.

[5] 龚小冬, 李飞, 张志敏, 等. 一种新颖的高分辨率SAR舰船检测方法[J]. 雷达科学与技术, 2015, 13(1):70-75.

[6] 殷雄, 王超, 张红, 等. 基于结构特征的高分辨率TerraSAR-X图像船舶识别方法研究[J]. 中国图象图形学报, 2012, 17(1):106-113.