陳國(guó)忠,趙 迪,王曉鵬,侯雨生,陳 亮,張德新

(1.上海衛(wèi)星工程研究所,上海 201109; 2.上海航天技術(shù)研究院,上海 201109;3.電子科技大學(xué) 機(jī)械電子工程學(xué)院,四川 成都 610054; 4.上海交通大學(xué) 航空航天學(xué)院,上海 200240)

干涉SAR多基線分析與設(shè)計(jì)

陳國(guó)忠1,趙 迪1,王曉鵬2,侯雨生1,陳 亮3,張德新4

(1.上海衛(wèi)星工程研究所,上海 201109; 2.上海航天技術(shù)研究院,上海 201109;3.電子科技大學(xué) 機(jī)械電子工程學(xué)院,四川 成都 610054; 4.上海交通大學(xué) 航空航天學(xué)院,上海 200240)

為對(duì)干涉合成孔徑雷達(dá)(InSAR)進(jìn)行優(yōu)化設(shè)計(jì)以提高系統(tǒng)的測(cè)高精度，對(duì)干涉SAR多基線設(shè)計(jì)進(jìn)行了分析和研究。用多基線最大似然相位估計(jì)方法，比較干涉相位誤差最小和相位解纏穩(wěn)健性最大兩種準(zhǔn)則，確定選擇后者作為多基線系統(tǒng)的基線比例設(shè)計(jì)準(zhǔn)則。給出了多基線長(zhǎng)度優(yōu)化設(shè)計(jì)方法：先設(shè)計(jì)多基線系統(tǒng)中的長(zhǎng)基線，基于單基線對(duì)測(cè)高性能的傳遞模型，按測(cè)高精度要求反推長(zhǎng)基線的長(zhǎng)度范圍；再以多基線相位解纏的穩(wěn)健性最大化為準(zhǔn)則進(jìn)行基線比例尋優(yōu)，獲得短基線長(zhǎng)度范圍，解決長(zhǎng)基線的相位解纏。給出了針對(duì)山區(qū)地形的一個(gè)多基線InSAR系統(tǒng)設(shè)計(jì)樣例，結(jié)果表明：采用最大似然法，通過(guò)短基線輔助長(zhǎng)基線相位解纏，獲得了平滑的解纏相位和數(shù)字高程模型(DEM)重建結(jié)果。該信號(hào)處理方法具提高DEM產(chǎn)品精度的潛力。

分布式衛(wèi)星； 合成孔徑雷達(dá)； 多基線干涉； 最大似然法； 長(zhǎng)短基線干涉； 相位解纏； 相位估計(jì)； 基線比例； DEM重建

0 引言

基線是干涉SAR技術(shù)中的重要參量，基線長(zhǎng)度直接影響圖像相干性、測(cè)高靈敏度(測(cè)高誤差的傳遞系數(shù))，進(jìn)而影響系統(tǒng)性能?；€越長(zhǎng)，測(cè)高誤差的傳遞系數(shù)越小，但干涉圖像對(duì)的相干性下降，干涉相位誤差也越大，影響高程的測(cè)量?；€越短，干涉圖像對(duì)的相干性增大，干涉相位誤差也越小，但測(cè)高誤差傳遞系數(shù)將越大，也影響高程測(cè)量精度。另外，基線的選擇還與地形有關(guān)，對(duì)陡峭區(qū)域，基線不能過(guò)長(zhǎng)。由于上述這些約束條件，使單基線干涉在使用中存在很大局限性。多基線SAR干涉測(cè)量有多個(gè)基線對(duì)應(yīng)的干涉相位圖，可優(yōu)勢(shì)互補(bǔ)，在保持相位解纏的穩(wěn)健性、減小干涉相位誤差，提高高程估計(jì)精度等方面有其優(yōu)勢(shì)，能克服傳統(tǒng)單基線SAR干涉測(cè)量方法的局限性[1-6]。在系統(tǒng)實(shí)現(xiàn)中，傳統(tǒng)多基線干涉主要通過(guò)單星多航過(guò)實(shí)現(xiàn)，滿足合適的多基線干涉圖像對(duì)概率較小。隨著分布式衛(wèi)星InSAR系統(tǒng)成為研究熱點(diǎn)，以及第一個(gè)實(shí)際系統(tǒng)TanDEM-X在軌成功運(yùn)行，對(duì)多基線干涉系統(tǒng)優(yōu)化設(shè)計(jì)和信號(hào)處理的研究有其實(shí)用意義[7-9]。

