主权项 |
一种基于极化合成孔径雷达的植被下地形方位角估计方法,其特征在于,包括以下步骤:步骤1,雷达影像噪声去除,包括以每个待恢复点为中心开一个大小为n×n的窗口,用窗口中的所有像素进行多视处理估计极化相干矩阵;然后,在n×n的窗口内,对极化相干矩阵每一个元素进行快速傅立叶变换,在频率域构建低通滤波器,将快速傅立叶变换后的极化相干矩阵输入低通滤波器,进行逆傅立叶变化获取低通滤波以后的低噪声极化相干矩阵,作为观测值矩阵Tobserver;步骤2,针对植被下地形有关的参数θg,构建植被层与植被下地表模型如下,Tobserver=Tv(fv,θv,σ)+Tg(fg,α,β,θg)其中,Tv(fv,θv,σ)与Tg(fg,α,β,θg)分别为植被层矩阵与植被下地表矩阵,Tv(fv,θv,σ)由以下植被层模型提供, <mrow> <msub> <mi>T</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>v</mi> </msub> <mo>,</mo> <msub> <mi>θ</mi> <mi>v</mi> </msub> <mo>,</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>16</mn> </mfrac> <msub> <mi>f</mi> <mi>v</mi> </msub> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>8</mn> </mtd> <mtd> <mo>-</mo> <mn>4</mn> <mi>q</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>4</mn> <mi>p</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mn>4</mn> <mi>p</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>4</mn> <mo>+</mo> <mn>4</mn> <mi>q</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>4</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>-</mo> <mn>4</mn> <mi>q</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>4</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>4</mn> <mi>p</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>-</mo> <mn>4</mn> <mi>q</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>4</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mn>4</mn> <mo>-</mo> <mn>4</mn> <mi>q</mi> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>4</mn> <mi>θ</mi> </mrow> <mi>v</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>Tg(fg,α,β,θg)由以下植被下地表模型提供, <mrow> <msub> <mi>T</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>g</mi> </msub> <mo>,</mo> <mi>α</mi> <mo>,</mo> <mi>β</mi> <mo>,</mo> <msub> <mi>θ</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mo>·</mo> <mfrac> <msub> <mi>f</mi> <mi>g</mi> </msub> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mi>α</mi> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>β</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msup> <mi>β</mi> <mo>*</mo> </msup> </mtd> <mtd> <msup> <mrow> <mo>|</mo> <mi>α</mi> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>·</mo> <mi>R</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> </mrow> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mrow> <mi>cos</mi> <mn>2</mn> <mi>θ</mi> </mrow> <mi>g</mi> </msub> </mtd> <mtd> <mi>sin</mi> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>g</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <mn>2</mn> <msub> <mi>θ</mi> <mi>g</mi> </msub> </mtd> <mtd> <mi>cos</mi> <msub> <mrow> <mn>2</mn> <mi>θ</mi> </mrow> <mi>g</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>其中,fv、θv、σ为植被层矩阵Tv(fv,θv,σ)包含的三个未知参数,分别为植被层散射强度、植被层方位偏移角、植被层随机度,p(σ)与q(σ)为以σ为输入量的多项式函数;fg、θg为地表散射强度、植被下地形方位偏移角,α与β为地表散射系数,β*为β的共轭系数;步骤3,根据步骤1所得观测值矩阵Tobserver,求解步骤2中植被层与植被下地表模型的未知参数θg。 |