主权项 |
1.一种基于空间相关特征超像素块的SAR图像分割方法,包含以下步骤:(1)预处理待分割SAR图像:(1a)输入一幅待分割SAR图像;(1b)提取当前待分割SAR图像中所有像素点的小波特征描述子;(2)计算超像素块:(2a)采用Turbopixels方法,将当前待分割SAR图像分割成超像素块;(2b)对当前待分割SAR图像的每个超像素块中所有像素点的小波特征描述子求平均值,将该平均值作为每个超像素块的小波特征描述子;(3)初试化加权隶属度矩阵;(4)计算修正隶属度;按照下式,得到关于待分割SAR图像的修正隶属度:<maths num="0001"><![CDATA[<math><mrow><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>u</mi><mi>ij</mi></msub><msub><mi>h</mi><mi>ij</mi></msub></mrow><mrow><msubsup><mi>Σ</mi><mrow><mi>c</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></msubsup><msub><mi>u</mi><mi>cj</mi></msub><msub><mi>h</mi><mi>cj</mi></msub></mrow></mfrac></mrow></math>]]></maths>其中,r<sub>ij</sub>表示第j块超像素属于第i类聚类中心v<sub>i</sub>的修正隶属度,h<sub>ij</sub>表示第j块超像素属于第i类聚类中心v<sub>i</sub>的加权隶属度,u<sub>ij</sub>表示第j块超像素属于第i类聚类中心v<sub>i</sub>的隶属度,∑表示求和操作;(5)计算位置信息隶属度和纹理信息隶属度;(5a)按照下式,计算位置信息隶属度:<maths num="0002"><![CDATA[<math><mrow><msub><mi>e</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></msubsup><msub><mi>r</mi><mi>ik</mi></msub><msub><mi>d</mi><mi>jk</mi></msub></mrow><mrow><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></msubsup><msub><mi>d</mi><mi>jk</mi></msub></mrow></mfrac></mrow></math>]]></maths>其中,e<sub>ij</sub>表示第j块超像素属于第i类聚类中心v<sub>i</sub>的位置信息隶属度,S表示与第j块超像素相关联的超像素块的总数,S=a+b,a表示与第j块超像素小波特征相似的超像素块的个数,b表示与第j块超像素位置相邻的超像素块的个数,∑表示求和操作,r<sub>ik</sub>表示第k块超像素属于第i类聚类中心v<sub>i</sub>的修正隶属度,d<sub>jk</sub>表示第j块超像素中心位置与第k块超像素中心位置之间的高斯距离;(5b)按照下式,计算纹理信息隶属度:<maths num="0003"><![CDATA[<math><mrow><msub><mi>f</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></msubsup><msub><mi>r</mi><mi>ik</mi></msub><msub><mi>X</mi><mi>jk</mi></msub></mrow><mrow><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></msubsup><msub><mi>X</mi><mi>jk</mi></msub></mrow></mfrac></mrow></math>]]></maths>其中,f<sub>ij</sub>表示第j块超像素属于第i类聚类中心v<sub>i</sub>的纹理信息隶属度,S表示与第j块超像素相关联的超像素块的总数,S=a+b,a表示与第j块超像素小波特征相似的超像素块的个数,b表示与第j块超像素位置相邻的超像素块的个数,∑表示求和操作,r<sub>ik</sub>表示第k块超像素属于第i类聚类中心v<sub>i</sub>的修正隶属度,X<sub>jk</sub>表示第j块超像素的小波特征与第k块超像素的小波特征之间的高斯距离;(6)计算隶属度相关量;(7)优化加权因子:采用粒子群优化方法,对加权隶属度中的加权因子α<sub>1</sub>和α<sub>2</sub>进行优化;(8)判断是否满足循环终止条件:判断相邻两次隶属度相关量之间的差值是否大于0.01,如果是,则表示满足循环终止条件,执行步骤(9);否则,执行步骤(4);(9)获得分割图像:(9a)将修正隶属度矩阵中,每列最大隶属度值对应的隶属度矩阵的行数作为超像素块的类别标签;(9b)依次遍历所有的超像素块,得到所有超像素块的类别标签;(9c)选取所有超像素块中类别标签相同的超像素块,对其着相同的颜色,完成SAR图像分割。 |