吳一戎*

?

多維度合成孔徑雷達(dá)成像概念

吳一戎

(中國科學(xué)院電子學(xué)研究所 北京 100190)

該文通過回顧合成孔徑雷達(dá)(SAR)成像技術(shù)的發(fā)展歷程，論述了在極化、頻率、角度和時相空間內(nèi)開展多樣本聯(lián)合觀測的內(nèi)涵、必要性及發(fā)展趨勢。在此基礎(chǔ)上，基于信息綜合利用的需求建立了多維度SAR成像的基本定義，描述了多維度SAR成像的特點，并提出了基于克羅內(nèi)克積分解技術(shù)的觀測對象散射機理描述方法。最后列出了多維度SAR信號處理框架和基本假設(shè)條件，建立了SAR工作方式的多維度表述及多維度數(shù)的定義方法。

多維度合成孔徑雷達(dá)(MSJosSAR)；散射機理；克羅內(nèi)克積分解

1 引言

合成孔徑雷達(dá)(Synthetic Aperture Radar, SAR)成像技術(shù)是以微波譜段的電磁波作為探測手段，利用成像傳感器獲取觀測對象散射特征和相關(guān)信息的信息獲取技術(shù)。與傳統(tǒng)的光學(xué)成像相比，SAR成像不受日照和天氣條件的限制，能夠全天時、全天候的對目標(biāo)或場景進行觀測，已發(fā)展成為資源勘察、環(huán)境監(jiān)測和災(zāi)害評估的重要手段，是21世紀(jì)最具發(fā)展?jié)摿Φ膽?zhàn)略高技術(shù)領(lǐng)域之一，是國家綜合競爭力的具體體現(xiàn)。

一方面，SAR系統(tǒng)自身持續(xù)、迅猛的發(fā)展不斷推動著SAR成像理論的探索和關(guān)鍵技術(shù)的研究。最初的SAR數(shù)據(jù)獲取方式較為單一，主要工作在單極化、單角度和單波段。隨著SAR成像技術(shù)的發(fā)展及其應(yīng)用需求的推動，數(shù)據(jù)獲取方式日臻多樣化，逐步由單極化、單角度和單波段發(fā)展到多極化、多頻率、多角度和多時相等不同觀測方式的組合。從SAR成像技術(shù)自身的發(fā)展而言，在聯(lián)合觀測條件下，被觀測對象的散射機理及其在SAR圖像中的表征呈現(xiàn)出特性上的差異，上述特性直接關(guān)系到對觀測對象的認(rèn)知和理解。原有的針對單一觀測量或幾個觀測量簡單組合的SAR成像理論和方法已經(jīng)不適于描述觀測對象的復(fù)雜散射特征，也不利于提高定量遙感精度。

另一方面，國家重大應(yīng)用需求對未來SAR成像理論和方法提出了更高的要求。從應(yīng)用需求的角度來說，現(xiàn)有的SAR成像技術(shù)也很難解決實際應(yīng)用中的一系列難題，特別是高精度地形測繪、海洋復(fù)雜環(huán)境觀測以及森林資源與生態(tài)質(zhì)量動態(tài)監(jiān)測等國家重大需求。例如在高精度測繪方面，目前SAR成像技術(shù)只能實現(xiàn)對地表平均高度的測量和方位-斜距平面成像等。迭掩、陰影和體散射去相干等因素是制約高精度地形測量(包括形變檢測)和地表覆蓋物制圖發(fā)展的瓶頸所在。研究結(jié)果初步表明，能夠獲取極化、頻率、角度和時相多個空間內(nèi)聯(lián)合觀測量的SAR成像技術(shù)，具有解決上述難題的潛力。在海洋復(fù)雜環(huán)境觀測方面，目前單一觀測量的信息難以滿足對高精度水下地形探測和海冰厚度、類型及冰間水道探測的要求。相關(guān)研究初步表明，利用多時相的SAR成像數(shù)據(jù)可提高水下地形的探測精度，利用多極化SAR成像數(shù)據(jù)可提高海冰的分類精度。此外，國家近年來逐步啟動了天然林資源保護工程、退耕還林工程、京津風(fēng)沙源治理工程、“三北”及長江流域等防護林體系建設(shè)工程等六大林業(yè)重點生態(tài)工程。對上述林業(yè)重點生態(tài)工程的進展和實施效果進行監(jiān)測和評價，迫切需要提高森林3維空間分布信息的遙感提取精度。多波段、多時相、極化和極化干涉SAR 對于森林垂直結(jié)構(gòu)信息的定量反演具有獨特的優(yōu)勢。

