層析SAR系統(tǒng)基線優(yōu)化設計

盧紅喜劉宏偉 羅 濤 索志勇 糾 博 保 錚

(西安電子科技大學雷達信號處理國家重點實驗室 西安 710071)

層析SAR系統(tǒng)基線優(yōu)化設計

盧紅喜*劉宏偉 羅 濤 索志勇 糾 博 保 錚

(西安電子科技大學雷達信號處理國家重點實驗室 西安 710071)

自然場景3維影像重構(gòu)是層析合成孔徑雷達(TomoSAR)的重要應用。傳統(tǒng)方法在沿高程維進行層析處理時，均通過對等距線陣模型進行加權(quán)的方式，以主瓣展寬為代價來抑制方向圖旁瓣水平。針對該問題，該文建立一種基于非等距線陣的峰值旁瓣比極小極大優(yōu)化模型，即在約束主瓣寬度一定的情況下，通過陣元位置的優(yōu)化配置來獲取觀測場景視角范圍內(nèi)任意指向的最優(yōu)旁瓣水平；提出一種目標函數(shù)離散柵格化方法，進而采用序列二次規(guī)劃(SQP)方法并結(jié)合差分進化算法特有的全局記憶能力，以獲取最優(yōu)陣元配置方案。對PolSARpro全極化層析SAR仿真數(shù)據(jù)的處理結(jié)果表明該方法能夠有效地應用于自然場景的3維影像重構(gòu)。

層析合成孔徑雷達；基線優(yōu)化設計；極小極大最優(yōu)化；序列二次規(guī)劃

1 引言

層析合成孔徑雷達(SAR Tomography, Tomo SAR)3維成像，具有傳統(tǒng)SAR全天時全天候、高分辨率及電磁穿透等優(yōu)勢，能夠?qū)ψ匀粓鼍斑M行3維重構(gòu)，是復雜空間結(jié)構(gòu)檢測與識別技術(shù)的基礎，在地球遙感、海洋監(jiān)測、軍事偵察、環(huán)境保護以及災情監(jiān)測等方面具有廣泛應用前景[1,2]。目前非等距線陣設計方法主要可分為兩大類，一是陣元間距為半波長整數(shù)倍的非均勻線陣，如最小冗余陣列[3]、最大連續(xù)延遲線陣[4]和最小間隙陣列[5]等，此類方法陣元間距要求嚴格使其實際應用受到了較大限制；二是在陣列孔徑和陣元數(shù)確定等特定應用背景下的陣列優(yōu)化方法，如利用智能仿生學優(yōu)化算法來實現(xiàn)優(yōu)化布陣[6?8]，此類方法求解結(jié)果受初始值影響較大，無法保證全局最優(yōu)；文獻[9]提出一種基于插值的大規(guī)模非線性陣方向圖快速優(yōu)化方法，然而該方法無法解決局部視角范圍內(nèi)的方向圖旁瓣優(yōu)化問題。針對這些問題，本文將層析SAR系統(tǒng)的基線設計構(gòu)造為非等距線陣峰值旁瓣比的極小極大最優(yōu)化問題，即在陣列孔徑和陣元數(shù)確定的情況下，考慮陣列在觀測視角范圍內(nèi)各個方向的聚焦性能，以陣元位置為自變量，以陣列指向觀測范圍內(nèi)各個方向目標點形成的方向圖最大峰值旁瓣比(Peak Side Lobe Ratio, PSLR)為代價函數(shù)，建立非線性數(shù)學優(yōu)化模型。同時，提出一種目標函數(shù)離散柵格化方法，將原問題轉(zhuǎn)換為極小極大優(yōu)化問題，再結(jié)合差分進化算法(Differential Evolution, DE)特有的記憶能力，采用序列二次規(guī)劃方法(Sequence Quadratic Programming, SQP)獲得最優(yōu)陣元位置配置方案。為了驗證該方法應用于3維重構(gòu)的有效性，針對機載單天線SAR重復航過[10]系統(tǒng)，本文采用歐空局PolSARpro軟件獲取的層析SAR仿真數(shù)據(jù)進行了試驗。仿真試驗結(jié)果表明，在陣列孔徑及陣元數(shù)一定的情況下，本文算法可以穩(wěn)健地獲取最優(yōu)陣元位置配置方案，能夠有效地抑制旁瓣水平并獲取觀測場景的3維影像。

