[發(fā)明專(zhuān)利]一種彈載SAR圖像校正方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201811218825.5
• 申請(qǐng)日: 2018-10-19
• 公開(kāi)（公告）號(hào): CN109471101B
• 公開(kāi)（公告）日: 2021-05-04
• 發(fā)明（設(shè)計(jì)）人: 李芬芬;黃龍;侍偉偉;周世平;胡哲;桂陽(yáng);靳永亮;韓建莉;石穩(wěn);李姣
• 申請(qǐng)（專(zhuān)利權(quán)）人: 湖北航天技術(shù)研究院總體設(shè)計(jì)所
• 主分類(lèi)號(hào): G01S13/90
• 分類(lèi)號(hào): G01S13/90;G01S7/41
• 代理公司: 武漢智權(quán)專(zhuān)利代理事務(wù)所(特殊普通合伙) 42225
• 代理人: 余浩

1.一種彈載SAR圖像校正方法，其特征在于，其包括：

根據(jù)彈體在孔徑中心時(shí)刻的參數(shù)建立地面坐標(biāo)系O-x-y-z；

以波束中心點(diǎn)P為原點(diǎn)，建立圖像坐標(biāo)系P-xp-yp；

讀入回波信號(hào)，得到聚焦后的斜平面圖像；

以所述圖像坐標(biāo)系建立圖像網(wǎng)格，根據(jù)地面實(shí)時(shí)圖像尺寸M×N和分辨率得到所述圖像網(wǎng)格上每一點(diǎn)的坐標(biāo)(xi，yi)；

建立旋轉(zhuǎn)矩陣，將所述圖像網(wǎng)格上每一點(diǎn)的坐標(biāo)(xi，yi)進(jìn)行旋轉(zhuǎn)變換，得到該點(diǎn)在地面坐標(biāo)系中的坐標(biāo)(xm，ym)；

根據(jù)所述坐標(biāo)(xm，ym)，計(jì)算所述圖像網(wǎng)格上的每一點(diǎn)在成像斜平面的距離信息和多普勒信息，并換算成對(duì)應(yīng)斜平面圖像的行信息和列信息；

對(duì)所述行信息和列信息以及對(duì)應(yīng)的斜平面圖像的灰度信息進(jìn)行二維sinc插值，得到圖像網(wǎng)格上所有點(diǎn)的幅度信息，從而得到校正后的SAR圖像；

所述圖像坐標(biāo)系P-xp-yp建立在沿波束方向，所述波束中心點(diǎn)P的坐標(biāo)為：

x0＝H/tan(β)*cos(α)

y0＝H/tan(β)*sin(α)

其中，H為彈體高度，α為波束方位角，β為波束俯仰角；

計(jì)算所述圖像網(wǎng)格上的每一點(diǎn)在成像斜平面的距離信息和多普勒信息，并換算成對(duì)應(yīng)斜平面圖像的行信息和列信息，具體包括：

根據(jù)所述坐標(biāo)(xm，ym)得到地面網(wǎng)格上任意點(diǎn)的斜距矢量根據(jù)慣導(dǎo)參數(shù)，得到速度矢量

根據(jù)所述和獲得該點(diǎn)對(duì)應(yīng)的斜視角θ_P，進(jìn)而計(jì)算出圖像網(wǎng)格上的該點(diǎn)的距離信息R_m，以及該點(diǎn)的多普勒信息f_a；

換算成對(duì)應(yīng)斜平面圖像的行信息row和列信息col為：

row＝round((R_m-R₀)/deltaR)

col＝round((f_a+prf/2)/delta_prf)

其中，deltaR＝c/2/fs，delta_prf＝prf/nan，round表示四舍五入取整，R₀為波門(mén)前沿對(duì)應(yīng)的距離，c為光速，nan為方位向采樣點(diǎn)數(shù)，prf為脈沖重復(fù)頻率，fs為采樣率。

2.如權(quán)利要求1所述的彈載SAR圖像校正方法，其特征在于：在所述地面坐標(biāo)系O-x-y-z中，Oxz平面在地平面上，Ox指向速度在地面投影方向，Oy垂直地面向上，Oz與Ox和Oy滿足右手法則。

3.如權(quán)利要求1所述的彈載SAR圖像校正方法，其特征在于：所述圖像坐標(biāo)系P-xp-yp，建立在沿波束方向、速度方向或地面正北向。