2 SAR成像技術(shù)發(fā)展歷程

伴隨著各種先進技術(shù)的出現(xiàn)、定量化應(yīng)用需求的增加以及科學(xué)研究的不斷深入，國際上SAR成像技術(shù)在極化、頻率、角度和時相等空間觀測能力的發(fā)展大致歷經(jīng)了以下幾個重要階段：

第一階段主要是從距離向壓縮到方位向合成孔徑，再發(fā)展到干涉高程測量。1951年美國的Carl Wiley首先提出頻率分析方法改善雷達(dá)角分辨率，從而使得對觀測區(qū)域的方位向合成孔徑成像成為可能。1969年，Rogers等人提出將雷達(dá)干涉技術(shù)應(yīng)用于金星表面觀測。1974年，美國航空航天局(NASA)利用機載SAR系統(tǒng)首次對地球地形進行了雷達(dá)干涉測量，從而使得微波成像技術(shù)開始向干涉高程測量發(fā)展。從1996年到2001年，美國啟動“航天飛機雷達(dá)地形測量任務(wù)”(SRTM)計劃，利用“奮進號”航天飛機完成了約占全球80%陸地面積的測繪，經(jīng)處理可制成地球表面的數(shù)字高程模型。

第二階段重點是在單一觀測空間內(nèi)的延伸。主要包括極化、角度、頻率等。波段由微波波段向低波段和太赫茲(THz)擴展；極化由單極化向多極化到全極化的發(fā)展，如日本的星載全極化SAR系統(tǒng)ALOS/PALSAR，德國全極化SAR系統(tǒng)TerraSAR- X和加拿大的RADARSAT-2，極化SAR技術(shù)已成為全球環(huán)境監(jiān)測和地表測繪的重要手段之一；觀測角度由單站SAR成像向雙/多基方向發(fā)展。

第三階段是對兩個觀測空間的綜合研究，是目前國外的研究熱點。在綜合極化和角度兩個觀測空間方面，極化干涉SAR已開展了大量的研究工作。最近，又提出了將頻率和角度結(jié)合的MIMO-SAR。1997年，Papathanassiou等人首先利用SIR-C/X- SAR數(shù)據(jù)研究了頻率、極化對干涉相干性的影響，這是極化信息和干涉信息相結(jié)合的極化干涉SAR的最初起源。極化干涉SAR將極化和角度兩個空間的信息綜合，充分利用了干涉SAR對地表散射體高程敏感的特性和極化SAR對體散射體的形狀、方向和介電特性敏感的特性，具有提高傳統(tǒng)干涉測量相干性的能力以及提取體散射體結(jié)構(gòu)信息的潛力，使得SAR應(yīng)用于體散射體結(jié)構(gòu)提取等成為可能。目前，國際上針對于極化干涉研制了不同種類的機載實驗平臺，主要有德國DLR E-SAR(P/L/S三波段全極化)、日本PiSAR(L/X雙波段全極化)、美國AIRSAR(P/L/C三波段全極化)以及法國RAMSES系統(tǒng)(P/L/S/C/X/Ku/Ka/W多波段全極化)等，都具有重復(fù)飛行極化干涉測量能力，可實現(xiàn)多頻、高空間分辨和干涉測量等多種組合模式的數(shù)據(jù)獲取。近年來，歐、美和加拿大等已經(jīng)開始星載極化干涉SAR系統(tǒng)的初步研究，已開展收發(fā)分置星載極化干涉SAR系統(tǒng)的可行性論證。星載極化和極化干涉SAR系統(tǒng)已成為目前國外研究和研制的重點，但這是建立在他們已經(jīng)研制了若干部機載極化SAR和極化干涉SAR系統(tǒng)，并開展大量實驗研究的基礎(chǔ)之上的。從2003年起，歐空局隔年組織一次全球范圍的POLinSAR研討會，到目前為止一共舉行了五屆，進一步使得該技術(shù)的應(yīng)用研究獲得了廣泛關(guān)注，利用E-SAR, AIRSAR, RAMSES等系統(tǒng)獲取的機載極化干涉SAR數(shù)據(jù)，以及ENVISAT ASAR, ALOS/PALSAR, TerraSAR-X和RADARSAT-2等星載重復(fù)飛行極化干涉數(shù)據(jù)開展大量的理論和方法的研究。

