买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：西安电子科技大学

摘要：本发明公开了一种基于图和Wishart距离的极化SAR图像分类方法，主要解决传统有监督分类方法中杂点过多，地物类别模糊的问题。其实现步骤是：1.读取极化SAR数据协方差矩阵并对其进行滤波；2.用滤波后的协方差矩阵构建权值图,并将其按权值递增排序；3.按序遍历权值图，判断每条权值边的权值连接的两个区域是否合并，若合并则更新两个区域的标签号，否则不变，遍历完后得到过分割标签；4.合并该过分割标签中像素点过小的区域，得到新的标签；5.将新标签中的每个区域用有监督分类法标注为标准图中与其最相似的那一类区域的标签号，得到最终分类标签。本发明抗干扰强，地物类别模糊小，分类精度高，可用于极化SAR目标识别。

主权项：一种基于图和Wishart距离的极化SAR图像分类方法，包括如下步骤：1读取极化SAR数据协方差矩阵m×n为像素个数，对协方差矩阵C0进行精致Lee滤波处理得到滤波后的协方差矩阵CC1,C2,...,Cm×n；2根据协方差矩阵CC1,C2,...,Cm×n构建权值图GV,E：2a将极化SAR图像的每个像素点都作为权值图中的一个顶点，设vi为第i个顶点，vj为顶点vi八邻域中的一个顶点，计算顶点vi和vj的不相似度evi,vj：evi,vj=retrCiCj-1+retrCi-1Cj2---1]]其中，i∈1,m×n，Ci为vi的协方差矩阵，Cj为顶点vj的协方差矩阵，是对Ci的求逆，是对Cj的求逆，tr是指矩阵的迹，re·表示取实部；2b利用步骤式1，计算顶点集Vv1,v2,...,vm×n中的所有顶点与其对应的八邻域的不相似度即图的边集Ee1,e2,...,ew，w为权值图中权值边的个数；2c用顶点集V中的各顶点和连接各顶点与其对应的八邻域中的每个顶点之间的权值边集合E组成全连图即权值图GV,E；3根据权值图GV,E对极化SAR数据进行过分割处理，得到过分割标签label0：3a对权值图GV,E按照权值边Ee1,e2,...,ew中元素递增的顺序排序；3b定义初始分割标签S0，S0中每个顶点分别为一类；3c第p个权值边为p＝1；3d如果Sp‑1分割结果中第p个权值边连接的两个顶点vi、vj不在同一个块中，令顶点vi、vj在Sp‑1中对应的标签号分别为I和J，令Sp‑1中标签号为I、J的区域为DI、DJ，计算DI、DJ的内部不相似度MDI,DJ：MDI,DJ=mindifDI+kcI,difDJ+kcJ---2]]其中，difDI为区域DI的最小生成树中的最大值，difDJ为区域DJ的最小生成树中的最大值，cI为区域DI中的像素个数，cJ为区域DJ中的像素个数，k为大于零的常数；3e判断epMDI,DJ是否成立，其中ep表示第p个权值边的权值，若成立，则将分割标签Sp‑1中DI、DJ这两个区域的标签号I、J中的较小值作为新的分割标签Sp中这两个区域的标签号，即合并区域DI、DJ，否则，Sp＝Sp‑1；3f令p＝p+1，返回到步骤3d，直到p＝w时停止，其中w为权值边的个数，得到分割结果label0＝Sw；4设定一个后处理阈值:m＝20，统计分割结果label0中每一个标签的像素点个数，如果该标签的像素点个数小于m，则将该标签的像素点与它邻域中不相似度最小的块进行合并，得到新的标签图label1；5将标签label1中的每个标签号对应的区域作为一个分类单元，用Wishart有监督分类方法得到label1中每个区域在标准图中对应的标签号，即最终的分类结果；步骤4所述的将该标签的像素点与它邻域中不相似度最小的块进行合并，其实现如下：4a统计分割结果label0中标签类数L以及每个标签为l的区域Dl的像素点个数nl，其中l＝1,2,…,L；4b令l＝1；4c将区域Dl中的像素点个数nl与阈值m进行比较，如果nlm，则统计标签号为l的区域Dl的邻域个数zl以及各邻域的标签号oc，其中c＝1,2,…,zl；4d根据邻域个数zl来判断Dl是否只有一个邻域，若zl＝1，则表示Dl只有一个邻域，并将该邻域的标签号oc作为Dl的新标签号；若zl1，则表示Dl有多个邻域，并计算Dl与每个邻域Doc的不相似度WDl,Doc：WDl,Doc=|trClCoc-1+trCl-1Coc|---3]]其中oc为区域Dl的第c个邻域的标签号，c＝1,2,…,zl，为区域Dl中所有像素点的协方差矩阵的平均值，为区域Doc中所有像素点的协方差矩阵的平均值，是的逆矩阵，是的逆矩阵，tr·是指矩阵的迹，|·|表示取绝对值，4e假设区域Dl与其第c0个邻域的不相似度WDl,Doc最小，则将标签号oc0作为区域Dl中像素点的新标签号，即合并区域Dl和区域4f令l＝l+1，返回到步骤4c，直到l＝L时停止，得到新的标签label1。

