基于表征转换机的SAR图像目标分割方法

赵晓辉 姜义成朱同宇

（哈尔滨工业大学 哈尔滨 150001）

基于表征转换机的SAR图像目标分割方法

赵晓辉 姜义成*朱同宇

（哈尔滨工业大学 哈尔滨 150001）

针对SAR（Synthetic Aperture Radar）图像中的目标分割问题，由于目标与杂波空间模式（像素强度和分布）不同，通过分析图像空间模式的方式可达到分辨目标和杂波并分割目标的目的。该文基于表征转换机理论提出一种有效的SAR图像目标分割方法，该算法分析SAR图像中的空间模式，计算其与参考杂波图像的相似程度，最后将与参考杂波相似程度较高的部分消除以达到分割目标的目的，并在衡量相似度部分使用基于累积直方图的自动阈值选取办法。仿真和实测数据的实验验证了此算法的有效性。

合成孔径雷达；目标分割；表征转换机；极限学习机

引用格式：赵晓辉， 姜义成， 朱同宇.基于表征转换机的SAR图像目标分割方法［J］.雷达学报， 2016， 5（4）: 402-409.

DOI: 10.12000/JR16066.

Reference format: Zhao Xiaohui， Jiang Yicheng， and Zhu Tongyu.Target segmentation method in SAR images based on appearance conversion machine［J］.Journal of Radars， 2016， 5（4）: 402-409.DOI: 10.12000/JR16066.

1 引言

合成孔径雷达（Synthetic Aperture Radar，SAR）具有全天时、全天候执行遥感监视任务的能力，其可以通过相干积累目标的多次回波得到目标的2维SAR图像［1，2］。一幅SAR图像通常由感兴趣目标和杂波组成。在实际环境中，由于较强的海杂波或地杂波的影响，使得感兴趣目标与杂波背景间的分界更不明显，这种弱分界增加了目标分割的难度。精确的目标分割结果有利于提高后续图像处理的有效性，如根据分割得到的目标轮廓可更精确地估计目标所处的姿态，进而提高姿态矫正的有效性。

图像分割算法主要包括两类：基于局部的分割算法和基于全局的分割算法。传统分割算法的通过计算阈值的方式，根据阈值将给定图像分割为前景和背景两部分，即大于阈值的一部分和小于阈值的另一部分。通过选取合适的阈值可得到感兴趣目标（前景图像）［3］。更具体地讲，对于SAR图像，前景分割的主要目的是选择能够分割感兴趣目标像素点和背景杂波像素点的合适阈值。

基于全局的方法相对容易实现且计算复杂度较低。文献［4］通过遍历的方法，寻找使图像前景和背景两部分方差最大的阈值以达到分割的目的，但需要前景和背景有较明显的不同。文献［5］提出了一种基于全局的简便二值化分割方法，其直接选取累积直方图二阶导数的第1个拐点作为分割阈值，但只能应用于背景像素值较低的场景。不难发现，基于全局的方法对待分割图像有较严格的要求，即待分割图像的背景应相对简单，且前景和背景间有明显对比。因此，学者们的研究重心逐渐由基于全局的分割方法转至基于局部的分割方法。

与基于全局的方法相比，基于局部的分割方法更能有效解决困难场景下的分割问题。如被广泛使用的恒虚警率（Const False Alarm Rate， CFAR）算法已被应用于解决SAR图像的前景分割问题［6-9］。更详细地讲，CFAR的本质是保证即使测试单元（Cell Under Test， CUT）内杂波/噪声水平有所改变时仍保持恒定的虚警率，这时的分割阈值（Threshold Value， TV）是通过分析CUT的近邻单元得到的。目前，研究CFAR的学者已提出较多改进版本的CFAR算法。如单元平均最大CFAR（Cell Averaging Greatest of CFAR）算法通过求单元平均值的方式计算TV，然后将CUT内大于TV的像素点归类为感兴趣目标的像素点。有序统计CFAR（Ordered Statistic CFAR）算法选择近邻单元内一个相对大的值作为TV，然后只保留CUT内大于TV的像素。然而，基于CFAR的算法都需要对近邻单元大小、乘积因子和累加因子等多个参数进行设置，并针对不同的环境仔细调节。

