本发明公开了一种基于曲线波滤波器和卷积结构学习的SAR图像分割方法，主要解决现有技术分割SAR图像不准确的问题。其实现步骤是：(1)SAR图像素描化，得到素描图；(2)根据SAR图像的区域图，划分SAR图像的像素子空间；(3)构建曲线波滤波器集合；(4)构造卷积结构学习模型；(5)采用基于曲线波滤波器和卷积结构学习模型的SAR图像分割方法，分割混合聚集结构地物像素子空间；(6)基于素描线聚拢特征的独立目标分割；(7)基于视觉语义规则的线目标分割；(8)采用基于多项式逻辑回归先验模型对匀质区域像素子空间进行分割；(9)合并分割结果，得到SAR图像分割结果。本发明获得了SAR图像的良好分割效果，可用于SAR图像的语义分割。
1.一种基于曲线波滤波器和卷积结构学习的SAR图像分割方法，包括如下步骤：(1)SAR图像素描化：
(1a)对输入的合成孔径雷达SAR图像，根据SAR图像的分布特点得到其素描模型；
(1b)利用SAR图像素描模型，对输入的合成孔径雷达SAR图像进行素描化处理，获得输入的合成孔径雷达SAR图像对应的素描图；
(2a)采用素描线区域化方法，对合成孔径雷达SAR图像的素描图进行区域化处理，得到包括聚集区域、无素描线区域和结构区域的合成孔径雷达SAR图像的区域图；
(2b)将包括聚集区域、无素描线区域和结构区域的区域图，映射到输入的合成孔径雷达SAR图像中，得到合成孔径雷达SAR图像中的混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间；
(3a)从合成孔径雷达SAR图像的区域图中提取混合聚集结构地物像素子空间的对应的聚集区域，在[0°，180°]区间内以间隔为10°划分为18个区间即18个方向，分别统计每个区间内该聚集区域中素描线段的线段条数；
(3b)对该聚集区域所有的素描线段条数，按照各个间隔内的线段条数的多少进行排序，得到方向的排序序列，将方向的排序序列中前6个方向的度数分别作为曲线波滤波器中的方向参数；
(3c)将曲线波滤波器的尺度参数在[-4,0]区间内，以间隔为0.2进行离散化，得到离散化后的尺度参数；将曲线波滤波器的位移参数在[0，9]区间内以间隔为0.5进行离散化，得到离散化后的位移参数；
其中，t表示曲线波滤波器的曲线波函数，b表示曲线波滤波器中的尺度参数，Db表示将尺度参数b作为参数的曲线波滤波器中的尺度算子，θ表示曲线波滤波器中的方向参数，Rθ表示将方向参数θ作为参数的曲线波滤波器中的旋转算子，(m,n)表示曲线波滤波器中像素点的坐标位置，m0表示水平方向的位移参数，n0表示垂直方向的位移参数， 表示将水平方向的位移参数m0和垂直方向的位移参数n0作为参数的曲线波滤波器中的的位移参数；
其中，c(t)表示以曲线波函数t作为参数的曲线波滤波器，cos表示余弦操作，exp表示以自然常数e为底的指数操作；
(3f)将计算得到的每一个曲线波滤波器组合成为曲线波滤波器集合；
(4a)对合成孔径雷达SAR图像中的混合聚集结构地物像素子空间的每个互不连通的区域，按31×31的窗口进行隔一滑窗采样，得到每个区域对应的多个图像块，将多个图像块依次输入到卷积结构学习模型中，得到卷积结构学习模型的输入层；
(4b)采用曲线波滤波器，对卷积结构学习模型的输入层中的图像块进行卷积操作，得到卷积结构学习模型的卷积层；
其中，E(c)表示数据保真项，c表示卷积结构学习模型卷积层中的滤波器，N表示卷积结构学习模型中要学习的每个互不连通的区域包含的图像块的总个数，||·||F表示做frobenius范数操作， 表示frobenius范数的平方操作，xi表示卷积结构学习模型中输入的第i个图像块， 表示提取xi中间的大小为n×n的特征图像块，Mi表示卷积结构学习模型中输入的第i个图像块对应的曲线波滤波器的总数，∑表示求和操作，*表示卷积操作， 表示卷积结构学习模型中输入的第i个图像块对应的第j个曲线波滤波器；
其中，G(c)表示结构保真项，R(·)表示求素描图中所有素描线总长度的操作，SM(·)表示提取与输入图像块一一对应的素描图块的操作；
其中，L(c)表示目标函数， 表示在目标函数L(c)值最小时，求取曲线波滤波器c的操作；
(4f)输出由目标函数指导学习得到的曲线波滤波器，得到卷积结构模型的输出层；
(5b)将步骤(4a)采样得到图像块依次输入到卷积结构学习模型中；
(5c)从曲线波滤波器集合中，随机选取六个滤波器，其方向参数由步骤(3b)中统计的6个方向得到，其位移参数和尺度参数随机初始化，将这些初始六个滤波器组成的滤波器集合作为所选取的曲线波滤波器集合；
(5d)将当前输入的图像块与所选取的曲线波滤波器集合中的每个曲线滤波器进行卷积操作,得到与每个曲线滤波器对应的特征图；
(5e)利用步骤(4d)中的结构保真项公式，计算输入图像块的结构保真项；
(5f)判断当前输入图像块的结构保真项是否小于结构误差阈值，若是，则执行步骤(5i)，否则，执行步骤(5g)；
(5g)利用尺度参数更新公式和位移参数更新公式，分别更新步骤(4c)数据保真项公式中曲线波滤波器的尺度参数和位移参数,得到更新后的曲线波滤波器；
(5h)将更新后的曲线波滤波器集合作为所选取的曲线波滤波器集合，返回步骤(5d)，对输入图像块重新进行训练；
(5i)将学习得到的曲线波滤波器保存至该输入图像块学习好的曲线波滤波器集合中，完成对该输入图像块的特征的学习，并输出该输入图像块学习好的曲线波滤波器集合；
(5j)判断所有图像块是否通过卷积结构学习模型完成了特征的学习，若是，结束程序，否则，输入下一个图像块并执行步骤(5c)；
(6a)将所有互不连通区域的特征集合拼接，将拼接后的特征集合作为码本；
(6b)对每个互不连通区域的所有特征，分别计算与码本中的每个特征的内积，得到每个区域所有特征在码本上的投影向量；
(6c)对每个互不连通区域的所有投影向量进行最大值汇聚，得到每个区域对应的一个结构特征向量；
(6d)利用近邻传播AP聚类算法，对所有互不连通区域的结构特征向量进行聚类，得到混合聚集结构地物像素子空间的分割结果；
(7c)对线目标和独立目标分割的结果进行合并，得到结构像素子空间的分割结果；
采用基于多项式逻辑回归先验模型的无素描线区域分割方法，对匀质区域像素子空间进行分割，得到匀质区域像素子空间的分割结果；
将混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间的分割结果进行合并，得到合成孔径雷达SAR图像的最终分割结果。
2.根据权利要求1所述的基于曲线波滤波器和卷积结构学习 的SAR图像分割方法，其特征在于，步骤(1b)中所述利用SAR图像素描模型，对输入的合成孔径雷达SAR图像进行素描化处理的具体步骤如下：第1步，在[100,150]范围内，任意选取一个数，作为模板的总数；
第2步，构造具有不同方向和尺度的由像素点组成的边、线的一个模板，利用模板的方向和尺度信息构造各向异性高斯函数，通过该高斯函数，计算模板中每个像素点的加权系数，统计模板中所有像素点的加权系数，其中，尺度个数取值为3～5，方向个数取值为18；
第3步，按照下式，计算与模板区域坐标相对应的合成孔径雷达SAR图像中像素点的均值：其中，μ表示与模板区域坐标相对应的合成孔径雷达SAR图像中所有像素点的均值，∑表示求和操作，g表示模板第Ω个区域中任意一个像素点对应的坐标，∈表示属于符号，wg表示模板第Ω个区域中像素点在坐标g处的权重系数，wg的取值范围为wg∈[0,1]，Ag表示与模板第Ω个区域中像素点在坐标g处对应的合成孔径雷达SAR图像中的像素点的值；
第4步，按照下式，计算与模板区域坐标相对应的合成孔径雷达SAR图像中像素点的方差值：其中，ν表示与模板区域坐标相对应的合成孔径雷达SAR图像中所有像素点的方差值；
第5步，按照下式，计算合成孔径雷达SAR图像中每个像素点针对比值算子的响应值：其中，R表示合成孔径雷达SAR图像中每个像素点针对比值算子的响应值，min{·}表示最小值操作，a和b分别表示模板中的两个不同的区域，μa表示模板区域a中所有像素点的均值，μb表示模板区域b中所有像素点的均值；
第6步，按照下式，计算合成孔径雷达SAR图像中每个像素针对相关性算子的响应值：其中，C表示合成孔径雷达SAR图像中每个像素针对相关性算子的响应值，表示平方根操作，a和b分别表示模板中两个不同区域，νa表示模板区域a中所有像素点的方差值，νb表示模板区域b中所有像素点的方差值，μa表示模板区域a中所有像素点的均值，μb表示模板区域b中所有像素点的均值；
第7步，按照下式，计算合成孔径雷达SAR图像中每个像素点针对各个模板的响应值：其中，F表示合成孔径雷达SAR图像中每个像素点针对各个模板的响应值， 表示平方根操作，R和C分别表示合成孔径雷达SAR图像中像素点针对比值算子和合成孔径雷达SAR图像中像素点针对相关性算子的响应值；
第8步，判断所构造的模板是否等于所选取模板的总数，若是，则执行第9步，否则，执行第2步；
第9步，从各个模板中选择具有最大响应值的模板，作为合成孔径雷达SAR图像的模板，并将该模板的最大响应值作为合成孔径雷达SAR图像中像素点的强度，将该模板的方向作为合成孔径雷达SAR图像中像素点的方向，获得合成孔径雷达SAR图像的边线响应图和梯度图；
第10步，按照下式，计算合成孔径雷达SAR图像强度图的强度值，得到强度图：其中，I表示合成孔径雷达SAR图像强度图的强度值，r表示合成孔径雷达SAR图像边线响应图中的值，t表示合成孔径雷达SAR图像梯度图中的值；
第11步，采用非极大值抑制方法，对强度图进行检测，得到建议草图；
第12步，选取建议草图中具有最大强度的像素点，将建议草图中与该最大强度的像素点连通的像素点连接形成建议线段，得到建议素描图；
第13步，按照下式，计算建议素描图中素描线的编码长度增益：
其中，CLG表示建议素描图中素描线的编码长度增益，∑表示求和操作,J表示当前素描线邻域中像素点的个数，Aj表示当前素描线邻域中第j个像素点的观测值，Aj,0表示在当前素描线不能表示结构信息的情况下，该素描线邻域中第j个像素点的估计值，ln(·)表示以e为底的对数操作，Aj,1表示在当前素描线能够表示结构信息的情况下，该素描线邻域中第j个像素点的估计值；
第14步，在[5,50]范围内，任意选取一个数，作为阈值T；
第15步，选出所有建议素描线中CLG>T的建议素描线，将其组合成合成孔径雷达SAR图像的素描图。
3.根据权利要求1所述的基于曲线波滤波器和卷积结构学习 的SAR图像分割方法，其特征在于，步骤(2a)中所述的素描线区域化方法的具体步骤如下：第1步，按照合成孔径雷达SAR图像的素描图中素描线段的聚集度，将素描线划分为表示聚集地物的聚集素描线和表示边界、线目标、孤立目标的边界素描线、线目标素描线、孤立目标素描线；
第2步，根据素描线段聚集度的直方图统计，选取聚集度等于最优聚集度的素描线段作为种子线段集{Ek,k＝1,2,...,m}，其中，Ek表示种子线段集中的任一条素描线段，k表示种子线段集中任一条素描线段的标号，m表示种子线段的总条数，{·}表示集合操作；
第3步，以未被选取添加进种子线段集和的线段作为基点，以此基点递归求解线段集合；
第4步，构造一个半径为最优聚集度区间上界的圆形基元，用该圆形基元对线段集合中的线段进行膨胀，对膨胀后的线段集合由外向内进行腐蚀，在素描图上得到以素描点为单位的聚集区域；
第5步，对表示边界、线目标以及孤立目标的素描线，以每个素描线的每个素描点为中心构造大小为5×5的几何结构窗，得到结构区域；
第6步，将素描图中除去聚集区域和结构区域以外的部分作为不可素描区域；
第7步，将素描图中的聚集区域、结构区域和不可素描区域合并，得到包括聚集区域、无素描线区域和结构区域的合成孔径雷达SAR图像的区域图。
4.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(5g)中所述的尺度参数的更新公式如下：其中，bt表示第t次迭代的曲线波滤波器的尺度参数，bt-1表示第t-1次迭代的尺度参数，η表示步长，取值范围为[0,1]，c(·)表示曲线波滤波器，sin表示余弦操作，X表示曲线波滤波器的水平方向参数，其值由公式X＝DbRθ(m-m0)得到，Y表示曲线波滤波器的垂直方向参数，其值由公式Y＝DbRθ(n-n0)得到， 表示曲线波滤波器对自变量X的偏导，其值由以下公式得到：表示曲线波滤波器水平方向参数X对尺度参数的偏导，其值由以下公式得到：表示曲线波滤波器对自变量Y的偏导，其值由以下公式得到：
表示曲线波滤波器垂直方向Y对尺度参数的偏导，其值由以下公式得到:
5.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(5g)中所述的位移参数的更新公式如下：其中， 表示第t次迭代的曲线波滤波器的水平方向位移参数， 表示第t-1次迭代的水平方向位移参数，其中， 表示第t次迭代的曲线波滤波器的垂直方向位移参数， 表示第t-1次迭代的垂直方向位移参数。
6.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(7a)中所述的视觉语义规则如下：设第i条素描线li与第j条素描线lj之间的距离为Dij，li的方向为Oi，lj的方向为Oj，i,j∈[1,2,...,S]，S为素描线的总条数；
将宽度大于3个像素的线目标用两条素描线li和lj表示，li和lj之间的距离Dij小于T1且方向差(Oi-Oj)小于10度，其中T1＝5；
设第s条素描线ls的几何结构窗ws内每一列的平均灰度为Ai，设相邻列的灰度差为ADi＝|Ai-Ai+1|，设zs＝[zs1,zs2,...,zs9]为相邻列的灰度差ADi的标记向量；
将宽度小于3个像素的线目标用单个素描线ls表示，ls的几何结构窗ws内，计算相邻列的灰度差ADi，如果ADi＞T2，则zsi＝1；否则zsi＝0，zs中有两个元素的值为1，其余为0，其中T2＝34；
设L1,L2是表示线目标的素描线的集合，如果Dij＜T1并且|Oi-Oj|＜10，则li,lj∈L1；如果sum(zs)＝2，则ls∈L2，其中sum(·)表示对向量的所有分量求和的操作。
7.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(7a)中所述的分割线目标的具体步骤如下：第1步，在结构像素子空间，根据线目标的素描线的集合L1，将第i条素描线li和第j条素描线lj之间的区域作为线目标；
第2步，在结构像素子空间，根据线目标的素描线的集合L2，将覆盖第s条素描线ls的区域作为线目标。
8.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(7b)中所述的分割独立目标的具体步骤如下：第1步，在区域图的结构区域中，将不表示线目标的所有素描线标记为候选素描线集合中的素描线；
第2步，从候选素描线集合中随机选取一条素描线，以所选取的素描线的一个端点为中心，构造大小为5×5的一个几何结构窗；
第3步，判断几何结构窗内是否存在其它素描线的端点，若存在，执行第4步；否则，执行第6步；
第4步，判断是否只存在一个端点，若是，将该端点所在素描线和当前素描线进行连接；
第5步，连接所选取素描线与各端点所在的素描线，从所有连接线中选取其中夹角最大的两条素描线作为连接完成的素描线；
第6步，判断素描线的另一个端点的几何结构窗内是否存在其他素描线的端点，若存在，执行第4步；否则，执行第7步；
第7步，对完成连接操作的素描线，选取包含两条及两条以上素描线段的素描线，统计所选取素描线中包含素描线段的条数n，其中n≥2；
第8步，判断素描线的条数n是否等于2，若是，则执行第9步；否则，执行第10步；
第9步，将素描线顶点的角度值在[10°,140°]的范围内的素描线作为具有聚拢特征的素描线；
第10步，选出素描线对应的n-1个顶点的角度值都在[10°,140°]范围内的素描线；
第一种情况，判断第i条素描线段相邻两条素描线段第i-1条、第i+1条是否在第i条素描线段所在直线的同一侧，2≤i≤n-1，若素描线上的所有素描线段与相邻线段都在同一侧，则标记该素描线为具有聚拢特征的素描线；
第二种情况，判断第i条素描线段相邻两条素描线段第i-1条、第i+1条是否在第i条素描线段所在直线的同一侧，2≤i≤n-1，若素描线上有n-1条素描线段与相邻线段在同一侧，而有一条素描线段与其相邻线段在非同一侧，也标记该素描线为具有聚拢特征的素描线；
第12步，在具有聚拢特征的素描线中任选一条素描线，由所选取素描线的两个端点坐标，确定两个端点间的距离，若该端点距离在[0,20]范围内，则将所选取素描线作为表示独立目标的素描线；
第13步，判断未处理的具有聚拢特征的素描线是否全部选完，若是，执行第14步；否则，执行第12步；
第14步，用超像素分割的方法，对合成孔径雷达SAR图像中表示独立目标的素描线周围的像素进行超像素分割，将分割后超像素的灰度值在[0,45]或[180,255]内的超像素作为独立目标超像素；
第15步，合并独立目标超像素，将合并后的独立目标超像素的边界作为独立目标的边界，得到独立目标的分割结果。
9.根据权利要求1所述的基于曲线波滤波器和卷积结构学习的SAR图像分割方法，其特征在于，步骤(8)中所述基于多项式逻辑回归先验模型的无素描线区域分割方法的具体步骤如下：第1步，从匀质区域像素子空间中任意选取一个像素点，以所选取的像素点为中心建立
第2步，将方形窗口的边长增加2，得到新的方形窗口，计算新方形窗口的标准差σ2；
第3步，设标准差阈值T3＝3，如果|σ1-σ2|＜T3，则将标准差为σ2的方形窗口作为最终的方形窗口，执行第4步；否则，执行第2步；
第4步，按照下式，计算方形窗口内中心像素的先验概率：
其中，p′1表示方形窗口内中心像素的先验概率，η'表示概率模型参数，η'取值为1，xk′′表示方形窗口内属于第k'类的像素个数，k'∈[1,...,i',...,K']，K'表示分割的类别数，K'取值为5，xi′′表示第3步得到的方形窗口内属于第i'类的像素个数；
第5步，将像素灰度的概率密度与纹理的概率密度相乘，得到似然概率p'2，其中，灰度的概率密度由衰落信道Nakagami分布得到，纹理的概率密度由t分布得到；
第6步，将先验概率p1'与似然概率p2'相乘，得到后验概率p12'；
第7步，判断匀质区域像素子空间中是否还有未处理的像素点，若有，执行第1步；否则，执行第8步；
第8步，根据最大后验概率准则，得到匀质区域像素子空间的分割结果。
基于曲线波滤波器和卷积结构学习的SAR图像分割方法
技术领域
[0001] 本发明属于图像处理技术领域，更进一步涉及图像分割技术领域中的一种基于曲线波滤波器和卷积结构学习的合成孔径雷达SAR(Synthetic Aperture Radar)图像分割方法。本发明能够对合成孔径雷达SAR图像的具有不同特征的区域准确地进行分割，并且可用于后续的合成孔径雷达SAR图像的目标检测与识别。
背景技术
[0002] 合成孔径雷达SAR由于其特殊的成像机制获得了在军事和民用多个领域的广泛应用。SAR具有全天时、全天候、多波段、多极化、可变侧视角以及高分辨率等优点，不仅可以详细、准确的测绘地形、地貌，获取地球表面的信息，还可以透过地表和自然植被收集地下信息，甚至在恶劣的环境下也能以较高的分辨率提供详细的地面测绘数据和图像。
[0003] 作为图像理解与解译中基本而关键的技术之一，图像分割是连接图像处理和图像分析理解的纽带，在图像工程中占据重要的地位，分割结果的优劣将影响随后图像分析、理解的质量和效率。但是SAR图像却有相干斑噪声明显，同类目标差异较大，异类目标非常接近，多尺度目标的同时存在，数据规模庞大，缺乏有效的语义表示等问题，这些问题都增加了SAR图像分割的难度。
[0004] 同时SAR图像分割需要引入高层语义知识，当人们面对一幅图像的时候，关注的并不是一个一个的像素，而是多个像素聚集而成的目标，这些有意义的目标才是人们理解图像的基本单元，而像素只不过是图像离散化的数字串而已，因此图像的语义单元与单个的像素是有差别的，这也就是著名的“语义鸿沟”问题。作为SAR图像理解前期工作的SAR图像分割，也包含有语义的要求。然而，与之相关的研究却非常少，现有的语义表示方法均是建立在自然图像之上的，如何针对SAR图像设计有效的语义规则，也是一个亟待解决的问题。
[0005] 中国科学院对地观测与数字地球科学中心的钟微宇和沈汀在其发表的论文“结合多特征和SVM的SAR图像分割”(《计算机应用研究》，2013，30(9):2846-2851.)中提出了一种结合非下采样轮廓变换(NSCT)和GLCM的纹理特征提取方法来进行SAR图像分割。该方法为实现灰度共生矩阵(GLCM)多尺度、多方向的纹理特征提取，提出了一种结合非下采样轮廓变换(NSCT)和GLCM的纹理特征提取方法。先用NSCT对合成孔径雷达(SAR)图像进行多尺度、多方向分解；再对得到的子带图像使用GLCM提取灰度共生量；然后对提取的灰度共生量进行相关性分析，去除冗余特征量，并将其与灰度特征构成多特征矢量；最后，充分利用支持向量机(SVM)在小样本数据库和泛化能力方面的优势，由SVM完成多特征矢量的划分，实现SAR图像分割。但是，该方法仍然存在的不足之处是，没有引入SAR图像的高层语义知识，只在像素级别对SAR图像进行了分割，导致了SAR图像分割结果不准确。
[0006] 武汉大学的林达等人在其发表的论文“一种有效的MSTAR SAR图像分割方法”(《武汉大学学报》，2015，40(10)：1377-1380.)中提出了一种MSTAR SAR图像分割方法。该方法首先对待处理图像进行过分割操作，得到过分割图像区域。其次对过分割后的图像进行图像区域级和像素级的特征提取，得到用于表示图像的特征向量，对MSTAR SAR图像使用空间隐含狄利克雷分配模型(sLDA)和马尔科夫随机场(MRF)建立本文所提出的模型，得到能量泛函。最后运用Graph-Cut算法和Branch-and-Bound算法对能量泛函进行优化，得到最终的分割结果。该方法存在的不足之处是，在求得SAR图像的特征向量时，没有学习SAR图像中由于像素之间的相关性而特有的结构特征。
[0007] 刘芳，段一平，李玲玲，焦李成等在其发表的论文“基于层次视觉语义和自适应邻域多项式隐模型的SAR图像分割”(IEEE Trancactions on Geoscience and Remote Sensing，2016,54(7):4287-4301.)中提出了一种基于层次视觉语义和自适应邻域多项式隐模型的SAR图像分割方法，该方法在SAR图像素描图的基础上，提出了SAR图像的层次视觉语义。该层次视觉语义把SAR图像划分为聚集区域，结构区域和匀质区域。基于该划分，对不同特性的区域采用了不同的分割方法。对于聚集区域，提取灰度共生矩阵特征，并采用局部线性约束编码的方法得到每个聚集区域的表示，进而采用层次聚类的方法进行分割。对结构区域，通过分析边模型和线模型，设计了视觉语义规则定位边界和线目标。另外，边界和线目标包含了强烈的方向信息，因此设计了基于几何结构窗的多项式隐模型进行分割。对匀质区域，为了能找到恰当的邻域去表示中心像素，设计了基于自适应窗口的多项式隐模型进行分割。这三个区域的分割结果被整合到一起得到最后的分割结果。该方法的不足之处是，对聚集区域的边界定位不够精确，对匀质区域类别数的确定不够合理，分割结果的区域一致性较差，而结构区域的分割中未对独立目标进行处理。
发明内容
[0008] 本发明的目的在于针对上述已有技术的不足，提出一种基于曲线波滤波器和卷积结构模型的SAR图像分割方法，以更加有效的完成SAR图像的分割。
[0009] 为实现上述目的，本发明的技术方案如下：
[0010] (1)SAR图像素描化：
[0011] (1a)对输入的合成孔径雷达SAR图像，根据SAR图像的分布特点得到其素描模型；
[0012] (1b)利用SAR图像素描模型，对输入的合成孔径雷达SAR图像进行素描化处理，获得输入的合成孔径雷达SAR图像对应的素描图；
[0013] (2)划分像素子空间：
[0014] (2a)采用素描线区域化方法，对合成孔径雷达SAR图像的素描图进行区域化处理，得到包括聚集区域、无素描线区域和结构区域的合成孔径雷达SAR图像的区域图；
[0015] (2b)将包括聚集区域、无素描线区域和结构区域的区域图，映射到输入的合成孔径雷达SAR图像中，得到合成孔径雷达SAR图像中的混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间；
[0016] (3)构建曲线波滤波器集合：
[0017] (3a)从合成孔径雷达SAR图像的区域图中提取混合聚集结构地物像素子空间的对应的聚集区域，在[0°，180°]区间内以间隔为10°划分为18个区间即18个方向，分别统计每个区间内该聚集区域中素描线段的线段条数；
[0018] (3b)对该聚集区域所有的素描线段条数，按照各个间隔内的线段条数的多少进行排序，得到方向的排序序列，将方向的排序序列中前6个方向的度数分别作为曲线波滤波器中的方向参数；
[0019] (3c)将曲线波滤波器的尺度参数在[-4,0]区间内，以间隔为0.2进行离散化，得到离散化后的尺度参数；将曲线波滤波器的位移参数在[0，9]区间内以间隔为0.5进行离散化，得到离散化后的位移参数；
[0020] (3d)按照下式，计算曲线波滤波器的曲线波函数：
[0021]
[0022] 其中，t表示曲线波滤波器的曲线波函数，b表示曲线波滤波器中的尺度参数，Db表示将尺度参数b作为参数的曲线波滤波器中的尺度算子，θ表示曲线波滤波器中的方向参数，Rθ表示将方向参数θ作为参数的曲线波滤波器中的旋转算子，(m,n)表示曲线波滤波器中像素点的坐标位置，m0表示水平方向的位移参数，n0表示垂直方向的位移参数， 表示将水平方向的位移参数m0和垂直方向的位移参数n0作为参数的曲线波滤波器中的的位移参数；
[0023] (3e)按照下式，计算每一个曲线波滤波器：
[0024]
[0025] 其中，c(t)表示以曲线波函数t作为参数的曲线波滤波器，cos表示余弦操作，exp表示以自然常数e为底的指数操作；
[0026] (3f)将计算得到的每一个曲线波滤波器组合成为曲线波滤波器集合；
[0027] (4)构造卷积结构学习模型：
[0028] (4a)对合成孔径雷达SAR图像中的混合聚集结构地物像素子空间的每个互不连通的区域，按31×31的窗口进行隔一滑窗采样，得到每个区域对应的多个图像块，将多个图像块依次输入到卷积结构学习模型中，得到卷积结构学习模型的输入层；
[0029] (4b)采用曲线波滤波器，对卷积结构学习模型的输入层中的图像块进行卷积操作，得到卷积结构学习模型的卷积层；
[0030] (4c)按照下式，计算数据保真项：
[0031]
[0032] 其中，E(c)表示数据保真项，c表示卷积结构学习模型卷积层中的滤波器，N表示卷积结构学习模型中要学习的每个互不连通的区域包含的图像块的总个数，||·||F表示做frobenius范数操作， 表示frobenius范数的平方操作，xi表示卷积结构学习模型中输入的第i个图像块， 表示提取xi中间的大小为n×n的特征图像块，Mi表示卷积结构学习模型中输入的第i个图像块对应的曲线波滤波器的总数，∑表示求和操作，*表示卷积操作， 表示卷积结构学习模型中输入的第i个图像块对应的第j个曲线波滤波器；
[0033] (4d)按照下式，计算结构保真项：
[0034]
[0035] 其中，G(c)表示结构保真项，R(·)表示求素描图中所有素描线总长度的操作，SM(·)表示提取与输入图像块一一对应的素描图块的操作；
[0036] (4e)按照下式，计算目标函数：
[0037]
[0038] 其中，L(c)表示目标函数， 表示在目标函数L(c)值最小时，求取曲线波滤波器c的操作；
[0039] (4f)输出由目标函数指导学习得到的曲线波滤波器，得到卷积结构模型的输出层；
[0040] (5)训练卷积结构学习模型：
[0041] (5a)将结构误差阈值设置为0.1；
[0042] (5b)将步骤(4a)采样得到图像块依次输入到卷积结构学习模型中；
[0043] (5c)从曲线波滤波器集合中，随机选取六个滤波器，其方向参数由步骤(3b)中统计的6个方向得到，其位移参数和尺度参数随机初始化，将这些初始六个滤波器组成的滤波器集合作为所选取的曲线波滤波器集合；
[0044] (5d)将当前输入的图像块与所选取的曲线波滤波器集合中的每个曲线滤波器进行卷积操作,得到与每个曲线滤波器对应的特征图；
[0045] (5e)利用步骤(4d)中的结构保真项公式，计算输入图像块的结构保真项；
[0046] (5f)判断当前输入图像块的结构保真项是否小于结构误差阈值，若是，则执行步骤(5i)，否则，执行步骤(5g)；
[0047] (5g)利用尺度参数更新公式和位移参数更新公式，分别更新步骤(4c)数据保真项公式中曲线波滤波器的尺度参数和位移参数,得到更新后的曲线波滤波器；
[0048] (5h)将更新后的曲线波滤波器集合作为所选取的曲线波滤波器集合，返回步骤(5d)，对输入图像块重新进行训练；
[0049] (5i)将学习得到的曲线波滤波器保存至该输入图像块学习好的曲线波滤波器集合中，完成对该输入图像块的特征的学习，并输出该输入图像块学习好的曲线波滤波器集合；
[0050] (5j)判断所有图像块是否通过卷积结构学习模型完成了特征的学习，若是，结束程序，否则，输入下一个图像块并执行步骤(5c)；
[0051] (6)分割SAR图像混合聚集结构地物像素子空间：
[0052] (6a)将所有互不连通区域的特征集合拼接，将拼接后的特征集合作为码本；
[0053] (6b)对每个互不连通区域的所有特征，分别计算与码本中的每个特征的内积，得到每个区域所有特征在码本上的投影向量；
[0054] (6c)对每个互不连通区域的所有投影向量进行最大值汇聚，得到每个区域对应的一个结构特征向量；
[0055] (6d)利用近邻传播AP聚类算法，对所有互不连通区域的结构特征向量进行聚类，得到混合聚集结构地物像素子空间的分割结果；
[0056] (7)分割结构像素子空间：
[0057] (7a)用视觉语义规则，分割线目标；
[0058] (7b)基于素描线的聚拢特征，分割独立目标；
[0059] (7c)对线目标和独立目标分割的结果进行合并，得到结构像素子空间的分割结果；
[0060] (8)分割匀质区域像素子空间：
[0061] 采用基于多项式逻辑回归先验模型的无素描线区域分割方法，对匀质区域像素子空间进行分割，得到匀质区域像素子空间的分割结果；
[0062] (9)合并分割结果：
[0063] 将混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间的分割结果进行合并，得到合成孔径雷达SAR图像的最终分割结果。
[0064] 本发明与现有技术相比具有以下优点：
[0065] 第一，本发明利用合成孔径雷达SAR图像素描图，采用素描线区域化方法，得到合成孔径雷达SAR图像的区域图，将区域图映射到输入的合成孔径雷达SAR图像，得到合成孔径雷达SAR图像中的混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间，在混合聚集结构地物的像素子空间中进行采样和特征学习，克服了现有技术在分割合成孔径雷达SAR图像时只在像素级别对合成孔径雷达SAR图像进行分割的不足，使得本发明具有引入了合成孔径雷达SAR图像的高层语义知识，而提升了合成孔径雷达SAR图像分割结果的准确性。
[0066] 第二，在本发明构造的卷积结构模型中，采用输入图像块的数据保真项和输入图像块的结构保真项来提取输入图像块的结构信息，克服了现有技术在求得合成孔径雷达SAR图像的特征向量时，只用到合成孔径雷达SAR图像的像素级特征分割合成孔径雷达SAR图像，而忽略了合成孔径雷达SAR图像的结构信息的的不足，使得本发明具有能够学习合成孔径雷达SAR图像中由于像素之间的相关性而特有结构特征的优点。
[0067] 第三，本发明对合成孔径雷达SAR图像中的混合聚集结构地物像素子空间的每个互不连通的区域，按32×32的窗口进行隔一滑窗采样，得到多个图像块，将得到的图像块输入到卷积结构模型中，克服了现有技术所用的自动提取图像特征的深度自编码器的输入为一维向量，会破坏图像的空间结构特征的不足，使得本发明具有可以提取图像的本质特征，提升合成孔径雷达SAR图像分割的精度的优点。
附图说明
[0068] 图1是本发明的流程图；
[0069] 图2是本发明的仿真图；
[0070] 图3是本发明独立目标分割结果图；
[0071] 图4是本发明仿真结果示意图。
具体实施方式
[0072] 下面结合附图对本发明做进一步的描述。
[0073] 参照附图1，本发明的具体步骤如下。
[0074] 步骤1，SAR图像素描化。
[0075] 对输入的合成孔径雷达SAR图像，根据SAR图像的分布特点得到其素描模型。
[0076] 按照以下步骤，利用SAR图像素描模型，对输入的合成孔径雷达SAR图像进行素描化处理，获得输入的合成孔径雷达SAR图像对应的素描图：
[0077] 第1步，在[100,150]范围内，任意选取一个数，作为模板的总数；
[0078] 第2步，构造具有不同方向和尺度的由像素点组成的边、线的一个模板，利用模板的方向和尺度信息构造各向异性高斯函数，通过该高斯函数，计算模板中每个像素点的加权系数，统计模板中所有像素点的加权系数，其中，尺度个数取值为3～5，方向个数取值为18；
[0079] 第3步，按照下式，计算与模板区域坐标相对应的合成孔径雷达SAR图像中像素点的均值：
[0080]
[0081] 其中，μ表示与模板区域坐标相对应的合成孔径雷达SAR图像中所有像素点的均值，∑表示求和操作，g表示模板第Ω个区域中任意一个像素点对应的坐标，∈表示属于符号，wg表示模板第Ω个区域中像素点在坐标g处的权重系数，wg的取值范围为wg∈[0,1]，Ag表示与模板第Ω个区域中像素点在坐标g处对应的合成孔径雷达SAR图像中的像素点的值；
[0082] 第4步，按照下式，计算与模板区域坐标相对应的合成孔径雷达SAR图像中像素点的方差值：
[0083]
[0084] 其中，ν表示与模板区域坐标相对应的合成孔径雷达SAR图像中所有像素点的方差值；
[0085] 第5步，按照下式，计算合成孔径雷达SAR图像中每个像素点针对比值算子的响应值：
[0086]
[0087] 其中，R表示合成孔径雷达SAR图像中每个像素点针对比值算子的响应值，min{·}表示最小值操作，a和b分别表示模板中的两个不同的区域，μa表示模板区域a中所有像素点的均值，μb表示模板区域b中所有像素点的均值；
[0088] 第6步，按照下式，计算合成孔径雷达SAR图像中每个像素针对相关性算子的响应值：
[0089]
[0090] 其中，C表示合成孔径雷达SAR图像中每个像素针对相关性算子的响应值， 表示平方根操作，a和b分别表示模板中两个不同区域，νa表示模板区域a中所有像素点的方差值，νb表示模板区域b中所有像素点的方差值，μa表示模板区域a中所有像素点的均值，μb表示模板区域b中所有像素点的均值；
[0091] 第7步，按照下式，计算合成孔径雷达SAR图像中每个像素点针对各个模板的响应值：
[0092]
[0093] 其中，F表示合成孔径雷达SAR图像中每个像素点针对各个模板的响应值， 表示平方根操作，R和C分别表示合成孔径雷达SAR图像中像素点针对比值算子和合成孔径雷达SAR图像中像素点针对相关性算子的响应值；
[0094] 第8步，判断所构造的模板是否等于所选取模板的总数，若是，则执行第9步，否则，执行第2步；
[0095] 第9步，从各个模板中选择具有最大响应值的模板，作为合成孔径雷达SAR图像的模板，并将该模板的最大响应值作为合成孔径雷达SAR图像中像素点的强度，将该模板的方向作为合成孔径雷达SAR图像中像素点的方向，获得合成孔径雷达SAR图像的边线响应图和梯度图；
[0096] 第10步，按照下式，计算合成孔径雷达SAR图像强度图的强度值，得到强度图：
[0097]
[0098] 其中，I表示合成孔径雷达SAR图像强度图的强度值，r表示合成孔径雷达SAR图像边线响应图中的值，t表示合成孔径雷达SAR图像梯度图中的值；
[0099] 第11步，采用非极大值抑制方法，对强度图进行检测，得到建议草图；
[0100] 第12步，选取建议草图中具有最大强度的像素点，将建议草图中与该最大强度的像素点连通的像素点连接形成建议线段，得到建议素描图；
[0101] 第13步，按照下式，计算建议素描图中素描线的编码长度增益：
[0102]
[0103] 其中，CLG表示建议素描图中素描线的编码长度增益，∑表示求和操作,J表示当前素描线邻域中像素点的个数，Aj表示当前素描线邻域中第j个像素点的观测值，Aj,0表示在当前素描线不能表示结构信息的情况下，该素描线邻域中第j个像素点的估计值，ln(·)表示以e为底的对数操作，Aj,1表示在当前素描线能够表示结构信息的情况下，该素描线邻域中第j个像素点的估计值；
[0104] 第14步，在[5,50]范围内，任意选取一个数，作为阈值T；
[0105] 第15步，选出所有建议素描线中CLG>T的建议素描线，将其组合成合成孔径雷达SAR图像的素描图。
[0106] 本发明使用的合成孔径雷达SAR图像素描模型是Jie-Wu等人于2014年发表在IEEE Transactions on Geoscience and Remote Sensing杂志上的文章《Local maximal homogenous region search for SAR speckle reduction with sketch-based geometrical kernel function》中所提出的模型。
[0107] 步骤2，划分像素子空间。
[0108] 采用素描线区域化方法，对合成孔径雷达SAR图像的素描图进行区域化处理，得到包括聚集区域、无素描线区域和结构区域的合成孔径雷达SAR图像的区域图。
[0109] 按照合成孔径雷达SAR图像的素描图中素描线段的聚集度，将素描线划分为表示聚集地物的聚集素描线和表示边界、线目标、孤立目标的边界素描线、线目标素描线、孤立目标素描线。
[0110] 根据素描线段聚集度的直方图统计，选取聚集度等于最优聚集度的素描线段作为种子线段集{Ek,k＝1,2,...,m}，其中，Ek表示种子线段集中的任一条素描线段，k表示种子线段集中任一条素描线段的标号，m表示种子线段的总条数，{·}表示集合操作。
[0111] 以未被选取添加进种子线段集和的线段作为基点，以此基点递归求解线段集合。
[0112] 构造一个半径为最优聚集度区间上界的圆形基元，用该圆形基元对线段集合中的线段进行膨胀，对膨胀后的线段集合由外向内进行腐蚀，在素描图上得到以素描点为单位的聚集区域。
[0113] 对表示边界、线目标以及孤立目标的素描线，以每个素描线的每个素描点为中心构造大小为5×5的几何结构窗，得到结构区域。
[0114] 将素描图中除去聚集区域和结构区域以外的部分作为不可素描区域。
[0115] 将素描图中的聚集区域、结构区域和不可素描区域合并，得到包括聚集区域、无素描线区域和结构区域的合成孔径雷达SAR图像的区域图。
[0116] 将包括聚集区域、无素描线区域和结构区域的区域图，映射到输入的合成孔径雷达SAR图像中，得到合成孔径雷达SAR图像中的混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间。
[0117] 步骤3，构建曲线波滤波器。
[0118] 从合成孔径雷达SAR图像的区域图中提取混合聚集结构地物像素子空间的对应的聚集区域，在[0°，180°]区间内以间隔为10°划分为18个区间即18个方向，分别统计每个区间内该聚集区域中素描线段的线段条数。
[0119] 对该聚集区域所有的素描线段条数，按照各个间隔内的线段条数的多少进行排序，得到方向的排序序列，将方向的排序序列中前6个方向的度数分别作为曲线波滤波器中的方向参数。
[0120] 将曲线波滤波器的尺度参数在[-4,0]区间内，以间隔为0.2进行离散化，得到离散化后的尺度参数；将曲线波滤波器的位移参数在[0，9]区间内以间隔为0.5进行离散化，得到离散化后的位移参数。
[0121] 按照下式，计算曲线波滤波器的曲线波函数：
[0122]
[0123] 其中，t表示曲线波滤波器的曲线波函数，b表示曲线波滤波器中的尺度参数，Db表示将尺度参数b作为参数的曲线波滤波器中的尺度算子，θ表示曲线波滤波器中的方向参数，Rθ表示将方向参数θ作为参数的曲线波滤波器中的旋转算子，(m,n)表示曲线波滤波器中像素点的坐标位置，m0表示水平方向的位移参数，n0表示垂直方向的位移参数， 表示将水平方向的位移参数m0和垂直方向的位移参数n0作为参数的曲线波滤波器中的的位移参数。
[0124] 按照下式，计算每一个曲线波滤波器：
[0125]
[0126] 其中，c(t)表示以曲线波函数t作为参数的曲线波滤波器，cos表示余弦操作，exp表示以自然常数e为底的指数操作。
[0127] 将计算得到的每一个曲线波滤波器组合成为曲线波滤波器集合。
[0128] 步骤4，构造卷积结构学习模型。
[0129] 对合成孔径雷达SAR图像中的混合聚集结构地物像素子空间的每个互不连通的区域，按31×31的窗口进行隔一滑窗采样，得到每个区域对应的多个图像块，将多个图像块依次输入到卷积结构学习模型中，得到卷积结构学习模型的输入层。
[0130] 采用曲线波滤波器，对卷积结构学习模型的输入层中的图像块进行卷积操作，得到卷积结构学习模型的卷积层。
[0131] 按照下式，计算数据保真项：
[0132]
[0133] 其中，E(c)表示数据保真项，c表示卷积结构学习模型卷积层中的滤波器，N表示卷积结构学习模型中要学习的每个互不连通的区域包含的图像块的总个数，||·||F表示做frobenius范数操作， 表示frobenius范数的平方操作，xi表示卷积结构学习模型中输入的第i个图像块，表示提取xi中间的大小为n×n的特征图像块，Mi表示卷积结构学习模型中输入的第i个图像块对应的曲线波滤波器的总数，∑表示求和操作，*表示卷积操作， 表示卷积结构学习模型中输入的第i个图像块对应的第j个曲线波滤波器。
[0134] 按照下式，计算结构保真项：
[0135]
[0136] 其中，G(c)表示结构保真项，R(·)表示求素描图中所有素描线总长度的操作，SM(·)表示提取与输入图像块一一对应的素描图块的操作。
[0137] 按照下式，计算目标函数：
[0138]
[0139] 其中，L(c)表示目标函数， 表示在目标函数L(c)值最小时，求取曲线波滤波器c的操作。
[0140] 输出由目标函数指导学习得到的曲线波滤波器，得到卷积结构模型的输出层。
[0141] 步骤5，按照以下步骤，训练卷积结构学习模型：
[0142] 第1步，将结构误差阈值设置为0.1；
[0143] 第2步，将步骤4采样得到图像块依次输入到卷积结构学习模型中；
[0144] 第3步，从曲线波滤波器集合中，随机选取六个滤波器，其方向参数由步骤(3b)中统计的6个方向得到，其位移参数和尺度参数随机初始化，将这些初始六个滤波器组成的滤波器集合作为所选取的曲线波滤波器集合；
[0145] 第4步，将当前输入的图像块与所选取的曲线波滤波器集合中的每个曲线滤波器进行卷积操作,得到与每个曲线滤波器对应的特征图；
[0146] 第5步，利用步骤4中的结构保真项公式，计算输入图像块的结构保真项；
[0147] 第6步，判断当前输入图像块的结构保真项是否小于结构误差阈值，若是，则执行第9步，否则，执行第7步；
[0148] 第7步，利用尺度参数更新公式和位移参数更新公式，分别更新步骤4数据保真项公式中曲线波滤波器的尺度参数和位移参数,得到更新后的曲线波滤波器；
[0149] 尺度参数更新公式如下：
[0150]
[0151] 其中，bt表示第t次迭代的曲线波滤波器的尺度参数，bt-1表示第t-1次迭代的尺度参数，η表示步长，取值范围为[0,1]，c(·)表示曲线波滤波器，sin表示余弦操作，X表示曲线波滤波器的水平方向参数，其值由公式X＝DbRθ(m-m0)得到，Y表示曲线波滤波器的垂直方向参数，其值由公式Y＝DbRθ(n-n0)得到， 表示曲线波滤波器对自变量X的偏导，其值由以下公式得到：
[0152]
[0153] 表示曲线波滤波器水平方向参数X对尺度参数的偏导，其值由以下公式得到：
[0154]
[0155] 表示曲线波滤波器对自变量Y的偏导，其值由以下公式得到：
[0156]
[0157] 表示曲线波滤波器垂直方向Y对尺度参数的偏导，其值由以下公式得到:
[0158]
[0159] 位移参数更新公式如下：
[0160]
[0161] 其中， 表示第t次迭代的曲线波滤波器的水平方向位移参数， 表示第t-1次迭代的水平方向位移参数，
[0162]
[0163] 其中， 表示第t次迭代的曲线波滤波器的垂直方向位移参数， 表示第t-1次迭代的垂直方向位移参数。
[0164] 第8步，将更新后的曲线波滤波器集合作为所选取的曲线波滤波器集合，返回步骤第4步，对输入图像块重新进行训练；
[0165] 第9步，将学习得到的曲线波滤波器保存至该输入图像块学习好的曲线波滤波器集合中，完成对该输入图像块的特征的学习，并输出该输入图像块学习好的曲线波滤波器集合；
[0166] 第10步，判断所有图像块是否通过卷积结构学习模型完成了特征的学习，若是，结束程序，否则，输入下一个图像块并执行第3步。
[0167] 步骤6，分割SAR图像混合聚集结构地物像素子空间。
[0168] 将所有互不连通区域的特征集合拼接，将拼接后的特征集合作为码本。
[0169] 对每个互不连通区域的所有特征，分别计算与码本中的每个特征的内积，得到每个区域所有特征在码本上的投影向量。
[0170] 对每个互不连通区域的所有投影向量进行最大值汇聚，得到每个区域对应的一个结构特征向量。
[0171] 利用近邻传播AP聚类算法，对所有互不连通区域的结构特征向量进行聚类，得到混合聚集结构地物像素子空间的分割结果。
[0172] 步骤7，分割结构像素子空间。
[0173] 用视觉语义规则，分割线目标。
[0174] 设第i条素描线li与第j条素描线lj之间的距离为Dij，li的方向为Oi，lj的方向为Oj，i,j∈[1,2,...,S]，S为素描线的总条数。
[0175] 将宽度大于3个像素的线目标用两条素描线li和lj表示，li和lj之间的距离Dij小于T1且方向差(Oi-Oj)小于10度，其中T1＝5。
[0176] 设第s条素描线ls的几何结构窗ws内每一列的平均灰度为Ai，设相邻列的灰度差为ADi＝|Ai-Ai+1|，设zs＝[zs1,zs2,...,zs9]为相邻列的灰度差ADi的标记向量。
[0177] 将宽度小于3个像素的线目标用单个素描线ls表示，ls的几何结构窗ws内，计算相邻列的灰度差ADi，如果ADi＞T2，则zsi＝1。否则zsi＝0，zs中有两个元素的值为1，其余为0，其中T2＝34。
[0178] 设L1,L2是表示线目标的素描线的集合，如果Dij＜T1并且|Oi-Oj|＜10，则li,lj∈L1；如果sum(zs)＝2，则ls∈L2，其中sum(·)表示对向量的所有分量求和的操作。
[0179] 在结构像素子空间，根据线目标的素描线的集合L1，将第i条素描线li与第j条素描线lj之间的区域作为线目标。
[0180] 在结构像素子空间，根据线目标的素描线的集合L2，将覆盖第s条素描线ls的区域作为线目标。
[0181] 按照以下步骤，基于素描线的聚拢特征，分割独立目标。
[0182] 第1步，在区域图的结构区域中，将不表示线目标的所有素描线标记为候选素描线集合中的素描线；
[0183] 第2步，从候选素描线集合中随机选取一条素描线，以所选取的素描线的一个端点为中心，构造大小为5×5的一个几何结构窗；
[0184] 第3步，判断几何结构窗内是否存在其它素描线的端点，若存在，执行第4步；否则，执行第6步；
[0185] 第4步，判断是否只存在一个端点，若是，将该端点所在素描线和当前素描线进行连接；否则，执行第5步；
[0186] 第5步，连接所选取素描线与各端点所在的素描线，从所有连接线中选取其中夹角最大的两条素描线作为连接完成的素描线；
[0187] 第6步，判断素描线的另一个端点的几何结构窗内是否存在其他素描线的端点，若存在，执行第4步；否则，执行第7步；
[0188] 第7步，对完成连接操作的素描线，选取包含两条及两条以上素描线段的素描线，统计所选取素描线中包含素描线段的条数n，其中n≥2；
[0189] 第8步，判断素描线的条数n是否等于2，若是，则执行第9步；否则，执行第10步；
[0190] 第9步，将素描线顶点的角度值在[10°,140°]的范围内的素描线作为具有聚拢特征的素描线；
[0191] 第10步，选出素描线对应的n-1个顶点的角度值都在[10°,140°]范围内的素描线；
[0192] 第11步，在所选出的素描线中，定义如下两种情况：
[0193] 第一种情况，判断第i条素描线段相邻两条素描线段第i-1条、第i+1条是否在第i条素描线段所在直线的同一侧，2≤i≤n-1，若素描线上的所有素描线段与相邻线段都在同一侧，则标记该素描线为具有聚拢特征的素描线；
[0194] 第二种情况，判断第i条素描线段相邻两条素描线段第i-1条、第i+1条是否在第i条素描线段所在直线的同一侧，2≤i≤n-1，若素描线上有n-1条素描线段与相邻线段在同一侧，而有一条素描线段与其相邻线段在非同一侧，也标记该素描线为具有聚拢特征的素描线；
[0195] 第12步，在具有聚拢特征的素描线中任选一条素描线，由所选取素描线的两个端点坐标，确定两个端点间的距离，若该端点距离在[0,20]范围内，则将所选取素描线作为表示独立目标的素描线；
[0196] 第13步，判断未处理的具有聚拢特征的素描线是否全部选完，若是，执行第14步；否则，执行第12步；
[0197] 第14步，用超像素分割的方法，对合成孔径雷达SAR图像中表示独立目标的素描线周围的像素进行超像素分割，将分割后超像素的灰度值在[0,45]或[180,255]内的超像素作为独立目标超像素；
[0198] 第15步，合并独立目标超像素，将合并后的独立目标超像素的边界作为独立目标的边界，得到独立目标的分割结果。
[0199] 对线目标和独立目标分割的结果进行合并，得到结构像素子空间的分割结果。
[0200] 步骤8，分割匀质区域像素子空间。
[0201] 按照以下步骤，采用基于多项式逻辑回归先验模型的匀质区域分割方法，对匀质区域像素子空间进行分割，得到匀质区域像素子空间的分割结果：
[0202] 第1步，从匀质区域像素子空间中任意选取一个像素点，以所选取的像素点为中心建立3×3的方形窗口，计算该窗口的标准差σ1；
[0203] 第2步，将方形窗口的边长增加2，得到新的方形窗口，计算新方形窗口的标准差σ2；
[0204] 第3步，设标准差阈值T3＝3，如果|σ1-σ2|＜T3，则将标准差为σ2的方形窗口作为最终的方形窗口，执行第4步；否则，执行第2步；
[0205] 第4步，按照下式，计算方形窗口内中心像素的先验概率：
[0206]
[0207] 其中，p1'表示方形窗口内中心像素的先验概率，exp(·)表示指数函数操作，η'表示概率模型参数，η'取值为1，xk′′表示方形窗口内属于第k'类的像素个数，k'∈[1,...,i',...,K']，K'表示分割的类别数，K'取值为5，xi′′表示第3步得到的方形窗口内属于第i'类的像素个数；
[0208] 第5步，将像素灰度的概率密度与纹理的概率密度相乘，得到似然概率p'2，其中，灰度的概率密度由衰落信道Nakagami分布得到，纹理的概率密度由t分布得到；
[0209] 第6步，将先验概率p1'与似然概率p2'相乘，得到后验概率p12'；
[0210] 第7步，判断匀质区域像素子空间中是否还有未处理的像素点，若有，执行第1步；否则，执行第8步；
[0211] 第8步，根据最大后验概率准则，得到匀质区域像素子空间的分割结果。
[0212] 步骤9，合并分割结果。
[0213] 将混合聚集结构地物像素子空间、匀质区域像素子空间和结构像素子空间的分割结果合并，得到合成孔径雷达SAR图像的最终分割结果。
[0214] 下面结合仿真图对本发明的效果做进一步的描述。
[0215] 1.仿真条件：
[0216] 本发明仿真的硬件条件为：智能感知与图像理解实验室图形工作站；本发明仿真所使用的合成孔径雷达SAR图像为：Ku波段分辨率为1米的Piperiver图。
[0217] 2.仿真内容：
[0218] 本发明的仿真实验是对SAR图像中的Piperiver图进行分割，如图2(a)所示的Piperiver图。该图来源于Ku波段分辨率为1米的合成孔径雷达SAR图像。
[0219] 采用本发明的SAR图像素描化步骤，对图2(a)所示的Piperiver图素描化，得到如图2(b)所示的素描图。
[0220] 采用本发明的划分像素子空间步骤，对图2(b)所示的素描图区域化，得到如图2(c)所示的区域图。图2(c)中的白色空间表示混合聚集结构语义空间，其他的为匀质结构语义空间和结构语义空间。将图2(c)所示区域图映射到图2(a)所示原SAR图像，得到如图2(d)所示的SAR图像混合聚集结构地物像素子空间图。
[0221] 采用本发明的基于素描线的聚拢特征分割独立目标步骤，提取图2(b)中非聚集区域的素描线，得到图3(a)，提取图3(a)中的表示线目标的素描线，得到图3(b)所示的结果。图3(b)中的黑色素描线为表示线目标的素描线。对图3(b)中不表示线目标的素描线，提取具有聚拢特征的素描线，得到图3(c)所示的结果，其中黑色素描线表示独立目标。对合成孔径雷达SAR图像，求取表示疑似独立目标素描线周围的超像素，得到图3(d)所示的结果。将分割后超像素的灰度值在[0,45]或[180,255]内的超像素作为独立目标超像素，并且合并独立目标超像素，将合并后的独立目标超像素的边界作为独立目标的边界，得到的独立目标分割结果如图3(e)所示。
[0222] 采用本发明的分割SAR图像混合聚集结构地物像素子空间步骤，对图2(d)所示的Piperiver图的混合聚集结构地物像素子空间图进行分割，得到图4(a)所示的混合地物像素子空间分割结果图，其中灰色区域表示未处理的地物空间，其余相同颜色的区域表示同一种地物，不同颜色的区域表示不同的地物。
[0223] 采用本发明的合并分割结果步骤，合并图4(a)所示的混合聚集结构地物像素子空间分割结果和匀质区域像素子空间分割结果和结构像素子空间分割结果，得到图4(b)，图4(b)是图2(a)Piperiver图像的最终分割结果,图4(c)是基于层次视觉语义和自适应邻域多项式隐模型的SAR图像分割方法对Piperiver图像的最终分割结果图。
[0224] 3.仿真效果分析：
[0225] 图4(b)是本发明方法对Piperiver图像的最终分割结果图，图4(c)是基于层次视觉语义和自适应邻域多项式隐模型的SAR图像分割方法对Piperiver图像的最终分割结果图，通过分割结果图的对比，可得出结论，本发明方法对于混合聚集结构地物像素子空间的边界确定更精确，对于匀质区域像素子空间的分割，区域一致性明显较好，类别数更加合理，且对结构像素子空间中的独立目标进行了较好的分割处理。使用本发明方法对合成孔径雷达SAR图像进行分割，有效的将SAR图像进行了分割，并提高了SAR图像分割的准确性。