從20世紀(jì)70年代開始，中國科學(xué)院電子學(xué)研究所開展SAR成像技術(shù)的研究，并于1979年成功研制第一部機載SAR原理樣機。中航607所、中電14所和中電38所等單位也相繼開展了SAR系統(tǒng)的研制工作。經(jīng)過近30年的發(fā)展，尤其是近10年來，我國的SAR成像技術(shù)得到了迅猛的發(fā)展，2維成像分辨率已由最初的幾十米提高到目前的亞米級，從而為目標(biāo)的檢測與識別等應(yīng)用奠定了堅實的基礎(chǔ)?！笆濉逼陂g，在國家863計劃支持下，成功研制了我國第一部機載雙天線干涉SAR系統(tǒng)。

可以預(yù)見，在一定約束條件下，獲取極化、頻率、角度和時相空間中多個空間的聯(lián)合觀測量，是未來SAR成像技術(shù)及其應(yīng)用發(fā)展的重要趨勢。

3 多維度SAR的定義和基本內(nèi)涵

多維度SAR (Multidimensional Space Joint- observation SAR, MSJosSAR)是在一定約束條件下，以SAR的基本觀測方式在極化、頻率、角度和時相空間中的至少兩個空間內(nèi)，分別獲得多個觀測量集合的聯(lián)合觀測技術(shù)手段。與之相對的是，在單一空間內(nèi)獲取單/多個觀測量集合的探測手段是單維度SAR。多維度SAR通過信號與信息綜合處理，有可能更準(zhǔn)確地區(qū)分被觀測對象的不同散射機理，進而更準(zhǔn)確地獲得其幾何特征和物理特征。針對特定應(yīng)用，分別從數(shù)據(jù)域、信號域和特征域?qū)τ^測對象進行客觀描述。以觀測對象與電磁波之間的相互作用規(guī)律為基礎(chǔ)，發(fā)展成像機制優(yōu)化方法、成像處理技術(shù)和信息提取方法，并最終實現(xiàn)對觀測對象各類特征的精細(xì)刻畫和定量反演。

在單色平面波假設(shè)下可得任意位置處觀測對象的散射場，表示為

式中為電磁波頻率；表示傳播時刻；描述了平面波傳遞方向；與位置的點乘運算代表了波矢量在該位置處的投影量；可以通過一對單位瓊斯矢量來表達(dá)電磁場在傳播過程 中的極化狀態(tài)；為處觀測對象對入射波的2×2復(fù)散射系數(shù)矩陣，也是SAR觀測的目標(biāo)反演量。該項不僅與觀測對象自身的特征有關(guān)，還與入射波的極化、頻率、角度和時相都有關(guān)。由式(1)可知，通過改變雷達(dá)工作的極化、頻率、角度和時相參數(shù)，就可以獲得對同一個觀測對象的多個像。

為了直觀地表示多維度SAR的觀測空間，需要建立多維度參數(shù)坐標(biāo)系。首先，極化空間可以利用龐加萊(Poincaré)球來完全表征，如圖1所示。

龐加萊球可以表征由極化基變換式定義的任意極化矢量：

(2b)

(2c)

式中為特殊酉變換；式(2a)，式(2b)通過極化方向角、極化橢圓率角及初始相位角給出了任意極化基的表達(dá)式；式(2c)定義了任意兩組正交極化基下場強的變換關(guān)系。由此可知，全極化合成孔徑雷達(dá)具備了在有限頻點和觀測角度內(nèi)獲取觀測對象完整的瞬態(tài)極化散射特征的能力，即在極化空間內(nèi)的完備性。正是由于這一完備性，極化SAR具有了對觀測對象極化散射特征的描述、分類及定量化比較能力。

圖1 龐加萊球

