曹 曄，師亞輝，閆海鵬，張劍琦

艦船目標(biāo)混合式SAR/ISAR成像算法研究

曹 曄，師亞輝，閆海鵬，張劍琦

（北京遙測(cè)技術(shù)研究所 北京 100076）

針對(duì)高動(dòng)態(tài)平臺(tái)和大斜視角成像條件下運(yùn)動(dòng)艦船成像模糊的問題，從提高成像處理各個(gè)環(huán)節(jié)的精度的角度出發(fā)，提出了一種基于加速度模型的混合式SAR/ISAR成像算法，實(shí)現(xiàn)了兩種算法的有效融合：首先利用運(yùn)動(dòng)估計(jì)結(jié)合距離徙動(dòng)校正完成回波包絡(luò)粗對(duì)齊，同時(shí)在距離多普勒域通過相位濾波去除調(diào)頻率空變，然后通過相關(guān)對(duì)準(zhǔn)實(shí)現(xiàn)包絡(luò)精對(duì)齊，補(bǔ)償由于運(yùn)動(dòng)估計(jì)不準(zhǔn)、模型失配以及慣導(dǎo)測(cè)量等引入的誤差，最后結(jié)合加權(quán)多特顯點(diǎn)自聚焦法對(duì)剩余非空變相位項(xiàng)進(jìn)行統(tǒng)一校正。仿真結(jié)果表明，所提算法能夠有效改善成像效果，得到聚焦良好的目標(biāo)圖像。

混合式SAR/ISAR；艦船目標(biāo)；平動(dòng)補(bǔ)償；高分辨成像

引 言

海上超大型艦船是海洋艦隊(duì)的重要組成部分，針對(duì)超大型艦船目標(biāo)的高分辨成像在當(dāng)代海戰(zhàn)中具有重要意義。而在實(shí)際工程應(yīng)用當(dāng)中，由于雷達(dá)平臺(tái)的高速機(jī)動(dòng)和前斜視角的成像條件，針對(duì)非合作艦船目標(biāo)的高分辨成像更加困難。為獲取更多的目標(biāo)信息，需要研究提高成像結(jié)果分辨率的算法。

目前針對(duì)艦船目標(biāo)的成像算法主要包括逆合成孔徑雷達(dá)ISAR（Inverse Synthetic Aperture Radar）成像算法、基于運(yùn)動(dòng)補(bǔ)償?shù)暮铣煽讖嚼走_(dá)SAR（Synthetic Aperture Radar）成像算法以及混合式SAR/ISAR成像算法。ISAR成像算法通過平動(dòng)補(bǔ)償和轉(zhuǎn)臺(tái)成像進(jìn)行成像處理，在雷達(dá)平臺(tái)高速運(yùn)動(dòng)和成像斜視角很大的情況下，目標(biāo)各散射點(diǎn)距離徙動(dòng)量增大且空變嚴(yán)重，相鄰回波相關(guān)性下降，包絡(luò)對(duì)齊難以精確對(duì)準(zhǔn)，而且，基于特顯點(diǎn)的針對(duì)所有距離單元的統(tǒng)一初相校正忽略了距離和方位空變，也存在誤差，導(dǎo)致ISAR圖像散焦，成像質(zhì)量下降?；谶\(yùn)動(dòng)補(bǔ)償?shù)腟AR成像算法利用回波數(shù)據(jù)對(duì)艦船平動(dòng)速度進(jìn)行估計(jì)和補(bǔ)償，但在艦船運(yùn)動(dòng)非平穩(wěn)的情況下，比如目標(biāo)存在加速和拐彎運(yùn)動(dòng)時(shí)，現(xiàn)有算法難以精確估計(jì)艦船的運(yùn)動(dòng)參數(shù)，估計(jì)誤差導(dǎo)致補(bǔ)償不準(zhǔn)確，造成圖像散焦，而且，SAR成像算法及回波模型依賴于成像幾何模型，在實(shí)際應(yīng)用當(dāng)中，難以實(shí)現(xiàn)對(duì)雷達(dá)運(yùn)動(dòng)軌跡的精確建模，模型誤差也會(huì)導(dǎo)致成像過程中的補(bǔ)償不準(zhǔn)確，最終影響成像質(zhì)量。