除选择TV进行图像分割的传统方法外，基于马尔科夫随机场（Markov Random Field， MRF）的SAR图像分割算法也被提出。基于MRF的分割算法为图像像素构建联合概率分布模型，然后计算像素点属于某一特定类别的概率，并根据概率完成图像分割［10，11］。文献［10］提出最小化Gibbs能量函数以避免MRF中计算密集的全局优化问题。文献［11］提出先使用k均值算法获得初始分割，然后使用MRF算法优化分割结果。尽管目前已有较多基于MRF方法的研究，但基于MRF算法的分割结果仍在很大程度上受初始分割结果及潜在分割类数的影响。

另一类在SAR图像分割领域被广泛研究的算法是基于聚类的算法［12，13］。此类算法能在特定条件下获得较好的性能，例如k均值算法针对球形簇才能获得最好的分割效果，谱聚类方法需要计算关联矩阵的特征向量来衡量成对像素间的相似度，而这需要较高的计算效能才能足够快地获得计算结果。

为实现更精确的SAR图像目标分割，该文认为目标与杂波在图像中的空间模式（像素能量强度和空间分布规律）并不相同且是相互可区分的，并基于表征转换机（Appearance Conversion Machine，ACM）理论提出一种有效的SAR图像目标分割算法，该算法的本质是分析给定图像中的像素空间模式，然后根据图像中像素空间模式的分布规律相应地改变像素值反映其与参考杂波图像的相似程度，最后根据衡量结果构建蒙版图像并完成分割。

本文组织结构如下：第2节介绍ACM理论；第3节提出本文的目标分割算法，并介绍算法的实现方法和原理；第4节给出本文算法对仿真和真实数据的实验处理结果；第5节总结全文。

2 表征转换机原理

对同一目标在不同采集条件下获得的不同表征图像，在没有获得足够信息的情况下完美地获取另外一种采集条件下的表征图像是非常困难的。即基于目标的任意单幅表征图像，难以提供足够的信息来获取目标在另一种采集条件下的相应表征图像。然而，若允许存在一定的误差，则可近似获取另一采集条件的对应表征图像［14，15］。可将式（1）改写为：

3 基于表征转换机的目标分割方法

在一幅SAR图像中，属于不同目标的像素点通常呈现明显不同的空间模式，如能量强度和像素分布方式等。因此，通过分析图像空间模式的方式可辨别像素的归属，进而完成目标分割。本文提出基于ACM理论的目标分割方法，该算法首先衡量目标图像与参考杂波图像的相似程度，然后基于相似程度构建蒙版图像并完成分割。

本文算法的基本原理如图1所示，图中的每个方块指代矩阵的一个元素，方块内的数字为相应元素的值。矩阵左上角的3×3区域与右下角区域具有明显不同的空间模式。若对此矩阵进行一定处理后，可不改变左上角3×3区域像素的值并令其余区域的元素值有较大提升，通过对比中间结果与原始图像各元素相对改变的幅度即可识别出各元素的归属，进而达到图像分割的目的。实际上，图1中的原始矩阵和中间结果分别为ACM的两个可相互转化的表征图像，基于ACM理论可构建一组回归函数，它们根据目标图像与参考杂波图像的相似程度改变目标图像的像素强度并构成一幅相似程度图。相似程度图内像素值的大小代表了相应区域与杂波空间模式的相似程度，基于相似程度图擦除与杂波相似程度较高的区域后即可获得分割后的目标图像。

3.1 增强非相干模式

本节讨论如何获得图1的中间结果（相似程度图）以衡量SAR图像与参考杂波图像的相似程度。给定一幅杂波模式为Pb，目标模式为Pt的SAR图像（M和N分别为方位维和距离维分辨单元的个数）。则可以定义由vec（I）到vec（I）的映射函数：