與單基線干涉相同，多基線干涉中各基線的絕對(duì)長(zhǎng)度對(duì)高程估計(jì)的精度有很大影響，而且各基線的相對(duì)比例對(duì)高程估計(jì)精度亦有較大影響。對(duì)多基線干涉的研究多集中于信號(hào)處理，對(duì)多基線的優(yōu)化設(shè)計(jì)研究較少。文獻(xiàn)[10]將多基線等效為雷達(dá)陣列間距，從陣列方向圖優(yōu)化角度設(shè)計(jì)基線長(zhǎng)度，實(shí)質(zhì)為三維成像應(yīng)用。文獻(xiàn)[11]從最大后驗(yàn)概率方法出發(fā)進(jìn)行多基線設(shè)計(jì)，并基于長(zhǎng)方體等理想幾何體場(chǎng)景仿真進(jìn)行驗(yàn)證，該法利用了像素間的先驗(yàn)信息。本文基于一種常用的多基線相位估計(jì)方法——最大似然法，無(wú)需像素間的先驗(yàn)信息，通過(guò)比較兩種優(yōu)化準(zhǔn)則：干涉相位估計(jì)方差的Cramer-Rao界最小準(zhǔn)則和多基線相位解纏穩(wěn)健性最大化準(zhǔn)則，優(yōu)選后者進(jìn)行多基線設(shè)計(jì)，并通過(guò)長(zhǎng)短基線仿真實(shí)例進(jìn)行驗(yàn)證。

1 最大似然法相位估計(jì)原理

最大似然法是一種常見(jiàn)的相位估計(jì)方法，它利用最大似然準(zhǔn)則，從多幅SAR圖像中估計(jì)解纏后的干涉相位值，可作為基線比例設(shè)計(jì)的理論依據(jù)[1]。

討論最大似然法相位估計(jì)原理。考慮接收星K顆，則有K幅進(jìn)行干涉的SAR圖像，以第1幅圖像為基準(zhǔn)，其矢量形式可表示為

(1)

Φ=

(2)

此處：b1k為第k顆衛(wèi)星與第1顆衛(wèi)星形成的干涉有效基線；φ1K為干涉相位；j為虛數(shù)單位。式(1)中各相位以第K顆衛(wèi)星與第1顆衛(wèi)星形成的干涉相位圖φ1K為基準(zhǔn)進(jìn)行表示。

用最大似然法進(jìn)行多基線相位估計(jì)，即可使

(3)

最小的φs(1X)值。此處：L為視數(shù)；Cin(m,n)為相干系數(shù)矩陣的逆在(m，n)處的值；上標(biāo)“*”表示共軛。估計(jì)的方差用Cramer-Rao界衡量，Cramer-Rao界可表示為

(4)

式中：γmn為相干系數(shù)矩陣在(m，n)處的值，即第m，n幅圖像對(duì)的相干系數(shù)；bmn為m，n星間基線相對(duì)基準(zhǔn)基線(1，K星基線)的比例，且

(5)

至此，得到了用最大似然法進(jìn)行多基線相位估計(jì)的最優(yōu)值及其對(duì)應(yīng)的基線比例。

以3幅圖為例，形成的兩根獨(dú)立基線的比例為0.7，歸一化E[f(φ)]隨φ的變化關(guān)系如圖1所示。此處：函數(shù)f定義見(jiàn)式(3)。由圖1可知：似然函數(shù)存在多個(gè)相近的極值，導(dǎo)致在噪聲較大時(shí)相位估計(jì)值存在模糊，直接影響相位解纏處理。

圖1 多基線系統(tǒng)的似然函數(shù)期望值(基線比0.7)Fig.1 Likelihood function of multi-baselinesystem with baseline 0.7

2 多基線系統(tǒng)的基線比例設(shè)計(jì)準(zhǔn)則

根據(jù)最大似然法原理，采用多基線相對(duì)單基線能提高干涉相位估計(jì)精度。多基線的處理性能與多基線配置有關(guān)，多基線配置的影響主要表現(xiàn)為相位估計(jì)和相位解纏兩方面。因此，多基線系統(tǒng)的基線比例設(shè)計(jì)準(zhǔn)則主要如下。

a)準(zhǔn)則1：按多基線融合后，干涉相位估計(jì)方差的Cramer-Rao界最小為準(zhǔn)則進(jìn)行基線比例尋優(yōu)。

