发明名称 基于柱面几何的三维微波成像方法
摘要 本发明公开了一种基于柱面几何的三维微波成像方法,涉及微波成像技术,该方法通过设计各成像圆柱面对应的聚焦函数,依次完成不同圆柱面上的二维成像,最后实现整个目标区域的三维成像。本发明方法的处理过程没有近似,并且主要操作在高度波数域和方位波数域内进行,无需在三维波数域进行插值,因此可实现精确、高效的三维成像。
申请公布号 CN102540183A 申请公布日期 2012.07.04
申请号 CN201010587882.8 申请日期 2010.12.14
申请人 中国科学院电子学研究所 发明人 洪文;王彦平;谭维贤;林贇;吴一戎
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 周国城
主权项 1.一种基于柱面几何的三维微波成像方法,其特征在于,包括步骤如下:步骤S1:根据发射信号的形式将原始回波变换至距离压缩域,如果发射信号是线性调频信号,对获得的回波信号进行距离向压缩,获得信号S<sub>1</sub>(r,φ′,z′),其中,r为距离空间域,φ′∈[0,2π)为天线位置方位角,φ′方向定义为方位向,z′为天线位置高度,z′方向定义为高度向;如果发射信号为步进频连续波信号,对获得的回波信号作距离向逆傅里叶变换,获得信号S<sub>1</sub>(r,φ′,z′);步骤S2:对步骤S1所得的结果S<sub>1</sub>(r,φ′,z′)进行波传播损耗补偿,获得信号S<sub>2</sub>(r,φ′,z′);步骤S3:对步骤S2所得的结果S<sub>2</sub>(r,φ′,z′)作距离向傅里叶变换,获得信号S<sub>3</sub>(K<sub>ω</sub>,φ′,z′),其中K<sub>ω</sub>=2πf/c为波数,f为信号的发射频率,c为光速;步骤S4:对步骤S3所得的结果S<sub>3</sub>(K<sub>ω</sub>,φ′,z′)作高度向傅里叶变换,获得信号S<sub>4</sub>(K<sub>ω</sub>,φ′,K<sub>z</sub>),其中,K<sub>z</sub>为高度向波数;步骤S5:对步骤S4所得的结果S<sub>4</sub>(K<sub>ω</sub>,φ′,K<sub>z</sub>)作方位向傅里叶变换,获得信号S<sub>5</sub>(K<sub>ω</sub>,K<sub>φ</sub>,K<sub>z</sub>),其中,K<sub>φ</sub>为方位向波数;步骤S6:对步骤S5所得的结果S<sub>5</sub>(K<sub>ω</sub>,K<sub>φ</sub>,K<sub>z</sub>)与聚焦函数<img file="FDA0000038059290000011.GIF" wi="791" he="78" />相乘,获得信号S<sub>6</sub>(K<sub>ω</sub>,K<sub>φ</sub>,K<sub>z</sub>),其中,符号*表示共轭,ρ<sub>n</sub>=nΔρ为给定的成像圆柱面半径,n=0,1,2,..,(N-1),Δρ=c/2B,<img file="FDA0000038059290000012.GIF" wi="293" he="62" /><img file="FDA0000038059290000013.GIF" wi="177" he="62" />表示取ρ′/Δρ的整数部分;步骤S7:对步骤S6所得的结果S<sub>6</sub>(K<sub>ω</sub>,K<sub>φ</sub>,K<sub>z</sub>)沿K<sub>ω</sub>积分,获得信号S<sub>7</sub>(K<sub>φ</sub>,K<sub>z</sub>);步骤S8:对步骤S7所得的结果S<sub>7</sub>(K<sub>φ</sub>,K<sub>z</sub>)作方位向和高度向的二维逆傅里叶变换,获得半径为ρ<sub>n</sub>的圆柱面上的重建图像<img file="FDA0000038059290000021.GIF" wi="224" he="63" /><maths num="0001"><![CDATA[<math><mrow><mover><mi>I</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>n</mi></msub><mo>,</mo><mi>&phi;</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mrow><msub><mi>K</mi><mi>&phi;</mi></msub><mo>,</mo><msub><mi>K</mi><mi>z</mi></msub></mrow></munder><msub><mi>S</mi><mn>7</mn></msub><mrow><mo>(</mo><msub><mi>K</mi><mi>&phi;</mi></msub><mo>,</mo><msub><mi>K</mi><mi>z</mi></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><msub><mi>K</mi><mi>&phi;</mi></msub><mi>&phi;</mi><mo>+</mo><mi>j</mi><msub><mi>K</mi><mi>z</mi></msub><mi>z</mi><mo>)</mo></mrow><mi>d</mi><msub><mi>K</mi><mi>&phi;</mi></msub><mi>d</mi><msub><mi>K</mi><mi>z</mi></msub></mrow></math>]]></maths>其中,φ为目标位置方位角,z为目标位置高度;步骤S9:完成目标区域内所有圆柱面的重建,最后获得目标区域的三维图像;步骤S10:图像显示,直接在圆柱坐标系下显示图像,或通过坐标变换,将圆柱坐标系下的图像变换到直角坐标系(x,y,z)显示;其坐标变换方程为:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>x</mi><mo>=</mo><mi>&rho;</mi><mi>cos</mi><mi>&phi;</mi></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><mi>&rho;</mi><mi>sin</mi><mi>&phi;</mi></mtd></mtr><mtr><mtd><mi>z</mi><mo>=</mo><mi>z</mi></mtd></mtr></mtable></mfenced></math>]]></maths>其中,ρ为目标位置半径。
地址 100190 北京市海淀区北四环西路19号
您可能感兴趣的专利