全文数据：基于图和Wishart距离的极化SAR图像分类方法技术领域[0001]本发明属于图像处理领域，特别涉及一种涉及极化SAR图像分类的方法，可应用于目标识别。背景技术[0002]SAR是一种全天时、全天候的高分辨率雷达体制，应用于军事，导航，农业，地理监视等众多领域。与光学成像系统和其它遥感成像系统相比，SAR被广泛的应用在军事目标识别中的目标检测领域，而相比单极化SAR来说，极化SAR通过测量每个分辨单元的不同的发收组合的散射特征，更加完备的记录了目标的后向散射信息，为更好的分析目标的散射特性提供依据。在过去的二十多年的科学研究里，很多学者证明了极化SAR图像分类是图像解译中的一个非常重要的步骤，分类结果图可作为一个中间步骤为图像的边缘提取、目标检测与识别提供帮助，或者其本身就可以作为用户需求的最终结果。对极化SAR分类的研究于探索目标散射特性以及提高极化SAR系统应用水平有着举足轻重的作用。[0003]因为极化SAR图像成像的特殊性，使得极化SAR图像不能使用自然图像中的特征如灰度值，RGB值等。极化SAR图像含有大量的相干斑噪声，对一般的分类方法来说，噪声敏感度很强。因此对于极化SAR而言，使用到的特征多数为散射矩阵及其协方差矩阵。现有的极化SAR分类方法基本都以像素为分类单元，现有的分类方法涉及基于散射机理和极化目标分解的极化SAR图像分类，也有运用现存的数学方法如支持向量机，马尔科夫随机场和基于统计分布的分类方法。[0004]极化SAR分类中用到的特征，既有图像信息也有极化信息，图像信息为灰度共生矩阵等，而对于极化SAR来说，更重要的、使用更广泛的信息为极化信息，极化信息利用的方式也有很多种，大致可分为三种，一是直接利用散射矩阵以及散射矢量，二是利用散射矩阵的协方差矩阵或者相干矩阵，三是利用极化SAR的极化分解信息或者其它的极化特征。[0005]现有的多数分类方法皆是以像素为分辨单元的方法，但是由于相干斑的影响，即使在去噪的情况下，仍因为其高敏感度而不能避免干扰，因此，有学者开始尝试先将极化SAR数据初始化成像素块，而后以像素块代替像素来进行分割、分类。[0006]几种经典的极化SAR分类算法有：[0007]1.非监督分类经典方法之一的Η-α方法，这种方法通过计算极化SAR数据的每个点的散射熵和散射角，根据每个散射熵与散射角的大小分别将数据分为3类，组合后的图像被分为9类，由于Η-α方法对所有的数据都分为9类，使得类别数存在误差，且由于相干斑噪声的无法抑制，易出现地物类别模糊问题。[0008]2·20世纪末Freeman、Durden等人提出了三分量分解也被称为Freeman-Durden分解，三分量分解是将散射矩阵的协方差矩阵分解成表面散射、偶次散射、体积散射。由于这个三分量模型不适合城区地物，Yamaguchi等人在这个算法上又作了改进，提出了四分量模型，在原有的三分量模型基础上添加了螺旋散射分量，也对Freeman-Durden分解中的体积散射做了改进，然而改进后的算法对于地物非常类似的同构地区区分度仍旧不高。[0009]21世纪初至现今，尽管很多学者和研究者们对这些无监督方法进行了改进，还引用了一些算法来对极化SAR图像进行分类，如支持向量机SVM，神经网络，小波分析，统计，以及分形等，但由于相干斑噪声的干扰使得这些算法的分类精度无法进一步提高。发明内容[0010]本发明的目的在于克服上述已有技术的不足，提出一种基于图和Wishart距离的极化SAR图像分类方法，以减小相干斑噪声的干扰，改善地物类别模糊，提高整体分类精度。[0011]实现本发明目的的技术思路是：用极化SAR数据的协方差矩阵特征代替现有的图方法中的灰度值特征，将极化SAR数据过分割成均匀像素块的集合，避免了相干斑噪声的干扰，以像素块作为Wishart有监督方法的分类单元，提高分类的精度。其技术步骤包括如下：[0012]1读取极化SAR数据协方差矩阵，mXn为像素个数，对协方差矩阵C*3进行精致Lee滤波处理得到滤波后的协方差矩阵：&，&，...，CmXn;[0013]2根据协方差矩阵：&，：2，...，^，计算顶点集VV1，v2,...，vmXn中的所有顶点与其对应的八邻域的不相似度即图的边集Eei，e2,...，ew，w为权值图中权值边的个数；[0018]2c用顶点集V中的各顶点和连接各顶点与其对应的八邻域中的每个顶点之间的权值边集合E组成全连图即权值图GV，E;[0019]3根据权值图GV，E对极化SAR数据进行过分割处理，得到过分割标签labelo:[0020]3a对权值图GV，E按照权值边集Eei，e2，...，ew中元素递增的顺序排序；[0021]3b定义初始分割标签S^S*3中每个顶点分别为一类；[0022]3c，第p个权值边为p=1;[0023]3d如果F1分割结果中第p个权值边连接的两个顶点Vl、Vj不在同一个块中，令顶点Vi、Vj在Sp-1中对应的标签号分别为I和J，令Sp-1中标签号为I、J的区域为D^Dj，计算D^Dj的内部不相似度M^^Dj:[0025]其中，difΦύ为区域D:的最小生成树中的最大值，difDj为区域Dj的最小生成树中的最大值，CI为区域出中的像素个数，cj为区域Dj中的像素个数，k为大于零的常数；[0026]3e判断eP〈MD:，Dj是否成立，其中eP表示第p个权值边的权值，若成立，则将分割标签f1中D^Dj这两个区域的标签号I、J中的较小值作为新的分割标签Sp中这两个区域的标签号，即合并区域D:、Dj，否则，Sp=f1;[0027]3f令p=p+l，返回到步骤3d，直到p=w时停止，其中w为权值边的个数，得到分割结果labelo=Sw;[0028]4设定一个后处理阈值:m=20,统计分割结果labelo中每一个标签的像素点个数，如果该标签的像素点个数小于m，则将该标签的像素点与它邻域中不相似度最小的块进行合并，得到新的标签图labeli。[0029]5将标签label!中的每个标签号对应的区域作为一个分类单元，用Wishart有监督分类方法得到labeh*每个区域在标准图中对应的标签号，即最终的分类结果。[0030]本发明与现有的技术相比具有以下优点：[0031]1.本发明采取先将极化SAR数据过分割成许多均匀的小区域，然后以小区域为分辨单元分类，避免了相干斑噪声的干扰，改善了传统分类结果中杂点过多的情况；[0032]2.本发明中对初步过分割结果中像素点个数过小的区域进行合并，减少了分类时的类别t吴糊；[0033]3.将现有图方法中的不相似度计算公式变为修正的Wishart距离，更能适应极化SAR数据，使得分类结果精确度明显提高。附图说明[0034]图1是本发明的实现流程不意图；[0035]图2是本发明在AIRSARFlevoland农田数据上的实验结果图；[0036]具体实现方式[0037]以下结合附图对本发明的方案和效果作进一步详细描述。[0038]参照图1，本发明的具体实现步骤如下：[0039]步骤一、读取极化SAR数据协方差矩阵，mXη为像素个数，对协方差矩阵C用7X7的窗口进行精致Lee滤波处理，得到滤波后的协方差矩阵...，CmXn；[0040]步骤二、根据协方差矩阵C&，C2,...，CmXn构建权值图GV，E。[0041]2a将极化SAR图像的每个像素点都作为权值图中的一个顶点，设Vi为第i个顶点，Vj为顶点Vi八邻域Λζ8中的一个顶点，计算顶点Vi和Vj的不相似度eVi，Vj:^为Vi的协方差矩阵，Cj为顶点Vj的协方差矩阵，Cf1是对匕的求逆，C；1是对Q的求逆，tr是指矩阵的迹，re·表示取实部；[0044]2b利用步骤式〈1，计算顶点集VV1，v2，...，vmXn中的所有顶点与其对应的八邻域的不相似度即图的边集Eei，e2,...，ew，w为权值图中权值边的个数；[0045]2c用顶点集V中的各顶点和连接各顶点与其对应的八邻域中的每个顶点之间的权值边集合E，组成全连图即权值图GV，E。[0046]步骤三、根据权值图GV，E对极化SAR数据进行过分割处理，得到过分割标签labelo〇[0047]3a对权值图GV，E按照权值边集Eei，e2，...，ew中元素递增的顺序排序；[0048]3b定义初始分割标签S'S*3中每个顶点分别为一类；[0049]3c设第p个权值边为p=1;[0050]3d判断第p个权值边连接的两个顶点在分割标签SM中对应的两个区域是否在同一个块中：如果在同一个块中，则跳到步骤3f;如果不在同一个块中，则令顶点vi、vj在Sp1，则表示1区域Di有多个邻域，并计算1区域Di与每个邻域D。c的不相似度WDi，D。c:[0061]其中oc为1区域Di的第c个邻域的标签号，c=l，2,…，ζ1，ζ为1区域Di中所有像素点的协方差矩阵的平均值，为邻域D^co中所有像素点的协方差矩阵的平均值，是G的逆矩阵，^&是^^的逆矩阵，tr·是指矩阵的迹，Μ·||表示取绝对值；[0062]4e假设1区域Di与其第c〇个邻域化⑷的不相似度WQhDofc最小，则将标签号〇co作为1区域D冲像素点的新标签号，即合并区域Di和区域化,。;[0063]4f令1=1+1，返回到步骤4c，直到1=L时停止，得到新的标签labelu[0064]步骤五、将标签labeli中的每个标签号对应的区域作为一个分类单元，用Wishart有监督分类方法得到labeh*每个区域在标准图中对应的标签号，即最终的分类结果。[0065]5a统计分割标签labeh*的标签类别数Q，新的标签labeh中标签号为q的区域表示为Dq，其中q=l，2，-_，Q，计算新的标签labeh*的区域Dq中所有像素点的协方差矩阵平均值0，统计标准图中标签类别数T，标准图中标签号为t的区域表示为Dt，其中t=l，2,…，T，计算标准图中区域Dt中所有像素点的协方差矩阵的平均值巧；[0066]5b定义分类标签为label，并将所有像素点的标签号初始化为0,令q=l;[0067]5c根据步骤5a中得到的所述协方差矩阵平均值0和G，计算新的标签labell中标签号为q的区域Dq与标准图中的每个区域Dt的不相似度WDq，Dt:[0069]其中，tr为矩阵的迹，t=l，2,…，Τ，（《Γ1表示矩阵取逆，||·||表示取绝对值；[0070]5d比较新的标签labeh中标签号为q的区域Dq与标准图中的每个区域Dt的不相似度WDq，Dt之间的大小，t=1，2，…，T，将值最小的不相似度WDq，所对应标准图中区域A。的标签号to作为区域Dq中所有像素点在分类标签label中的标签号；[0071]5e令q=q+l，返回到步骤5c，当q=Q时停止，得到Q个区域每个区域的分类标签号，即最终分类标签label。[0072]本发明的效果通过以下实验可以进一步说明[0073]1.仿真内容[0074]用本方法和传统的以像素点为分类单元的Wishart有监督方法分别对AIRSARFlevoland农田极化SAR数据进行分割分类，农田数据的大小为380X420像素，实验结果如图2所示，其中：[0075]图2a表示农田极化SAR数据Pauli分解伪彩图；[0076]图2b为分类结果标准图；[0077]图2c表示的是用传统的Wishart有监督方法对农田极化SAR数据分类的结果图；[0078]图2d是对2c所示分类结果去除非测试地区后的示意图；[0079]图2e表示用本发明方法对农田极化SAR数据初步过分割的结果示意图；[0080]图2f表示用本发明方法对初步过分割结果中像素点个数较小的区域合并后的示意图；[0081]图2g表示用本发明方法分类的结果示意图；[0082]图2h是对图2g所示分类结果去除非测试区域后的示意图；[0083]2.仿真结果分析[0084]从图2e和图2f可以看出，本发明合并小区域后使得区域与区域之间的边缘更清晰。[0085]从图2c和图2g可以看出，传统的分类结果中，由于相干斑噪声的影响，每一类中都存在很多杂点，且两类地物的分界处出现地物类别模糊的问题，[0086]可见，本发明很好的避免了分类中相干斑噪声的干扰，改善了类别的模糊，精度明显提尚。

