基于随机相似圆图的全极化合成孔径雷达图像分类方法
【技术领域】
[0001] 本发明设及图像分类领域，特别设及一种基于随机相似圆图的全极化合成孔径雷 达图像分类方法。
【背景技术】
[0002] 传统单极化雷达测量无法全面获取目标散射信息，因此基于单极化测量无法实现 对复杂分布雷达目标的精细理解。全极化散射测量的提出旨在通过在一组正交极化基下交 替发射同时接收的方式获取更多的目标散射信息，但也进一步增加了信号处理难度，因为 此时获得的目标散射不再为一个带幅度和相位的复数（即后向散射系数），而是由六个复数 和Ξ个实数组成的抽象的3 X 3目标相干矩阵。全极化雷达信号处理的核屯、就是要设计算法 从该矩阵中提取有用的、可W理解的目标几何物理信息，实现对目标的有效识别和分类。
[0003] 一个常用的方法就是将目标相干矩阵进行极化分解，运一概念由化ynen于1970年 首次提出（参见参考文献 l:''J.R.Huynen,"Phenomenological Theory of Radar 了日巧613,"陆.0.013361'1日1：；[0]1,了6油.1]]1；[￥.061的，061的，化0化1：1161'1日]1(1,1970")，引领了 该领域40多年的蓬勃发展。此外，基于散射相似性的目标识别和分类方法近期也受到了广 泛应用和关注，我国科学家在此做了重要贡献。清华大学杨健等首先定义了两个单目标散 射矩阵间的相似性（参见参考文献2: "J. Yang, Y.N. Peng, and S. Μ丄in/Similarity between two scattering matrices ,"Electron . Lett.,vol.37,no .3,pp.193-194, Feb. 200Γ );国防科大陈强等将该参数进行拓展，W描述分布目标与标准单目标散射体间 的相似性（参见参考文献3: "Q. Chen, Y.Μ. Jiang,L. J. Zhao,and G. Y.Kuang,"化larimel：!·iC scattering similarity between a random scatterer and a canonical scatterer," IE邸 Geosci .Remote Sens丄ett. ,vol .7,no.4,pp.866-869,Oct.2010");在杨和陈工作的 基础上，中科院空间中屯、李东和张云华等近期进一步发展出一个随机散射相似性参数(参 见参考文南犬 4 : "D . Li and Υ . Zhan 邑，"Random similarity between two mixed scatterers, "IEEE Geosci .Remote Sens 丄ett. , vol. 12,no. 12,pp. 1-5,Dec. 2015")，用从 描述任意两个散射体，特别是两个分布/混合散射体间的相似性。该参数的一个优势在于， 当待描述的两个散射体为同一散射体时，所得到的自随机散射相似性参数能直接刻画该散 射体的随机性，故可实现对目标相似性和随机性的直接描述。
[0004] 在现有技术中使用最广泛、最成功、且最认可的全极化合成孔径雷达图像分类法 为赌/alpha角分类法，但该方法存在错分现象，在准确性上还有提升空间。
【发明内容】
[0005] 本发明的目的在于克服已有的全极化合成孔径雷达图像分类方法存在错分现象 的缺陷，从而提供一种准确性更高的分类方法。
[0006] 为了实现上述目的，本发明提供了一种基于随机相似圆图的全极化合成孔径雷达 图像分类方法，包括：
[0007] 步骤1)、读入待分类图像的目标相干矩阵并进行去取向处理；
[0008] 步骤2)、由步骤1)所得到的结果分别求取目标的自随机相似参数rrrs，W及目标 对标准表面散射的相似性参数rrru、目标对二面散射的相似性参数rrrsd和目标对体散射 的相似性参数rrrcv;
[0009] 步骤3)、基于步骤2)所得到的自随机相似参数rrrs构建随机相似圆图，所述随机 相似圆图具有立层环状结构，从内到外分别代表高、中、低立种散射随机情形;其中，所述随 机相似圆图从外环到内环的划分方法及其物理意义如下：
[0010]
[0011]步骤4)、基于目标对标准表面散射的相似性参数rrru、目标对二面散射的相似性 参数rrrcd和目标对体散射的相似性参数rrrcv的排列组合，进一步将步骤3)所得到的随机 相似圆图划分为十块，分别表示十种不同的散射机制;其中，划分为十块后的随机相似圆图 表不为：
[001引其中，rrrcs > {rrrcd, rrrcv}表示;rrrcs大于rrrcd和rrrcv;随机相似圆图的外环被 划分为Ξ块，对应S、D、V^种散射机制，S表示目标散射"相似于表面散射"、D表示目标散射 "相似于二面散射"、V表示目标散射"相似于体散射"；随机相似圆图的中环被划分为六块， 对应50、5￥、05、0￥、￥5和￥0六种散射机制，50表示目标散射"相似于表面散射和二面散射"， SV表示目标散射"相似于表面散射和体散射"，DS表示目标散射"相似于二面散射和表面散 射"，DV表示目标散射"相似于二面散射和体散射"，VS表示目标散射"相似于体散射和表面 散射"，VD表示目标散射"相似于体散射和二面散射"；随机相似圆图的内环代表高随机性散 射R;
[0016] 步骤5)、用不同的标记将待分类图像根据之前所做的划分加 W标识，形成最终分 类图。
[0017] 上述技术方案中，所述步骤1)具体包括:对待分类图像中每一个像素位置的目标 相干矩阵做去取向操作;其中，
[0018] 若原始的目标相干矩阵τ为：
[0019]
[0020] 取向角9通过下式估计：
[0021]
[0022] 则去取向后的目标相干矩阵Τ'表示为：
[0023] 其中上标"Η"表示矩阵 , 共辆转置操作。
[0024] 上述技术方案中，在步骤2)中，目标对标准表面散射的相似性参数rrrcs、目标对标 准二面散射的相似性参数rrrcd、W及目标对标准体散射的相似性参数rrrcv的计算公式分 别为：
[0025]
[0026] 目标的自随机相似参数rr。的计算公式为：
[0027]
[00%]上述技术方案中，所述步骤5)进一步包括：
[0029] 步骤5-1)、对待分类图像中的某一像素，基于该像素位置的目标T矩阵计算相似性 参数口_。3、1'1'。(1、1'1'。巾^及自随机相似参数1·!·。；
[0030] 步骤5-2)、然后利用随机相似圆图的模型判断该像素位于随机相似圆图的哪一环 W及该环的哪一类；
[0031 ]步骤5-3)、用赋予该类别的区分标记标识该像素位置；
[0032] 步骤5-4)、对待分类图像中的所有像素均做步骤5-1)-步骤5-3)操作，得到最终分 类图。
[0033] 本发明的优点在于：
[0034] 本发明的相似圆图分类方法有效地利用随机散射相似性参数在相似性和随机性 刻画方面的双重优势实现对全极化合成孔径雷达图像的精细分类，该方法实现简单且易于 理解，与当前最流行的分类方法的对比实验表明，其具有优异的雷达目标区分性能。
【附图说明】
[0035] 图1是本发明的基于随机相似圆图的全极化合成孔径雷达图像分类方法的总体流 程图；
[0036] 图2是本发明的基于随机相似圆图的全极化合成孔径雷达图像分类方法的具体流 程图；
[0037] 图3是一个实施例中所采用的待分类全极化合成孔径雷达图像的示意图；
[0038] 图4(a)是本发明的方法所实现的随机相似圆图模型的示意图；
[0039] 图4(b)是现有技术中的赌/alpha角分类模型的示意图；
[0040] 图5是图3经本发明的方法做分类后的结果示意图；
[0041] 图6是图3采用现有技术中的赌/alpha角分类方法做分类后的结果示意图。
【具体实施方式】
[0042] 现结合附图对本发明作进一步的描述。
[0043] 参考图1和图2,本发明的基于随机相似圆图的全极化合成孔径雷达图像分类方法 包括W下步骤：
[0044] 步骤1)、读入待分类图像的目标相干矩阵并进行去取向处理；
[0045] 步骤2)、由步骤1)所得到的结果分别求取目标的自随机相似参数rrrs，W及目标 对标准表面散射、二面散射和体散射的相似性参数rrres、rrred和rrrcv;
[0046] 步骤3)、基于步骤2)所得到的rrrs构建随机相似圆图，其具有Ξ层环状结构，从内 到外分别代表高、中、低Ξ种散射随机情形；
[0047] 步骤4)、基于参数；TTTcsjrrcd和rrrcv的排列组合，进一步将步骤3)所得到的相似 圆图划分为十块，分别表示十种不同的散射机制；
[0048] 步骤5)、用不同的标记将待分类图像根据步骤4)所做的划分加 W标识，形成最终 分类图。
[0049] 下面对本发明方法中的步骤做进一步描述。
[0050] 在步骤1)中，读入待分类的全极化合成孔径雷达图像。在一个实施例中，所读入的 图像如图3所示，该图像由NASA/肝L研制的机载AIRSAR系统获取于美国旧金山地区，该地区 的地面验证信息易于获得，因此经常被用于各种分解和分类算法的有效性验证中。
[0051 ]在读入待分类的全极化合成孔径雷达图像后，对图像每一个像素位置的目标相干 矩阵(描述了对应像素的散射信息)进行去取向操作。若原始目标相干矩阵T为：
[0化2] U)
[0化3] 则取向角^^可通过下式估计：
[0055]去取向后的目标相干矩阵Τ'可表示为：
[0化4] α):
[0化6]
(3J
[0057] 上式中，上标"Η"表示矩阵共辆转置操作。目标若在垂直于雷达视线的平面内存在 非零倾角，将会使同极化散射能量分解至交叉极化上，引起交叉极化散射能量增加，最终可 导致建筑物被错识别为森林。经过本步骤中的上述去取向操作，目标物理取向产生的散射 贡献得W补偿，保证了后续分类结果不受影响。
[0058] 本发明的分类方法是基于随机散射相似性参数实现的，对于步骤1)所得到的去取 向后的分布目标相干矩阵Τ'，在步骤2)中，其与任意一个标准分布雷达目标相干矩阵Tc间 的随机散射相似性参数rrr定义为：
[0化9] (4)
[0060]其中，Tr( ·)表示矩阵的迹，即对矩阵对角线元素做求和操作。若T。退化为一个单 目标，那么Tc可表示为：
[006。
(5)
[0062] 其中，k。称之为单目标T。的化uli矢量。将式巧)带入式(4)有：
[0063] (6)
[0064] 其中，I I . I U表示二范数。常见的标准表面散射、标准二面散射和标准体散射皆 存在单目标散射形式，其化U1 i矢量kts、ktd、ktv分别为：
[00 化] (7)
[0066] 将式(7)带入式(5)，我们即可得到