仍以3幅圖為例，取基線比例k為0.5，1.0，相位標(biāo)準(zhǔn)差隨相干系數(shù)的變化規(guī)律如圖3所示。圖中：SIMO表示一發(fā)多收；SISO表示單發(fā)單收。由圖2可知：雙基線時(shí)用多干涉相位圖的融合可降低相位噪聲，進(jìn)而提高測(cè)高精度；基線比例對(duì)干涉相位估計(jì)精度有影響但并不顯著。因此，準(zhǔn)則1不適于作為基線比例設(shè)計(jì)的主要約束。

圖2 相位標(biāo)準(zhǔn)差隨相干系數(shù)的變化關(guān)系Fig.2 Root phase plot as function of correlation

b)準(zhǔn)則2：按多基線相位解纏的穩(wěn)健性最大化為準(zhǔn)則進(jìn)行基線比例尋優(yōu)，優(yōu)化的目標(biāo)函數(shù)是似然函數(shù)旁瓣峰值比。

由圖1可知：似然函數(shù)存在多個(gè)相近的極值，相位估計(jì)值存在模糊?？紤]似然函數(shù)與基線比例密切相關(guān)，基線比例直接影響似然函數(shù)的形狀，為此可用似然函數(shù)的旁瓣與主瓣比值表示解纏性能，旁瓣相對(duì)主瓣的比值越低，解纏性能就越穩(wěn)健。根據(jù)此準(zhǔn)則優(yōu)化基線比例。

因基線較短時(shí)，誤差傳遞系數(shù)偏大，為滿足高精度高程測(cè)量需求，必須選擇較長(zhǎng)基線。但長(zhǎng)基線又導(dǎo)致相位解纏難題，這可通過(guò)短基線輔助多基線相位解纏處理解決。2010年7月22日，德國(guó)TanDEM-X雙星系統(tǒng)利用長(zhǎng)短基線配合(有效基線長(zhǎng)度分別約2 600，200 m)獲得了DEM精度達(dá)10 cm左右，從原理和實(shí)踐證明了短基線輔助長(zhǎng)基線干涉提高測(cè)高精度的可行性[7]。

綜上所述，多基線長(zhǎng)度優(yōu)化設(shè)計(jì)的方法是：

a)先設(shè)計(jì)多基線系統(tǒng)中的長(zhǎng)基線?；趩位€對(duì)測(cè)高性能的傳遞模型，按測(cè)高精度要求反推長(zhǎng)基線的長(zhǎng)度范圍。

b)再用多基線相位解纏的穩(wěn)健性最大化為準(zhǔn)則進(jìn)行基線比例尋優(yōu)，獲得短基線長(zhǎng)度范圍，解決長(zhǎng)基線相位解纏難的問(wèn)題。

3 長(zhǎng)短基線分析與設(shè)計(jì)仿真

3.1 長(zhǎng)基線選擇

根據(jù)單基線與測(cè)高精度的約束關(guān)系，先進(jìn)行長(zhǎng)基線的設(shè)計(jì)，得到長(zhǎng)基線的長(zhǎng)度范圍[12]。用仿真分析與設(shè)計(jì)的方法，確定InSAR性能仿真參數(shù)為：衛(wèi)星高度500 km；信號(hào)帶寬400 MHz；相位同步剩余誤差8°；波束空間重疊率90%；主星定軌精度[2 2 2] m；等效斜距誤差4 m；干涉基線測(cè)量誤差[0.007 0.007 0.007] m；方位向/距離向配準(zhǔn)誤差0.1像素；沿航跡基線去相干0.95。

考慮長(zhǎng)短基線配合主要解決山區(qū)的高精度測(cè)高，本文針對(duì)山地地形(地面坡度25°)進(jìn)行仿真，給出不同的基線長(zhǎng)度及多視處理視數(shù)(山地區(qū)域分別為9，4視)下，對(duì)系統(tǒng)的測(cè)高性能進(jìn)行預(yù)估，仿真結(jié)果見(jiàn)表1。

表1 測(cè)高性能仿真結(jié)果

由表1可得初步結(jié)論：從相干性約束來(lái)說(shuō)，當(dāng)相干系數(shù)小于0.5時(shí)，干涉相位處理將很難進(jìn)行。從1∶25 000比例尺測(cè)繪標(biāo)準(zhǔn)的要求來(lái)說(shuō)，山地的高程精度優(yōu)于4 m。依據(jù)這兩個(gè)約束，長(zhǎng)基線干涉的基線長(zhǎng)度初步范圍取為1 000～1 800 m。

3.2 基線比例設(shè)計(jì)