為了解決艦船成像模糊的問題，早在1989年，Ward K D[1]就首次提出了混合式SAR/ISAR的概念和設(shè)想；隨后幾年，Ward K D、Tough R J A和Porter N J等一同致力于這方面的研究，并于1993年聯(lián)合發(fā)表文章統(tǒng)一了SAR、ISAR、混合式SAR/ISAR在雷達(dá)成像中的應(yīng)用方法[2,3]?；旌鲜絊AR/ISAR的實(shí)現(xiàn)主要有兩種方式[4]，一是在SAR圖像中重新聚焦運(yùn)動(dòng)目標(biāo)：首先對(duì)SAR的子孔徑數(shù)據(jù)進(jìn)行ISAR處理，然后將處理后的子孔徑數(shù)據(jù)進(jìn)行相干組合，得到新的SAR數(shù)據(jù)，最后用SAR重新處理新的數(shù)據(jù)便可以得到聚焦的運(yùn)動(dòng)目標(biāo)圖像，這種方法屬于利用ISAR的優(yōu)勢(shì)來重新聚焦SAR中的運(yùn)動(dòng)目標(biāo)；二是結(jié)合多個(gè)ISAR的數(shù)據(jù)來等效為SAR的數(shù)據(jù)：SAR的處理用于組合數(shù)據(jù)可以得到高的距離橫向分辨率，并提取目標(biāo)三維特性，這種處理方法屬于利用SAR的優(yōu)勢(shì)來處理ISAR數(shù)據(jù)。

國內(nèi)對(duì)混合式SAR/ISAR成像技術(shù)的研究起步較晚，但近年來也有一些成果，主要利用上述方式一，即利用ISAR來處理SAR的數(shù)據(jù)[5-7]。典型的如文獻(xiàn)[6]中利用SAR粗成像和ISAR成像進(jìn)行序貫處理。本文在此基礎(chǔ)上進(jìn)行了改進(jìn)，采用加速度模型下SAR成像結(jié)合運(yùn)動(dòng)估計(jì)，并利用ISAR成像步驟對(duì)SAR成像誤差進(jìn)行補(bǔ)償?shù)幕旌鲜匠上穹椒ǎ瑢?duì)成像處理的各個(gè)環(huán)節(jié)進(jìn)行更精確的補(bǔ)償處理，最終實(shí)現(xiàn)了運(yùn)動(dòng)艦船目標(biāo)的兩維高分辨成像。

1 場(chǎng)景建模及回波仿真

圖1 成像幾何關(guān)系

圖2 艦船目標(biāo)模型

表1 仿真參數(shù)

2 混合式SAR/ISAR成像過程

2.1 問題分析

若采用單一的ISAR成像算法進(jìn)行處理，雖然不用考慮成像模型的影響，但由于雷達(dá)平臺(tái)運(yùn)動(dòng)速度快，以及成像斜視角大，導(dǎo)致同一散射點(diǎn)在相鄰回波中的距離徙動(dòng)量較大（大于距離分辨單元）且存在空變性，無論采用非參數(shù)化的方法還是參數(shù)化的方法都難以實(shí)現(xiàn)回波包絡(luò)的精確對(duì)齊；同時(shí)目前廣泛應(yīng)用的初相校正方法也是通過選擇特顯點(diǎn)單元進(jìn)行估計(jì)進(jìn)而對(duì)所有距離單元進(jìn)行統(tǒng)一校正，當(dāng)存在較嚴(yán)重的空變時(shí)必然會(huì)產(chǎn)生殘留誤差，從而影響方位聚焦，導(dǎo)致成像質(zhì)量下降。

若采用單一的SAR成像算法對(duì)運(yùn)動(dòng)艦船目標(biāo)進(jìn)行成像處理，需要依靠基于回波的多普勒參數(shù)估計(jì)和復(fù)雜的運(yùn)動(dòng)補(bǔ)償，由于雷達(dá)平臺(tái)的高速機(jī)動(dòng)和艦船目標(biāo)的非合作運(yùn)動(dòng)，實(shí)際成像幾何模型十分復(fù)雜。仿真模型與實(shí)際情況存在偏差、艦船運(yùn)動(dòng)的復(fù)雜性及估計(jì)不精確、慣性測(cè)量誤差等因素會(huì)導(dǎo)致運(yùn)動(dòng)補(bǔ)償不夠精確，最終影響成像效果。

圖3 原始回波數(shù)據(jù)

表2 艦船在5級(jí)海況下的各維晃動(dòng)參數(shù)

