孟智超, 張 磊,*, 盧景月, 李 軍

(1. 中山大學(xué)深圳電子與通信工程學(xué)院, 深圳 518107; 2. 西安電子科技大學(xué)計(jì)算機(jī)科學(xué)與 技術(shù)學(xué)院, 陜西 西安 710071; 3. 北京無(wú)線電測(cè)量研究所, 北京 100854)

0 引 言

合成孔徑雷達(dá)(synthetic aperture radar, SAR)[1-2]自問(wèn)世以來(lái)便得到了廣泛的關(guān)注。在理想的平臺(tái)飛行條件下,SAR通過(guò)對(duì)均勻采樣的慢時(shí)間回波進(jìn)行相干處理來(lái)獲得方位高分辨的聚焦圖像。然而,對(duì)于空基平臺(tái),尤其是安裝在小型有人飛機(jī)以及更小的無(wú)人機(jī)平臺(tái)上的SAR系統(tǒng)[3-5]而言,受駕駛員控制以及外部大氣擾流的影響,載機(jī)平臺(tái)往往無(wú)法進(jìn)行理想的勻速直線飛行。換言之,實(shí)際的機(jī)載平臺(tái)將不可避免地產(chǎn)生運(yùn)動(dòng)誤差。為了解決這一問(wèn)題,實(shí)際中的SAR系統(tǒng)均配備有記錄其天線相位中心(antenna phase center, APC)實(shí)時(shí)位置的慣性測(cè)量單元(inertial measurement unit, IMU)設(shè)備,而利用高精度的IMU數(shù)據(jù)便可對(duì)運(yùn)動(dòng)誤差實(shí)現(xiàn)精確補(bǔ)償。由于運(yùn)動(dòng)誤差會(huì)影響回波數(shù)據(jù)的多普勒調(diào)制特性,因此IMU設(shè)備的精度必須與SAR系統(tǒng)的工作頻段相匹配。即工作頻段越高,所需的IMU設(shè)備精度要求也越高。一方面,小型無(wú)人機(jī)平臺(tái)受載荷體積與重量限制,往往搭載的是設(shè)備更小的高頻段的(如X波段、Ka波段等)SAR系統(tǒng)。另一方面,為了降低成本,中等精度甚至低精度的IMU設(shè)備往往是這類(lèi)SAR系統(tǒng)的首要選擇。這就使得在利用IMU數(shù)據(jù)進(jìn)行運(yùn)動(dòng)補(bǔ)償后依然存在一部分殘余運(yùn)動(dòng)誤差,而在高分辨SAR成像處理中,該部分運(yùn)動(dòng)誤差仍然會(huì)引起圖像的方位向散焦。因此,需要利用數(shù)據(jù)本身進(jìn)行自聚焦處理[6-9]。