2 層析SAR系統(tǒng)基本原理與數(shù)據(jù)獲取體制

層析SAR系統(tǒng)利用合成孔徑原理獲得沿航向高分辨能力，利用脈沖壓縮技術(shù)實現(xiàn)距離向高分辨能力，利用波束形成技術(shù)獲取垂直距離向(cross range，即斜平面法向)的高分辨能力。圖1所示為層析SAR系統(tǒng)雷達幾何示意圖，雷達平臺高度為H，雷達下視角為θ，層析SAR系統(tǒng)N個天線沿垂直航向呈直線分布，形成有效孔徑長度為L的線性陣列，線陣傾角為α(以下考慮有效基線的情況，即α=θ)。層析SAR系統(tǒng)通過距離向脈沖壓縮和方位向合成孔徑處理分別對各副天線接收的信號進行成像，對于各“方位-斜距”分辨單元，聯(lián)合N副天線獲取的SAR圖像數(shù)據(jù)，通過層析處理技術(shù)(如數(shù)字波束形成)可區(qū)分該分辨單元中沿垂直距離向分布的各個散射目標，從而獲取“方位-斜距-垂直距離向”分辨單元對應的目標散射特征，實現(xiàn)3維空間的高分辨成像。

層析處理需要大量不同視角的空間采樣，其數(shù)據(jù)獲取體制主要包括單天線SAR系統(tǒng)重復航過、多天線InSAR系統(tǒng)重復航過和天線陣列SAR系統(tǒng)單航過等。天線陣列SAR系統(tǒng)采用天線陣列單次航過進行觀測即可獲取層析SAR數(shù)據(jù)，但系統(tǒng)設計難度和成本較大，對于機載情況天線陣列孔徑也難以提高。因此，為了降低系統(tǒng)設計的成本和復雜性，空間采樣的增加主要通過重軌觀測實現(xiàn)，即采用單天線SAR系統(tǒng)重復航過或多天線InSAR系統(tǒng)重復航過體制。

3 層析SAR系統(tǒng)基線優(yōu)化設計

層析SAR系統(tǒng)的基線設計主要考慮垂直距離向分辨率和不模糊觀測范圍兩方面的問題。為了獲得垂直距離向的高分辨能力，必須保證足夠大的垂直航向孔徑長度；為了保證觀測場景內(nèi)不產(chǎn)生模糊，又要求陣元間距足夠小。在滿足所需分辨力要求的情況下，配置過多的陣元會導致數(shù)據(jù)冗余并增加不必要的系統(tǒng)載荷或運行成本。因此，應盡量減少陣元數(shù)量(觀測次數(shù))，但這會引起所需的觀測角度內(nèi)出現(xiàn)柵瓣和高旁瓣。針對這些問題，本文在選定天線陣列長度L(即保證所需的空間角分辨率，約為λ/2L )的前提下，研究了層析SAR系統(tǒng)的基線優(yōu)化設計方法，在陣元數(shù)量確定的情況下，通過陣元位置的優(yōu)化設計抑制觀測視角范圍內(nèi)的柵瓣和高旁瓣，以得到更好的陣列指向性及聚焦性能。

