一种极化SAR图像的超像素快速生成方法技术

本发明专利技术提供一种极化SAR图像的超像素快速生成方法。技术方案包括下述步骤：第一步：初始化；设不稳定点集的元素为极化SAR图像中的所有像素点，将不稳定点集里的元素称为不稳定点，设不稳定点标签为不稳定点对应像素点的标签值；第二步：更新不稳定点标签，第三步：更新图像区域中心、不稳定点集、迭代次数；第四步：基于不相似度度量的后处理。通过上述四步所得到的一个图像区域就是一个超像素。本发明专利技术可使得面积小的且不同于其周边邻域图像区域的超像素得到保留。

A super pixel fast generation method for polarimetric SAR images

The invention provides a super pixel fast generation method for polarized SAR images. Technical scheme comprises the following steps: the first step: initialization; unstable set point set elements for all pixels in polarimetric SAR images, the unstable points in the set of elements called unstable points, with unstable labels for the unstable points corresponding to the pixel label value; the second step: the update is not stable the third step: update the point label, image region center, unstable point set, the number of iterations; the fourth step: Based on similarity measure postprocessing. An image region obtained by the above four steps is a super pixel. The present invention enables reserved pixels with small areas and different image areas from their neighboring neighbors.

【技术实现步骤摘要】
一种极化SAR图像的超像素快速生成方法
本专利技术涉及图像分割
，更具体地说，涉及一种极化SAR(SyntheticApertureRadar，合成孔径雷达)图像的超像素快速生成方法。
技术介绍
极化SAR有四个通道：HH(水平发送水平接收)极化方式、HV(水平发送垂直接收)极化方式、VH(垂直发送水平接收)极化方式、VV(垂直发送垂直接收)极化方式，因此能提供更多的信息。极化SAR图像超像素生成是一种通过对SAR图像像素进行聚类，得到一系列大小相近的具有特定意义的图像区域的方法。极化SAR图像超像素生成本质上是一种图像分割方法，伴随着极化SAR图像的兴起，基于超像素的极化SAR图像处理也得到越来越多的关注。现有的超像素生成方法，包括Normalized-Cut方法(见文献Learningaclassificationmodelforsegmentation.NinthIEEEInt.Conf.onComputerVision:2003,pp.10-17.)、QuickShift方法(Quickshiftandkernelmethodsformodeseeking.EuropeanConf.onComputerVision:2008；pp.705-718.)和SLIC(Slicsuperpixels.Dept.SchoolComput.Commun.Sci.,EPFL,Lausanne,Switzerland,Tech.Rep,2012；149300.)方法等，其中SLIC方法在边缘贴合度和效率上更具有优势。SLIC方法包括两个主要步骤：(1)利用局部k-均值算法对图像像素聚类；(2)利用连通分量算法消除孤立的图像小区域。由于极化SAR图像受相干斑噪声影响，直接采用SLIC方法时，生成的超像素边缘贴合度较低。因此，学者们已经提出了一些针对极化SAR图像超像素生成的SLIC改进算法(参见文献SuperpixelsegmentationforpolarimetricSARimageryusinglocaliterativeclustering.IEEEGeosci.RemoteSens.Lett.2015,12,13-17)，用修正Wishart距离替代SLIC方法中的CIELAB颜色空间的欧氏距离，此方法生成超像素的边缘贴合度较高，但是计算效率太低。因此，研究极化SAR图像的超像素快速生成方法具有重要的理论意义和应用价值。
技术实现思路
本专利技术提出了一种极化SAR图像的超像素快速生成方法。本方法能够有效地克服极化SAR图像中相干斑噪声的影响，从而生成边缘贴合度较高的超像素，同时还能够极大地提高计算效率。本专利技术的基本思路是：首先，将极化SAR图像中所有的像素点都纳入不稳定点集，使得初始不稳定点集包含极化SAR图像中所有真实的边缘点；然后每次迭代中，基于提出的数据距离为不稳定点集中的每个不稳定点进行重贴标签，再通过判断准则更新不稳定点集；最后，基于不相似度度量进行后处理，移除生成的孤立小面积超像素。本专利技术的技术方案是：一种极化SAR图像的超像素快速生成方法，具体包括下述步骤：第一步：初始化。将极化SAR图像划分为一系列相邻并且不相交的大小均为L×L像素的矩形区域，每个矩形区域赋予一个区域序号，L的取值根据实际需要确定。对极化SAR图像的每个像素赋予一个标签，标签值为该像素所在矩形区域的区域序号。图像区域定义为具有相同标签值的像素点组成的区域，图像区域的序号为该图像区域所含像素点的标签值。由图像区域的定义可知图像区域的初始化值即上述的矩形区域。对极化SAR图像的每个像素赋予一个标签变化标识，标签变化标识取值为0。对任意图像区域计算以下两个参数：计算每个图像区域中所有像素点的平均相干矩阵作为其数据中心；计算每个图像区域中所有像素点的平均坐标值作为其几何中心。设不稳定点集的元素为极化SAR图像中的所有像素点，将不稳定点集里的元素称为不稳定点。设不稳定点标签为不稳定点对应像素点的标签值。设置迭代次数v＝0。第二步：更新不稳定点标签。首先，对任意不稳定点进行下述操作：设不稳定点p的搜索范围Sp为：以不稳定点p为中心的大小为2L×2L的窗口。计算不稳定点p和图像区域Rk的综合距离D(p,Rk)：其中：Rk∈Sp表示图像区域Rk位于不稳定点p的搜索范围Sp内，即图像区域Rk的几何中心在不稳定点p的搜索范围Sp内；d1是不稳定点p与图像区域Rk之间的空间欧式距离；m是可调的紧致度因子。利用下式计算数据距离d2：其中，Gk表示图像区域Rk的数据中心，Tp表示不稳定点p的相干矩阵，上标T表示矩阵的转置。tp＝f(Tp)，wk＝f(((Gk)-1)T)，其中f(T)＝[T11,T21,T31,T12,T22,T32,T13,T23,T33]，Τij表示相干矩阵T的第i行第j列的元素，其中i＝1,2,3，j＝1,2,3，T＝Tp或((Gk)-1)T，上标-1表示矩阵的逆。然后，利用下述规则更新极化SAR图像的不稳定点标签：如果不稳定点p与其搜索范围Sp内某一图像区域Rb的综合距离小于等于p与其搜索范围Sp内其他图像区域的综合距离，则：如果不稳定点p的标签值不等于图像区域Rb的序号，则令不稳定点p对应的像素标签变化标识取值为1，并且将不稳定点p的标签值更新为图像区域Rb的序号。第三步：更新图像区域中心、不稳定点集、迭代次数。首先，计算每个图像区域的数据中心和几何中心；其次，利用下述规则重新确定不稳定点集：如果在不稳定点p对应的像素的4邻域中：存在某个像素点的标签变化标识取值为1，并且该像素点的标签值与不稳定点p的标签值不同，则在不稳定点集中保留该不稳定点p，否则，从不稳定点集中剔除该不稳定点p。然后，令迭代次数v加1。如果不稳定点集为空或者重复迭代次数v达到设定的最大值vmax，则执行第四步，否则，返回第二步。第四步：基于不相似度度量的后处理。对极化SAR图像中的每个图像区域进行下述操作：如果图像区域Rn中的像素数目小于ThN，则利用下式计算图像区域Rn与其周边邻域的图像区域Rg间的不相似度F(Rn,Rg)，图像区域Rn的周边邻域的图像区域Rg是指图像区域Rg的边缘像素点与图像区域Rn的边缘像素点互为8邻域的图像区域：其中，Rg表示图像区域Rn周边邻域的第g个图像区域，M表示图像区域Rn周边邻域的图像区域数目，Gn表示图像区域Rn的数据中心，Gg表示图像区域Rg的数据中心，上标diag表示矩阵的对角元素构成的向量。如果图像区域Rn与其周边邻域的某个图像区域Rc的不相似度F(Rn,Rc)小于等于图像区域Rn与其周边邻域的其他图像区域的不相似度，并且F(Rn,Rc)＜Cth，Gth是一个取值范围为[0,1]的预设门限值，则图像区域Rn中所有元素所对应的像素点的标签值取值为图像区域Rc的序号。通过上述四步所得到的一个图像区域就是一个超像素。本专利技术的有益效果是：本专利技术在初始化时将所有像素点归类为不稳定点，使得初始不稳定点集包含极化SAR图像中所有真实的边缘点，然后通过迭代过程对不稳定点标签进行调整，在调整时采用提出的数据距离计算综合距离，这些步骤均使极化SAR图像生成的超像素边缘贴合度高。因为在迭代过程中只对不稳定点的标签进行调本文档来自技高网...

