发明名称 合成孔径激光成像雷达两维卷积成像方法
摘要 一种合成孔径激光成像雷达两维卷积成像方法,基本思想在于将复数化后的合成孔径激光成像雷达目标回波信号首先在交轨向附加二次项相位因子,然后进行两维卷积同时实现交轨向,顺轨向聚焦成像,并对成像结果强度信号进行输出显示。由于交轨向,顺轨向均采用卷积进行聚焦成像,降低了成像处理对计算机性能的要求,另外,两维卷积直接得到空域成像结果,无需进行两维频域-空域转换,直接进行输出显示,减少成像步骤,是合成孔径激光成像雷达回波信号成像处理的重要技术改进。
申请公布号 CN103885066B 申请公布日期 2016.05.11
申请号 CN201410108847.1 申请日期 2014.03.21
申请人 中国科学院上海光学精密机械研究所 发明人 孙志伟;职亚楠;孙建锋;周煜;侯培培;刘立人
分类号 G01S17/89(2006.01)I;G01S7/48(2006.01)I 主分类号 G01S17/89(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯
主权项 一种合成孔径激光成像雷达两维卷积成像方法,其特征在于,包括如下步骤:①合成孔径激光成像雷达发射系统向目标照射激光光波,照射到目标后被目标反射,反射光波经过合成孔径激光成像雷达接收系统进行光电接收、数字化与复数化处理后存储在计算机中,计算机存储的目标回波信号为i(t<sub>f</sub>,vt<sub>s</sub>):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>i</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>,</mo><msub><mi>vt</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mi>k</mi></munder><mi>A</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>S</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>v</mi><mi>y</mi></msub><msub><mi>t</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>W</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000011.GIF" wi="932" he="117" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>&times;</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mrow><mi>e</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mi>t</mi><mi>f</mi></msub><mo>+</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><mi>&lambda;</mi><mi>F</mi></mrow></mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>v</mi><mi>y</mi></msub><msub><mi>t</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000012.GIF" wi="1268" he="141" /></maths>式中,有如下设定:交轨向为x方向,顺轨向为y方向,A(x<sub>k</sub>,y<sub>k</sub>)为与SAIL结构,系统安排和目标点特性有关的接收光场因子,S(x<sub>k</sub>,y<sub>k</sub>‑v<sub>y</sub>t<sub>s</sub>)为由发射光斑和接收方向性特性决定的光学足趾函数,侧视SAIL是由光学天线孔径的衍射产生,直视SAIL是由发射机内口径成像产生,t<sub>f</sub>为交轨向快时间,t<sub>s</sub>为顺轨向慢时间,v<sub>y</sub>为雷达搭载平台顺轨向的运动速度,交轨向与快时间有关的窗函数为W(t<sub>f</sub>),f<sub>eq</sub>(x<sub>k</sub>)为目标在交轨向快时间上产生的等效拍频频率,λ为雷达发射激光中心波长,F为雷达光学足趾在目标面上的等效曲率半径;②对所述的目标回波信号交轨向附加二次项相位因子,所述的二次项相位因子为i<sub>a</sub>(t<sub>f</sub>):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>i</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><msub><mi>&lambda;f</mi><mi>a</mi></msub></mrow></mfrac><msubsup><mi>t</mi><mi>f</mi><mn>2</mn></msubsup><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000013.GIF" wi="1052" he="158" /></maths>式中,f<sub>a</sub>为所述的相位二次项的等效曲率半径,将所述的复数化后的目标回波信号i(t<sub>f</sub>,vt<sub>s</sub>)与i<sub>a</sub>(t<sub>f</sub>)式按下列(3)式相乘,得交轨向附加二次项相位因子后的回波信号为i<sub>b</sub>(t<sub>f</sub>,t<sub>s</sub>):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>i</mi><mi>b</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>,</mo><msub><mi>t</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>i</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>,</mo><msub><mi>vt</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>i</mi><mi>a</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000014.GIF" wi="554" he="94" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>=</mo><munder><mo>&Sigma;</mo><mi>k</mi></munder><mi>A</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>S</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>v</mi><mi>y</mi></msub><msub><mi>t</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>W</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>f</mi></msub><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mrow><mi>e</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mi>t</mi><mi>f</mi></msub><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000015.GIF" wi="1336" he="119" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>&times;</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><msub><mi>&lambda;f</mi><mi>a</mi></msub></mrow></mfrac><msubsup><mi>t</mi><mi>f</mi><mn>2</mn></msubsup><mo>&rsqb;</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><mi>&lambda;</mi><mi>F</mi></mrow></mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>v</mi><mi>y</mi></msub><msub><mi>t</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000882442990000016.GIF" wi="797" he="158" /></maths>③对交轨向附加二次项相位因子后的回波信号进行两维卷积,同时实现目标交轨向,顺轨向聚焦:将所述的交轨向附加二次项相位因子后的回波信号为i<sub>b</sub>(t<sub>f</sub>,t<sub>s</sub>)与二次项相位因子<img file="FDA0000882442990000021.GIF" wi="738" he="159" />按下列(4)式进行两维卷积,得到目标成像结果:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>i</mi><mrow><mi>i</mi><mi>m</mi><mi>a</mi><mi>g</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></munderover><munderover><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><mo>+</mo><mi>&infin;</mi></mrow></munderover><munder><mo>&Sigma;</mo><mi>k</mi></munder><mi>A</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>S</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>,</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><mi>&eta;</mi><mo>)</mo></mrow><mi>W</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mrow><mi>e</mi><mi>q</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>&xi;</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000882442990000022.GIF" wi="1387" he="142" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mo>&times;</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><msub><mi>&lambda;f</mi><mi>a</mi></msub></mrow></mfrac><msup><mi>&xi;</mi><mn>2</mn></msup><mo>&rsqb;</mo><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><mi>&lambda;</mi><mi>F</mi></mrow></mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><mi>&eta;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><msub><mi>&lambda;f</mi><mi>a</mi></msub></mrow></mfrac><msup><mrow><mo>(</mo><mi>&xi;</mi><mo>-</mo><mi>u</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000882442990000023.GIF" wi="1202" he="158" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mo>&times;</mo><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mi>j</mi><mfrac><mi>&pi;</mi><mrow><mi>&lambda;</mi><mi>F</mi></mrow></mfrac><msup><mrow><mo>(</mo><mi>&eta;</mi><mo>-</mo><mi>v</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mi>d</mi><mi>&xi;</mi><mi>d</mi><mi>&eta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000882442990000024.GIF" wi="1469" he="143" /></maths><img file="FDA0000882442990000025.GIF" wi="1372" he="149" /><img file="FDA0000882442990000026.GIF" wi="1212" he="159" />式中,u,v分别为成像结果交轨向,顺轨向空间坐标,ξ,η分别为卷积中交轨向,顺轨向坐标,<img file="FDA0000882442990000029.GIF" wi="74" he="55" />代表傅里叶变换,符号*代表卷积,交轨向成像位置为u=‑λf<sub>a</sub>f<sub>eq</sub>(x<sub>k</sub>),<img file="FDA0000882442990000027.GIF" wi="350" he="87" />的零值全宽为交轨向成像分辨率,顺轨向成像位置为v=y<sub>k</sub>,顺轨向成像分辨率为<img file="FDA0000882442990000028.GIF" wi="478" he="86" />的零值全宽;④将所述的目标成像结果的强度信号|i<sub>image</sub>(u,v)|<sup>2</sup>进行输出,并由显示器进行显示。
地址 201800 上海市嘉定区800-211邮政信箱