因此，單一的SAR或者ISAR成像算法均無法滿足高分辨成像的要求?；谏鲜鲈?，本文提出了一種基于加速度模型的混合式SAR/ISAR成像算法，能夠?qū)崿F(xiàn)回波包絡(luò)的精確對(duì)齊以及回波相位的精確補(bǔ)償：首先利用更精確的加速度模型和艦船目標(biāo)運(yùn)動(dòng)參數(shù)估計(jì)，通過距離徙動(dòng)校正等完成包絡(luò)粗對(duì)齊和空變消除；其次在距離-多普勒域利用高次相位濾波來去除調(diào)頻率的空變性；最后通過ISAR處理中的相關(guān)對(duì)準(zhǔn)進(jìn)一步實(shí)現(xiàn)回波包絡(luò)的精確對(duì)齊，并結(jié)合加權(quán)多特顯點(diǎn)自聚焦法對(duì)剩余非空變項(xiàng)進(jìn)行統(tǒng)一校正。圖4是成像算法流程。

圖4 算法流程

2.2 混合式SAR/ISAR成像過程

2.2.1 條件假設(shè)

為便于分析，首先進(jìn)行如下假設(shè)：

2.2.2 成像處理過程

下面詳細(xì)介紹混合式SAR/ISAR成像處理的過程，具體如下：

①距離壓縮

首先將回波信號(hào)沿距離維做FFT變換，在距離頻域完成距離向壓縮，對(duì)應(yīng)的壓縮函數(shù)為

②基于距離徙動(dòng)校正的包絡(luò)粗對(duì)齊

由于雷達(dá)平臺(tái)的高速運(yùn)動(dòng)，同一散射點(diǎn)在相鄰回波間的距離徙動(dòng)量較大，通常大于距離分辨單元，且由于成像斜視角較大，各散射點(diǎn)對(duì)應(yīng)的距離方位空變不能忽略，直接采用相關(guān)函數(shù)法進(jìn)行包絡(luò)對(duì)齊必然存在誤差，效果不夠理想。下面首先通過距離徙動(dòng)校正對(duì)回波包絡(luò)進(jìn)行粗對(duì)齊，從而將距離徙動(dòng)量減小到一定范圍內(nèi)，增加相鄰回波的相關(guān)性，后續(xù)通過相關(guān)對(duì)準(zhǔn)能夠得到精度更高的包絡(luò)對(duì)齊結(jié)果。

其中

③相位濾波

將信號(hào)變換到距離-多普勒域，由于成像斜視角較大，回波信號(hào)方位調(diào)頻率空變不能忽略，這意味著無法對(duì)方位進(jìn)行統(tǒng)一的補(bǔ)償聚焦處理，下面通過相位濾波來去除方位調(diào)頻率隨方位位置的變化。

④包絡(luò)精對(duì)齊

步驟5 返回步驟2對(duì)第3次回波進(jìn)行包絡(luò)對(duì)齊，直至所有回波對(duì)齊。

⑤加權(quán)多特顯點(diǎn)自聚焦

經(jīng)過上述步驟的處理，已經(jīng)完成了回波包絡(luò)的精確對(duì)齊，通過步驟③的相位濾波處理，沿方位向回波信號(hào)相位中關(guān)于方位慢時(shí)間的各項(xiàng)系數(shù)是非方位空變的，可以進(jìn)行統(tǒng)一補(bǔ)償。在SAR成像處理中，可以利用高階多項(xiàng)式模型統(tǒng)一補(bǔ)償，由于采用模型的方法會(huì)不可避免地引入誤差，下面采用基于數(shù)據(jù)估計(jì)的加權(quán)多特顯點(diǎn)自聚焦法對(duì)非方位空變項(xiàng)進(jìn)行補(bǔ)償，具體步驟如下：

⑥方位壓縮

經(jīng)過上述處理，實(shí)現(xiàn)了回波包絡(luò)的精確對(duì)齊以及回波相位的精確補(bǔ)償。此外，由于本文選擇的是大型艦船作為目標(biāo)進(jìn)行成像，其晃動(dòng)較為平緩，此時(shí)回波信號(hào)在方位向近似為單頻信號(hào)，通過FFT即可得到聚焦良好的成像結(jié)果。

3 仿真結(jié)果及成像效果分析

3.1 仿真結(jié)果