層析SAR系統(tǒng)垂直距離向聚焦性能主要由沿垂直航向分布的天線陣列方向圖特性決定。當主瓣不展寬或展寬系數(shù)較小時，天線陣列方向圖在觀測視角范圍內(nèi)峰值旁瓣比越低，陣列系統(tǒng)指向性以及聚焦性能也就越好。對于確定的天線陣列構(gòu)型，陣列方向圖具有確定的空間分布形式。根據(jù)波束空間傳播理論，如圖2所示，斜距分辨單元ri中的目標點Di處與斜距分辨單元r1中的目標點D1處具有相同的角分辨能力，即其天線陣列方向圖的PSLR一致。因而，斜距分辨單元ri中各方向的聚焦性能可由斜距分辨單元r1(稱為參考距離單元)中相應方向的聚焦性能進行表征。如圖2所示，參考距離單元中AC段是考慮了天線陣列系統(tǒng)波束同步后，測繪帶場景范圍AB段沿垂直距離向的投影(目標點Di的投影為D1)。在實際處理中，可選取場景中心的斜距單元作為參考。沿垂直距離向的層析處理在2維合成孔徑成像處理基礎之上進行，因而只需抽取各陣元所錄取的對應于同一“方位-距離”分辨單元的信號進行高程維聚焦即可。基于這一前提，對于如圖2中所示的聚焦方向，綜合考慮雷達視線跨越的各個距離分辨單元，進行層析處理所需考慮的方向圖柵瓣模糊及旁瓣抑制范圍僅在圖中標示的矩形區(qū)域，即認為陣列接收的數(shù)據(jù)僅僅來源于該矩形區(qū)域范圍，而該區(qū)域以外沒有回波，圖中ε=hsin θ，h為觀測場景深度。根據(jù)以上分析，非等距線陣近場陣列基線優(yōu)化設計的雷達幾何等效模型可簡化為如圖3所示。

圖1 層析SAR系統(tǒng)雷達幾何示意圖

圖2 非均勻線陣近場陣 列方向圖優(yōu)化幾何模型

圖3 層析SAR系統(tǒng)基線優(yōu)化 設計雷達幾何等效模型

層析SAR系統(tǒng)的最優(yōu)基線配置可通過對指向參考距離單元中各方向目標點的陣列方向圖的聯(lián)合優(yōu)化進行實現(xiàn)。如圖3所示，參考距離單元中考慮陣列系統(tǒng)波束同步的觀測范圍為AC段，L為孔徑長度，xn為陣元位置，Hs為陣列到場景中心處的垂直距離。基于這一雷達幾何等效模型，本文將層析SAR系統(tǒng)的基線優(yōu)化設計構(gòu)造為極小極大最優(yōu)化問題，即在陣列有效孔徑長度和陣元數(shù)量確定的情況下，綜合考慮觀測場景范圍內(nèi)各個方向的聚焦性能，以各陣元的相對位置為自變量，以陣列指向參考距離單元中各方向目標點形成的天線方向圖最大PSLR為代價函數(shù)，建立多元非線性數(shù)學優(yōu)化模型，并通過非線性最優(yōu)化方法求解最優(yōu)的基線配置方案。若不考慮時間去相干、航跡控制偏差以及航跡擾動等影響因素，單天線SAR系統(tǒng)重復航過和天線陣列SAR系統(tǒng)(可等效為各陣元自發(fā)自收)單航過在層析數(shù)據(jù)獲取上并沒有本質(zhì)的區(qū)別。然而，實際中天線陣列SAR系統(tǒng)設計難度大、成本高，且陣列長度難以提高，無法保證所需的空間分辨率。因此，下文主要針對機載SAR系統(tǒng)重復航過體制探討和分析基線優(yōu)化設計方法。

下面以機載單天線SAR重軌系統(tǒng)為例進行分析，暫不考慮因重軌引起的時間去相干對各幅圖像間相干性的影響，且所形成的有效孔徑長度小于極限基線長度[11]。如圖3所示的雷達幾何等效模型，坐標系3軸分別為方位向(垂直紙面向外)、垂直距離向和雷達視線(反)方向，Hs為雷達平臺等效高度(沿雷達中心視線方向)，L為陣列長度，陣元個數(shù)為N，其坐標分別記為(xn,Hs), n=0, 1,…, N?1。以陣列中心點為參考，第1個陣元坐標為x0=?L /2，第N個陣元的坐標為xN?1=L/2，陣元坐標向量記為x=[x0, x1,…, xN?1]T。圖中AC段為波束同步觀測的場景視角范圍，記?0=[ζa,ζc]。考慮觀測場景中某一方向點目標D，坐標記為(ζd, 0)，則點目標D對應的數(shù)據(jù)接收導向矢量為

