干涉合成孔徑雷達仿真策略研究
摘要： 將干涉合成孑L徑雷達(InSAR-Interferometric Synthetic perture Radar)仿真劃分為3個(gè)層次，以 仿真層次的劃分為出發(fā)點(diǎn)，從干涉三維高程重建模型的建模策略和誤差傳遞模型的建模策略?xún)蓚�(gè)方面詳細 對InSAR仿真策略進(jìn)行論述，提出分階段進(jìn)行誤差傳遞模型建模的建模策略．對于具體的仿真模型約束條 件和仿真所要實(shí)現的目的，在InSAR仿真中應該采用相應的仿真策略． 仿真技術(shù)在遙感測繪領(lǐng)域中的應用可以追溯到最初的衛星編程，隨著(zhù)成像傳感技術(shù)、圖像處理技術(shù) 和制圖技術(shù)等的發(fā)展而不斷發(fā)展．校正器| 轉換器| 傳送器| 變送器| 傳感器| 記錄儀| 有紙記錄儀| 無(wú)紙記錄儀|
仿真的目的是為了預測或可視化顯示地形在給定數據獲取設備和數據 獲取機理下所呈現的特征，進(jìn)而利用仿真技術(shù)進(jìn)行數據獲取設備的參數優(yōu)化設計、系統定標，并開(kāi)展遙 感應用技術(shù)或數據獲取方式方面的研究和評估． 仿真技術(shù)在合成孔徑雷達(SAR—Synthetic Aperture Radar)成像中也得到了廣泛應用[】q]，干涉 SAR(InSAR)是SAR與干涉技術(shù)相結合的產(chǎn)物，其突出特點(diǎn)在于以相位信息作為處理和分析的焦點(diǎn)， 能夠提供地形表面的高程信息，仿真技術(shù)同樣也是InSAR系統優(yōu)化設計的強有力工具．由于InSAR不 再局限于實(shí)現常規SAR獲取灰度圖的目的，而是提取地形表面的高程信息，因此InSAR仿真和SAR 仿真的側重點(diǎn)自然也就發(fā)生變化，不再側重于考察各種因素對灰度圖質(zhì)量的影響，其重點(diǎn)在于評估各種 因素對高程測量精度的影響，這種評估需要借助SAR仿真模型(包括地面目標模型、雷達與地面目標相 互作用建模等)、InSAR三維高程重建模型和誤差傳遞模型來(lái)進(jìn)行E‘ ]，而SAR仿真模型的建立可以借 助相應的SAR仿真技術(shù)進(jìn)行．因此，如何建立InSAR三維高程重建模型和相應的誤差傳遞模型便成為 干涉SAR仿真的關(guān)鍵． 本文以InSAR仿真層次劃分為出發(fā)點(diǎn)，從三維高程重建模型建模和誤差傳遞模型建模兩個(gè)方面詳 細論述InSAR仿真所應采取的仿真策略．