图1 本文提出算法的工作原理Fig.1 Basic concept of the proposed method

其中，vec（I）代表向量化的I，即将I中各行数据首尾相连获得单个行向量。

类似地，对只含杂波背景且杂波模式为Pb的 SAR图像，可定义一个由vec（B）到vec（B）的映射函数：

式中，E为全1矩阵。

结合式（5）和式（7）后可构建一个函数f （·），输入向量化的杂波图像vec（B）时，输出为向量化的全1矩阵vec（E），输入向量化的目标图像vec（I）时输出为向量化的矩阵vec（I）。则有映射函数f （·）可表示为：

式中，t（xi， yi）为I中子图像中的第i个像素点；Ns=（2d+1）2为子图像中的像素个数；d为尺度因子。

令q（m， n）为ROI（m， n）的因变像素点，即输入为ROI（m， n）时的相应输出。针对这组输入和输出可构建一组回归函数smn为：

式中，qi（m， n）为围绕Qi第m行第n列的像素，ROIi（m， n）为围绕Pi第m行第n列像素的子图像，i为输入输出图像的序号。

至此，式（8）求解的单一回归函数问题已经转化为求解一组回归函数的问题：

3.2 构建蒙版图像

将原图像I输入计算得到的回归函数后，可计算得到新图像I'，新图像如图2（b）所示。不难发现，处理后的图像中杂波背景的像素值被显著提高了，而待分割目标区域像素的值却没有显著改变。PIP处理后的图像I'减去原始图像I后可得到像素值的变化幅度图（相似程度图），如图2（c）所示。分析变化幅度图可知各区域与杂波背景的相似程度，基于此可以生成一幅蒙版图像，将此蒙版图像与原图像按位与操作即可完成目标分割。

图2（c）的变化幅度图反映了图像的各个区域与杂波的相似程度。不难理解，通过统计变化幅度图分布在各个灰度区间的像素个数，可以计算得到存在大部分杂波像素的最小相似阈值。实际上，这可以通过统计变化幅度图的灰度累积直方图来实现。灰度累积直方图反映了灰度区间相对前一个灰度区间的像素个数的改变程度，因此像素数量急剧增加的灰度区间（即灰度直方图曲率最大的点）将被选做分割阈值。

为获得灰度累积直方图，首先将图像的像素值按灰度区间分为R份，则对于第r个区间有式（12）成立

式中，n为图像的像素总数，nr为第r个灰度区间内的像素个数。在直角坐标系中做出的r与pr的关系图形即为灰度累积直方图。

图3是基于图2（c）计算获得的累积直方图，这100份灰度区间内第76个灰度区间的曲率最大。

为满足更精细的分割需求，也可用基于微调参数的方式使用具体公式计算。此时，计算蒙版图像的数学表达式为：

3.3 完整实现

本文基于ELM算法求解算法中的回归问题，这是因为ELM能以较低的计算复杂度获得有效的图像特征，这主要归因于ELM使用随机核和更少的优化参数，因此更易实现且拥有更强的扩展性［16］。

根据ELM学习理论［17-19］，正则化后的图像像素值位于欧几里德空间的一个有界可测紧凑子集上，对于任意给定的连续目标函数s（x）都存在较多映射函数和相应的输出权重使

式中，L为ELM特征空间的维数，有大部分非线性分段连续函数都可用于特征映射［18］

图2 计算变化幅度图Fig.2 The calculation of the variation scale map

图3 累积直方图和分割阈值Fig.3 Cumulative histogram and segmentation threshold