在高分辨SAR成像中自聚焦技術(shù)是保證成像性能的重要技術(shù)之一。目前,大多數(shù)自聚焦技術(shù)都是針對(duì)空不變誤差提出的。如基于代價(jià)函數(shù)優(yōu)化類(lèi)的自聚焦算法,其主要代表有最小熵自聚焦(minimum entropy autofocus, MEA)[10-11],最大對(duì)比度自聚焦(maximum contract autofocus, MCA)[12]等算法。這類(lèi)算法以某一代價(jià)函數(shù)評(píng)價(jià)圖像的聚焦質(zhì)量,并通過(guò)對(duì)該代價(jià)函數(shù)的迭代優(yōu)化來(lái)求解誤差相位。其特點(diǎn)是穩(wěn)健性高,能適應(yīng)不同場(chǎng)景下的自聚焦任務(wù),且不需要選取特顯點(diǎn)。但面臨的問(wèn)題是計(jì)算復(fù)雜度較高,無(wú)法適應(yīng)數(shù)據(jù)量大的自聚焦任務(wù)。相比于代價(jià)函數(shù)優(yōu)化類(lèi)算法,相位梯度類(lèi)算法在實(shí)際應(yīng)用中較為廣泛,其典型代表為相位梯度自聚焦(phase gradient autofocus, PGA)類(lèi)算法[13-14]。PGA類(lèi)算法通過(guò)對(duì)場(chǎng)景中的特顯點(diǎn)進(jìn)行加窗提取,并利用相位差分技術(shù)獲得誤差相位的差分值,而后利用相位誤差的相干特性提高估計(jì)精度,最終通過(guò)對(duì)誤差相位梯度的積分實(shí)現(xiàn)誤差相位的估計(jì)。由于PGA需要選擇一定數(shù)量的強(qiáng)特顯點(diǎn),因而對(duì)特顯點(diǎn)數(shù)量以及質(zhì)量的要求較高。質(zhì)量PGA(qualification PGA, QPGA)[15]通過(guò)二維選取特顯點(diǎn)緩解了對(duì)特顯點(diǎn)數(shù)量的要求,特顯點(diǎn)數(shù)量的增加提高了相位梯度估計(jì)的精度以及算法穩(wěn)健性。而加權(quán)PGA(weighted PGA, WPGA)[16]和相位加權(quán)估計(jì)PGA(phase weighted estimation PGA, PWE-PGA)[17]通過(guò)對(duì)不同特顯點(diǎn)施加不同的權(quán)重,弱化信雜比(signal clutter ratio, SCR)較小的特顯點(diǎn)在相位梯度估計(jì)中的貢獻(xiàn)度而增強(qiáng)高質(zhì)量特顯點(diǎn)的貢獻(xiàn),實(shí)現(xiàn)了更加穩(wěn)健的相位梯度估計(jì)性能。然而,上述提到的各類(lèi)算法,大多是為補(bǔ)償空不變運(yùn)動(dòng)誤差而提出的。實(shí)際中,對(duì)于飛行高度較低而俯仰向波束較寬的高分辨SAR系統(tǒng)而言,運(yùn)動(dòng)誤差往往是距離依賴的[17]。如果利用上述方法,對(duì)距離向分塊進(jìn)行空不變運(yùn)動(dòng)誤差補(bǔ)償,則存在算法低精度以及低效率的問(wèn)題。