下圖6（a）、圖6（b）分別仿真了基于艦船平動(dòng)補(bǔ)償?shù)膯我籗AR成像算法和單一ISAR成像算法的成像結(jié)果。文獻(xiàn)6中的混合式SAR/ISAR成像算法是目前應(yīng)用較為廣泛的一種混合式成像算法，其算法流程如圖5所示：首先通過距離壓縮、方位壓縮得到SAR的粗成像結(jié)果（忽略了回波的空變），再通過目標(biāo)提取、方位向反壓縮，得到距離向壓縮-方位向未壓縮的兩維時(shí)域結(jié)果，最后通過平動(dòng)補(bǔ)償和轉(zhuǎn)臺(tái)成像完成ISAR成像部分的處理，得到目標(biāo)圖像。這種算法涉及兩個(gè)完整成像過程（SAR成像和ISAR成像）的序貫進(jìn)行，沒有從本質(zhì)上進(jìn)行算法的融合，其優(yōu)點(diǎn)是能夠從中間步驟獲得目標(biāo)的粗成像結(jié)果，能夠預(yù)先對(duì)最終高分辨成像結(jié)果有基本掌握，也便于現(xiàn)有成像系統(tǒng)的升級(jí)改造，缺點(diǎn)是存在冗余計(jì)算，降低了成像效率。圖6（c）為針對(duì)文獻(xiàn)6成像過程的仿真結(jié)果。圖6（d）為本文算法的仿真結(jié)果。仿真具體參數(shù)見表1，艦船晃動(dòng)參數(shù)見表2。

圖5 文獻(xiàn)6算法流程

圖6 斜視角70°下不同算法的成像結(jié)果

對(duì)比圖6（a）和圖6（b）可以看出，后者的仿真結(jié)果總體上優(yōu)于前者：圖中紅色方框標(biāo)出了艦船桅桿的兩條輪廓線，兩個(gè)紅色橢圓標(biāo)出了船頭和船側(cè)的部分散射點(diǎn)，從這些位置可以明顯看出兩種算法在此處成像效果的差異，特別是在船頭部分，在SAR成像結(jié)果當(dāng)中，4個(gè)散射點(diǎn)存在模糊黏連，而在ISAR成像結(jié)果中，4個(gè)散射點(diǎn)清晰可辨；對(duì)比圖6（c）和圖6（d）兩種混合式SAR/ISAR成像算法的成像結(jié)果發(fā)現(xiàn)，兩種方法的成像效果均較好，大部分的散射點(diǎn)成像均十分清楚，桅桿的四條豎棱清晰可見，圖中兩個(gè)綠色橢圓標(biāo)出了船尾兩個(gè)拐角處的部分散射點(diǎn)，對(duì)比可以看出本文所提混合式SAR/ISAR成像算法的成像效果更勝一籌。

圖7 斜視角80°下不同算法的成像結(jié)果

3.2 成像效果分析

為對(duì)上述各種算法的成像效果進(jìn)行定量的分析，下面選取艦船目標(biāo)的5個(gè)邊緣點(diǎn)和1個(gè)中心點(diǎn)來進(jìn)行觀察分析。以峰值點(diǎn)下降3dB對(duì)應(yīng)的距離為分辨率，依次仿真6個(gè)散射點(diǎn)的兩維幅度譜，并通過局部插值來提高讀取圖像分辨率的準(zhǔn)確性。表3列出了斜視角70°情況下對(duì)應(yīng)圖6各種算法的兩維分辨率結(jié)果。

表3 斜視角70°時(shí)的分辨率結(jié)果

從表3可以看出，在70°的大斜視角情況下，各成像算法距離向分辨率差別不大，均達(dá)到1.5米以內(nèi)，方位向分辨率屬本文所提的混合式SAR/ISAR成像算法最優(yōu)，分辨率在3m以內(nèi)；文獻(xiàn)6中的混合式SAR/ISAR成像算法和ISAR成像算法次之，兩種算法成像效果差別不大，方位向分辨率在4m以內(nèi)；基于平動(dòng)補(bǔ)償?shù)腟AR成像算法最差，方位向分辨率在5m以內(nèi)。

4 結(jié)束語