a (ζd)

式中，k=2πλ為電磁波波數(shù)，λ為波長。Δrn(ζd)=為目標到第n個陣元和陣列中心點的單程波程差。由于各天線為自發(fā)自收模式，導向矢量相位對應為雙程波程差。

將天線陣列波束形成指向點目標D，即以點目標D的陣列接收數(shù)據(jù)導向矢量作為權(quán)值，則天線陣列方向圖可表述為

式中，ζ∈?0, ?ζd∈?0。其中，x0=?L2, xN?1= L/2, xn∈(?L/2, L/2), n=1, 2,…, N ?2。利用泰勒展開式1+x2≈1+x2/2, 0＜x?1，則式(2)可表示為

通過典型參數(shù)仿真分析驗證，式(3)的近似誤差約為-30 dB，滿足工程精度要求。

由式(2)可見，對于確定的雷達工作頻率和平臺高度Hs，指向位置ζd處點目標D的陣列波束形成與陣元位置x密切相關(guān)，陣列方向圖PSLR是各個陣元位置x的函數(shù)。在等距線陣情況下，陣元間距d必須滿足奈奎斯特(Nyquist)準則，以抑制高旁瓣和柵瓣模糊對層析圖像的影響。對于非等距線陣，在觀測場景高度范圍內(nèi)不出現(xiàn)柵瓣模糊的有效(即沿垂直雷達視線方向投影)最大陣元間距必須滿足[12]

式中，r0為雷達觀測距離，h為場景深度，θ為雷達視角，hn為待觀測場景沿垂直距離向的高度。因此所需的最少觀測次數(shù)為

對于確定的陣列長度和陣元數(shù)量，在約束主瓣寬度一定的情況下，通過陣元的非等距配置可在最大程度上保證陣列的指向性，并減小其旁瓣對垂直距離向聚焦質(zhì)量的影響。因此，本文綜合考慮指向觀測場景內(nèi)各個目標單元方向的陣列方向圖，構(gòu)建以陣列指向各方向目標點形成的天線方向圖最大PSLR為代價函數(shù)，以陣元位置x為自變量的多元非線性數(shù)學優(yōu)化模型，再利用非線性優(yōu)化方法求解得到最優(yōu)的基線配置方案。結(jié)合圖3所示的雷達幾何等效模型，對于觀測場景內(nèi)任意的目標單元位置ζd，要求在目標單元的鄰域[ζd?ε, ζd+ε]內(nèi)相應的方向圖PSLR最低，則該數(shù)學模型可描述為

其中，?0=[ζa,ζc]為測繪帶觀測區(qū)間，?1=[ζd?ε, ζd+ε]為待聚焦層析單元所需考慮的實際觀測范圍，ε=h/sin θ。ζ3dB為方向圖3 dB主瓣寬度，模型中約束其低于系統(tǒng)垂直距離向分辨率δhn。目標函數(shù)F(x)為指向場景內(nèi)各方向目標單元陣列方向圖的最大峰值旁瓣比，自變量x為陣元位置向量。PSLR(·)表示在給定的觀測區(qū)間?1上的峰值旁瓣比算子，其它參數(shù)含義同上。針對式(6)進行變量替換，上述模型可簡化為

式中，ζ'=ζ?ζd, ?'=[?ε,ε]。