1 InSAR仿真層次劃分 根據InSAR 數據的獲取原理和InSAR處理流程，得到如圖1所示InSAR仿真框圖．圖中地面后 向散射特性模型的建立可以通過(guò)兩種途徑進(jìn)行：① 用概率模型近似表示；② 從真實(shí)的SAR 圖像中提取 地面的后向散射特性；SAR 數據獲取模型包括傳感器模型和運動(dòng)平臺模型以及成像幾何模型3部分；土 地模型包括數字高程模型和土地的使用狀況信息；虛線(xiàn)框內部分為InSAR仿真器． InSAR仿真的主要目的是考察各種因素對高程測量精度的影響，對于不同的因素，可以采用不同的 仿真流程，為便于問(wèn)題的分析，我們將InSAR仿真劃分為3個(gè)層次． 層次一的仿真為全過(guò)程仿真，仿真流程為(1)～(2)～(3)～(4)～(5)～(6)～ (7)，如圖1所示．該層 次仿真利用地形的后向散射特性模型、SAR 數據獲取模型、土地模型結合SAR數據獲取機理，仿真得到 SAR原始回波信號，后經(jīng)SAR成像處理和干涉處理得到干涉圖．該層次仿真的特點(diǎn)在于：采用的SAR處 理器為實(shí)際應用的SAR 處理器，沒(méi)有對SAR 處理器的有效性施加任何假設條件，便于對波長(cháng)級影響 InSAR 高程測量精度的因素進(jìn)行考察．該層次的仿真可以利用G．Alberti提出的Coherent仿真器[6]實(shí)現， 該仿真器從DEM 和地形的后向散射特性出發(fā)，考慮傳感器參數、成像幾何關(guān)系以及軌道參數和運動(dòng)平臺 的姿態(tài)，在時(shí)域中通過(guò)脈沖到脈沖綜合方法進(jìn)行干涉原始數據的仿真，仿真框圖如圖2所示． ⋯ 孽 數 SAR SAR SAR 復 據 干 原 原 復 圖 干 涉 獲 始 始 像 涉 圖 取 圖 處 生 機 數 數 配 理 制 據 據 像 準 成 (1) (2) (3) (4) (6) _ (5) (7) 圖1 干涉合成 L徑雷達仿真框圖 n舀1 Block diagram of InSAR simulation I型和翠平窖臺萎姿篡態(tài)H{ {仿真器卜l SAR 信 成 號 像 合 處 成 理 器 Il涉l 干 I處理I 一涉圖 I．．．．．．．_．J 圖2 層次一仿真策略下的InSAR 仿真框圖 n舀2 Block diagr am of InSAR simulation bac,~ on level 1 層次二的仿真避開(kāi)原始數據的仿真，直接利用地形的后向散射特性模型、SAR 數據獲取模型、土地 模型結合SAR 數據獲取機制仿真得到復圖像，后經(jīng)復圖像配準、干涉處理得到干涉圖，仿真流程為 (1)～ (4)～(5)～(6)～(7)，該層次的仿真需要SAR 系統沖激響應函數的支持，文獻[7]中的InSAR 仿 真器即為該層次的仿真器，該仿真器直接利用地形的后向散射特性、SAR 系統的沖激響應和實(shí)際的干涉 處理設備產(chǎn)生干涉圖．應用該層次仿真的前提條件是已知SAR系統的沖激響應或者沖激響應可以通過(guò) 獨立的仿真過(guò)程得到． 層次三的仿真直接根據地面后向散射特性模型、土地模型和干涉成像幾何模型以及成像機理直接得 到干涉高程圖，仿真流程為(1)～(7)．在該層次的仿真中，建立正確的InSAR成像幾何模型和地面后向 散射特性模型是該層次仿真的關(guān)鍵，該層次的仿真需要SAR系統沖激響應函數和干涉處理系統函數的 支持．L．S．Wray等提出的Incoherent仿真器L8 便是該層次仿真的典型代表，該仿真器在給定SAR圖像 的條件下，利用圖像的統計特性，結合成像幾何關(guān)系和地形的后向散射特性直接獲取干涉圖．
在InSAR仿真中，針對不同的仿真目的，需要采用不同層次的仿真．如要考察地形特性和SAR參 數對InSAR高程測量性能的影響，InSAR 的仿真應該采用層次一的仿真．對于復圖像配準和干涉處理 過(guò)程的仿真，層次二的仿真就足夠了．這種仿真從幾何的角度來(lái)說(shuō)，必須是嚴格的和準確的，但是可以 按照簡(jiǎn)單的方式進(jìn)行散射特性和去相關(guān)現象建模，甚至可以不考慮土地使用因素的影響．干涉圖可以通 過(guò)幾何關(guān)系直接得到，也可以用層次一或層次二仿真得到的復數據按照常規處理方法得到．由前者仿真 得到的相位差來(lái)源于斜矩的差值，可成為條紋估計，可以將失配或時(shí)間去相關(guān)引入的噪聲加入到干涉圖 中逼近真實(shí)效果，但是前提是有噪聲模型供利用．考察地形對相位差進(jìn)而對干涉高程精度影響的唯一方法是對影響相位差的其它因素進(jìn)行仿真，然后從相位差中剔除這些因素的影響，僅保留地形的影響，這 種仿真過(guò)程必須是基于層次一的仿真．但是，通常情況下，這些因素非常復雜，難以進(jìn)行預測(如大氣和 植被)，這些因素對相位的影響仍然存在于DEM 中，因此，InSAR仿真結果與真實(shí)結果間存在差異主要 歸于物理模型的準確性而不是仿真算法的有效性．

2 InSAR仿真策略 InSAR仿真算法的設計主要基于兩個(gè)約束條件[1 ]：① 應用方面的約束條件(高程測量精度、特定的 成像區域、高程測量的影響因素)；② 算法運行中的約束條件(軟、硬件條件和模型的可用性)．通盤(pán)考慮 這些約束條件，就能夠得到一個(gè)適當的、可行的仿真設計，本文稱(chēng)之為InSAR仿真策略． InSAR仿真建立在地形模型、干涉處理模型(包括復圖像配準、相位展開(kāi))和三維重建模型等一系列 仿真模型的基礎之上，應用誤差傳遞模型可以考察各模型輸入參數對高程精度的影響．其中，地形建模 可以參考SAR仿真技術(shù)的相關(guān)內容進(jìn)行，干涉處理模型依賴(lài)于一定的復圖像配準算法和相位展開(kāi)算法． 因此，我們從三維重建模型的建模策略和誤差傳遞模型的建模策略?xún)蓚�(gè)方面論述InSAR的仿真策略．

2．1 InSAR 目標的三維重建模型建模策略 利用SAR位置信息、SAR參數以及InSAR成像幾何模型，確定圖像上每一點(diǎn)的三維位置的過(guò)程稱(chēng) 為目標的三維重建．目標三維重建模型的建模策略有幾何建模和解析建模兩種． 幾何建模策略從雷達測距、多普勒頻率測量以及干涉相位測量的基本原理出發(fā)，利用距離球、多普 勒錐和相位雙曲面三表面的交點(diǎn)確定目標的三維位置，如圖3所示．幾何建模給出了InSAR 目標三維 重建的直觀(guān)解釋?zhuān)�但是幾何建模無(wú)法用于干涉高程精度的定量分析． 解析建模從InSAR 的基本方程入手，通過(guò)聯(lián) 立方程組的求解獲取目標的三維位置，實(shí)現三位重 建．基于這一建模策略的干涉仿真可以考察各干涉 參數對三維重建精度的影響．為了簡(jiǎn)化解析建模過(guò) 程，目前的解析建模主要建立在兩個(gè)假設基礎上： 電磁波波前的平面波假設和地球模型的平地假設． 最簡(jiǎn)單的三維重建算法是同時(shí)引入平地假設和平面 波假設[9]，本文作者對該算法進(jìn)行了仿真研究，其 定位誤差是不可接受的[1 ． Madsen[¨]和 Wilkinson L1。]的三維重建算法是在地心笛卡兒坐標 系中描述InSAR成像幾何關(guān)系，在保留平面波假 設的同時(shí)揚棄了平地假設，定位精度較文獻[9]有 了提高．另外，D．Small等人采用WGS84橢球模 型作為地球模型，將干涉得到的地形高度附加到橢 球半徑上形成數字高程模型，利用距離球、多普勒 目標位置 圖3 InSAR 目標三維重建的幾何解釋 Fig．3 Geometric Interpretation of Target 3-D Reconstruct|0n in InSAR 錐和生成的數字高程模型對目標進(jìn)行三維重建[1引．該算法通過(guò)迭代方式求解，計算非常復雜．考慮到現 有三維重建算法的缺陷，本文作者從InSAR的基本方程入手，充分利用InSAR的基本測量值，剔除現 有算法中的平面波假設和平地假設，提出了一種建立在球形地球和球面電磁波的普遍模型基礎上的 InSAR 目標三維重建算法．解析建模用解析函數的形式明確了目標高程位置和干涉參數間的關(guān)系，基于 這一建模策略的InSAR仿真可以定量考察各干涉參數對高程精度的影響．

2．2 誤差傳遞模型建模策略 從統計的觀(guān)點(diǎn)來(lái)看，干涉處理過(guò)程可以用一個(gè)依賴(lài)于一系列輸入參數的復雜的數學(xué)函數／( )進(jìn)行 表述，這一系列的參數在干涉處理的不同階段引入．輸入參數和表征干涉處理過(guò)程的廠(chǎng)( )函數均存在不確定性，這些不確定因素在計算過(guò)程中不可避免的引入誤差，并且這種計算誤差會(huì )傳遞到后續處理的 每一步．在已知輸入參數值及其精度的條件下，應用誤差 傳遞模型(EPM-Error Propagation Mode1)可以估計在特 定階段引入的誤差對函數f(x)的影響．由于InSAR 與 SAR的最大區別在于它能夠提供地形的高程信息，而這種 信息的獲得是通過(guò)干涉相位的測量實(shí)現的，因此，考察各 種因素對干涉相位的影響就成為InSAR仿真的核心內容， 而這種考察可以通過(guò)誤差傳遞模型有效實(shí)現 ，EPM 的典 型框圖如圖4所示． 誤差傳遞模型的建立通常有以下幾種方法[s 14]： FOSM (First一0rder Second—Moment)法、SOSM (Second— Order Second Moment)法、Rosenblueth法和Monte Carlo 圖4 誤差傳遞模型框圖 Fig．4 Block diagram of error propagation model 仿真法．FOSM 法和SOSM 法統稱(chēng)為臺勞級數法，在系統函數容易獲得而且較為簡(jiǎn)單的情況下，將系統 函數對各輸入參數求偏導數可以得到干涉靈敏度方程，根據干涉靈敏度方程和方差傳遞模型，采用 FOSM 法或SOSM 法就可以很容易的建立EPM；在已知系統輸入變量概率密度函數和已知系統函數的 情況下，Rosenblueth法是建立EPM 的有效方法，該方法將輸入變量的連續概率密度函數用若干離散點(diǎn) 的概率密度函數近似表示，利用系統函數的傳遞作用得到輸出變量的離散概率密度函數，進(jìn)而利用離散 概率密度函數的低階矩確定輸出及輸出變量的精度；而對于復雜過(guò)程來(lái)說(shuō)，很難獲得函數明確的解析表 達式，也無(wú)法得到輸入參量的概率密度函數，在這種條件下，FOSM 法、SOSM 法和Rosenblueth法均 不適用于EPM 的建立，通常通過(guò)Monte Carlo仿真的方法建立EPM，Monte Carlo仿真通過(guò)對具有一 定精度的確定性的輸入參數進(jìn)行隨機樣，將確定性問(wèn)題轉化為隨機問(wèn)題，通過(guò)對函數廠(chǎng)( )的重復計算， 近似得到輸出參數的精度，如果計算次數足夠大，就能夠得到輸出參數的完整分布和可靠的測量精度， 從而確立系統的誤差傳遞模型．
具體到干涉合成孔徑雷達EPM 的建立，在干涉的不同階段，根據問(wèn)題的繁易程度，采取分階段 EPM 建模的策略．下面從SAR原始信號仿真、干涉處理、InSAR成像幾何模型建模幾個(gè)階段分別闡述 EPM 建模策略． 在InSAR原始數據仿真中，地形特征和大氣效應均會(huì )對干涉的高程精度產(chǎn)生影響，此時(shí)的InSAR 仿真應該是基于層次一的仿真．在考察地形特征對干涉高程精度的影響中，地形特征通常用一定的概率 模型近似表示，其中聯(lián)合圓高斯概率分布模型[12,15 是一種通常采用的地面散射特性概率模型，系統函數 可以參照常規SAR的系統函數建立方法建立，此時(shí)EPM 宜采用Rosenblueth建模策略進(jìn)行建模，仿真 應該是基于層次一的仿真，G．Alberti提出的Coherent仿真器 即是這一建模策略的仿真器．如果干涉 仿真建模中要考慮大氣效應在內，由于大氣折射系數對斜距的作用是在波長(cháng)級而不是在象素級，此時(shí)的 仿真也應該是基于層次一的仿真，近似情況下，大氣折射系數模型可以采用分層建模策略 進(jìn)行，該模 型可以用解析表達式表達，因此大氣折射系數模型的EPM 可以采用FOSM 或者是SOSM 建模策略進(jìn) 行建模，因此，如考慮大氣效應在內的InSAR原始數據的仿真中，EPM 的建模策略應選用混合建模策 略：大氣效應模型的EPM 建模采用FOSM(或SOSM)建模策略，而地形模型的EPM 則采用 Rosenblueth建模策略． InSAR數據處理包括復圖像配準、相位展開(kāi)等，處理過(guò)程非常復雜，而且對于不同的處理算法，輸 入參數精度對輸出參數精度的影響各異，處理過(guò)程的系統函數難以用解析表達式明確表示，臺勞展開(kāi)法 和點(diǎn)估計法不適合該類(lèi)EPM 的建立，由于仿真建模只需知道輸入參數名義值和輸入參數精度，無(wú)需了 解干涉處理中與處理算法有關(guān)的細節內容，EPM 仿真建模是該階段EPM 的最佳建模策略．其中Monte Carlo仿真建�！�4]是最常用的EPM 仿真建模策略．仿真建模中，輸出參數的精度估計有時(shí)較為困難．
針 對參數精度估計問(wèn)題，已經(jīng)進(jìn)行了大量研究，如1994年Solaas對基線(xiàn)問(wèn)題的研究 ]． 在InSAR成像幾何關(guān)系建模中，成像的幾何模型通�？梢杂媒馕霰磉_式表示，如文獻[16]中二維InSAR和三維InSAR成像幾何模型的解析建模方法，以及文獻[17]中解決干涉直線(xiàn)軌道約束問(wèn)題而提 出的任意軌道建模策略． 由于成像幾何模型可以用明確的解析表達式表示，此階段的EPM 建模應用FOSM 法或者SOSM 法即可容易得到，這種建模策略的基本思想如下：利用成像幾何模型的基本關(guān)系式，確立目標高程對參 數變量的靈敏度方程，然后利用FOSM 法和SOSM 法建立輸出參數精度、輸入參數精度和靈敏度方程 三者間的聯(lián)系，通過(guò)方差傳遞模型，得到成像幾何模型的EPM．

3 結 論 InSAR是合成孔徑雷達技術(shù)與干涉技術(shù)相結合的產(chǎn)物，其獨特優(yōu)勢使之成為遙感測繪領(lǐng)域中一顆閃 亮的新星，但是由于InSAR對地形的測繪過(guò)多的依賴(lài)于SAR傳感器和地物反射特性，且成像過(guò)程具有 瞬時(shí)性，單純依靠硬件測試和室外飛行試驗難以完成InSAR系統的優(yōu)化設計，應用仿真技術(shù)可以取代 絕大部分室外試驗，節省系統研制的時(shí)間和經(jīng)費開(kāi)銷(xiāo)，同時(shí)可以避免錯誤決策和不當操作帶來(lái)的風(fēng)險． 本文給出的InSAR仿真策略將為順利開(kāi)展InSAR仿真提供一些借鑒意義．