本文通過分析ISAR成像算法與SAR成像算法對(duì)運(yùn)動(dòng)艦船目標(biāo)進(jìn)行高分辨成像處理時(shí)面臨的問題，綜合考慮雷達(dá)平臺(tái)高速機(jī)動(dòng)、艦船目標(biāo)非合作，以及成像斜視角大的特點(diǎn)，提出了一種基于加速度模型的混合式SAR/ISAR成像算法。本文從提高成像過程中各個(gè)步驟補(bǔ)償精度的角度出發(fā)，將已有的兩種算法進(jìn)行了有效融合，將回波包絡(luò)對(duì)齊分步進(jìn)行，提高了包絡(luò)對(duì)齊精度，利用SAR處理中的相位濾波去除調(diào)頻率的空變性，結(jié)合ISAR處理中的相位校正方法對(duì)剩余相位項(xiàng)統(tǒng)一校正，最終得到了兩維高分辨的結(jié)果。

[1] WARD K D. Radar imaging of ships at sea[C]. The Royal Signals and Radar Establishment, Gt Malvern, Worcs WR14 3PS, UK, 1989:1–5.

[2] PORTER N J, TOUGH R J A, WARD K D. SAR, ISAR and hybrid SAR/ISAR-a unified treatment of radar imaging[C]. IEEE National Radar Conference, UK, 1993:134–139.

[3] TOUGH R J A, PORTER N J. The hybrid SAR/ISAR imaging of a rocking ship target[C]. Defence Research Agency, Malvern, 1993:4793.

[4] CHEN V C, LIU B. Hybrid SAR/ISAR for distributed ISAR imaging of moving targets[C]. IEEE Radar Conference, 2015:658–663.

[5] 栗君. 基于混合式SAR/ISAR處理的艦船目標(biāo)聚束SAR成像方法研究[D].哈爾濱: 哈爾濱工業(yè)大學(xué), 2015.LI Jun. Study on hybrid SAR/ISAR imaging of ship target based on spotlight SAR[D]. Harbin: Harbin Institute of Technology, 2015.

[6] 李增. 機(jī)載 SAR/ISAR對(duì)艦船成像算法的研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2016.LI Zeng. Study on airborne SAR/ISAR imaging of ship[D]. Harbin: Harbin Institute of Technology, 2016.

[7] 闞學(xué)超. 基于混合式SAR/ISAR的艦船目標(biāo)成像技術(shù)研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2017.KAN Xuechao. The imaging technology research of the ship target based on the hybrid SAR and ISAR[D]. Harbin: Harbin Institute of Technology, 2017.

[8] WHENER D R. High resolution radar[J]. Boston: Artech House, 1995:1–6.

[9] 大斜視高海情下艦船目標(biāo)成像算法研究[D]. 北京: 北京理工大學(xué), 2016. HU Guangli. Research on imaging algorithms of ship targets in large squint mode and atrocious sea condition[D]. Beijing: Beijing Institute of Technology, 2016.

Research on hybrid SAR/ISAR imaging algorithm of ship target

CAO Ye, SHI Yahui, YAN Haipeng, ZHANG Jianqi

（Beijing Research Institute of Telemetry, Beijing 100076, China）

Aiming at the problem of fuzzy imaging of moving ship target under high dynamic platform and large squint angle, this paper proposes a hybrid SAR/ISAR imaging algorithm based on acceleration model from the perspective of improving the accuracy of each part of imaging processing. The effective fusion of the two algorithms is achieved. Firstly, the motion estimation is combined with the distance migration correction to complete the roughly alignment of the echo envelopes, and at the same time, the phase modulation is eliminated by phase filtering in Range-Doppler domain. Then the accurate envelope alignment is achieved by correlation alignment to compensate the errors caused by the motion estimation inaccuracy, the model mismatch and the inertial measurement. Finally, the residual non-space-variant phase term is uniformly corrected by the weighted multi-special point autofocus method. The simulation results show that the proposed algorithm can effectively improve the imaging effect and obtain a well-focused target image.

Hybrid SAR/ISAR; Ship target; Translation compensation; High-resolution imaging

TN958

A

CN11-1780(2019)05-0034-08

Email:ycyk704@163.com TEL:010-68382327 010-68382557

2019-04-15

2019-05-25

曹 曄 1994年生，在讀研究生，主要研究方向?yàn)槔走_(dá)信號(hào)處理。

師亞輝 1977年生，研究員，主要研究方向?yàn)镾AR信號(hào)處理算法及系統(tǒng)設(shè)計(jì)仿真。

閆海鵬 1988年生，碩士，主要研究方向?yàn)槔走_(dá)信號(hào)處理。

張劍琦 1985年生，高級(jí)工程師，主要研究方向?yàn)槔走_(dá)總體技術(shù)。