由式(7)可見，式(6)所描述的場景范圍?0內(nèi)陣列方向圖峰值旁瓣比的極小極大問題，最終可等效為將波束形成指向目標點ζ'=0時(而與ζd無關(guān))的陣列方向圖在觀測區(qū)間[?ε, ε]內(nèi)的峰值旁瓣比最小化問題。根據(jù)式(7)陣列方向圖的對稱性，峰值旁瓣比最小化的尋優(yōu)過程可在觀測區(qū)間[0, ε]內(nèi)進行。然而，由于目標函數(shù)具有較復雜的非線性形式，PSLR的表達形式無法解析求取。因此，本文首先根據(jù)系統(tǒng)分辨率δhn要求所需的主瓣寬度ζ3dB計算出方向圖的主瓣第1零陷位置ζ0，則峰值旁瓣比最小化的尋優(yōu)過程可在觀測區(qū)間[ζ0, ε]內(nèi)進行，即目標函數(shù)可表示為式中，?=[ζ0,ε]。由于第1零陷位置ζ0由系統(tǒng)分辨率δhn限定，式(8)等效于在目標函數(shù)中對方向圖的主瓣寬度進行了約束，即ζ3dB≤δhn。

由于方向圖具有非凸性質(zhì)，旁瓣Gx(ζ)峰值無法通過常規(guī)凸優(yōu)化方法求解。因此，本文提出了一種方向圖PSLR的離散解析表達方法，即以主瓣寬度的十分之一為間隔對ζ進行采樣，經(jīng)離散柵格化得到Gx(ζ)的離散形式Gx(ζm)，則模型式(7)可轉(zhuǎn)換為極小極大問題

其中，ζm∈?=[ζ0,ε], m =1, 2, …,M。ζ0為主瓣第1零陷位置，M為離散采樣點數(shù)。

對于有約束的極小極大問題，目前已有成熟算法進行求解，如SQP方法[13]，文獻[14]也給出了一種有效的超線性收斂方法。然而，這些較為成熟的算法依然無法避免收斂到局部最小的情況，因此本文在運用SQP算法求解式(9)的同時，還結(jié)合了差分進化算法的全局搜索能力，以獲取精確的全局最優(yōu)解。DE算法基于差分的簡單變異操作和一對一的競爭生存策略，降低了遺傳操作的復雜性并保留了基于種群的全局搜索策略。同時，DE特有的記憶能力使其可以動態(tài)跟蹤當前的搜索情況，以調(diào)整其搜索策略，具有較強的全局收斂能力和魯棒性[15]。本文首先隨機初始化若干不同的初始值，并通過SQP方法解出由各初始值得到的極小點。然后，以這些極小點作為粒子，再增加若干隨機粒子以增強其多樣性，作為DE算法的粒子初始集，利用DE算法基于種群的全局搜索策略近似逼近以得到式(7)的最優(yōu)解。

以第4節(jié)表1的仿真參數(shù)為例，結(jié)合式(5)可知，當陣列指向觀測場景范圍內(nèi)任意目標單元ζd時，其鄰域[ζd?ε,ζd+ε]內(nèi)不出現(xiàn)柵瓣所需的最少觀測次數(shù)為N=12，此時由模型式(9)得到的非等距配置最優(yōu)解也具有高達-15.2 dB的旁瓣水平。通過增加觀測次數(shù)N并結(jié)合模型式(9)可在約束主瓣寬度不變的情況下進一步降低旁瓣水平。圖4給出了表1中雷達幾何參數(shù)下，場景深度分別為h=10 m, h= 30 m和h=50 m時，場景觀測范圍內(nèi)陣列方向圖最優(yōu)PSLR隨航過次數(shù)(即觀測次數(shù)N)的變化曲線。如圖4所示，當觀測次數(shù)增加到時，對應的最優(yōu)PSLR不再有明顯的改善。根據(jù)表1參數(shù)的仿真結(jié)果，隨著觀測次數(shù)的增加，最優(yōu)基線配置中陣元空間位置發(fā)生重合的陣元數(shù)也隨之增加，而不同的空間位置(以下稱為有效空間觀測位置)的陣元數(shù)量Neff也趨于穩(wěn)定，有Neff≈h/λr0sin θ圖5給出了基線配置最優(yōu)時有效空間觀測位置數(shù)隨著實際觀測次數(shù)的變化曲線。當陣元空間位置發(fā)生重合時，由于陣元數(shù)目確定，相當于對陣列進行整數(shù)加權(quán)約束。因此，雖然有效空間觀測位置數(shù)趨于穩(wěn)定，但有效空間觀測位置會有所不同。對于位置發(fā)生重合的陣元，其接收信號以加權(quán)的形式疊加到陣列波束形成的輸出信號中，相應的權(quán)值即為其位置重合的數(shù)量。考慮陣元加權(quán)的陣列基線配置最優(yōu)化模型可直接根據(jù)模型式(7)進行擴展，因此本文不作深入探討。