权利要求：1.一种基于图和Wishart距离的极化SAR图像分类方法，包括如下步骤：1读取极化SAR数据协方差矩阵，mXn为像素个数，对协方差矩阵C*3进行精致Lee滤波处理得到滤波后的协方差矩阵C&，C2,...，CmXn;2根据协方差矩阵C&，C2，...，CmXn构建权值图GV，E:2a将极化SAR图像的每个像素点都作为权值图中的一个顶点，设vi为第i个顶点，vj为顶点Vi八邻域，计算顶点集Vvi，v2,...，vmXn中的所有顶点与其对应的八邻域的不相似度即图的边集Eei，e2,...，ew，w为权值图中权值边的个数；2c用顶点集V中的各顶点和连接各顶点与其对应的八邻域中的每个顶点之间的权值边集合E组成全连图即权值图GV，E;3根据权值图GV，E对极化SAR数据进行过分割处理，得到过分割标签labelo:3a对权值图GV，E按照权值边Eei，e2，...，ew中元素递增的顺序排序；3b定义初始分割标签中每个顶点分别为一类；3c第p个权值边为p=1;3d如果f1分割结果中第p个权值边连接的两个顶点Vl、Vj不在同一个块中，令顶点Vl、Vj在S1"1中对应的标签号分别为I和J，令S1"1中标签号为I、J的区域为D^Dj，计算D^Dj的内部不相似度MDi，Dj:其中，difD:为区域D:的最小生成树中的最大值，difDj为区域Dj的最小生成树中的最大值，CI为区域出中的像素个数，cj为区域Dj中的像素个数，k为大于零的常数；3e判断eP〈MD:，Dj是否成立，其中eP表示第p个权值边的权值，若成立，则将分割标签f1中D^Dj这两个区域的标签号I、J中的较小值作为新的分割标签Sp中这两个区域的标签号，即合并区域D^Dj，否则，SP=SP、3f令p=p+l，返回到步骤3d，直到p=w时停止，其中w为权值边的个数，得到分割结果labelo=Sw；4设定一个后处理阈值:m=20,统计分割结果labelo中每一个标签的像素点个数，如果该标签的像素点个数小于m，则将该标签的像素点与它邻域中不相似度最小的块进行合并，得到新的标签图labeli;5将标签labeli中的每个标签号对应的区域作为一个分类单元，用Wishart有监督分类方法得到labeh*每个区域在标准图中对应的标签号，即最终的分类结果；步骤4所述的将该标签的像素点与它邻域中不相似度最小的块进行合并，其实现如下：4a统计分割结果labelo中标签类数L以及每个标签为1的区域Di的像素点个数n1，其中1=1，2,…，L;4b令1=1;4c将区域0冲的像素点个数n⑴与阈值m进行比较，如果n⑴1，则表示Di有多个邻域，并计算口丄与每个邻域Dofe的不相似度WDi：其中〇c为区域Di的第c个邻域的标签号，c=1，2，…，z1，C:为区域Di中所有像素点的协方差矩阵的平均值，为区域Dw。中所有像素点的协方差矩阵的平均值，是G的逆矩阵是^^的逆矩阵，tr〇是指矩阵的迹，|•|表示取绝对值，4e假设区域口:与其第c〇个邻域的不相似度WDl，DciW最小，则将标签号〇CQ作为区域Di中像素点的新标签号，即合并区域m和区域;4f令1=1+1，返回到步骤4c，直到1=L时停止，得到新的标签labelu2.根据权利要求1所述的方法，其中步骤5所述的用Wishart有监督分类方法得到labeh*每个区域在标准图中对应的标签号，按如下步骤进行：5a统计分割标签labeh*的标签类别数Q，标签labeh中标签号为q的区域表示为Dq，其中q=l，2,…，Q，计算label:*的区域Dq中所有像素点的协方差矩阵的平均值句，统计标准图中标签类别数T，标准图中标签号为t的区域表示为Dt，其中t=l，2,…，T，计算标准图中区域Dt中所有像素点的协方差矩阵的平均值巧.;5b令q=l;5c根据步骤5a中得到label:*的区域Dq中所有像素点的协方差矩阵的平均值泛和标准图中的区域Dt中所有像素点的协方差矩阵的平均值计算label:*的区域Dq与标准图中的每个区域Dt的不相似度WDq，Dt:其中，tr为矩阵的迹，t=l，2,…，T，。）^1表示矩阵取逆，|•|表示取绝对值；5d假设label冲的区域Dq与标准图中的第to个区域马。的不相似度f最小，则将to作为labeh*的区域Dq的最终标签号；5e令q=q+l，返回到步骤5c，当q=Q时停止，得到最终的分类标签label。

百度查询： 西安电子科技大学 基于图和Wishart距离的极化SAR图像分类方法