总地来讲，对任意给定SAR图像，可相应获得一组参考图像。为保证空间模式的有效分析，参考区域的最小尺寸应大于式（11）中的ROI尺寸。由于杂波背景以相似的分布模式遍布几乎整幅图像，即处于不同空间位置的杂波也可具有相似的分布模式。因此，参考背景区域空间位置的选择对算法的影响较小，对目标图像中任意只含有杂波背景的区域，通过平移和复制的方式即可生成与目标图像等尺寸的多幅参考杂波图像。基于这组图像数据，可相应构建一组回归函数计算给定目标SAR图像与杂波背景的相似程度，进而构建掩模图像并完成目标分割。本文算法的流程图如图4所示，具体实现步骤如下：

图4 本文算法的流程图Fig.4 Flowchart of the proposed method

输入：ELM的成本参数C、核参数P和数据集

其中，Ns为中的变量个数（图像像素数）；NC为输入的图像数（参考杂波图像的数量加1张原SAR图像）。

步骤1 为M×N个回归函数分别计算隐藏层输出矩阵H

步骤3 将原始图像输入得到的函数，计算变化幅度图；

步骤4 根据变化幅度图选择分割阈值，构建蒙版图像完成图像分割。

4 实验结果

本节测试算法的性能，并将本文算法与多个具有代表性的方法对比，如k-means算法和马尔科夫随机场算法［10］。实验环节使用多种实验数据测试算法在不同场景下的目标分割效果，实验数据包括MSTAR公开数据集中的高分辨SAR图像，MSTAR数据集数据的分辨率为0.3 m×0.3 m，均由X波段的聚束式SAR雷达采集得到。此外，还基于OKTAL获得目标在不同SCR环境下的仿真SAR图像。OKTAL基于射线追踪技术针对目标的CAD模型计算生成相应目标的仿真SAR图像，仿真参数载频35 GHz，斜距23.6 km，生成SAR图像的分辨率为0.3 m×0.3 m。

在算法衡量方面使用模版图像来衡量本文算法对目标细节的分割精度，感兴趣目标和杂波区域均基于原始图像人工标记获得［20］。在衡量参数方面，被正确分割的前景像素数（正确率）、被误判为背景的前景像素数（虚警率）和被正确分割的总像素数（精度）分别被定义为Ptp， Pfa和Pc。在实验参数方面，基于MRF模型方法的潜在分割类别数为2，期望值为0.4，最大循环数为50。本文涉及的所有算法均基于Matlab R2013a编程实现，测试电脑主机的内存为4 GB， CPU主频为1.8 GHz。

4.1 仿真SAR图像实验对比

本节使用OKTAL生成不同SCR环境下的仿真SAR图像，以测试本文算法在各种SCR环境下的目标分割性能。测试时共使用了4种SCR环境下的仿真SAR图像，它们的SCR分别为-2.3861 dB，-1.6118 dB，-0.7058 dB和0.3323 dB并分别用图A，图B，图C和图D表示，如图5所示。处理结果如表1所示，分析对比后不难发现，即使在原SAR图像SCR较低，目标与杂波间分界不明显的情况下也能得到较好的目标分割效果。但由于在SCR较低的图像中，目标的部分区域已被湮没在杂波中（如图5（a）中目标的上下边缘），因此分割后的图像中的目标也在一定程度上残缺了。随着待分割SAR图像SCR的提高则可相应地获得更好的分割效果，如图C和图D中的SAR图像均在分割后可得到保存较完好的目标。

图5 使用各SCR环境的OKTAL仿真数据时的算法处理效果对比Fig.5 Segmentation results of OKTAL simulation data in different SCR scenarios

表1 各SCR环境下算法的分割性能对比Tab.1 Performance comparison of the proposed method in different SCR scenarios

4.2 真实SAR图像实验对比

本节使用MSTAR数据集内的高分辨SAR图像更进一步地测试算法的分割性能，并与其它算法对比分割结果。高分辨SAR图像详尽地展现了目标的细节信息，因此适合衡量算法对图像前景和背景的精确分割能力。原SAR图像和处理结果见图6。