【技术保护点】
一种极化SAR图像的超像素快速生成方法，SAR是指合成孔径雷达，其特征在于，具体包括下述步骤：第一步、初始化：将极化SAR图像划分为一系列相邻并且不相交的大小均为L×L像素的矩形区域，每个矩形区域赋予一个区域序号，L的取值根据实际需要确定；对极化SAR图像的每个像素赋予一个标签，标签值为该像素所在矩形区域的区域序号；定义具有相同标签值的像素点组成的区域为图像区域，图像区域的序号为该图像区域所含像素点的标签值；图像区域的初始化值即上述的矩形区域；对极化SAR图像的每个像素赋予一个标签变化标识，标签变化标识取值为0；对任意图像区域计算以下两个参数：计算每个图像区域中所有像素点的平均相干矩阵作为其数据中心；计算每个图像区域中所有像素点的平均坐标值作为其几何中心；设不稳定点集的元素为极化SAR图像中的所有像素点，将不稳定点集里的元素称为不稳定点；设不稳定点标签为不稳定点对应像素点的标签值；设置迭代次数v＝0；第二步、更新不稳定点标签：首先，对任意不稳定点进行下述操作：设不稳定点p的搜索范围Sp为：以不稳定点p为中心的大小为2L×2L的窗口；计算不稳定点p和图像区域Rk的综合距离D(p,Rk)：

