主权项 |
一种基于旋转域极化相干特征全极化SAR建筑物倒损率估计方法,利用配准的灾前和灾后包含建筑物区域的全极化SAR图像,分别获得建筑物区域倒损前的极化散射矩阵S<sub>i</sub><sup>pre‑event</sup>,i=1,2,…,W,和建筑物区域倒损后的极化散射矩阵S<sub>i</sub><sup>post‑event</sup>,i=1,2,…,W,i表示全极化SAR图像第i个像素,W为全极化SAR图像的像素点个数;令极化散射矩阵S=S<sub>i</sub><sup>pre</sup><sup>‑event</sup>或S=S<sub>i</sub><sup>post‑event</sup>,其特征在于:先对上述每一个极化散射矩阵进行下面第一步至第三步的处理:第一步,极化散射矩阵绕全极化雷达视线作旋转处理;令极化散射矩阵<img file="FDA0001115367080000011.GIF" wi="363" he="148" />S<sub>HH</sub>为在水平极化H发射和水平极化H接收条件下获取的复后向散射系数;S<sub>VH</sub>为在水平极化H发射和垂直极化V接收条件下获取的复后向散射系数;S<sub>HV</sub>为在垂直极化V发射和水平极化H接收条件下获取的复后向散射系数;S<sub>VV</sub>为在垂直极化V发射和垂直极化V接收条件下获取的复后向散射系数;在绕全极化雷达视线方向,对极化散射矩阵S进行旋转处理,对旋转域中的旋转角θ,θ∈[‑π,π],计算经旋转处理后的旋转极化散射矩阵S(θ)的表达式:<maths num="0001"><math><![CDATA[<mrow><mi>S</mi><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><msubsup><mi>SR</mi><mn>2</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001115367080000012.GIF" wi="444" he="77" /></maths>其中,旋转极化散射矩阵<img file="FDA0001115367080000013.GIF" wi="579" he="158" />旋转矩阵<img file="FDA0001115367080000014.GIF" wi="525" he="146" />上标T为转置处理;由上式可得,旋转极化散射矩阵S(θ)各元素在旋转域中的表达式为:S<sub>HH</sub>(θ)=S<sub>HH</sub>cos<sup>2</sup>θ+S<sub>HV</sub>cosθsinθ+S<sub>VH</sub>cosθsinθ+S<sub>VV</sub> sin<sup>2</sup>θS<sub>HV</sub>(θ)=‑S<sub>HH</sub>cosθsinθ+S<sub>HV</sub> cos<sup>2</sup>θ‑S<sub>VH</sub>sin<sup>2</sup>θ+S<sub>VV</sub> cosθsinθS<sub>VH</sub>(θ)=‑S<sub>HH</sub>cosθsinθ‑S<sub>HV</sub>sin<sup>2</sup>θ+S<sub>VH</sub>cos<sup>2</sup>θ+S<sub>VV</sub>cosθsinθS<sub>VV</sub>(θ)=S<sub>HH</sub>sin<sup>2</sup>θ‑S<sub>HV</sub>cosθsinθ‑S<sub>VH</sub>cosθsinθ+S<sub>VV</sub> cos<sup>2</sup>θ;第二步,计算旋转域中的极化相干特征序列;将旋转域中的旋转角θ进行离散化处理,得到离散化的旋转角序列θ<sub>j</sub>,<img file="FDA0001115367080000015.GIF" wi="236" he="117" />j=0,1,…,N,j表示第j个旋转角,N+1为旋转角θ的个数,其取值根据实际情况确定;对极化通道HH和VV,计算旋转角序列θ<sub>j</sub>对应的极化相干特征序列|γ<sub>HH‑VV</sub>(θ<sub>j</sub>)|:<maths num="0002"><math><![CDATA[<mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>=</mo><mfrac><mrow><mo>|</mo><mo><</mo><msub><mi>S</mi><mrow><mi>H</mi><mi>H</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>·</mo><msubsup><mi>S</mi><mrow><mi>V</mi><mi>V</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>></mo><mo>|</mo></mrow><mrow><msqrt><mrow><mo><</mo><msub><mi>S</mi><mrow><mi>H</mi><mi>H</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>·</mo><msubsup><mi>S</mi><mrow><mi>H</mi><mi>H</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>></mo></mrow></msqrt><mo>·</mo><msqrt><mrow><mo><</mo><msub><mi>S</mi><mrow><mi>V</mi><mi>V</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>·</mo><msubsup><mi>S</mi><mrow><mi>V</mi><mi>V</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>θ</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>></mo></mrow></msqrt></mrow></mfrac></mrow>]]></math><img file="FDA0001115367080000016.GIF" wi="1186" he="206" /></maths>其中,<·>为集合平均处理,|·|为取绝对值处理,上标*为共轭处理;第三步,计算极化相干起伏度;计算极化相干特征序列|γ<sub>HH‑VV</sub>(θ<sub>j</sub>)|的极化相干起伏度(|γ<sub>HH‑VV</sub>|<sub>std</sub>)<sub>i</sub>:(|γ<sub>HH‑VV</sub>|<sub>std</sub>)<sub>i</sub>=std{|γ<sub>HH‑VV</sub>(θ<sub>j</sub>)}其中,std{·}为求序列的极化相干特征标准差;利用上述第一步至第三步,分别求得建筑物倒损前和倒损后全极化SAR图像中建筑物区域的极化相干起伏度(|γ<sub>HH‑VV</sub>|<sub>std</sub>)<sub>i</sub>,分别记为<img file="FDA0001115367080000021.GIF" wi="362" he="87" />和<img file="FDA0001115367080000022.GIF" wi="389" he="86" />第四步,计算建筑物极化倒损因子;以倒损后和倒损前建筑物区域的极化相干起伏度的比值作为建筑物极化倒损因子<img file="FDA0001115367080000023.GIF" wi="271" he="77" />即<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo>=</mo><mfrac><msubsup><mrow><mo>(</mo><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><mi>i</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mo>-</mo><mi>e</mi><mi>v</mi><mi>e</mi><mi>n</mi><mi>t</mi></mrow></msubsup><msubsup><mrow><mo>(</mo><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><mi>i</mi><mrow><mi>p</mi><mi>o</mi><mi>s</mi><mi>t</mi><mo>-</mo><mi>e</mi><mi>v</mi><mi>e</mi><mi>n</mi><mi>t</mi></mrow></msubsup></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001115367080000024.GIF" wi="694" he="182" /></maths>第五步,建筑物倒损率估计;建筑物极化倒损因子<img file="FDA0001115367080000025.GIF" wi="243" he="78" />与建筑物倒损率DL<sub>i</sub>的关系为:<img file="FDA0001115367080000026.GIF" wi="245" he="78" />取值越大,建筑物倒损率DL<sub>i</sub>越大,这样,建立的建筑物倒损率DL<sub>i</sub>的反演表达式为:<maths num="0004"><math><![CDATA[<mrow><msub><mi>DL</mi><mi>i</mi></msub><mi>=</mi><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>·</mo><msup><mrow><mo>(</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo>)</mo></mrow><mn>4</mn></msup><mo>+</mo><msub><mi>p</mi><mn>2</mn></msub><mo>·</mo><msup><mrow><mo>(</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo>)</mo></mrow><mn>3</mn></msup><mo>+</mo><msub><mi>p</mi><mn>3</mn></msub><mo>·</mo><msup><mrow><mo>(</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>p</mi><mn>4</mn></msub><mo>·</mo><mrow><mo>(</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>p</mi><mn>5</mn></msub></mrow></mtd><mtd><mrow><mn>1</mn><mo>≤</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mn>0</mn><mo><</mo><msubsup><mi>Ratio</mi><mi>i</mi><mrow><mo>|</mo><msub><mi>γ</mi><mrow><mi>H</mi><mi>H</mi><mo>-</mo><mi>V</mi><mi>V</mi></mrow></msub><msub><mo>|</mo><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow></msub></mrow></msubsup><mo><</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001115367080000027.GIF" wi="1942" he="194" /></maths>其中,p<sub>1</sub>、p<sub>2</sub>、p<sub>3</sub>、p<sub>4</sub>和p<sub>5</sub>为拟合参数,基于建筑物区域内真实倒损率的样本像素,利用该真实倒损率和估计倒损率,通过多项式拟合得到;遍历极化SAR图像中的所有建筑物区域,即可得到建筑物倒损率估计图。 |