在實際中，單天線SAR系統(tǒng)重復航過受時間去相干影響(尤其是星載情況重訪周期較長)嚴重，且重復飛行成本較高(如機載情況)。多天線InSAR系統(tǒng)能夠很好地彌補這些缺陷，單次航過多副天線既能獲取相干性較好的數(shù)據(jù)，還能通過較少的重復航過實現(xiàn)天線陣列孔徑提升，其系統(tǒng)設計載荷輕、成本小、難度低(例如機載三天線InSAR系統(tǒng))，也大大地降低了飛行成本。相比較于單天線SAR系統(tǒng)情況，多天線InSAR重復航過系統(tǒng)增加了對陣元位置的約束，即多天線InSAR系統(tǒng)各天線的相對位置確定，而各次航過觀測基線長度可變。此外，由于多天線InSAR系統(tǒng)通常為一發(fā)多收模式，單次飛行各陣元間接收信號的波程差為單程，而各次航過間各參考陣元接收信號波程差為雙程。因此，在模型設計時還需要將一發(fā)多收模式通過等效相位中心原理[16]等效為自發(fā)自收模式，即各陣元間的波程差為雙程，而各基線長度也對應為等效之后的有效基線長度。

4 實驗與結(jié)果分析

本節(jié)利用歐空局研發(fā)的PolSARpro仿真軟件驗證了本文提出的層析SAR系統(tǒng)基線優(yōu)化設計方法，并對仿真得到的L波段全極化層析SAR數(shù)據(jù)進行成像處理[12]，驗證了所提算法應用于層析SAR系統(tǒng)3維成像的有效性。根據(jù)表1給出的層析SAR系統(tǒng)及觀測場景基本參數(shù)，求解的最優(yōu)基線配置結(jié)果由表2給出，其中各次航過對應的觀測位置坐標系統(tǒng)對應于圖3所示的雷達幾何坐標系。

表1 雷達幾何及場景仿真參數(shù)

表2 最優(yōu)基線配置下空間觀測相對位置

圖4 不同場景深度時最優(yōu)峰值旁瓣比隨陣元數(shù)變化曲線

圖5 不同場景深度時有效空間觀測位置數(shù)隨實際陣元數(shù)變化曲線

根據(jù)表1和表2所提供的層析SAR系統(tǒng)及仿真場景參數(shù)，即可利用歐空局研發(fā)的PolSARpro仿真軟件計算全極化層析SAR數(shù)據(jù)。本試驗中設置的仿真自然場景為包含164棵樹的灌木叢林，樹高Hv為21 m、種植密度約為48棵/公頃。圖8為仿真場景的下視投影圖像，圖中給出了仿真場景中的灌木樹高度、樹冠大小及其生長分布，其中灰度值表示每棵樹木的實際高度，灰度條為場景內(nèi)樹木的高度值范圍。圖9為VV極化通道的層析成像結(jié)果方位向剖面圖，其灰度條表示層析圖像的像素灰度值范圍，該方位向剖面位置由圖8中的虛線切剖線所示。

圖6 最優(yōu)基線配置陣列分別對準不同方向的方向圖

圖7 非等距陣列方向圖優(yōu)化算法性能對比

圖8 試驗仿真場景下視投影圖像及剖面線位置

圖9 SAR層析成像結(jié)果方位向剖面圖