在長(zhǎng)基線分別為1 000，1 800 m時(shí)，用似然函數(shù)旁瓣峰值比最小準(zhǔn)則求解短基線與長(zhǎng)基線的比例系數(shù)，得到似然函數(shù)旁瓣峰值比與基線比例的關(guān)系曲線如圖3所示。由圖3可知：對(duì)兩種長(zhǎng)基線，短基線長(zhǎng)度約為長(zhǎng)基線的0.1～0.2倍時(shí)，多基線解纏的性能最優(yōu)。按上述多基線解纏性能最優(yōu)時(shí)基線比例的仿真結(jié)果，長(zhǎng)基線1 800 m時(shí)對(duì)應(yīng)的短基線長(zhǎng)度應(yīng)選取180～360 m。

圖3 似然函數(shù)旁瓣峰值比與基線比例的關(guān)系Fig.3 PLSR of likelihood function asfunction of baseline ratio

3.3 仿真實(shí)驗(yàn)

基于3星編隊(duì)針對(duì)某山區(qū)場(chǎng)景的一組DEM數(shù)據(jù)進(jìn)行了信號(hào)仿真，基線1長(zhǎng)度1 800 m，基線2長(zhǎng)度340 m，基線2與基線1的比例約0.19(滿足本文的基線設(shè)計(jì)要求)，得到SAR圖像后進(jìn)行干涉處理，結(jié)果如圖4所示。

由仿真結(jié)果可知：基線1為長(zhǎng)基線，其干涉相位對(duì)地形起伏更敏感，干涉條紋(圖4(a))也更密集，存在條紋斷裂，相位殘差點(diǎn)多的問(wèn)題，此時(shí)采用單基線相位處理，導(dǎo)致解纏結(jié)果出現(xiàn)跳變的現(xiàn)象，DEM重建出現(xiàn)明顯偏差(如圖4(b)、(c)所示)。基線2為短基線，其干涉條紋圖則十分稀疏，能順利獲得相位解纏結(jié)果(如圖4(d)、(e)所示)，但其DEM重建結(jié)果較粗糙，表明測(cè)高精度低(如圖4(f)所示)。結(jié)合長(zhǎng)短基線干涉相位圖，采用最大似然法，通過(guò)短基線輔助長(zhǎng)基線相位解纏，獲得了平滑的解纏相位和DEM重建結(jié)果(如圖4(g)、(h)所示)?？梢?jiàn)，基于最大似然估計(jì)的多基線InSAR信號(hào)處理方法具提高DEM產(chǎn)品精度的潛力，仿真驗(yàn)證了用本文基線設(shè)計(jì)方法得到的短基線能輔助長(zhǎng)基線進(jìn)行相位解纏，進(jìn)而提高系統(tǒng)的測(cè)高精度。

4 結(jié)束語(yǔ)

本文對(duì)干涉SAR多基線進(jìn)行了分析并對(duì)設(shè)計(jì)進(jìn)行了研究。根據(jù)最大似然法估計(jì)相位的原理，多基線主要從提高相位解纏的穩(wěn)健性、減小干涉相位誤差兩個(gè)方面發(fā)揮作用，因此有兩種基線優(yōu)化的準(zhǔn)則：干涉相位估計(jì)方差的Cramer-Rao界最小準(zhǔn)則和多基線相位解纏穩(wěn)健性最大化準(zhǔn)則。通過(guò)比較，本文優(yōu)選后一準(zhǔn)則進(jìn)行多基線設(shè)計(jì)，得到多基線長(zhǎng)度優(yōu)化設(shè)計(jì)的方法為：首先按傳統(tǒng)方法得到單基線的長(zhǎng)度范圍，其次用多基線相位解纏穩(wěn)健性最大化準(zhǔn)則進(jìn)行基線比例尋優(yōu)，得到短基線長(zhǎng)度范圍，解決長(zhǎng)基線難以相位解纏的問(wèn)題。長(zhǎng)短基線仿真實(shí)例驗(yàn)證結(jié)果表明：所設(shè)計(jì)的短基線能很好地輔助長(zhǎng)基線進(jìn)行相位解纏，驗(yàn)證了本文方法的有效性。本文方法無(wú)需像素間的先驗(yàn)信息，適用性強(qiáng)，但后續(xù)還需解決基線比例多值性的問(wèn)題。

圖4 仿真數(shù)據(jù)多基線聯(lián)合處理結(jié)果Fig.4 Results of multi-baseline interfermetric data

[1] LOMBARDINI F. Optimum absolute phase retrieval in a three-element SAR interferometer[J]. Electron Let, 1998, 34: 1522-1524.

[2] ROBERTSON E. Multi-baseline interferometric SAR for iterative height estimation[D]. Provo: Brigham Young University, 1998.