為了解決這一問(wèn)題,針對(duì)距離依賴的自聚焦算法越來(lái)越多地受到科研人員的關(guān)注。文獻(xiàn)[18]對(duì)傳統(tǒng)的圖像偏移(map drift, MD)算法[19]進(jìn)行了改進(jìn),作者將MD算法逐距離子塊應(yīng)用,來(lái)估計(jì)線性距離依賴的運(yùn)動(dòng)誤差模型參數(shù),實(shí)現(xiàn)了距離依賴的MD(range-dependent MD, RDMD)算法。另外,在文獻(xiàn)[18]中,作者將PWE-PGA算法應(yīng)用到了距離依賴的誤差相位估計(jì)任務(wù)中,得到了良好的估計(jì)效果。然而,由于距離依賴的誤差相位估計(jì)需要盡可能多的距離單元參與估計(jì)以提高模型參數(shù)的估計(jì)精度。一方面,若在每個(gè)距離單元均提取特顯點(diǎn),傳統(tǒng)PGA算法中的相干疊加優(yōu)勢(shì)將不復(fù)存在,則低SCR樣本可能會(huì)導(dǎo)致模型參數(shù)估計(jì)出現(xiàn)偏差。另一方面,若僅提取固定門(mén)限以上的特顯點(diǎn)樣本則可能存在相當(dāng)一部分距離單元未被利用的情況(這種情況在同時(shí)存在建筑群以及田地的場(chǎng)景中較為常見(jiàn)),導(dǎo)致估計(jì)的模型在特顯點(diǎn)集中的區(qū)域性能好,而在特顯點(diǎn)稀疏的區(qū)域性能差。為了解決這些問(wèn)題,文獻(xiàn)[20]在PWE-PGA的基礎(chǔ)上結(jié)合WPGA加權(quán)策略提出了局部最大似然WPGA(local maximum likelihood WPGA, LML-WPGA)算法。該算法對(duì)待處理數(shù)據(jù)進(jìn)行距離向分塊,并假設(shè)每個(gè)子塊內(nèi)的距離依賴效應(yīng)可以被忽略。在子塊內(nèi),該算法保證每個(gè)距離塊中均有一定數(shù)量的特顯點(diǎn)被選中以提供該距離段的誤差相位特性,并通過(guò)應(yīng)用PGA的相干疊加特性提高子塊內(nèi)的相位梯度估計(jì)精度。在子塊間,該算法通過(guò)對(duì)各個(gè)子塊獲得的相位梯度利用SCR加權(quán)獲得更加穩(wěn)健的相位梯度估計(jì)性能。通過(guò)以上的改進(jìn),LML-WPGA算法實(shí)現(xiàn)了更加精確且穩(wěn)健的自聚焦性能。然而,即使這樣,LML-WPGA算法選擇特顯點(diǎn)的方法依然存在一定的問(wèn)題。在每次迭代過(guò)程中,利用一定的閾值在距離子塊內(nèi)僅能夠選取固定的距離單元樣本。在模型參數(shù)估計(jì)過(guò)程中,存在一部分始終未被選擇的距離單元。在此情況下,子塊相位梯度估計(jì)結(jié)果也會(huì)因樣本不充分的原因產(chǎn)生估計(jì)偏差。因此,為了保證算法的估計(jì)精度以及迭代收斂的穩(wěn)健性,在迭代估計(jì)過(guò)程中,應(yīng)當(dāng)保證盡可能多的距離單元樣本至少有一次被選中以提供該距離單元的誤差相位特性。

為了實(shí)現(xiàn)這一目的,本文針對(duì)LML-WPGA算法在樣本選擇方面的不足進(jìn)行了改進(jìn)。在距離子塊內(nèi)選擇樣本時(shí),不再以固定閾值對(duì)樣本進(jìn)行硬剔除,而是利用樣本的信噪比特性建立對(duì)應(yīng)的選擇概率密度函數(shù),并利用該概率密度函數(shù)以及所需要的樣本數(shù)量,在每次迭代過(guò)程中進(jìn)行樣本隨機(jī)選擇。因此,在迭代過(guò)程中,每個(gè)距離單元樣本均有一定的概率被選擇進(jìn)入模型參數(shù)估計(jì)中,保證了模型參數(shù)的全局最優(yōu)性能。與此同時(shí),在迭代過(guò)程中,信噪比較高的高質(zhì)量樣本總能被多次選擇,以此保證了算法收斂的穩(wěn)健性。由于每次迭代過(guò)程中均是隨機(jī)選擇樣本并估計(jì)誤差相位梯度值,因而本文將這種算法稱(chēng)為局部最大似然加權(quán)隨機(jī)PGA(local maximum like-hood weighted stochastic PGA, LML-WSPGA)方法。最后,本文通過(guò)對(duì)實(shí)測(cè)SAR數(shù)據(jù)的處理結(jié)果驗(yàn)證了所提算法的可行性以及有效性。

1 距離依賴的運(yùn)動(dòng)誤差模型

距離依賴的運(yùn)動(dòng)誤差模型的建立是空變自聚焦的前提,本節(jié)從SAR實(shí)際觀測(cè)幾何出發(fā)建立較為精確的運(yùn)動(dòng)誤差模型,為后續(xù)的求解提供理論基礎(chǔ)。

