一種改進的高分辨率sar成像自聚焦方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于雷達成像技術(shù)領(lǐng)域，特別涉及一種改進的高分辨率SAR成像自聚焦方 法，適用于機載SAR成像的運動誤差補償。
【背景技術(shù)】
[0002] 自聚焦方法是從數(shù)據(jù)本身出發(fā)，估計出對成像結(jié)果有影響的非空變相位誤差數(shù) 據(jù)，并對該非空變相位誤差數(shù)據(jù)進行補償，從而獲得良好成像結(jié)果。
[0003] 現(xiàn)有的自聚焦方法大體上可以分為三類，具體是：基于圖像偏移的方法，相位梯度 自聚焦（Phase Gradient Autofocus，PGA)的方法和基于代價函數(shù)的方法?；趫D像偏移 的方法雖易于實現(xiàn)，但在高分辨率SAR成像中估計形式復雜的運動誤差時，得到的運動誤 差估計精度達不到要求；相位梯度自聚焦方法是SAR成像中的運動補償過程最常采用的一 種自聚焦方法，該方法對于含有孤立散射點的數(shù)據(jù)，可以精確估計出其任意形式的運動誤 差，但對于信噪比較低、場景中無明顯強散射點的數(shù)據(jù)，如在沙漠、海面情況下得到的數(shù)據(jù)， 該方法的運動誤差估計性能也較差；基于代價函數(shù)的方法計算量一般比較大，并且該方法 很容易受到代價函數(shù)本身特性的影響，使用該方法也得不到很好的聚焦效果。

【發(fā)明內(nèi)容】