【技术特征摘要】
1.一种极化SAR图像的超像素快速生成方法，SAR是指合成孔径雷达，其特征在于，具体包括下述步骤：第一步、初始化：将极化SAR图像划分为一系列相邻并且不相交的大小均为L×L像素的矩形区域，每个矩形区域赋予一个区域序号，L的取值根据实际需要确定；对极化SAR图像的每个像素赋予一个标签，标签值为该像素所在矩形区域的区域序号；定义具有相同标签值的像素点组成的区域为图像区域，图像区域的序号为该图像区域所含像素点的标签值；图像区域的初始化值即上述的矩形区域；对极化SAR图像的每个像素赋予一个标签变化标识，标签变化标识取值为0；对任意图像区域计算以下两个参数：计算每个图像区域中所有像素点的平均相干矩阵作为其数据中心；计算每个图像区域中所有像素点的平均坐标值作为其几何中心；设不稳定点集的元素为极化SAR图像中的所有像素点，将不稳定点集里的元素称为不稳定点；设不稳定点标签为不稳定点对应像素点的标签值；设置迭代次数v＝0；第二步、更新不稳定点标签：首先，对任意不稳定点进行下述操作：设不稳定点p的搜索范围Sp为：以不稳定点p为中心的大小为2L×2L的窗口；计算不稳定点p和图像区域Rk的综合距离D(p,Rk)：其中：Rk∈Sp表示图像区域Rk位于不稳定点p的搜索范围Sp内，即图像区域Rk的几何中心在不稳定点p的搜索范围Sp内；d1是不稳定点p与图像区域Rk之间的空间欧式距离；m是可调的紧致度因子，根据实际情况确定；利用下式计算数据距离d2：其中，Gk表示图像区域Rk的数据中心，Tp表示不稳定点p的相干矩阵，上标T表示矩阵的转置；tp＝f(Tp)，wk＝f(((Gk)-1)T)，其中f(T)＝[T11,T21,T31,T12,T22,T32,T13,T23,T33]，Τij表示相干矩阵T的第i行第j列的元素，其中i＝1,2,3，j＝1,2,3，T＝Tp或((Gk)-1)T，上标-1表示矩阵的逆；然后，利用下...

【专利技术属性】
技术研发人员：邹焕新，张月，周石琳，孙浩，计科峰，雷琳，杜春，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南,43