與極化空間不同，在現(xiàn)有技術(shù)水平及傳感器使用條件下，式(1)中的頻率、角度和時相都無法提供和極化空間相比擬的完備集。具體來說，頻率的完備集可以描述為一條由從低到高不同頻點組成的直線，反映了觀測對象散射特征的頻譜響應(yīng)曲線。而角度的完備集則可以通過觀測波矢在目標(biāo)坐標(biāo)系的俯仰角及方位角來表征，即空間視角，反映了觀測對象散射特征的視角變化響應(yīng)曲面。復(fù)波數(shù)

定義了電磁波傳播的頻率和方向。與龐加萊球相類似，可以將平面波復(fù)波數(shù)集合直觀地表示為描述3維波數(shù)域的球，如圖2所示。

圖2 球

圖3多維度球簇

通過在極化、頻率、角度和時相4個空間對發(fā)射/接收電磁波的描述，并結(jié)合現(xiàn)有合成孔徑雷達(dá)的觀測能力、使用模式及獲取信息的完備性分析，可以將多維度合成孔徑雷達(dá)定性地描述為以極化空間為核心，能夠獲取可以反映觀測對象散射特征隨頻率、角度和時相變化的離散多維度觀測量子集，結(jié)合多維度信號處理方法來擬合目標(biāo)在極化、頻率、角度和時相4個空間坐標(biāo)上的完整散射特征像，用于直接定量化反演觀測對象的生物物理參數(shù)。

4 多維度SAR信號模型

有關(guān)文獻指出，極化層析SAR信號可以在滿足一定條件的前提下表達(dá)為觀測對象散射特征之和的形式。本文將其拓展到多維度SAR觀測空間，即若觀測對象的散射特征響應(yīng)隨極化、頻率、角度和時相的變化滿足：

H1：每種散射機制相互統(tǒng)計獨立；

H2：每種散射機制的極化特征穩(wěn)定；

H3：每種散射機制在有限帶寬內(nèi)的頻率響應(yīng)特征穩(wěn)定；

H4：每種散射機制在有限角度內(nèi)的角度響應(yīng)特征穩(wěn)定；

H5：每種散射機制在有限時間內(nèi)的時相響應(yīng)特征穩(wěn)定。

則有，觀測量的復(fù)相關(guān)的期望是多種獨立散射機制下的克羅內(nèi)克積之和。如下式所示：

其中

(5)

(7)

(8)

(1) 極化復(fù)散射系數(shù)模型 例如F-D極化目標(biāo)分解模型可以將式(9)退化為單頻，單視角，單一時刻觀測下的極化SAR情況，

(2) 極化層析(Tomo)成像模型 可以將式(9)退化為單頻，多視角，單一時刻觀測下的極化SAR情況，

(3) 極化相干散射模型 特例為RVoG極化干涉模型

多維度SAR中的維度可以定義為由極化、頻率、角度和時相空間所張成雷達(dá)觀測空間，用表示。由不同雷達(dá)工作方式所獲取的觀測量構(gòu)成雷達(dá)觀測空間中的一個或多個樣本。例如，以單極化、單頻、單視角在單一時刻獲取的SAR觀測量對應(yīng)于空間中的一個樣本，獲取樣本數(shù)是1；多維度SAR通過個極化、個頻、個視角、個時相的觀測量得到觀測空間的個樣本。需要注意的是，判斷是否屬于多維度SAR觀測的標(biāo)準(zhǔn)并不僅依據(jù)獲取樣本數(shù)量，還需要考慮是否在多個空間內(nèi)分別獲取了多個樣本。例如表1中的差分干涉SAR雖然可以在時相空間內(nèi)獲取個樣本，但由于在其它3個空間內(nèi)的樣本數(shù)都是1，因此也在單維度SAR范疇之內(nèi)。與之相對，極化干涉SAR在極化空間內(nèi)的樣本數(shù)是4，同時在角度空間內(nèi)的樣本數(shù)是2，因此屬于多維度SAR。表1列出了部分常見的SAR工作方式，表中有底色的工作方式都屬于多維度SAR觀測。由此，將多維度數(shù)定義為包含多樣本的空間數(shù)量。

5 總結(jié)與展望

本文總結(jié)了合成孔徑雷達(dá)(SAR)成像技術(shù)的發(fā)展歷程，論述了多維度SAR成像的必要性及發(fā)展趨勢。在此基礎(chǔ)上，基于信息綜合利用的需求解釋了多維度SAR成像的基本內(nèi)涵；建立了多維度SAR成像的定義；提出了基于觀測對象特征分解的描述方法、多維度信號處理框架和基本假設(shè)條件；建立了SAR工作方式的多維度表述及多維度數(shù)的定義方法。

