发明名称 极化SAR多成分目标分解方法
摘要 本发明提供了一种极化SAR多成分目标分解方法。该方法在以往四成分分解算法的基础上,增加了可能来自地形或粗糙表面的额外漫散射机制作为第五分解成分,即:将全极化SAR相干矩阵T或协方差矩阵C分解为五个散射成分(奇次散射、偶次散射、体散射、螺旋体散射和漫散射)的线性加权和。该算法采用改进的NNED方法求解体散射功率,并直接采用Pauli矩阵相应的奇次散射、偶次散射和漫散射模型,利用广义相似性参数GSP结合特征值分解直接求解奇次散射、偶次散射和漫散射的贡献值。本发明提供的新方法有利于提高极化SAR数据的地物分类精度。
申请公布号 CN105445712A 申请公布日期 2016.03.30
申请号 CN201510777393.1 申请日期 2015.11.11
申请人 中国测绘科学研究院 发明人 张继贤;魏钜杰;赵争;黄国满;余小萍
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 北京汇信合知识产权代理有限公司 11335 代理人 王杰
主权项 一种极化SAR多成分目标分解方法,其特征在于,包括如下步骤:步骤10:获得全极化SAR相干矩阵T,或者,获得全极化SAR协方差矩阵C并且将所述全极化SAR协方差矩阵C转化为全极化SAR相干矩阵T;对所述全极化SAR相干矩阵T进行极化方位角补偿,获得补偿后的相干矩阵T(θ);步骤20:对补偿后的相干矩阵T(θ)进行矩阵变换,获得协方差矩阵C;步骤30:按照公式1将协方差矩阵C分解为奇次散射、偶次散射、体散射、螺旋体散射和漫散射的线性加权和;C=f<sub>odd</sub>C<sub>odd</sub>+f<sub>dbl</sub>C<sub>dbl</sub>+f<sub>diff</sub>C<sub>diff</sub>+f<sub>vol</sub>C<sub>vol</sub>+f<sub>hlx</sub>C<sub>hlx</sub>         公式1式中,f<sub>odd</sub>、f<sub>dbl</sub>、f<sub>diff</sub>、f<sub>vol</sub>和f<sub>hlx</sub>分别表示待求的奇次散射、偶次散射、漫散射、体散射和螺旋体散射的贡献值;C<sub>odd</sub>、C<sub>dbl</sub>、C<sub>diff</sub>、C<sub>vol</sub>和C<sub>hlx</sub>分别表示已知的奇次散射模型、偶次散射模型、漫散射模型、体散射模型和螺旋体散射模型;步骤40:将协方差矩阵C和已知的五个散射模型代入公式1,可推导出螺旋体散射贡献值f<sub>hlx</sub>的计算表达式,见公式2:<math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>h</mi><mi>l</mi><mi>x</mi></mrow></msub><mo>=</mo><mn>2</mn><mo>|</mo><mi>Im</mi><mo>&lt;</mo><msubsup><mi>S</mi><mrow><mi>H</mi><mi>V</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>H</mi><mi>H</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>V</mi><mi>V</mi></mrow></msub><mo>)</mo></mrow><mo>&gt;</mo><mo>|</mo></mrow>]]></math><img file="FDA0000844435390000011.GIF" wi="582" he="102" /></maths>          公式2式中,Im(·)表示取复数的虚部;上标*表示复数共轭;运算符&lt;·&gt;表示窗口统计平均;运算符|·|表示复数求模;步骤50:确定排除螺旋体散射以外的第一残余协方差矩阵<img file="FDA0000844435390000012.GIF" wi="511" he="95" />步骤60:选择体散射模型,并利用第一残余协方差矩阵<img file="FDA0000844435390000013.GIF" wi="214" he="95" />和改进的NNED方法求解体散射功率f<sub>vol</sub>;步骤70:确定排除体散射以外的第二残余协方差矩阵<img file="FDA0000844435390000014.GIF" wi="692" he="94" />步骤80:对第二残余协方差矩阵<img file="FDA0000844435390000015.GIF" wi="215" he="95" />进行特征值分解;步骤90:利用广义相似性参数GSP,结合特征值分解结果,计算奇次散射、偶次散射和漫散射的贡献值f<sub>odd</sub>、f<sub>dbl</sub>、f<sub>diff</sub>。
地址 100830 北京市海淀区莲花池西路28号