分割结果方面，基于MRF模型的算法和本文算法均得到了较好的目标分割结果。然而，对于能量相对强的背景杂波，基于MRF模型的算法有一定背景像素被误判。分析表2中数据可知，本文算法的虚警率则相对较低，这意味着其在防止误判背景像素这方面有较好的表现。值得注意的是，尽管k-means算法的虚警率为最低的0%，但这是以擦除部分感兴趣目标为代价才得到的。此外，本文算法获得了最高的正确率和分割精度。这是因为本算法擅长分析图像像素的空间模式，即使对能量较强的杂波背景也可通过分析空间模式的方式得到更好的处理结果。在计算复杂度方面，k-means算法需要处理0.04 s，基于MRF模型的方法需要处理0.31 s，本文算法需要处理4.85 s。

图6 真实SAR图像的分割结果对比Fig.6 Comparison of the segmentation results of different methods on the real SAR image

表2 本文算法与其它算法的目标分割性能对比Tab.2 Performance comparison between the proposed method and other methods

总地来讲，本文算法处理后图像的正确率为93.08%，虚警率为0.62%，证明了本文算法在分析SAR图像空间模式方面的有效性及在SAR图像目标分割任务中的优势。

5 结论

（1） 本文提出一种有效的目标分割算法，此算法分割时也不会削弱与参考杂波不同的模式（空间分布规律）。

（2） 实验结果表明本文方法达到了预期的效果，在不同的信杂比下均有较好的目标分割效果和可靠性。

（3） 使用OKTAL软件获得了多种信杂比条件下的仿真数据，用于验证算法的有效性，为未来研究提供了更多可用的实验数据。

［1］ 程江华， 高贵， 库锡树， 等.高分辨率SAR图像道路交叉口检测与识别新方法［J］.雷达学报， 2012， 1（1）: 100-108.Cheng Jiang-hua， Gao Gui， Ku Xi-shu， et al..A novel method for detecting and identifying road junctions from high resolution SAR images［J］.Journal of Radars， 2012，1（1）: 100-108.

［2］ 李光廷， 杨亮， 黄平平， 等.SAR图像相干斑抑制中的像素相关性测量［J］.雷达学报， 2012， 1（3）: 301-308.Li Guang-ting， Yang Liang， Huang Ping-ping， et al..The pixel-similarity measurement in SAR image despeckling［J］.Journal of Radars， 2012， 1（3）: 301-308.

［3］Sauvola J and Pietikäinen M.Adaptive document image binarization［J］.Pattern Recognition， 2000， 33（2）: 225-236.

［4］Otsu N.A threshold selection method from gray-level histograms［J］.IEEE Transactions on Systems， Man， and Cybernetics， 1979， 9（1）: 62-66.

［5］Zhao X， Jiang Y， and Zhang Y.Automatic binarization method in ISAR image［C］.IEEE International Geoscience and Remote Sensing Symposium， Milan， 2015: 5415-5418.

［6］Stagliano D， Lupidi A， Berizzi F， et al..Exploitation of COSMO-SkyMed system for detection of ships responsible for oil spills［C］.2012 IEEE International Geoscience and Remote Sensing Symposium， Munich， 2012: 915-918.

［7］Leng X， Ji K， Yang K， et al..A bilateral CFAR algorithm for ship detection in SAR images［J］.IEEE Geoscience and Remote Sensing Letters， 2015， 12（7）: 1536-1540.

［8］Liao M， Wang C， Wang Y， et al..Using SAR images to detect ships from sea clutter［J］.IEEE Geoscience and Remote Sensing Letters， 2008， 5（2）: 194-198.

［9］Mcconnell A I and Oliver C J.Comparison of segmentation methods with standard CFAR for point target detection［C］.Proceedings SPIE 3497， SAR Image Analysis， Modeling， and Techniques， 1998.doi: 10.1117/12.331364.

［10］Lankoande O， Hayat M M， and Santhanam B.Segmentation of SAR images based on Markov random field model［C］.IEEE International Conference on Systems， Man and Cybernetics， 2005， 3: 2956-2961.