為便于3維成像結(jié)果對比分析，圖10給出了VV極化通道的SAR主圖像(灰度條表示該圖像像素的灰度值)，圖11為試驗仿真自然場景在高程值為Hv=21 m處的實際剖面情況(灰度值為該樹木的實際高度)，圖12為全極化層析SAR 3維成像結(jié)果在高程值為Hv=21 m處的層析偽彩色合成圖像剖面圖，圖中分別標出了3排灌木的真實空間分布和相應的成像結(jié)果。由此可見，本文提出的層析SAR系統(tǒng)基線優(yōu)化設計方法能夠有效地應用于自然場景3維影像重構(gòu)。

5 結(jié)束語

本文針對層析SAR系統(tǒng)的基線設計問題，建立了機載平臺重復航過系統(tǒng)基線優(yōu)化配置的數(shù)學模型，結(jié)合數(shù)學優(yōu)化方法給出了層析系統(tǒng)的基線最優(yōu)配置方案。在系統(tǒng)基本參數(shù)給定的情況下，本文算法首先通過構(gòu)建以陣列指向各方向目標點形成的天線方向圖最大PSLR為代價函數(shù)，在約束主瓣寬度一定的情況下，以陣元位置為自變量的多元非線性數(shù)學優(yōu)化模型，再利用非線性優(yōu)化方法求解得到最優(yōu)的基線配置方案。仿真試驗處理結(jié)果表明，在陣列孔徑及陣元數(shù)一定的情況下，本文算法通過對陣元位置進行優(yōu)化配置，能夠穩(wěn)健地抑制旁瓣水平的影響，有效地獲取觀測場景的3維影像。

圖10 SAR主圖像(VV極化通道)

圖11 試驗仿真場景高程(Hv=21 m)剖面圖

圖12 全極化層析SAR 3維成像高程(Hv=21 m)剖面圖