合成孔徑雷達(dá)的幾何模型如圖1所示,其中雷達(dá)飛行平均高度為H,平均飛行速度為v,波束下視角為θc,參考最短斜距為Rc。圖1中,黑色虛線為理想情況下的雷達(dá)平臺(tái)航跡,而紅色實(shí)線為包含了運(yùn)動(dòng)誤差的實(shí)際雷達(dá)平臺(tái)航跡。在某一時(shí)刻,理想的雷達(dá)APC位置為A(x,0,H),而其實(shí)際位置則為A′(x+Δx,Δy,H+Δz)。其中,Δx,Δy和Δz分別為平臺(tái)的X軸,Y軸和Z軸的運(yùn)動(dòng)誤差,它們是時(shí)間的高階復(fù)雜函數(shù),可由搭載的IMU設(shè)備以一定的精度實(shí)時(shí)測(cè)量獲得。

圖1 SAR成像運(yùn)動(dòng)誤差幾何示意圖Fig.1 Geometry schematic diagram of SAR imaging motion errors

假設(shè)場(chǎng)景中存在一點(diǎn)P(x0,y0),雷達(dá)APC到該點(diǎn)的最短斜距為R0。則在遠(yuǎn)場(chǎng)假設(shè)下,某一時(shí)刻的平臺(tái)三軸運(yùn)動(dòng)誤差所引起的斜距誤差可近似為

ΔR(t)≈Δx(t)·sinθsq(t)+Δy(t)·sinθ0+Δz(t)·cosθ0

(1)

式中:t為慢時(shí)間變量;θ0為雷達(dá)APC相對(duì)于點(diǎn)目標(biāo)P的入射角,且有

(2)

θsq(t)為雷達(dá)APC相對(duì)于點(diǎn)目標(biāo)P的實(shí)時(shí)斜視角,且有

(3)

從式(1)～式(3)中可以看到,運(yùn)動(dòng)誤差ΔR(t)不僅是慢時(shí)間的復(fù)雜高階函數(shù),還與目標(biāo)的斜距R0有關(guān),也即運(yùn)動(dòng)誤差是距離依賴的,這一點(diǎn)在SAR距離向波束較寬且飛行高度較低時(shí)尤為顯著。雖然,在實(shí)際的SAR成像信號(hào)處理過(guò)程中,沿航向的X軸運(yùn)動(dòng)誤差Δx(t)總能通過(guò)插值或者相位補(bǔ)償?shù)确椒ū蝗コ?。但由于IMU設(shè)備測(cè)量精度限制,其殘余分量依然存在。因此,在本文的算法中依然將殘余的X軸運(yùn)動(dòng)誤差Δx(t)視為距離空變誤差的一部分,同時(shí)忽略斜視角中與目標(biāo)方位位置有關(guān)的誤差項(xiàng)。

此時(shí),式(1)可以被重寫(xiě)為

ΔR(t,R0)≈Δx(t,R0)+Δy(t)·sinθ0+Δz(t)·cosθ0

(4)

式中:

(5)

為殘余X軸運(yùn)動(dòng)誤差的距離空變項(xiàng),且有

(6)

式中:θsq為孔徑中心的天線斜視角?？梢钥吹绞?4)在每一個(gè)慢時(shí)間點(diǎn)上都是目標(biāo)斜距R0的函數(shù),因此式(4)總能在R0=Rc處被泰勒展開(kāi)為

ΔR(t,Δr)≈α0Δy(t)+β0Δz(t)+γ0Δx(t,R0)+ (α1Δy(t)+β1Δz(t)+γ1Δx(t,R0))·Δr+ (α2Δy(t)+β2Δz(t)+γ2Δx(t,R0))·Δr2

(7)