表1 SAR工作方式的多維度表述

在多維度SAR應(yīng)用的第一個層面上，通過融合單維度復(fù)散射系數(shù)圖便具有了一定的觀測對象散射特征反演能力，不存在維度間的相干獲取限制；在多維度SAR應(yīng)用的第二個層面上，通過層析或3維成像等信號處理方法在維度間直接形成多維度復(fù)散射系數(shù)圖像，此類多維度SAR應(yīng)用對相干獲取有一定要求，反演精度中等；在多維度SAR應(yīng)用的第三個層面上，通過提取和分離相干散射形成的相位信息來直接反演觀測對象的生物物理參數(shù)，這是在多維度SAR應(yīng)用中對相干獲取的要求最高但反演也最精確的類型。在未來的多維度散射機理研究中，應(yīng)根據(jù)觀測對象的散射類型特征以及待反演量來建立場的具體表達(dá)形式，利用極化空間維度的完備性來作為分析驗證工具，利用張量等數(shù)學(xué)工具來解決多維度遙感定量化所面臨的獲取機制優(yōu)化、信號處理及參數(shù)反演等問題。

致謝 本文的工作是在國家自然科學(xué)基金重大項目“多維度微波成像基礎(chǔ)理論與關(guān)鍵技術(shù)”的資助下完成的。重大項目課題一“多維度微波成像散射機理與成像機制”的核心成員洪文研究員、王彥平研究員、王衛(wèi)延研究員、李洋、尹嬙、譚維賢、趙曜博士等對本文有重要貢獻。在項目研究和概念完善過程中得到了基金委信息學(xué)部及項目專家組的指導(dǎo)以及項目承研單位中國科學(xué)院電子學(xué)研究所、西安電子科技大學(xué)、北京理工大學(xué)、中國林業(yè)科學(xué)研究院、武漢大學(xué)和國家海洋局海洋一所等參研團隊的寶貴建議，在此一并表示感謝。

[1] Sherwin C W, Ruina J P, and Rawclisse R D. Some early developments in synthetic aperture radar systems[J]., 1962, 6(2): 111-115.

[2] Brown W M. Synthetic aperture radar[J]., 1967, 3(2): 217-229.

[3] Wiley C A. Synthetic aperture radars[J]., 1985, 21(3): 440-443.

[4] Currie A and Brown M A. Wide-swath SAR[J]., 1992, 139(2): 122-135.

[5] Moreira J, Schwabisch M, and Fornaro G. X-SAR interferometry: first results[J]., 1995, 33(4): 950-956.

[6] Gabriel A K and Goldstein R M. Crossed orbit interferometry: theory and experimental results from SIR-B[J]., 1988, 9(5): 857-872.

點評：不僅要綠化荒山，還要在村莊道路兩側(cè)及村莊周邊栽植綠化樹和景觀樹，在公路沿線、村集體荒灘、荒山、閑置地種植果樹、綠化樹和風(fēng)景樹，不僅動員村民們干，而且全體黨員以“示范林”的形式干，鄉(xiāng)村自然生態(tài)的改變指日可待！

[7] Goldstein R M, Zebker H A, and Werner C L. Satellite radar interferometry: two-dimension phase unwrapping[J]., 1988, 23(5): 713-720.

[8] Hoen E W and Zebker H A. Penetration depths in inferred from interferometric volume decorrelation observed over the Greenland Ice Sheet[J]., 2000, 38(6): 2571-2583.

[9] Bamler R and Eineder M. ScanSAR processing using standard high precision SAR algorithms[J]., 1996, 34(1): 212-218.

[10] Dubois-Fernandez P, Cantalloube H, and Dupuis X. Clutter statistical analysis for high resolution SAR data[C]. Proceedings of IGARSS, Anchorage, 2004, 3: 1781-1783.

[11] Cantalloube H and Dubois-Fernandez P. Airborne X-band SAR imaging with 10 cm resolution: technical challenge and preliminary results[J].-,, 2006, 153(2): 163-176.

