一种面向全极化合成孔径雷达图像的散射金字塔分类方法技术

本发明专利技术涉及一种面向全极化合成孔径雷达图像的散射金字塔分类方法，包括：读入待分类图像的相干矩阵并进行去取向处理；分别求取目标对表面散射、二面散射和体散射的偏好度参数SDoPs、SDoPd和SDoPv以及平均偏好度参数SDoP3；基于SDoP3构建三层金字塔模型，该三层金字塔模型从底到顶分别代表高、中、低三种散射随机情形；基于参数SDoPs、SDoPd和SDoPv的排列组合，将三层金字塔模型划分为十块，分别表示十种不同的散射机制；用不同的标记将待分类图像加以标识，形成最终分类图。

【技术实现步骤摘要】
一种面向全极化合成孔径雷达图像的散射金字塔分类方法
本专利技术涉及图像分类方法，特别涉及一种面向全极化合成孔径雷达图像的散射金字塔分类方法。
技术介绍
全极化合成孔径雷达图像分类是微波遥感和地球科学领域的研究热点。全极化雷达通过交替正交接收可获得关于目标物理散射和几何结构等信息，这些信息包含在获取的目标散射二阶统计矩阵，诸如Kennaugh矩阵、相干矩阵和协方差矩阵当中。要想从这些矩阵中解译出目标信息，需对其进行一定的处理。Huynen在1970年发现(参见参考文献1：“J.R.Huynen,"PhenomenologicalTheoryofRadarTargets,"Ph.D.Dissertation,Tech.Univ.Delft,Delft,TheNetherland,1970”)：目标分解是一种有效地理解目标散射信息的方法，其通过从目标中提取出主要散射信息或将其在不同标准散射机制上进行展开，可在一定程度上将特定目标散射信息揭露无遗。Huynen分解是第一个正式提出的极化目标分解技术，其奠定了该领域的基础。该领域在过去四十多年里受到了持续的广泛关注，大批不同极化分解方法被提出，并已广泛用于机载和星载极化雷达图像处理领域。这其中使用最广泛、最成功的莫属Cloude-Pottier分解，基于该分解发展的熵/alpha角分类在极化雷达图像处理领域具有广泛的影响力(参见参考文献2：“S.R.CloudeandE.Pottier,"AnentropybasedclassificationschemeforlandapplicationsofpolarimetricSAR,"IEEETrans.Geosci.RemoteSens.,vol.35,no.1,pp.68-78,Jan.1997”)。尽管Huynen分解具有重要理论意义，但并没有得到广泛的实际应用，一个主要原因是其只能处理规则/对称目标。而实际上许多自然目标和人造目标并非如此理想，具有复杂的非规则性和非对称性，例如茂密的森林和拥挤的城市建筑群。对于这类目标，Huynen分解将无能为力。标准Huynen分解的提出意欲克服Huynen分解在这方面的不足。从目标散射描述的物理可实现条件出发，Li和Zhang近期发现仍存在另外两种分解形式，分别偏好于散射非对称性和非规则性，其与原始Huynen分解一起构成了标准Huynen分解，可提供关于目标散射的全面描述，这在自适应雷达目标提取方面已得到了成功应用(参见参考文献3：“D.LiandY.Zhang,"Huynendichotomy-basedradartargetadaptiveextraction,"inProc.PIERS,Stockholm,Sweden,2013,pp.881-85”)。然而关于该分解的分类应用，尚没有一个有效的实现方案。分类应用考查的是一个分解算法能否用于对随机分布目标的描述，这比自适应目标提取更具挑战性，因为后者通常只需处理单目标即可。按照当前极化分解领域的发展趋势，一个优秀的分解算法不仅须能提取单目标，更须对广泛存在的分布目标具有适用性。因此标准Huynen分解能否有效用于目标分类，是影响其未来应用的关键指标。
技术实现思路
本专利技术的目的在于克服现有技术中尚缺乏对Huynen分解的有效实现方案，从而提供一种简单快速、易于实现的图像分类方法。为了实现上述目的，本专利技术提供了一种面向全极化合成孔径雷达图像的散射金字塔分类方法，包括：步骤1)、读入待分类图像的目标相干矩阵并进行去取向处理；步骤2)、由步骤1)所得到的结果分别求取目标对表面散射、二面散射和体散射的偏好度参数SDoPs、SDoPd和SDoPv以及平均偏好度参数SDoP3；步骤3)、基于步骤2)所得到的SDoP3构建三层金字塔模型，该三层金字塔模型从底到顶分别代表高、中、低三种散射随机情形；步骤4)、基于参数SDoPs、SDoPd和SDoPv的排列组合，进一步将步骤3)所得到的三层金字塔模型划分为十块，分别表示十种不同的散射机制；步骤5)、用不同的标记将待分类图像根据步骤4)所做的划分加以标识，形成最终分类图。上述技术方案中，所述步骤1)具体包括：对待分类图像中每一个像素位置的目标相干矩阵做去取向操作；其中，若原始的目标相干矩阵T为：取向角通过下式估计：则去取向后的目标相干矩阵T'表示为：上述技术方案中，在步骤2)中，目标对表面散射的偏好度参数SDoPs、目标对二面散射的偏好度参数SDoPd、目标对体散射的偏好度参数SDoPv的计算公式分别为：平均偏好度参数SDoP3的计算公式为：上述技术方案中，在步骤3)中，基于SDoP3所构建的三层金字塔模型为：上述技术方案中，在步骤4)中，将三层散射金字塔划分为十块包括：其中，金字塔顶层被划分为三块，其中的S表示偏好表面散射、D表示二面散射、V表示体散射；金字塔中层被划分为六块，分别用SD、SV、DS、DV、VS和VD表示，SD表示目标散射“偏好表面散射和二面散射”，SV表示目标散射“偏好表面散射和体散射”，DS表示目标散射“偏好二面散射和表面散射”，DV表示目标散射“偏好二面散射和体散射”，VS表示目标散射“偏好体散射和表面散射”，VD表示目标散射“偏好体散射和二面散射”。上述技术方案中，所述步骤5)进一步包括：步骤5-1)、对待分类图像中的某一像素基于该像素位置的T矩阵计算参数SDoPs、SDoPd、SDoPv以及它们的平均SDoP3，步骤5-2)、然后利用三层金字塔模型判断其位于金字塔哪一层以及该层的哪一类；步骤5-3)、用赋予该类别的区分标记标识该像素位置；步骤5-4)、对待分类图像中的所有像素均做步骤5-1)-步骤5-3)操作，得到最终分类图。本专利技术的优点在于：本专利技术的散射金字塔分类方法有效地将标准Huynen分解用于目标分类，该方法实现简单且易于理解，与当前最流行的分类方法的对比实验表明，其具有优异的雷达目标区分性能。附图说明图1是本专利技术的面向全极化合成孔径雷达图像的散射金字塔分类方法的总体流程图；图2是本专利技术的面向全极化合成孔径雷达图像的散射金字塔分类方法的具体流程图；图3是一个实施例中所采用的待分类全极化合成孔径雷达图像的示意图；图4(a)是本专利技术的方法所实现的散射金字塔模型的示意图；图4(b)是现有技术中的熵/alpha角分类模型的示意图；图5是图3经本专利技术的方法做分类后的结果示意图；图6是图3采用现有技术中的熵/alpha角分类方法做分类后的结果示意图。具体实施方式现结合附图对本专利技术作进一步的描述。参考图1和图2，本专利技术的面向全极化合成孔径雷达图像的散射金字塔分类方法包括以下步骤：步骤1)、读入待分类图像的相干矩阵并进行去取向处理；步骤2)、由步骤1)所得到的结果分别求取目标对表面散射、二面散射和体散射的偏好度(SDoP，ScatteringDegreeofPreference)参数SDoPs、SDoPd和SDoPv，以及平均偏好度参数SDoP3；步骤3)、基于步骤2)所得到的SDoP3构建三层金字塔模型，从底到顶分别代表高、中、低三种散射随机情形；步骤4)、基于参数SDoPs、SDoPd和SDoPv的排列组合，进一步将步骤3)所得到的三层金字塔模型划分为十本文档来自技高网...