式中:Δr=R0-Rc為目標(biāo)斜距相對(duì)于參考斜距的斜距差;系數(shù)(αi,βi)(i=1,2,…,n)可由sinθ0及cosθ0分別在Rc處求導(dǎo)獲得,而γi則與Δx(t,R0)的具體形式相關(guān)。當(dāng)這些系數(shù)被精確已知時(shí),通過(guò)對(duì)誤差相位的估計(jì)便可以利用式(7)獲得三軸運(yùn)動(dòng)誤差的精確估計(jì)值,而這些估計(jì)值可被用于精確的后向投影(back projection, BP)[21-23]成像算法中。然而,在實(shí)際中,Δx(t,R0)的具體形式未知,因而系數(shù)γi較難準(zhǔn)確獲得。但在頻域算法中,我們僅需要獲得式(7)中多項(xiàng)式各個(gè)指數(shù)項(xiàng)的系數(shù),即可對(duì)運(yùn)動(dòng)誤差進(jìn)行精確補(bǔ)償。也即式(7)的形式可以被簡(jiǎn)化為

ΔR(t,Δr)≈ΔR0(t)+ΔR1(t)·Δr+ΔR2(t)·Δr2+o(Δr)

(8)

為了能夠精確地估計(jì)隨距離空變的誤差相位,本文對(duì)式(8)進(jìn)行三階近似。則由式(8)可得,隨距離空變的誤差相位可表示為

(9)

式中:λ為SAR系統(tǒng)工作波長(zhǎng)。

為了清楚地表示式(9),我們將離散化的表達(dá)式用矩陣表示。這里假設(shè)待估計(jì)數(shù)據(jù)的維度為N×M,即認(rèn)為該數(shù)據(jù)具有N個(gè)距離單元和M個(gè)回波脈沖:

(10)

(11)

式中:AN×4為誤差相位模型系數(shù)矩陣,當(dāng)距離單元給定時(shí),矩陣A為一常系數(shù)矩陣,可以被預(yù)先計(jì)算并存儲(chǔ);Δφ4×M為待估計(jì)的誤差相位系數(shù)矩陣。對(duì)每一慢時(shí)間采樣時(shí)刻而言,待求解分量有4個(gè),而方程數(shù)則由所選取的不同距離單元樣本數(shù)決定,且一般情況下遠(yuǎn)大于待求解分量個(gè)數(shù)。因而,對(duì)每一慢時(shí)間采樣時(shí)刻的誤差相位系數(shù)求解都可以等同于求解一個(gè)超定方程組。因此,可以將極小范數(shù)最小二乘解作為每一慢時(shí)間采樣時(shí)刻的誤差相位系數(shù)估計(jì)值,即

(12)

至此,對(duì)第m時(shí)刻的誤差相位系數(shù)求解已完成。為了使求解更加穩(wěn)健,對(duì)該過(guò)程進(jìn)行多次迭代,便可以獲得精確估計(jì)的誤差模型。進(jìn)一步,距離空變的相位誤差補(bǔ)償函數(shù)H(Δr,t)可表示為

(13)

2 LML-WSPGA算法

(14)

式中:假設(shè)距離分塊數(shù)為K,且有

(15)

(16)

式中:假定第k個(gè)距離子塊中所選樣本數(shù)為I個(gè),則Δri,k為第k個(gè)距離子塊中第i個(gè)樣本所在的距離單元差值。而mi,k為該子塊內(nèi)對(duì)應(yīng)樣本的加權(quán)值,該值由SCR決定,具體可見(jiàn)文獻(xiàn)[16]。相應(yīng)地,最終的模型參數(shù)估計(jì)式變?yōu)?/p>

(17)

雖然從式(14)、式(16) 以及式(17)分析可以看到,LML-WPGA算法通過(guò)距離分塊,相位加權(quán)提高了傳統(tǒng)RDPGA算法的求解穩(wěn)健性,但該算法在距離子塊內(nèi)的樣本選擇方式上依然存在問(wèn)題。其僅利用距離單元能量或者對(duì)比度方法選擇給定門(mén)限以上的特顯點(diǎn)樣本的方式,會(huì)導(dǎo)致發(fā)生一部分潛在的特顯點(diǎn)樣本所在距離單元始終未被選擇的情況。這些不被選擇的距離單元始終未參與求解過(guò)程,將使得模型求解產(chǎn)生偏差進(jìn)而導(dǎo)致距離空變誤差相位補(bǔ)償?shù)氖　?/p>

