INSAR提取地面高程信息的實驗研究

馮 強

北京同創達勘測有限公司，北京 100048

利用合成孔徑雷達干涉測量（INSAR）技術快速獲取高精度數字高程模型（DEM）是目前INSAR技術的廣泛應用之一。它在地信測繪、地表變形監測、冰川運動研究等方面都表現出了全天候、全天時、不受大氣和氣候影響、高精度高效率等突出優勢。

目前，世界上有許多科研及應用單位致力于SAR干涉技術及其應用實踐和軟件的開發，比較著名的事Earthview、Gamma等，本文立足國內研發的軟件，利用中國測科院研究開發的軟件進行INSAR技術提取DEM實驗，發現該方法易于掌握且效果較好。

1 INSAR的基本干涉原理

INSAR獲取數字地面DEM技術的基本激光干涉原理，主要是利用同區域，具有激光條紋干涉幾何構成模型成像能力，兩幅或者多幅SAR天線，獲取同一地區具有一定視差的兩幅圖像，該圖像是有干涉信息，同時也是單視復數主輔圖像，并根據其干涉相位信息進行分析計算，用最小范數反演出地面的高程數值，從而建立該地區的DEM[1]。

本實驗采用的是雙天線INSAR系統，即在SAR平臺上（飛機）安置兩（多）幅SAR天線，其中一幅天線向地面發射雷達激光波，另外的兩（多）幅天線，在同時相內，接收雷達波到達地面后的向散射反射輻射回波，從而得到相應同區域的兩（多）幅單視復數主輔輻射圖像。采用該方式可以同時獲取同一區域的主、輔含干涉波的數字輻射圖像，這樣的回波圖像具有的相干性很好的特性，可以利用干涉的波譜原理進行處理、分析和應用。

2 INSAR數據圖像處理一般流程

INSAR獲取數字DEM工作流程主要包括：地理配準復數主輔圖像、生成雷達干涉主輔影像圖匹配、評價雷達干涉圖像質量因子、去平地雷達效應、地雷干涉圖像濾波、相位干涉解纏、雙基線計算與干涉參數確定、數字DEM采集和重建等。

2.1 復數圖像配準

對于機載雙天線INSAR主、輔圖像的配準，由于干涉主輔圖像之間的之間存在相對偏（平）移，但偏移量較小，可以通過軌道參數進行計算減弱，可采用數值一級像元級匹配、子像元級二級灰度匹配和地理配準模型計算三步完成配準工作。

2.2 干涉圖生成

利用回波干涉影像計算出的基準地理影像配準后的模型，對主、輔雷達圖像進行重采集采樣，并將主、輔圖像的同名相應像元的復數值進行向量共軛相乘，加入非平行軌道情況的多高程點基線估計方法，計算干涉圖。

2.3 評價干涉影像質量因子

為了評價雷達影像干涉質量因子好壞并為干涉濾波和后續激光相位解纏提供參考依據和數值參數，從而計算出干涉質量圖文件。

以上三步可通過軟件輸入原始圖像、參數文件，同時進行計算，計算結果見圖1～3所示。

圖1 原始機載SAR圖像

圖2 干涉圖

圖3 干涉質量圖

2.4 去平地效應

由于采用雙天線的情況，圖4干涉影像圖的條紋，明顯能看出表現為基本豎直，且基本平行的條紋紋理，對這些干涉條紋進行基礎處理，經一步降低了雷達干涉圖像的條紋幾何頻率，減小干涉圖像濾波和相位解纏的難度系數，利用計算的參數文件對干涉圖進行去平地效應工作，并計算平地干涉相位。

圖4 去平地干涉圖

2.5 干涉圖濾波

為減少干涉圖像里的相位噪聲，干涉圖像采用濾波前和后干涉相位的均方差來評價干涉圖濾波相位保持精度的策略，從而降低（減弱）相位解纏的難度，可通過人工輸入濾波指數進行調整，實現精度較高、效果較好的濾波效果。

2.6 相位解纏

從干涉條紋圖像中，只能得到該相位的主值，它實際是以2π為模的不足一個周期內的相位差，但每一次的單個測量值的整數相位周期卻不能從信號中得到（整周期）。為了重建數字DEM，需要恢復丟失的整數相位周期（整周模糊恢復），這一過程就是相位解纏[4]。

相位解纏作為INSAR技術中關鍵步驟，一直是INSAR數據處理中的難點和熱點問題，常用的方法有枝切法、最小費用流等方法，相位解纏處理后的相位圖如圖5所示。

圖5 解纏后的干涉圖

2.7 基線估計與干涉參數定標

為了精確地進行基線估計或干涉參數定標，本次實驗是利用的已有的地面控制點文件，計算相應地面的高程信息。

2.8 DEM重建

通過對干涉圖像的配準、干涉圖生成、去平地效應、干涉圖濾波、相位解纏和相位高程轉換等處理獲得了地面點的高程信息后，根據主（輔）影像的像點坐標和定標參數，經過投影變換，把高程模型投影到相應坐標系下，最后根據要求進行格網化間距重采樣得到標準的DEM，DEM成果。

3 結語

（1）通過實驗分析表明，利用INSAR技術進行DEM的快速提取，已成為解決困難地區快速獲取地形圖和更新的重要手段之一，該方法可靠，有效。

（2）INSAR技術應用尚面臨諸多困難：例如高山區域和城市地區的干涉處理與精確高程反演困難。例如疊掩、陰影的存在，容易導致干涉處理困難、高程精度低等情況。

（3）基線估計是難點。本次實驗是利用的已知的控制點進行了基線估計或干涉參數定標。對于控制點布設困難地區，干涉測量精度可能會受限。

[1] 肖國超，朱彩英.雷達攝影測量[M].北京：地震出版社，2001.

[2] 靳國，徐青，張紅敏.合成孔徑雷達干涉測量[M]. 國防工業出版社，2014.

[3] 靳國，徐青，何鈺.機載雙天線干涉SAR圖像的自動匹配[J].儀器儀表學報 （zl），2006：794-795.

[4] 廖明生，林琿.雷達干涉測量——原理與信號處理基礎[M].北京：測繪出版社，2003.