[12] Prati C and Rocca F. Improving slant-range resolution with multiple SAR surveys[J]., 1993, 29(1): 135-143.

[13] Curlander C and McDonough R N. Synthetic Aperture Radar: Systems and Signal Processing[M]. New York: Wiley, 1991.

[14] Langenberg K J, Brandfa? M, Fellinger P,.. Radar Target Imaging[M]. Berlin, Germany: Springer-Verlag, Eds. H. Uberall and W. M. Boerner, 1994: 113-151.

[15] Huynen J R. Phenomenological theory of radar targets[D]. [Ph.D. dissertation], Delft: Delft University of Technology, The Nether-lands, 1970.

[16] Boener W M. Need for developing multi-band single and multiple pass polInSAR platforms in air and space[C]. IEEE International Geoscience and Remote Sensing Symposium, Toulouse, France, July 21-25, 2003: 422-424.

[17] Cloude S R. Concept of polarisation entropy in optical scattering[J]., 1995, 34(6): 1599-1610.

[18] Cloude S R and Pottier E. A review of target decomposition theorems in radar polarimetry[J]., 1996, 34(2): 498-518.

[19] Cloude S R and Pottier E. An entropy based classification scheme for land applications of polarimetric SAR[J]., 1997, 35(1): 68-78.

[20] Cloude S R and Papathanassiou K P. Polarimetric SAR interferometry[J]., 1998, 36(5): 1551-1565.

[21] Papathanassiou K P and Cloude S R. Single baseline polarimetric SAR interferometry[J]., 2001, 39(11): 2352-2363.

[22] Knaell K K and P Cardillo G. Radar tomography for the generation of three-dimensional images[J].-,, 1995, 142(2): 54-60.

[23] Cafforio C, Prati C, and Rocca F. SAR data focusing using seismic migration techniques[J]., 1991, 27(2): 194-207.

[24] Jakowatz C V Jr, Wahl D E, Eichel P H,.. Spotlight- Mode Synthetic Aperture Radar: A Signal Processing Approach[M]. Dordrecht, The Netherlands: Kluwer Academic Publishers, 1996.

[25] Ishimaru A, Chan T-K, and Kuga Y. An imaging technique using confocal circular synthetic aperture radar[J]., 1998, 36(5): 1524-1530.

[26] She Zhishun, Gray D A, Bogner R E,.. Three- demensional SAR Imaging via multiple pass processing[C]. Proceedings of IGARSS, Hamburg, Germany, 1999.

[27] Reigber A and Moreira A. First demonstration of airborne SAR tomography using multibaseline L-band data[J]., 2000, 38(5): 2142-2152.

[28] Fornaro G, Lombardini F, and Serafino F. Three-dimensional multipass SAR focusing: experiments with long-term spaceborne data[J]., 2005, 43(4): 702-714.

[29] Massonnet D,.. The displacement field of the Landers earthquake mapped by radar interferometry[J]., 1993, 364: 138-142.

[30] Ferreti A, Prati C, and Rocca F. Permanent scatterers in SAR interferometry[J]., 2001, 39(1): 8-20.

[31] Noferini L, Pieraccini M, Mecatti D,.. Permanent scatterers analysis for atmospheric correction in ground-based SAR inter- ferometry[J]., 2005, 43(7): 1459-1471.

[32] Colesanti C, Ferreti A, Novali F,.. SAR monitoring of progressive and seasonal ground deformation using the permanent scatterers technique[J]., 2003, 41(7): 1685-1701.

[33] Lombardini F. Differential tomography: a new framework for SAR interferometry[J]., 2005, 43(1): 37-44.

[34] Singhroy V and Saint-Jean R. Geological uses of multi- incidence RADARSAT images[C]. Proceedings of IGARSS, Seattle, WA, USA, 1998.

[35] Raney R K, Runge H, Bamler R,.. Precision SAR processing using chirp scaling[J]., 1994, 32(4): 786-799.

[36] Moreira A, Mittermayer J, and Scheiber R. Extended chirp scaling algorithm for air- and spaceborne SAR data processing in stripmap and ScanSAR imaging modes[J]., 1996, 34(5): 1123-1136.

[37] Mittermayer J, Moreira A, and Loffeld O. Spotlight SAR data processing using the frequency scaling algorithm[J]., 1999, 37(5): 2198-2214.