[1] Dinh Ho Tong M, Thuy Le T, Rocca F, et al.. Relating P-band synthetic aperture radar tomography to tropical forest biomass[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(2): 967-979.

[2] 王金峰, 皮亦鳴, 曹宗杰. 一種機載SAR層析三維成像算法[J]. 電子與信息學報, 2010, 32(5): 1029-1033.

Wang Jin-feng, Pi Yi-ming, and Cao Zong-jie. An algorithm for airborne SAR tomography 3D imaging[J]. Journal of Electronics & Information Technology, 2010, 32(5): 1029-1033.

[3] Moffet A T. Minimum-redundancy linear arrays[J]. IEEE Transactions on Antennas and Propagation, 1968, 16(2): 172-175.

[4] Chambers C, Tozer T C, and Sharman K C. Temporal and spatial sampling influence on the estimates of superimposed narrowband signals: when less can mean more[J]. IEEE Transactions on Signal Processing, 1996, 44(12): 3085-3098.

[5] Abramovich Y I, Spencer N K, and Gorokhov A Y. Positive-definite Toeplitz completion in DOA estimation for nonuniform linear antenna arrays[J]. IEEE Transactions on Signal Processing, 1999, 47(6): 1502-1521.

[6] Ahmad A, Behera A K, and Mandal S K. Artificial bee colony algorithm to reduce the side lobe level of uniformly excited linear antenna arrays through optimized element spacing[C]. Proceedings of the IEEE Conference on Information & Communication Technologies 2013 (ICT 2013), Casablanca, Morocco, 2013: 1029-1032.

[7] Yang S H and Kiang J F. Optimization of asymmetrical difference pattern with memetic algorithm[J]. IEEE Transactions on Antennas and Propagation, 2014, 62(4): 2297-2302.

[8] Joshi P and Dubey R. Optimization of linear antenna array using genetic algorithm for reduction in Side lobs levels and improving directivity based on modulating parameter M[J]. Optimization, 2013, 1(7): 1475-1482.

[9] Yang K, Zhao Z, and Liu Q H. Fast pencil beam pattern synthesis of large unequally spaced antenna arrays[J]. IEEE Transactions on Antennas and Propagation, 2013, 61(2): 627-634.

[10] 鐘雪蓮, 向茂生, 郭華東, 等. 機載重軌干涉合成孔徑雷達的發(fā)展[J]. 雷達學報. 2013, 2(3): 367-381.

Zhong Xue-lian, Xiang Mao-sheng, Guo Hua-dong, et al.. Current development in airborne repeat-pass interferometric synthetic aperture radar[J]. Journal of Radars, 2013, 2(3): 367-381.

[11] Gatelli F, Guamieri A M, Parizzi F, et al.. The wavenumber shift in SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(4): 855-865.

[12] Reigber A. Airborne polarimetric SAR tomography[D]. [Ph.D. dissertation], University of Stuttgart, 2001.

[13] Han S P. A globally convergent method for nonlinear programming[J]. Journal of Optimization Theory and Applications, 1977, 22(3): 297-309.

[14] 薛毅. 求解Minimax優(yōu)化問題的SQP方法[J]. 系統(tǒng)科學與數(shù)學, 2002, 22(7): 355-364.

Xue Yi. An SQP method for minimax optimazation problem[J]. Journal of Systems Science and Mathematical Sciences, 2002, 22(7): 355-364.

[15] Storn R and Price K. Differential evolution - a simple and efficient adaptive scheme for global optimization over continuous spaces[R]. Berkeley: University of California, 2006.

[16] 李真芳. 分布式小衛(wèi)星SAR-InSAR-GMTI的處理方法[D]. [博士論文], 西安電子科技大學, 2006.

Li Zhen-fang. Approaches to SAR-InSAR-GMTI for distributed small satellite SAR systems[D]. [Ph.D. dissertation], Xi'dian University, 2006.

盧紅喜： 男，1987年生，博士生，研究方向為極化InSAR信號處理、多基線InSAR信號處理、陣列優(yōu)化設計與極化層析成像技術(shù)研究.

劉宏偉： 男，1971年生，教授，博士生導師，雷達信號處理國家重點實驗室主任，研究方向為自適應信號處理、雷達信號處理、雷達目標識別等.

羅 濤： 男，1982年生，博士生，研究方向為波形設計與MIMO雷達.

索志勇： 男，1979年生，副教授，研究方向為SAR成像和InSAR數(shù)據(jù)處理技術(shù).

Optimal Baseline Design for SAR Tomography System

Lu Hong-xi Liu Hong-wei Luo Tao Suo Zhi-yong Jiu Bo Bao Zheng
(National Key Laboratory of Radar Signal Processing, Xidian University, Xi'an 710071, China)

Nowadays 3-D reconstruction for natural scene is an important aspect of the Earth observation with SAR Tomography (TomoSAR). The general method for pattern sidelobe suppression, during the tomographic processing, is usually implemented with the weighting of uniform linear array, however at the expense of mainlobe broadening. In this paper, a minimax optimization model for cells configuration is constructed based on non-uniform linear array, to achieve the optimal peak sidelobe ratio with a fixed mainlobe width for any beam direction in the range of perspectives. For this, an objective function rasterisation is proposed and then the optimal solution can be figured out by the Sequence Quadratic Programming (SQP) with differential evolution for its unique ability of global memory. Finally, experimental results with PolSARpro polarimetric TomoSAR data validate the effectiveness of the proposal for natural scene 3-D image reconstruction.

SAR Tomography (TomoSAR); Optimal baseline design; Minimax optimization; Sequence Quadratic Programming (SQP)

TN957.51

： A

：1009-5896(2015)04-0919-07

10.11999/JEIT140710

2014-05-27收到，2014-07-31改回

國家自然科學基金(61271291, 61201285, 61271024, 61201292, 61201283)，新世紀優(yōu)秀人才支持計劃(NCET-09-0630)，全國優(yōu)秀博士學位論文作者專項資金(FANEDD-201156)和中央高校基本科研業(yè)務費專項資助課題

*通信作者：盧紅喜 xdkd@163.com