【技术保护点】
一种面向全极化合成孔径雷达图像的散射金字塔分类方法，包括：步骤1)、读入待分类图像的目标相干矩阵并进行去取向处理；步骤2)、由步骤1)所得到的结果分别求取目标对表面散射、二面散射和体散射的偏好度参数SDoPs、SDoPd和SDoPv以及平均偏好度参数SDoP3；步骤3)、基于步骤2)所得到的SDoP3构建三层金字塔模型，该三层金字塔模型从底到顶分别代表高、中、低三种散射随机情形；步骤4)、基于参数SDoPs、SDoPd和SDoPv的排列组合，进一步将步骤3)所得到的三层金字塔模型划分为十块，分别表示十种不同的散射机制；步骤5)、用不同的标记将待分类图像根据步骤4)所做的划分加以标识，形成最终分类图。

【技术特征摘要】
1.一种面向全极化合成孔径雷达图像的散射金字塔分类方法，包括：步骤1)、读入待分类图像的目标相干矩阵并进行去取向处理；步骤2)、由步骤1)所得到的结果分别求取目标对表面散射、二面散射和体散射的偏好度参数SDoPs、SDoPd和SDoPv以及平均偏好度参数SDoP3；步骤3)、基于步骤2)所得到的SDoP3构建三层金字塔模型，该三层金字塔模型从底到顶分别代表高、中、低三种散射随机情形；步骤4)、基于参数SDoPs、SDoPd和SDoPv的排列组合，进一步将步骤3)所得到的三层金字塔模型划分为十块，分别表示十种不同的散射机制；步骤5)、用不同的标记将待分类图像根据步骤4)所做的划分加以标识，形成最终分类图；所述步骤1)具体包括：对待分类图像中每一个像素位置的目标相干矩阵做去取向操作；其中，若原始的目标相干矩阵T为：取向角通过下式估计：则去取向后的目标相干矩阵T'表示为：在步骤2)中，目标对表面散射的偏好度参数SDoPs、目标对二面散射的偏好度参数SDoPd、目标对体散射的偏好度参数SDoPv的计算公式分别为：

【专利技术属性】
技术研发人员：李东，张云华，
申请(专利权)人：中国科学院空间科学与应用研究中心，
类型：发明
国别省市：北京;11