发明名称 一种SAR大范围森林场景遥感数据的快速模拟方法
摘要 本发明涉及一种SAR大范围森林场景遥感数据的快速模拟方法,属于合成孔径雷达森林遥感技术领域。步骤包括参数定义、场景划分、获取虚拟散射单元的中心位置和等效雷达后向散射系数和最终计算得出SAR森林遥感模拟数据。本发明通过等效散射处理的方法将森林场景多个真实散射单元的SAR回波信号响应用一个虚拟散射单元的SAR回波信号响应进行等效,降低了森林场景离散化散射单元的数量,达到了提高SAR原始回波数据模拟效率的目的,实现了在高分辨率SAR系统下大范围森林遥感原始数据快速、高精度模拟。
申请公布号 CN102162844A 申请公布日期 2011.08.24
申请号 CN201010588049.5 申请日期 2010.12.07
申请人 北京理工大学 发明人 曾涛;胡程;孙晗伟
分类号 G01S7/41(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 代理人
主权项 1.一种SAR大范围森林场景遥感数据的快速模拟方法,其特征在于步骤为:1)对相关参数进行定义,包括坐标系定义和森林场景定义:1.1定义三维坐标系,其中Y轴为SAR运动方向,又被称为方位向,X轴和Z轴分别为地距向和高度向,SAR观测场景的中心下视角为θ<sub>c</sub>,中心斜距为r<sub>0</sub>;1.2定义森林场景,森林场景为长方体形状,其方位向范围为W<sub>y</sub>,地距向范围为W<sub>x</sub>,高度向范围为W<sub>z</sub>;2)将整个森林场景沿方位向等分为宽度为D<sub>y</sub>的子场景,最后一个宽度不足D<sub>y</sub>的子场景取其剩余宽度;3)对于每一个子场景沿地距向和高度向等分为宽度为D<sub>x</sub>、高度为D<sub>z</sub>的场景块,最后一个宽度不足D<sub>x</sub>或者高度不足D<sub>z</sub>的场景块取其剩余宽度或者剩余高度;4)在SAR的孔径中心时刻,将雷达对每个“场景块”中心的视线方向定义为<img file="FSA00000385642900011.GIF" wi="28" he="41" />方向,将每个场景块沿在XZ平面内与<img file="FSA00000385642900012.GIF" wi="29" he="43" />方向垂直的方向等分为宽度为D<sub>r</sub>的场景带,最后一个宽度不足D<sub>r</sub>的场景带取其剩余宽度,将每个场景带定义为一个虚拟散射单元;5)取每个场景带内沿<img file="FSA00000385642900013.GIF" wi="30" he="44" />方向的中点作为每个虚拟散射单元的中心位置;6)对场景带内所有真实散射单元的后向散射补偿斜距误差引入多普勒项后再进行相干叠加得到每个场景带的雷达后向散射系数,其中第n个场景带的等效雷达后向散射系数为<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>n</mi><mi>pq</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>n</mi></msub></munderover><msubsup><mi>&sigma;</mi><mi>i</mi><mi>pq</mi></msubsup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;&Delta;</mi><msub><mi>R</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub></mrow><mi>&lambda;</mi></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>式中,M<sub>n</sub>为第n个场景带内真实散射单元的个数,则<img file="FSA00000385642900015.GIF" wi="81" he="59" />第n个场景带的等效雷达后向散射系数,<img file="FSA00000385642900021.GIF" wi="80" he="59" />为第n个场景带内第i个真实散射单元的雷达后向散射系数,ΔR<sub>i0</sub>表示真实散射单元与虚拟散射单元在孔径中心时刻的斜距差,λ为雷达信号波长,p、q分别表示接收极化方式和发射极化方式;7)对所有的子场景按照步骤3)至6)进行处理,就可以得到全部虚拟散射单元的中心位置和等效雷达后向散射系数;8)将全部虚拟散射单元的中心位置和等效雷达后向散射系数代入SAR原始信号产生公式,得到SAR森林遥感原始回波信号;SAR原始信号产生公式为<maths num="0002"><![CDATA[<math><mrow><msup><mi>S</mi><mi>pq</mi></msup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msubsup><mi>S</mi><mi>k</mi><mi>pq</mi></msubsup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中S<sup>pq</sup>(u,t)为SAR森林遥感原始回波信号,K为全部虚拟散射单元的个数,u和t分别代表雷达方位位置和电磁波传播时间,<img file="FSA00000385642900023.GIF" wi="205" he="71" />为第k个虚拟散射单元的SAR原始信号。
地址 100081 北京市海淀区中关村南大街5号