［11］黄宇， 付琨， 吴一戎.基于Markov随机场K-Means图像分割算法［J］.电子学报， 2015， 37（12）: 2700-2704.HUANG Yu， FU Kun， and WU Yi-Rong.Image segmentation method using K-means based on Markov random field［J］.Acta Electronica Sinica， 2009， 37（12）: 2700-2704.

［12］Fowlkes C， Belongie S， Fan C， et al..Spectral grouping using the Nyström method［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence， 2004， 26（2）: 214-225.

［13］Zhang X， Hao L， Liu F， et al..Spectral clustering ensemble applied to SAR image segmentation［J］.IEEE Transactions on Geoscience and Remote Sensing， 2008， 46（7）: 2126-2136.

［14］Kusakunniran W， Wu Q， Zhang J， et al..A new viewinvariant feature for cross-view gait recognition［J］.IEEE Transactions on Information Forensics and Security， 2013，8（10）: 1642-1653.

［15］Zhao X， Jiang Y， Stathaki T， et al..Gait recognition method for arbitrary straight walking paths using appearance conversion machine［J］.Neurocomputing， 2015，173（3）: 530-540.

［16］Huang G B， Zhou H， Ding X， et al..Extreme learning machine for regression and multiclass classification［J］.IEEE Transactions on Systems， Man， and Cybernetics， Part B: Cybernetics， 2012， 42（2）: 513-529.

［17］Huang G B， Chen L， and Siew C K.Universal approximation using incremental constructive feedforward networks with random hidden nodes［J］.IEEE Transactionson Neural Networks， 2006， 17（4）: 879-892.

［18］Huang G B and Chen L.Enhanced random search based incremental extreme learning machine［J］.Neurocomputing，2008， 71（16/18）: 3460-3468.

［19］Huang G B， Zhu Q Y， and Siew C K.Extreme learning machine: theory and applications［J］.Neurocomputing， 2006，70（1/3）: 489-501.

［20］Gao G， Zhao L， Zhang J， et al..A segmentation algorithm for SAR images based on the anisotropic heat diffusion equation［J］.Pattern recognition， 2008， 41（10）: 3035-3043.

Target Segmentation Method in SAR Images Based on Appearance Conversion Machine

Zhao Xiaohui Jiang Yicheng Zhu Tongyu
（Harbin Institute of Technology， Harbin 150001， China）

Differences between the spatial pattern （pixel intensity and distribution） of targets and clutter allow target segmentation to be achieved by analyzing spatial patterns in Synthetic Aperture Radar （SAR） images.This paper thus proposes a target segmentation method for SAR images based on the appearance conversion machine theory.The proposed method analyses the spatial patterns in SAR images and calculates the degree of similarity between the SAR image and the reference clutter images.Subsequently， regions that show high similarity to reference clutter images are erased so that segmentation can be achieved.To evaluate the degree of similarity， we also use an automatic threshold selection method based on the cumulative histogram of the similarity imge.Experimental results using simulation and real data verify the effectiveness of the proposed method.

Synthetic Aperture Radar （SAR）； Target segmentation； Appearance Conversion Machine （ACM）；Extreme Learning Machine （ELM）

The National Natural Science Foundation of China （201306120111）

TN958.3

A

2095-283X（2016）04-0402-08

10.12000/JR16066

赵晓辉（1988-），男，内蒙古赤峰人，哈尔滨工业大学博士研究生，主要研究方向为机器学习和图像目标识别。

姜义成（1964-），男，黑龙江哈尔滨人，教授，哈尔滨工业大学电子与信息学院电子工程系主任，博士生导师，主要研究方向为雷达信号处理。

朱同宇（1992-），男，黑龙江哈尔滨人，哈尔滨工业大学硕士研究生，主要研究方向为雷达成像与目标识别。

2016-04-05；改回日期：2016-06-20；网络出版：2016-07-08

姜义成 jiangyc@hit.edu.cn

国家自然科学基金资助项目（201306120111）