主权项 |
一种基于广义反射率模型的高分辨雷达快速成像方法,其特征在于,所述成像方法包括以下步骤:1)建立雷达成像系统,获取雷达散射数据:雷达成像系统包括T个发射机,R个接收机,发射信号的频率个数为F,发射机依次向目标成像区域发射信号,并由全部接收机接收回波信号,则当第t个发射机发射信号后,各个接收机处接收到的回波信号依次为y<sub>F,t</sub>=[y<sub>F,t,1</sub>;y<sub>F,t,2</sub>;…y<sub>F,t,R</sub>],其中,T和R分别为≥2的自然数,t=1,2,…,T,F为≥2的自然数;2)建立广义反射率模型:广义反射率模型是建立在波恩近似模型基础上,并包含成像目标的各向异性特性与频率特性,即不同角度发射机下成像目标的反射率不同,不同工作频率下成像目标的反射率不同的特性;3)数据变换和数据整合:a)接收机接收频域回波信号;b)根据雷达系统函数的结构性特征,依据广义反射率模型,将雷达成像系统划分为K个子孔径或子频带,并将目标成像区域划分为B个子成像区域,根据上述划分结果,对频域回波信号进行组合整理,得到关于第b个子成像区域与第k个子孔径或子频带的回波数据函数y<sub>(k,b)</sub>,其中,k=1,2,…,k,b=1,2,…,B,根据自由空间并矢格林函数构造雷达成像系统函数A<sub>(k,b)</sub>,则广义反射率x<sub>(k,b)</sub>满足方程(1):y<sub>(k</sub>,<sub>b)</sub>=A<sub>(k</sub>,<sub>b)</sub>x<sub>(k</sub>,<sub>b)</sub>+n<sub>(k</sub>,<sub>b)</sub> (1)其中,n<sub>(k,b)</sub>表示模型误差,则雷达成像问题转化为求解方程(1)中未知数x<sub>(k,b)</sub>的逆问题;c)在方程(1)两边分别乘以系统函数A<sub>(k,b)</sub>的共轭转置矩阵<img file="FDA0000952503700000011.GIF" wi="139" he="69" />即对方程(1)进行对偶变换得到方程(2):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></msub><mo>=</mo><msub><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></msub><msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>n</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952503700000012.GIF" wi="918" he="71" /></maths>其中,<img file="FDA0000952503700000013.GIF" wi="414" he="71" />是表示图像处理的系统函数,<img file="FDA0000952503700000014.GIF" wi="405" he="71" />是表示后向投影成像结果,n<sub>(k,b)</sub>是模型误差,上述对偶变换将雷达成像问题转化为基于物理机制的雷达图像处理问题;4)对子成像区域进行并行成像:基于广义反射率模型,针对方程(2)根据梯度迭代算法对各个子成像区域进行并行成像,第b个子成像区域下的第k个子孔径或子频带的第m次迭代得到广义反射率<img file="FDA0000952503700000015.GIF" wi="112" he="71" />满足方程(3):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mi>m</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>α</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mi>d</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952503700000021.GIF" wi="845" he="71" /></maths>其中,<img file="FDA0000952503700000022.GIF" wi="451" he="174" />为第b个子成像区域下的第k个子孔径或子频带的第m次迭代的步长因子,<img file="FDA0000952503700000023.GIF" wi="778" he="78" />为第b个子成像区域下的第k个子孔径或子频带的第m次迭代的梯度函数;利用广义反射率模型将第b个子成像区域下的所有K个子孔径的成像结果x<sub>(k,b)</sub>进行图像融合,得到第b个子成像区域的成像结果x<sub>(b)</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>x</mi><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow></msub><mo>=</mo><mo>|</mo><mo>|</mo><msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msup><mrow><mo>(</mo><mrow><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msup><mrow><mo>|</mo><msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mi>b</mi><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo></mrow><mi>p</mi></msup></mrow><mo>)</mo></mrow><mrow><mi>q</mi><mo>/</mo><mi>p</mi></mrow></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>q</mi></mrow></msup><mn>2</mn><mo>≤</mo><mi>p</mi><mo>≤</mo><mi>∞</mi><mo>,</mo><mn>0</mn><mo>≤</mo><mi>q</mi><mo>≤</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000952503700000024.GIF" wi="1717" he="119" /></maths>其中,N表示x<sub>(k,b)</sub>的元素个数,n代表第n个元素,n=1,2,…,N,p和q是范数指数;5)图像融合:将B个子成像区域的成像结果x<sub>(b)</sub>进行图像融合可得到一幅完整目标成像区域的高分辨雷达成像图像X。 |