基于全极化干涉合成孔径雷达数据的建筑物提取方法技术

本发明专利技术公开了基于全极化干涉合成孔径雷达数据的建筑物提取方法，所述方法包括：步骤1)依次计算全极化干涉合成孔径雷达图像中每个点的散射总功率SPAN；步骤2)确定功率阈值α1，依次比较图像中每个点的功率SPAN是否大于α1，如果判定结果是否定的，则该点不是建筑物，否则，对于SPAN比阈值α1大的点，转入步骤3)；步骤3)基于全极化干涉合成孔径雷数据，求解相干最优化问题，计算平均相干参数；步骤4)确定平均相干参数的阈值α2，判定值比阈值α2小的点不是建筑物，反之则判定为建筑物。本发明专利技术的方法可以提高雷达图像中建筑物的识别准确率。

【技术实现步骤摘要】

本专利技术涉及重复轨道全极化干涉合成孔径雷达的信息处理领域，特别涉及基于全极化干涉合成孔径雷达数据进行建筑物提取方法。
技术介绍
利用雷达图像提取或识别建筑物对于土地资源监管是非常重要的手段。目前，使用全极化合成孔径雷达数据进行建筑物提取是常用的方法。一般采用目标分解的方法来解译目标散射体的散射类型。当雷达的发射的电磁波照射到建筑物上，建筑物一般会产生很强的双次散射。如果目标散射体的回波中包含很强的双次散射，则认为目标很可能是建筑物。但是当建筑物的方向与雷达航迹有一定的夹角的时候，建筑物的双次散射减弱，并且会产生额外的交叉极化的后向散射。当建筑物的方向与雷达航迹的夹角很大时(45度左右)，额外产生的交叉极化很强，双次散射很弱，导致这样的建筑物被错误判定为森林。虽然目前一些目标分解方法是旋转不变的(比如H/Alpha分解、Touzi分解等)，但仍然也存在这个问题。目前，全极化干涉合成孔径雷达技术被广泛应用。对于重复轨道的全极化干涉合成孔径雷达，建筑物对雷达的时间基线不敏感，而自然目标则相反。所以可以利用极化干涉相干参数来区别建筑物和自然目标，并且和建筑物与雷达的航迹之间的夹角无关。
技术实现思路
本专利技术的目的在于克服了全极化合成孔径雷达数据受建筑物与雷达航迹之间的夹角调制的缺点，提出了基于全极化干涉合成孔径雷达数据进行建筑物提取方法，该方法利用同一幅场景的不同时间的全极化数据之间的干涉特性，弥补单一全极化数据受建筑物与雷达航迹之间夹角调制的缺点。为实现上述目的，本专利技术提出了基于全极化干涉合成孔径雷达数据的建筑物提取方法，所述方法包括：步骤1)依次计算全极化干涉合成孔径雷达图像中每个点的散射总功率SPAN；步骤2)确定功率阈值α1，依次比较图像中每个点的功率SPAN是否大于α1，如果判定结果是否定的，则该点不是建筑物，否则，对于SPAN比阈值α1大的点，转入步骤3)；步骤3)基于全极化干涉合成孔径雷达数据，求解相干最优化问题，计算平均相干参数；步骤4)确定平均相干参数的阈值α2，判定值比阈值α2小的点不是建筑物，反之则判定为建筑物。上述技术方案中，所述步骤1)具体为：首先读入全极化干涉合成孔径雷达图像中每个点的散射矩阵形式的数据，得到两个三维的散射矢量k1和k2；利用k1和k2构建6×6的半正定艾米特矩阵T6：其中，上标*T表示矩阵的共轭转置操作；对T6进行斑点滤波，降低斑点噪声；然后从T6中得到3×3矩阵T11、T22和Ω12数据：和最后得到总功率SPAN：SPAN＝trace(T11)+trace(T22)SPAN是两个通道的功率之和。上述技术方案中，所述步骤3)具体为：对3×3矩阵进行特征分解，得到特征矢量ω11、ω12和ω13，特征值υ11、υ12和υ13；对3×3矩阵进行特征分解，得到特征矢量ω21、ω22和ω23，特征值υ21、υ22和υ23；由此得到三个最优的相干参数γ1、γ2和γ3：利用三个特征值υ11、υ12和υ13定义三个伪概率p1、p2和p3：pi＝υ1i/(υ11+υ12+υ13)p1、p2和p3表示相应的特征矢量出现的概率，由于三个最优的相干参数是由特征矢量计算得到的，所以相干参数γ1、γ2和γ3出现的概率也是p1、p2和p3；那么，定义一个平均相干参数该参数包含了γ1、γ2和γ3的信息，也表示两个干涉通道的相干性。本专利技术的优点在于：1、本专利技术的方法提出了平均相干参数的概念，该参数包含了三个最优相干参数的全部信息，同时表示重复轨道全极化干涉合成孔径雷达的两次测量数据之间的相干性。建筑物的重复轨道全极化干涉合成孔径雷达的两次测量的数据是相干的，而且与建筑物与雷达航迹之间的夹角无关，克服了全极化合成孔径雷达数据受建筑物与雷达航迹之间的夹角调制的缺点；2、由于一些固定的自然目标(比如没有植被覆盖的土地)的两次测量数据也有很大的相干性，所以本专利技术的方法引入总功率SPAN来去除这些目标可能带来的误检测可能。附图说明图1是本专利技术的基于重复轨道全极化干涉合成孔径雷达数据的建筑物提取方法的流程图。图2是Oberpfaffenhofen区域的GoogleEarth光学图像。图3是本专利技术方法得到的总功率SPAN的图像。图4是仅利用SPAN信息的建筑物提取结果。图5是本专利技术方法提出的平均相干参数的图像。图6是仅利用平均相干参数得到的建筑物提取结果。图7是本专利技术方法得到的建筑物提取结果。具体实施方式下面结合附图和具体实验数据对本专利技术的方法进行详细说明。如图1所示，本专利技术的基于全极化干涉合成孔径雷达数据进行建筑物提取方法，包括下述步骤：步骤1)依次计算全极化干涉合成孔径雷达图像中每个点的散射总功率SPAN；首先读入重复轨道全极化干涉合成孔径雷达图像中每个点的散射矩阵形式的数据，得到两个三维的散射矢量k1和k2。利用k1和k2构建6×6的半正定艾米特矩阵T6其中，上标*T表示矩阵的共轭转置操作。对T6形式的数据进行斑点滤波，降低斑点噪声。然后从T6形式的数据中可以得到3×3矩阵T11、T22和Ω12数据：和最后得到总功率SPAN：SPAN＝trace(T11)+trace(T22)SPAN是两个通道的功率之和。步骤2)确定功率阈值α1，依次比较图像中点的功率SPAN是否大于α1，如果判定结果是否定的，则该点不是建筑物，否则，对于SPAN比阈值α1大的点，转入步骤3)；步骤3)基于重复轨道全极化干涉合成孔径雷达数据，求解S.R.Cloude和K.P.Papathanassiou提出的相干最优化问题，计算平均相干参数；对3×3矩阵进行特征分解，得到特征矢量ω11、ω12和ω13，特征值υ11、υ12和υ13；对3×3矩阵进行特征分解，得到特征矢量ω21、ω22和ω23，特征值υ21、υ22和υ23。可以得到三个最优的相干参数γ1、γ2和γ3：利用三个特征值υ11、υ12和υ13定义三个伪概率p1、p2和p3：pi＝υ1i/(υ11+υ12+υ13)表示相应的特征矢量出现的概率，由于三个最优的相干参数是由特征矢量计算得到的，所以相干参数γ1、γ2和γ3出现的概率也是p1、p2和p3。那么，可以定义一个平均相干参数该参数包含了γ1、γ2和γ3的信息，也表示两个干涉通道的相干性。步骤4)确定平均相干参数的阈值α2，判定值比阈值α2小的点不是建筑物，反之则判定为建筑物。下面结合实际的重复轨道全极化干涉合成孔径雷达数据验证本专利技术提出的平均相干参数以及建筑物提取方法的有效性。使用的数据是德国L波段机载全极化合成孔径雷达ESAR在Oberpfaffenhofen机场区域的沿同一航迹的两次飞行获得的数据，时间基线小于一天。该区域有大块的森林、机场建筑、机场草地、公路、耕地、居民区，但是居民区的房屋周围有树木。作为对照，图2为谷歌地球采集的该区域的光学图像。图3为总功率SPAN的图像，白色表示SPAN值较大。图4为SPAN检测结果，白色代表SPAN大于阈值1的点，黑色代表SPAN小于阈值1的点，可以看出建筑物和森林区域的功率较大。为了除去利用SPAN信息检测到的结果中的森林区域，可以进一步使用本专利技术提出的平均相干参数的信息，因为森林区域在时间上是去相干的，建筑物在时间上是相干的，即建筑物的平均相干参数的值较大，而森本文档来自技高网...