為了避免這一情況的出現(xiàn),并進(jìn)一步提高LML-WPGA算法的穩(wěn)健性,本文對(duì)該算法在距離子塊內(nèi)的樣本選擇方式進(jìn)行了改進(jìn)。通過(guò)引入基于SCR準(zhǔn)則的隨機(jī)樣本挑選方法,本文所提算法不再固定地選擇一部分距離樣本,所有樣本將依由其SCR決定的概率被選擇。在這種情況下,一方面,每個(gè)距離單元在迭代過(guò)程中均有機(jī)會(huì)被選中并參與模型參數(shù)求解,以此來(lái)保證模型參數(shù)求解的精度;另一方面,在迭代過(guò)程中,高SCR樣本將大概率被選擇多次,保證了算法將以較高的穩(wěn)健性收斂。為了區(qū)別于LML-WPGA算法,將本文所提算法命名為L(zhǎng)ML-WSPGA算法,這是因?yàn)檫x擇隨機(jī)樣本計(jì)算誤差相位梯度的方法類(lèi)似于深度學(xué)習(xí)中的隨機(jī)梯度訓(xùn)練方法[24]。

現(xiàn)將本文所提隨機(jī)樣本挑選方法的具體步驟闡述如下。

首先,假設(shè)第k個(gè)距離子塊內(nèi)的距離點(diǎn)數(shù)為J,則以該子塊內(nèi)各距離單元的SCR作為參數(shù)建立樣本的概率密度函數(shù)P(σk)可表示為

(18)

式中:σk,j為第k個(gè)距離子塊內(nèi)第j個(gè)距離單元的SCR。

從式(18)可以看到P(σk)通過(guò)每個(gè)距離單元的SCR決定了該距離單元樣本在模型參數(shù)求解過(guò)程中被選擇的概率大小。當(dāng)給定樣本選擇數(shù)量時(shí),各個(gè)距離單元樣本均有一定的概率被選擇,該方法顯然不同于原始LML-WPGA算法中的固定門(mén)限選擇方法。實(shí)際上可將該方法理解為軟閾值門(mén)限法,其軟閾值的優(yōu)勢(shì)在于保證高質(zhì)量樣本被大概率選擇的同時(shí)使得各個(gè)距離單元的樣本均有機(jī)會(huì)參與模型參數(shù)求解過(guò)程,既保證了算法的穩(wěn)健性又不失精確性。

其次,為了能夠利用式(18)選擇樣本,需要求得概率密度P(σk)的累積分布函數(shù)F(σk),即

(19)

(20)

最后,將所選樣本以及其對(duì)應(yīng)的信雜比、距離單元序號(hào)進(jìn)行記錄并用于該子塊的相位誤差估計(jì),進(jìn)而進(jìn)行模型參數(shù)的估計(jì),便完成了一次迭代估計(jì)過(guò)程。經(jīng)過(guò)數(shù)次迭代該算法便可以得到精確的誤差相位模型參數(shù)估計(jì)值,最終將估計(jì)得到的模型用于逐距離單元的相位補(bǔ)償并進(jìn)行方位聚焦即可得到聚焦的SAR圖像。

為了更加清晰地描述所提算法,將基于所提算法的高分辨SAR成像算法詳細(xì)框圖描述如圖2所示。

圖2 LML-WSPGA算法流程圖Fig.2 Flowchart of LML-WSPGA

3 實(shí)測(cè)數(shù)據(jù)驗(yàn)證

本文利用實(shí)測(cè)機(jī)載SAR數(shù)據(jù)對(duì)所提算法進(jìn)行驗(yàn)證,并在相同迭代次數(shù)下,將其性能與PGA算法以及LML-WPGA算法進(jìn)行比較。錄取數(shù)據(jù)的機(jī)載SAR系統(tǒng)基本參數(shù)如表1所示。

表1 SAR實(shí)測(cè)數(shù)據(jù)基本參數(shù)