[38] Bamler R and Eineder M. Accuracy of differential shift estimation by correlation and split-bandwidth inter- ferometry for wideband and delta-k SAR systems[J]., 2005, 2(2): 151-155.

[39] Cumming I G and Wong F H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation[M]. Artech House, Inc., 2005.

[40] Reigber A, Alivizatos E, Potsis A,.. Extended wavenumber-domain synthetic aperture radar focusing with integrated motion compensation[J].-,, 2006, 153(3): 301-310.

[41] Wu Yirong and Hong Wen. From Synthetic Aperture Radar Imaging to Multi-Parameter Multi-Dimensional Microwave Imaging[C]. Proceedings of APSAR, Huangshan, China, November 2007.

[42] Krieger G, Moreira A, Fiedler H,.. TanDEM-X: a satellite formation for high-resolution SAR interferometry[J]., 2007, 45(11): 3317-3340.

[43] Tebaldini S and Rocca F. Polarimetric options for SAR tomography of forested areas[C].Proceedings of the Fourth International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interferometry, Netherlands, 2009: 45-48.

[44] Lin Yun, Hong Wen, Tan Weixian,.. Interferometric circular SAR method for three-dimensional imaging[J]., 2011, 8(6): 1026-1030.

[45] Zhou Yongsheng, Hong Wen, Wang Yanping,.. Maximal effective baseline for polarimetric interferometric SAR forest height estimation[J].:, 2012, 55(4): 867-876.

[46] Peng X, Tan W, Wang Y,.. Convolution back-projection imaging algorithm for downward-looking sparse linear array three dimensional synthetic aperture radar[J]., 2012, 129: 287-313.

[47] Tebaldini S and Rocca F. Multibaselinepolarimetric SAR tomography of a borealforestat P- and L-bands[J]., 2012, 50(1): 232-246.

[48] Ardeshir G, Raj Mittra, and Philip J M. Electromagnetic Wave Interactions[M]. World Scientific Pub Co Inc., 1997: 169-172.

[49] Boerner W M, Wei-ling Y, Xi A-Q,.. On the basic principles of radar polarimetry: the target characteristic polarization state theory of Kennaugh, Huynen’s polarization fork concept, and its extension to the partially polarized case[J]., 1991, 79(10): 1538-1550.

[50] Michiel Hazewinkel. Tensor Product[M]. Encyclopedia of Mathematics, German: Springer, 2001.

[51] Freeman A and Durden S L. A three-component scattering model for polarimetric SAR data[J]., 1998, 36(3): 963-973.

[52] Rudin Walter. Real and Complex Analysis[M]. New York: McGraw-Hill Science Education, 1986.

Concept of Multidimensional Space Joint-observation SAR

Wu Yi-rong

(Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China)

The objective of this study is to investigate the connotation, necessity, and tendency of Synthetic Aperture Radar (SAR) imaging technology within the framework of multidimensional space joint-observation, which includes polarimetry, frequency, angle, and time spaces, on the basis of the key evolvement phases of the SAR imaging technology. Furthermore, the definition of Multidimensional Space Joint-observation SAR (MSJosSAR) is given on the basis of the demand for the information fusion of MSJosSAR images. Then, the advantage of MSJosSAR is revealed by using the Kronecker product decomposition technology for gaining a better understanding of the target scattering mechanisms. Moreover, the hypothesis and the basic framework of the MSJosSAR signal processing are listed and illustrated. Finally, the number of joint observation spaces in typical SAR configurations is enumerated.

Multidimensional Space Joint-observation SAR (MSJosSAR); Scattering mechanisms; Kronecker product decomposition

TN95

A

2095-283X(2013)02-0135-08

10.3724/SP.J.1300.2013.13047

吳一戎(1963-)，男，中國科學(xué)院院士，研究員，博士生導(dǎo)師，中國科學(xué)院電子學(xué)研究所所長，主要研究方向為高分辨率和干涉SAR系統(tǒng)技術(shù)、遙感衛(wèi)星地面處理與應(yīng)用系統(tǒng)技術(shù)等。

2013-05-02收到，2013-05-22改回；2013-05-24網(wǎng)絡(luò)優(yōu)先出版

國家自然科學(xué)基金(60890070, 60890071)資助課題

吳一戎 wyr@mail.ie.ac.cn