【技术保护点】
基于全极化干涉合成孔径雷达数据的建筑物提取方法，所述方法包括：步骤1)依次计算全极化干涉合成孔径雷达图像中每个点的散射总功率SPAN；步骤2)确定功率阈值α1，依次比较图像中每个点的功率SPAN是否大于α1，如果判定结果是否定的，则该点不是建筑物，否则，对于SPAN比阈值α1大的点，转入步骤3)；步骤3)基于全极化干涉合成孔径雷达数据，求解相干最优化问题，计算平均相干参数；步骤4)确定平均相干参数的阈值α2，判定值比阈值α2小的点不是建筑物，反之则判定为建筑物。

【技术特征摘要】
1.基于全极化干涉合成孔径雷达数据的建筑物提取方法，所述方法包括：步骤1)依次计算全极化干涉合成孔径雷达图像中每个点的散射总功率SPAN；步骤2)确定功率阈值α1，依次比较图像中每个点的功率SPAN是否大于α1，如果判定结果是否定的，则该点不是建筑物，否则，对于SPAN比阈值α1大的点，转入步骤3)；步骤3)基于全极化干涉合成孔径雷达数据，求解相干最优化问题，计算平均相干参数；步骤4)确定平均相干参数的阈值α2，判定值比阈值α2小的点不是建筑物，反之则判定为建筑物。2.根据权利要求1所述的基于全极化干涉合成孔径雷达数据的建筑物提取方法，其特征在于，所述步骤1)具体为：首先读入全极化干涉合成孔径雷达图像中每个点的散射矩阵形式的数据，得到两个三维的散射矢量k1和k2；利用k1和k2构建6×6的半正定艾米特矩阵T6：T6=<k1k2[k1*Tk2*T]>=T11Ω12Ω12*TT22]]>其中，*T表示矩阵的共轭转置操作；对T6进行斑点滤波，降低斑点噪声；然后从T6中得到3×3矩阵T11、T22和Ω12数据：和最后得到总功率SPAN：SPAN＝trace(T11)+trace(T22)SPAN是两个通道的功率之和。3.根据权利要...

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