一种基于双曲线模型的大斜视SAR成像算法的制作方法

文档序号:12359723阅读:来源:国知局

技术特征:

1.一种基于双曲线模型的大斜视SAR成像算法,其特征在于步骤如下:

步骤1:根据距离等式采用双曲线模型,确定雷达与目标之间的瞬时斜距为

<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <msub> <mi>R</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>V</mi> <mi>t</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>R</mi> <mn>0</mn> </msub> <mi>V</mi> <mi>t</mi> <mi> </mi> <msub> <mi>sin&theta;</mi> <mn>0</mn> </msub> </mrow> </msqrt> </mrow>

其中,R(t,R0)表示雷达与目标之间的瞬时斜距,t表示载机飞行的慢时间,R0表示载机航线与场景中心线之间的最短距离,V表示载机飞行的速度,θ0表示雷达发射信号的斜视角度;

步骤2:根据瞬时斜距和载机雷达工作的参数,确定经解调的二维时域中的回波信号

<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>s</mi> <mi>s</mi> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>,</mo> <mi>&tau;</mi> <mo>;</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mi>G</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>m</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>m</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>&tau;</mi> <mo>-</mo> <mn>2</mn> <mfrac> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>,</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&CenterDot;</mo> <mi>exp</mi> <mrow> <mo>{</mo> <mrow> <mo>-</mo> <msub> <mi>j&pi;K</mi> <mi>r</mi> </msub> <msup> <mrow> <mo>(</mo> <mrow> <mi>&tau;</mi> <mo>-</mo> <mn>2</mn> <mfrac> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>,</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </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> <mi>R</mi> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>,</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mi>&lambda;</mi> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>

其中,τ是距离向时间,G(t)表示二维时域中回波信号的振幅,ma(t)表示方位向窗函数,mr(·)表示距离向窗函数,c表示光速,Kr表示发射信号的调频率,λ表示发射信号的波长;

步骤3:对回波信号进行距离向FFT,将信号变换到距离频域中,在距离频域中对信号乘以距离走动校正函数,所述的距离走动校正函数:

<mrow> <msub> <mi>H</mi> <mi>w</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mo>+</mo> <msub> <mi>f</mi> <mi>c</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mi>V</mi> <mi>t</mi> <mi> </mi> <msub> <mi>sin&theta;</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow>

其中,fτ表示距离向频率,fc表示对发射信号进行调制的载频;

步骤4:对距离走动校正后的信号进行方位向FFT,将信号变换到二维频域中;对信号依次进行距离压缩、方位频率三次项补偿和残余距离徙动校正:

所述的距离压缩函数为:

<mrow> <msub> <mi>H</mi> <mrow> <mi>a</mi> <mi>r</mi> <mi>c</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>f</mi> <mo>,</mo> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mo>;</mo> <mi>r</mi> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>&lsqb;</mo> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <mrow> <mfrac> <mi>&pi;</mi> <msub> <mi>K</mi> <mi>r</mi> </msub> </mfrac> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&pi;R</mi> <mn>0</mn> </msub> </mrow> <mrow> <msub> <mi>cf</mi> <mi>c</mi> </msub> </mrow> </mfrac> <mfrac> <msup> <mrow> <mo>(</mo> <mrow> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mn>3</mn> </msup> </mfrac> </mrow> <mo>)</mo> </mrow> <msup> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mn>2</mn> </msup> </mrow> <mo>&rsqb;</mo> </mrow> </mrow>

其中,f为方位向频率,V′=Vcosθ0

对距离压缩后的信号进行方位频率三次项的补偿,这一步补偿了因距离走动校正而引起的压缩后目标方位向旁瓣的不对称畸变问题,所述的方位频率三次项补偿函数为:

<mrow> <msub> <mi>H</mi> <mrow> <mi>t</mi> <mi>p</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>f</mi> <mo>,</mo> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mo>;</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>&lsqb;</mo> <mrow> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <msub> <mi>&pi;R</mi> <mn>0</mn> </msub> <msub> <mi>sin&theta;</mi> <mn>0</mn> </msub> <msup> <mi>&lambda;</mi> <mn>3</mn> </msup> <msup> <mi>f</mi> <mn>3</mn> </msup> <mo>/</mo> <mn>4</mn> <msup> <mi>V</mi> <mrow> <mo>&prime;</mo> <mn>3</mn> </mrow> </msup> <msub> <mi>cos&theta;</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&lambda;</mi> <msup> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mn>3</mn> </msup> </mrow> </mfrac> </mrow> <mo>&rsqb;</mo> </mrow> </mrow>

所述的残余距离徙动校正函数为:

<mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>c</mi> <mi>m</mi> <mi>c</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>,</mo> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> <mo>;</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mfrac> <mn>2</mn> <mi>c</mi> </mfrac> <mrow> <mo>(</mo> <mfrac> <msub> <mi>R</mi> <mn>0</mn> </msub> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <msub> <mi>f</mi> <mi>&tau;</mi> </msub> </mrow>

步骤5:对步骤4得到的信号进行距离向IFFT,将信号变换到距离多普勒域中,然后沿着方位向对信号进行方位向压缩;所述的方位向压缩函数为:

<mrow> <msub> <mi>H</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>;</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <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> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <msub> <mi>&pi;R</mi> <mn>0</mn> </msub> <msub> <mi>sin&theta;</mi> <mn>0</mn> </msub> <msup> <mi>&lambda;</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mn>3</mn> </msup> </mrow> <mrow> <mn>4</mn> <mi>V</mi> <mi> </mi> <msup> <mi>cos</mi> <mn>4</mn> </msup> <msub> <mi>&theta;</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mfrac> <mn>1</mn> <msup> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>f</mi> <mi>&lambda;</mi> <mo>/</mo> <mn>2</mn> <msup> <mi>V</mi> <mo>&prime;</mo> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mn>3</mn> </msup> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>

步骤6:对压缩后的信号再进行方位向IFFT,得到聚焦目标图像。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1