[3] VINOGRADOV M V. Phase unwrapping method for the multifrequency and multibaseline interferometry[C]// Proceedings of IEEE Int Geoscience and Remote Sensing Conference (IGARSS 1998). Washington， DC: IEEE, 1998: 1103-1105.

[4] LI Zhen-fang, BAO Zheng, LI Hai, et al. Image auto-coregistration and InSAR interferogram estimation using joint subspace projection[J]. IEEE Trans on GRS, 2006, 44(2): 288-297.

[5] LACHAISE M, BAMLER R. Minimum cost flow phase unwrapping supported by multibaseline unwrapped gradient[C]// Proceedings of European Conference on Synthetic Aperture Radar (EUSAR 2010). Aachen: IEEE, 2010: 16-19.

[6] HUANG Hai-feng. A fast multi-baseline and multi-frequency band phase-unwrapping algorithm[J]. Measurement, 2014, 49: 401-406.

[8] MOREIRA A, KRIEGER G, HAJNSEK I, et al. TanDEM-L: a highly innovative bistatic SAR mission for global observation of dynamic processes on the earth’s surface[J]. IEEE Geoscience and Remote Sensing Magazine, 2015, 3(2): 8-23.

[9] BORLA TRIDON D B, BACHMANN M, MARTONE M, et al. The future of TanDEM-X: final DEM and beyond[C]// European Conference on Synthetic Aperture Radar (EUSAR 2016). Hamburg: IEEE, 2016: 124-126.

[10] 徐江, 唐勁松, 張春華, 等. 多基線InSAR信號(hào)處理中優(yōu)化陣列的應(yīng)用[J]. 信號(hào)處理, 2002, 18(5): 410-413.

[11] FERRAIUOLO G, MEGLIO F, PASCAZIO V. DEM reconstruction accuracy in multichannel SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(1): 191-201.

[12] 黃海風(fēng), 鄧泳, 梁甸農(nóng). 基于測(cè)高精度最優(yōu)值的星載分布式InSAR編隊(duì)構(gòu)形設(shè)計(jì)方法[J]. 國(guó)防科技大學(xué)學(xué)報(bào), 2005, 27(4): 85-90.

Multi-Baseline Analysis and Design of Interferometric Synthetic Aperture Radar

CHEN Guo-zhong1, ZHAO Di1, WANG Xiao-peng2,HOU Yu-sheng1, CHEN Liang3, ZHANG De-xin4

(1. Shanghai Institute of Satellite Engineering, Shanghai 201109, China; 2. Shanghai Academy of Spaceflight Technology, Shanghai 201109, China; 3. School of Electronic Engineering, University of Electronic Science and Technology, Chengdu 610054, Sichuan, China; 4. School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China)

To optimize the design of interfeometric synthetic aperture radar (InSAR) to improve altimetry data accuracy of the system, multi-baseline InSAR was analyzed and researched in this paper. The two criterions which were the minimum squared-error of interferometric phase and the robust of phase unwrapping were compared by maximum likelihood phase estimation method. The second one was selected as the criterion for designing baseline proportion. The optimal design method for multi-baseline length was given. First, the long baseline was determined. The scope of long baseline length was inversed according to the altimetry requirement based on the transfer model altimetry performance for single baseline. The baseline ratio was optimized by the robust of phase unwrapping criterion for multi baseline. The scope of the short baseline length was obtained to solve the phase unwrapping. The design example of multi baseline InSAR system for mountain terrain was presented. The results showed that the smooth phase unwrapping and digital elevation models (DEM) reconstruction were acquired through long baseline unwrapping assisted by short baseline using maximum likelihood method. The signal processing method proposed has the potential ability to improve the precision of DEM products.

Distributed satellite; Synthetic aperture radar; Multi-baseline interferometry; Maximum likelihood method; Long and short baseline interferomety; Phase unwrapping; Phase estimation; Baseline ratio; DEM reconstructed

1006-1630(2016)06-0026-05

2016-07-20；

2016-08-09

總裝備部“十二五”高分辨率對(duì)地觀測(cè)系統(tǒng)重大專項(xiàng)資助(GFZX0XXXXXXX);中國(guó)航天科技創(chuàng)新基金資助

陳國(guó)忠(1976—)，博士，高級(jí)工程師，主要研究方向?yàn)樾l(wèi)星總體設(shè)計(jì)、雷達(dá)信號(hào)和圖像處理。

TN958

A

10.19328/j.cnki.1006-1630.2016.06.003