[0004] 針對現(xiàn)有方法存在的不足，本發(fā)明的目的在于提出一種改進的高分辨率SAR成像 自聚焦方法，該方法在正側(cè)視模式下對SAR成像的包絡(luò)誤差和相位誤差均進行校正，能夠 得到聚焦良好的高分辨率SAR成像，并且對相位誤差的類型和復雜度沒有限制。
[0005] 為達到上述技術(shù)目的，本發(fā)明采用如下技術(shù)方案予以實現(xiàn)。
[0006] -種改進的高分辨SAR成像自聚焦方法，包括以下步驟：
[0007] 步驟1，對SAR雷達天線接收到的雷達回波信號進行距離脈沖壓縮后，得到距離脈 壓雷達回波信號4心,,,）；其中，丨表示方位快時間，t m表示方位慢時間。
[0008] 步驟2,對距離脈壓雷達回波信號進行運動補償，并做距離向快速傅立葉變 換（FFT)，依次得到運動補償雷達回波信號s(f；，t m)和該運動補償雷達回波信號的簡化形 式巧./；々）；其中，f表示方位快時間，k表示方位慢時間，表示脈沖重復頻率。
[0009] 步驟3，將運動補償雷達回波信號的簡化形式乘以方位Dechirp函數(shù) 札(仁，〇，得到去除方位調(diào)頻率的運動補償雷達回波信號；其中，1表示方位慢時間，匕表 示脈沖重復頻率。
[0010] 步驟4,對去除方位調(diào)頻率的運動補償雷達回波信號沿距離向和方位向分別進 行二維插值后，得到二維插值雷達回波信號Sjk x;ky)，然后對該二維插值雷達回波信號 Sr(kx;k y)沿距離向作逆快速傅里葉變換（IFFT)，得到相位歷史域數(shù)據(jù)；其中，匕表示方位 向插值函數(shù)，ky表示距離向插值函數(shù)。
[0011] 步驟5,對相位歷史域數(shù)據(jù)沿距離向進行（J-j)倍降采樣操作，然后采用圖像偏移 的SAR自聚焦算法估計（J-j)倍降采樣操作后的相位歷史域數(shù)據(jù)的第（j+1)相位誤差，根 據(jù)該第（j+1)相位誤差得到（J_j)倍降采樣操作后的相位歷史域數(shù)據(jù)的第（j+1)包絡(luò)誤 差，并根據(jù)該第（j+1)包絡(luò)誤差對相位歷史域數(shù)據(jù)進行包絡(luò)補償后，得到完成包絡(luò)誤差補 償?shù)南辔粴v史域數(shù)據(jù)；其中，J表示人為設(shè)定的降采樣倍數(shù)，j e {〇, 1，2,…，J-l，J}，J為正 整數(shù)。
[0012] 步驟6,利用共軛梯度算法求解完成包絡(luò)誤差補償?shù)南辔粴v史域數(shù)據(jù)的相位誤差 估計值，并用該相位誤差估計值的相反數(shù)對完成包絡(luò)誤差補償?shù)南辔粴v史域數(shù)據(jù)進行相位 補償，得到完成包絡(luò)誤差補償和相位誤差補償?shù)南辔粴v史域數(shù)據(jù)。
[0013] 步驟7,對完成包絡(luò)誤差和相位誤差補償?shù)南辔粴v史域數(shù)據(jù)進行方位向逆快速傅 里葉變換（IFFT)后，得到聚焦良好的高分辨率SAR圖像。
[0014] 本發(fā)明的有益效果為：本發(fā)明方法對SAR成像中運動誤差的估計具有很高的穩(wěn)健 性和精確性，對SAR成像中運動誤差的類型和復雜度也沒有限制，是一種高精度的自聚焦 算法；并且，在沒有明顯孤立散射點、低信噪比、低對比度場景的SAR成像中，使用本發(fā)明方 法能夠獲得很好的聚焦效果。
【附圖說明】
[0015] 下面結(jié)合附圖和【具體實施方式】對本發(fā)明作進一步詳細說明。
[0016] 圖1是本發(fā)明的一種改進的高分辨SAR成像自聚焦方法流程示意圖；
[0017] 圖2是SAR平臺在正側(cè)視聚束模式中的幾何模型示意圖；
[0018] 其中，以場景中心點0為坐標原點，建立平面直角坐標系，場景中任意點目標 P(x，y)的航跡方向為沿X軸正方向，場景中心點0到合成孔徑雷達（SAR)中心的方向為Y 軸正方向，P(x，y)為場景中任意點目標，合成孔徑雷達（SAR)平臺的理想航跡是以速度v沿 X軸正方向、并距場景中心點0高度為H處的一條直線，場景中任意點目標P(x，y)在方位慢 時間t m時刻的理想位置為N，N的坐標為（vtm，Rb)，場景中任意點目標P(x，y)在方位慢時間 t m時刻的實際位置為A，A的坐標為（vt m+ A X，Rb+ A y)，Rb表示場景中任意點目標P (X，y)到 場景中心線的最短斜距；其中，AR(tJ表示合成孔徑雷達（SAR)平臺的非理想軌跡引入的 運動誤差，A x (tm)表示點目標P (X，y)在x軸的運動誤差，A y (tm)表示點目標P (X，y)在y 軸的運動誤差；
[0019] 圖3(a)是使用本發(fā)明方法得到的距離脈壓雷達回波信號包絡(luò)示意圖，
[0020] 圖3(b)是使用本發(fā)明方法得到的二維插值雷達回波信號包絡(luò)示意圖，
[0021] 圖3(c)是使用本發(fā)明得到完成包絡(luò)誤差補償?shù)南辔粴v史域數(shù)據(jù)的信號包絡(luò)示意 圖；
[0022] 圖4(a)是使用本發(fā)明得到完成包絡(luò)誤差補償?shù)南辔粴v史域數(shù)據(jù)中，沿方位向加 入相位誤差后利用共軛梯度算法估計相位誤差的結(jié)果示意圖，
[0023] 圖4(b)是采用傳統(tǒng)相位梯度自聚焦（PGA)方法估計相位誤差的結(jié)果示意圖，
[0024] 圖4(c)是采用傳統(tǒng)相位梯度自聚焦（PGA)算法進行相位補償后的點目標方位向 脈沖響應(yīng)函數(shù)剖面圖；
[0025] 圖5(a)是得到完成包絡(luò)誤差補償?shù)南辔粴v史域數(shù)據(jù)中，沿方位向加入相位誤差 后采用本發(fā)明方法的圖像偏移的SAR成像自聚焦算法估計相位誤差的結(jié)果示意圖，
[0026] 圖5(b)為采用本發(fā)明方法的圖像偏移的SAR自聚焦算法進行相位補償后的點目 標方位向脈沖響應(yīng)函數(shù)剖面圖。
【具體實施方式】
[0027] 參照圖1，為本發(fā)明的一種改進的高分辨率SAR成像自聚焦方法流程示意圖，該種 改進的高分辨SAR成像自聚焦方法，包括以下步驟：
[0028] 步驟1，對SAR雷達天線接收到的雷達回波信號進行距離脈沖壓縮后，得到距離脈 壓雷達回波信號4心《);其中，丨表示方位快時間，t m表示方位慢時間。
[0029] 步驟1的具體過程為：
[0030] 參照圖2,為SAR雷達在正側(cè)視聚束模式中的幾何模型，選取場景中心點0為坐標 原點，建立平面直角坐標系，SAR雷達航跡方向為沿X軸正方向，場景中心到合成孔徑雷達 (SAR)中心的方向為Y軸正方向，P為場景中任意一個點目標，其坐標為（x，y)，SAR雷達理 想航跡是以速度v沿X軸正方向、且距場景中心點0高度為H的一條直線，SAR雷達在方位 慢時間t m時刻的理想位置為N，N的坐標為（vt m，Rb)，SAR雷達在方位慢時間乜時刻的實 際位置為A，A的坐標為（vt m+ A x，Rb+ A y)，Rb表示SAR雷達到場景中心點0的最短斜距， AR(tJ表示SAR雷達非理想軌跡的運動誤差，Ax(tJ表示點目標P在x軸的運動誤差， A y (tm)表示點目標P在y軸的運動誤差。
[0031] 在方位慢時間tm時刻，SAR雷達到點目標P的理想瞬時斜距為RP(tJ，SAR雷達到 場景中心點〇的理想瞬時斜距為％(〇,其表達式分別如下所示：
[0034] SAR雷達到點目標P的實際瞬時斜距為R(tJ，且該實際瞬時斜距R(tm)表達式為：
[0036]因此，SAR雷達天線接收到的雷達回波信號^可以表示為：
[0038]其中，丨表示方位快時間，tm表示方位慢時間，c表示光速，Tp表示雷達回波 信號的脈沖寬度，y表示雷達回波信號的調(diào)頻率，f。表示雷達回波信號的載波頻率，
.范圍內(nèi)幅度為1的矩形窗函數(shù)。
[0039] 接著，對SAR雷達天線接收到的雷達回波信號作距離向快速傅里葉變換 (FFT)，將該SAR雷達