基于反射对称性的极化SAR舰船目标CFAR检测

　摘 要： 反射对称性是极化散射中十分重要的一种性质。反射对称性检测量刻画了共极化通道和交叉极化通道间的相关性，在极化合成孔径雷达（PolSAR）目标检测中有重要应用。在复Wishart分布的基础上，结合乘积模型，推导出反射对称性检测量统计模型服从分布，利用实测数据证实该模型对杂波数据具有很好的拟合效果。最后基于分布模型，用恒虚警率（CFAR）方法进行了舰船目标的检测，实验结果验证了该方法的有效性。
　　关键词： 反射对称性； 舰船目标检测； 恒虚警率； 极化合成孔径雷达
　　中图分类号： TN957.52?34 文献标识码： A 文章编号： 1004?373X（2015）09?0019?06
　　Abstract： Reflection symmetry is a very importance property in polarization scattering. The correlation between the co?pola?rization channel and cross?polarization channel is described by reflection symmetry detection， it has important application in polarization synthetic aperture radar （PolSAR） target detection. Based on complex Wishart distribution and combined with product model， it is inferred that the statistic model of the reflection symmetry detection quantity obeys K?distribution. The measured data verified that good fitting effect of the clutter data is achieved by using the proposed model. Ship target detection is implemen?ted by the method of constant false alarm rate （CFAR） based on K?distribution model， and the experiment results verify the efficiency of the proposed detection.
　　Keywords： reflection symmetry； ship target detection； CFAR； PolSAR
　　0 引 言
　　合成孔径雷达（Synthetic Aperture Radar，SAR）具有独特的大范围、全天候、全天时的对地监视能力。利用SAR进行舰船目标检测，在海洋监视和海上救助方面具有十分关键的作用。随着极化SAR（Polarization Synthetic Aperture Radar，PolSAR）系统的快速发展，利用极化数据中丰富的极化信息进行舰船目标检测为SAR舰船目标检测开辟了新思路。
　　人造舰船和自然海洋表面具有不同的极化散射机理，基于这些极化散射差异有可能在极化SAR图像上把舰船目标从海洋中分离出来。Cameron于1996年提出了相干目标分解（Cameron Target Decomposition，CTD）理论。1999年，Ringrose利用SIR?C 的全极化数据，首次将Cameron分解方法应用于舰船目标检测。2000年，加拿大的Touzi等利用极化熵方法对Convair?580的极化数据在不同入射角下进行了舰船检测研究，实验结果证明入射角小于60°时，极化熵可增强目标与背景的对比度，但是不同的海况对该方法的影响较大。2001年，意大利学者Sciotti针对SIR?C全极化SAR数据提出了中心化法和非中心化法舰船检测方法，该类方法检测效果较好，但运算时间较长。
　　1981年Van de首次提出了散射体及与其所对称的响应散射体的极化散射矩阵形式。1992年首次研究了极化SAR遥感的反射对称性问题，2012年Nunziata与Migliaccio基于极化SAR数据中金属目标和海洋的反射对称性差异，首次利用极化协方差矩阵中的共极化（HH，VV）和交叉极化通道（HV）之间的相关性构建了一种新的检测量——反射对称检测量，用于观测海洋中人造金属目标，该方法在理论上表现出较好的检测性能，利用反射对称性构造反射对称检测量进行舰船目标检测已经成为近几年极化SAR舰船目标检测的研究热点。然而在已有文献中该检测量在应用于目标检测时都是根据经验选取检测阈值，当SAR成像系统、海洋背景等条件发生变化时经验阈值往往发生变化，不具有很好的适应性，造成了基于反射对称检测量的检测方法在实际应用中困难重重，迫切需要研究相应的自适应检测方法。这其中，恒虚警率（Constant False Alarm Rate，CFAR）检测是雷达目标检测领域的经典和常用方法，具有广泛的应用前景，其首要和关键问题是研究相应检测量的杂波统计模型。但遗憾的是目前还没有文献给出基于极化协方差矩阵的反射对称性检测量的杂波统计模型，这也制约了CFAR检测在基于反射对称性舰船目标检测中的应用。
　　本文基于多变量复高斯统计分布理论，从复Wishart分布出发，结合乘积模型推导了反射对称检测量的统计分布，并采用对数累计量方法给出了参数估计方法，实测数据实验表明推导的理论分布能较好地拟合检测量分布直方图。最后采用CFAR检测，在不同的数据源、不同虚警率下对反射对称性检测和极化交叉熵检测的性能进行了对比，实验证明反射对称性检测量具有更优良的检测性能。
　　1 基于极化协方差矩阵的反射对称性检测量
　　在实际的遥感应用中，不存在完全的确定性散射体，因此为了处理统计散射效应和分析局部散射体，通常引入极化协方差矩阵，它由lexicographic基散射矢量导出，描述了散射体的不确定性散射过程。利用lexicographic基表示的散射矢量直接包含了散射矩阵的复幅度信息，与系统的测量值直接对应，能够很好地解释物理散射机制。海面舰船目标由于其金属属性与海杂波之间的幅度对比度十分强，因此对海面目标进行极化检测可以强调利用幅度信息。因此，文献在进行极化舰船检测时均是从极化协方差矩阵出发构造反射对称性检测量。
　　表示共极化通道与交叉极化通道之间相关性的模。海平面满足反射对称性，所以趋于0。舰船目标不满足对称反射性，所以值较自由海面情况下会明显大于0。理论上，反射对称检测量同样可以构造为，但是研究表明，舰船目标在HH和HV极化通道中较为明显，而尾迹等海洋特征则在VV通道中较为明显，因此文中只对的统计分布进行研究。
　　2 基于极化协方差矩阵的反射对称性检测量的
　　统计建模
　　3 参数估计
　　K分布的参数估计方法常见的有最大似然估计（Maximum Likelihood，ML）、矩估计。ML估计方法是最优估计方法，但是通过对数似然函数很难得到形状参数和尺度参数。广泛采用的是基于矩估计的K分布参数估计法，其特点是形式简单、估计方便，但是该方法却存在诸多不足，突出表现为要单独估计有效视数而 有效视数的精度极大地影响了后续形状参数和尺度参数的估计精度，甚至使K分布失效。
　　文献给出了一种对数累计量参数估计方法，该方法把有效视数同形状参数尺度参数一样视为待估计参数。而且基于对数累计量的参数估计方法，相当于对SAR图像数据做了对数变换后进行参数估计，有效缩小了低阶样本对数累计量的变化范围，也就是说不论杂波均匀度如何，都不会出现K分布失效的问题。因此本文也采用对数累计量方法进行K分布参数估计。
　　4 基于反射对称性的舰船目标CFAR检测流程
　　检测流程主要分为三步：
　　（1） 要从极化协方差矩阵中水平极化数据垂直极化数据和水平极化与交叉极化的相关性数据构造反射对称性归一化检测量；
　　（2） 估计分布参数，建立精准的理论统计模型；
　　（3） 执行CFAR检测输出二值化检测结果图。
　　5 实验结果及分析
　　实验采用同一SAR平台，不同地区、不同分辨率、不同波段、不同海况下的极化数据充分验证K分布对反射对称检测量直方图的良好拟合性能以及反射对称检测方法的有效性。数据来源为美国NASA的AIR?SAR极化数据，日本东京湾极化数据和我国台湾高雄地区极化数据，数据示意图（Pauli基合成图）如图2所示。东京湾海域C波段水平极化（HH）数据如图2（a）所示，距离向和方位向理论分辨率为5 m×3 m，高雄港L波段水平极化数据如图2（b）所示，距离向和方位向理论分辨率为2.8 m×5 m。
　　5.1 检测量统计分布拟合
　　在水平极化（HH）、垂直极化（VV）和交叉极化（HV）通道数据中的相同位置割取一块海杂波切片，如图2中1号框标定区域所示，验证K分布对反射对称检测量直方图的拟合性能。首先按照式（8）对检测量进行归一化处理，并统计归一化后的检测量的直方图，然后根据文中第3节所提参数估计方法，对参数进行估计，拟合结果如图3所示。
　　另外，在高海况下，极化交叉熵检测结果中虚警显著增多，整体检测性能下降明显；而反射对称性检测方法检测性能无明显变化，也说明了该方法对海况的变化有较好的适应性。分析其原因，极化交叉熵方法是假定目标和杂波有着不同的散射机理，反映了相邻分辨单元极化散射的随机性，因此目标和杂波之间的极化熵值可区分，但是如果杂波呈现非布拉格散射特性（高海况时），那么这个假设并不是完全成立的，这也正是为什么海况复杂后极化交叉熵检测结果明显下降的原因。相反，反射对称性方法和杂波、目标的后向散射机理无明显关系，主要反映了人造反射物体与自然反射物体在散射对称性上的区别。结果对比如表1所示。
　　6 结 语
　　本文基于极化SAR数据反射对称性原理，在复Wishart分布的基础上，结合乘积模型，推导出反射对称检测量服从K分布模型，并利用不同地区、不同波段的实测数据进行了实验验证。实验结果表明K分布能够很好地对实测区域的反射对称检测量进行拟合。同时，基于K分布的CFAR舰船检测结果表明，该检测方法相比文献提出的极化交叉熵检测方法能够更好地实现对舰船目标的检测，且具有虚警率低的优点。
　　参考文献
　　 CAMERON W L， YOUSSEF N N， LEUNG L L. Simulated polarimetric signatures of primitive geometrical shapes . IEEE Transactions on Geoscience and Remote Sensing， 1996， 34（3）： 793?803.
　　 RINGROSE R， HARRIS N. Ship detection using polarimetric SAR data // Proceedings of Committee on Earth Observation Satellites. Toulouse： SAR Workshop， 1999， 450： 687?691.
　　 TOUZI R， CHARBONNEAU F， HAWKINS R K， et al. Ship?sea contrast optimization when using polarimetric SARs // IEEE 2001 International Conference on Geoscience and Remote Sensing Symposium. Sydney： IEEE， 2001， 1： 426?428.

本文链接：http://www.qk112.com/lwfw/jsjlw/zhinengkeji/242264.html