另外,我們將IMU測(cè)得的載機(jī)三軸運(yùn)動(dòng)誤差展示在圖3中。一方面,從表1的系統(tǒng)參數(shù)可以看到,平臺(tái)斜視角小且飛行高度低,容易引起運(yùn)動(dòng)誤差的距離依賴現(xiàn)象。另一方面,從圖3的載機(jī)運(yùn)動(dòng)誤差曲線可以看到,載機(jī)的Y軸和Z軸運(yùn)動(dòng)誤差遠(yuǎn)大于X軸,而該機(jī)載IMU設(shè)備的精度較低。因而,可以預(yù)見(jiàn)的是即使經(jīng)過(guò)了基于IMU數(shù)據(jù)的運(yùn)動(dòng)補(bǔ)償,殘余運(yùn)動(dòng)誤差依然會(huì)引起距離依賴的誤差相位,而在該數(shù)據(jù)的實(shí)際處理過(guò)程中也確實(shí)存在這一現(xiàn)象。綜合以上原因,本文利用該數(shù)據(jù)作為驗(yàn)證所提方法的實(shí)驗(yàn)數(shù)據(jù)。

進(jìn)一步,分別利用PGA算法、LML-WPGA算法以及本文所提LML-WSPGA算法對(duì)該數(shù)據(jù)進(jìn)行自聚焦處理,其結(jié)果如圖4所示。首先,從3幅圖像可以看到,PGA算法的處理結(jié)果存在明顯的距離空變問(wèn)題。距離遠(yuǎn)端被良好聚焦而距離近端則出現(xiàn)明顯的模糊現(xiàn)象。不同于PGA算法,LML-WPGA算法及LML-WSPGA均能夠?qū)嚯x依賴的運(yùn)動(dòng)誤差進(jìn)行估計(jì)。從圖4(b)和圖4(c)也可以明顯看出,距離近端和遠(yuǎn)端都能夠良好地聚焦,說(shuō)明了算法的有效性。另外,該3幅圖的圖像熵值分別為15.27、15.14和15.11。從圖像熵值可以看到,本文所提算法較LML-WPGA算法在補(bǔ)償精度上具有更大優(yōu)勢(shì)。為了更加清晰地看到所提算法相比于其他兩種算法的優(yōu)勢(shì),我們對(duì)場(chǎng)景中的3個(gè)典型區(qū)域進(jìn)行了放大比較。這3個(gè)區(qū)域分別位于距離近端、中間以及遠(yuǎn)端,且每個(gè)區(qū)域中均包含特殊結(jié)構(gòu)以更加容易地顯示聚焦差異。

圖4 不同自聚焦算法的處理結(jié)果Fig.4 Results performed by different autofows algorithms

對(duì)子區(qū)域1、子區(qū)域2和子區(qū)域3的放大圖被分別展示在圖5中并由左至右依次排列。從圖5可以看到,由PGA算法處理后的子區(qū)域1和子區(qū)域2依然存在嚴(yán)重散焦,而LML-WPGA算法和LML-WSPGA算法對(duì)3個(gè)子區(qū)域均能良好聚焦。另外,對(duì)比從子區(qū)域1的聚焦情況可以明顯地看出,本文所提LML-WSPGA算法的聚焦性能優(yōu)于LML-WPGA算法。

圖5 不同自聚焦算法的處理子圖結(jié)果比較 (由左至右依次為子區(qū)域1～區(qū)域3)Fig.5 Comparison of sub-images result performed by different autofocus algorithms (from left to right are Areas 1～3)

