发明名称 基于压缩感知的滑动聚束SAR三维成像方法
摘要 基于压缩感知的滑动聚束SAR三维成像方法,本发明涉及基于压缩感知的滑动聚束SAR三维成像方法。本发明是要解决现技术成像区域小并且成像精度低的问题,而提供了基于压缩感知的滑动聚束SAR三维成像方法。一、机载雷达根据场景辐照区接收到的原始回波信号建立滑动聚束SAR回波模型;二、对模型中的单个雷达的二维滑动聚束SAR成像;三、将步骤二成像后的多个雷达的SAR图像形成一个三维矩阵;四、对三维矩阵切片进行压缩感知处理,即完成了基于压缩感知的滑动聚束SAR三维成像方法。本发明涉及三维成像领域。
申请公布号 CN104062659B 申请公布日期 2016.08.24
申请号 CN201410330876.2 申请日期 2014.07.11
申请人 哈尔滨工业大学 发明人 刘梅;孙志鹏;王鹏飞;王志贵
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 基于压缩感知的滑动聚束SAR三维成像方法,其特征在于基于压缩感知的滑动聚束SAR三维成像方法按以下步骤实现:一、机载雷达根据场景辐照区接收到的原始回波信号建立滑动聚束SAR回波模型;二、对模型中的单个雷达的二维滑动聚束SAR成像;所述步骤二中对模型中的单个雷达的二维滑动聚束SAR成像:(一)、根据SAR信号域的原始回波信号进行方位傅里叶变换得到距离多普勒域对原始回波信号进行方位向傅里叶变换,将SAR信号域变到距离多普勒域,根据式(1‑1),得到原始回波信号的精确的距离多普勒域的表达式为:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>r</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><mi>&tau;</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mn>1</mn></msub><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>j&pi;K</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>&tau;</mi><mo>-</mo><mfrac><mn>2</mn><mi>c</mi></mfrac><msub><mi>R</mi><mi>f</mi></msub><msup><mrow><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><msub><mi>&pi;R</mi><mn>0</mn></msub></mrow><mi>&lambda;</mi></mfrac><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mn>2</mn><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000011.GIF" wi="1438" he="270" /></maths>其中,C<sub>s</sub>是一个系数表示多普勒频率对信号轨迹的影响,f<sub>s</sub>表示方位向频率,R<sub>f</sub>表示在距离多普勒域中的距离徙动,K<sub>s</sub>为等效调频斜率,R<sub>0</sub>为飞行航迹上飞行平台到场景中心最近的距离;R<sub>f</sub>(f<sub>s</sub>,R<sub>0</sub>)=R<sub>0</sub>(1+C<sub>s</sub>(f<sub>s</sub>))           (2‑2)<maths num="0002"><math><![CDATA[<mrow><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><mrow><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow></mrow></msqrt></mfrac><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000012.GIF" wi="1061" he="214" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>K</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>K</mi><mi>r</mi></msub><mrow><mn>1</mn><mo>-</mo><msub><mi>K</mi><mi>r</mi></msub><msub><mi>R</mi><mn>0</mn></msub><mfrac><mrow><mn>2</mn><mi>&lambda;</mi></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mfrac><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000892802690000013.GIF" wi="1326" he="335" /></maths>(二)采用第一步相位相乘对距离多普勒域完成补余RCMC中的Chirp Scaling校正,将信号变换到二维频域第一步相位相乘的相位函数为:<maths num="0004"><math><![CDATA[<mrow><msub><mi>H</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><msub><mi>j&pi;K</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&tau;</mi><mo>-</mo><mover><mi>&tau;</mi><mo>^</mo></mover><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000014.GIF" wi="1449" he="78" /></maths>其中<img file="FDA0000892802690000015.GIF" wi="117" he="71" />为相对回波时间;<maths num="0005"><math><![CDATA[<mrow><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>2</mn><mi>c</mi></mfrac><msub><mi>R</mi><mn>0</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>C</mi><mi>s</mi></msub><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000016.GIF" wi="1262" he="125" /></maths>对距离多普勒域信号进行距离向傅里叶变换,信号变换到二维频域,变换后信号为距离向的脉冲信号:<maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>S</mi><mrow><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&tau;</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mn>2</mn></msub><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>R</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow><mo>}</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>{</mo><mrow><mi>j</mi><mi>&pi;</mi><mfrac><msubsup><mi>f</mi><mi>&tau;</mi><mn>2</mn></msubsup><mrow><msub><mi>K</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mfrac></mrow><mo>}</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>c</mi></mfrac><msub><mi>f</mi><mi>&tau;</mi></msub><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>}</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;</mi></mrow><mi>c</mi></mfrac><msub><mi>f</mi><mi>&tau;</mi></msub><msub><mi>R</mi><mn>0</mn></msub><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mrow><mn>2</mn><mo>-</mo><mn>7</mn></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000021.GIF" wi="1621" he="247" /></maths>其中,C<sub>2</sub>为常数,f<sub>s</sub>表示方位向频率,f<sub>τ</sub>表示距离向频率;(三)采用第二步相位相乘完成距离压缩、二次距离压缩与补余RCMC校正,将距离向的脉冲信号压缩为距离向的一个点第二步相位相乘的相位函数为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&tau;</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mi>&pi;</mi><mfrac><msubsup><mi>f</mi><mi>&tau;</mi><mn>2</mn></msubsup><mrow><msub><mi>K</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>C</mi><mi>S</mi></msub><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>}</mo><mi>exp</mi><mo>{</mo><mi>j</mi><mfrac><mrow><mn>4</mn><msub><mi>&pi;f</mi><mi>&tau;</mi></msub><msub><mi>R</mi><mn>0</mn></msub><msub><mi>C</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000022.GIF" wi="1582" he="151" /></maths>对距离向的脉冲信号进行方位向傅里叶反变换,得到距离多普勒域的一个点;(四)采用第三步相位相乘完成方位压缩与相位校正在距离多普勒域对方位向进行方位向压缩,然后将目标校正到零多普勒位置,在方位向将目标的回波信号压缩为一个点;第三步相位相乘的相位函数为:<maths num="0008"><math><![CDATA[<mrow><msub><mi>H</mi><mn>3</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><msub><mi>&pi;R</mi><mn>0</mn></msub></mrow><mi>&lambda;</mi></mfrac><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&lambda;f</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892802690000023.GIF" wi="1325" he="142" /></maths>对信号进行方位向的傅里叶逆变换,信号返回SAR图像域,输出图像;三、将步骤二成像后的多个雷达的SAR图像形成一个三维矩阵;四、对三维矩阵切片进行压缩感知处理,即完成了基于压缩感知的滑动聚束SAR三维成像方法。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号