進(jìn)一步,將場(chǎng)景中的角反射器目標(biāo)方位點(diǎn)擴(kuò)展函數(shù)(point spread function, PSF)以及相應(yīng)的性能指標(biāo)展示在圖6與圖7中。從圖6可以看到,僅利用PGA算法無(wú)法對(duì)角反射器目標(biāo)有效聚焦。而利用本文算法不僅能夠有效聚焦角反射器目標(biāo),且得到的目標(biāo)方位點(diǎn)擴(kuò)展函數(shù)的峰值旁瓣比(peak sidelobe ratio, PSR)明顯低于LML-WPGA算法,兩者分別為-8.51 dB和-6.51 dB,說(shuō)明本文所提算法的聚焦性能優(yōu)于LML-WPGA。另外,從圖7中可以看到,對(duì)于脈沖響應(yīng)寬度(impulse response width, IRW)也即方位分辨率和積分旁瓣比(integrated sidelobe ratio, ISR)而言,本文所提算法均可達(dá)到與LML-WPGA相當(dāng)?shù)某潭?。測(cè)得LML-WSPGA和LML-WPGA的IRW及ISR分別為1.06 m,-1.53 dB和1.06 m,-1.98 dB。

圖6 角反射器目標(biāo)方位點(diǎn)擴(kuò)展函數(shù)Fig.6 PSF of target position of corner reflector

圖7 角反射器目標(biāo)方位聚焦性能指標(biāo)比較Fig.7 Comparison of azimuth focusing performance indicator for corner reflector

此外,為了說(shuō)明所提算法在不同信噪比下的穩(wěn)健性,本文驗(yàn)證了不同信噪比下的算法相位估計(jì)均方誤差值(mean square error, MSE)。

該實(shí)驗(yàn)以一個(gè)運(yùn)動(dòng)誤差被完全補(bǔ)償后的實(shí)測(cè)機(jī)載SAR數(shù)據(jù)為標(biāo)準(zhǔn)數(shù)據(jù),其大小為4 096×4 096點(diǎn)。通過(guò)人為加入距離空變的相位誤差,利用PGA,LML-WPGA及LML-WSPGA算法分別進(jìn)行誤差估計(jì)補(bǔ)償,并分別記錄第129、2 049以及3 969個(gè)距離單元的相位誤差估計(jì)MSE隨(signal to noise ratio, SNR)變化的曲線。實(shí)驗(yàn)中的SNR定義為方位壓縮前的平均SNR,利用AWGN函數(shù)實(shí)現(xiàn)。實(shí)驗(yàn)結(jié)果如圖8所示。可以看到,相比于PGA及LML-WPGA算法,所提算法具有較高的相位估計(jì)精度,且在不同SNR情況下性能較為穩(wěn)健。

圖8 不同算法相位估計(jì)MSE曲線Fig.8 Phase estimation MSE curves of different algorithms

綜上,在相同迭代次數(shù)下,本文所提算法相比于LML-WPGA算法具有更加優(yōu)異的全局聚焦性能,能夠?qū)嚯x依賴的運(yùn)動(dòng)誤差進(jìn)行較好的自聚焦補(bǔ)償。在不同SNR情況下具有穩(wěn)健的估計(jì)性能。因此,所提算法有望在低飛行高度、寬俯仰波束的SAR系統(tǒng)中得到實(shí)際應(yīng)用。

4 結(jié) 論

一直以來(lái),在低飛行高度、寬俯仰波束合成孔徑雷達(dá)成像中,距離依賴的運(yùn)動(dòng)誤差自聚焦問(wèn)題始終存在著算法補(bǔ)償精度有限與穩(wěn)健性不高的問(wèn)題。針對(duì)該問(wèn)題,本文通過(guò)對(duì)原有LML-WPGA算法的樣本選擇方法進(jìn)行了改進(jìn),提出了一種隨機(jī)樣本選擇的LML-WSPGA算法,其本質(zhì)相當(dāng)于一種軟閾值樣本選擇方法。該算法通過(guò)對(duì)特顯點(diǎn)樣本的隨機(jī)選擇,保證了相位梯度計(jì)算中的距離樣本完備性。在確保算法補(bǔ)償性能優(yōu)異的同時(shí),也提高了算法的穩(wěn)健性。