[發明專利]一種干涉相位仿真方法

在審

• 申請號: 201410539149.7
• 申請日: 2014-09-27
• 公開（公告）號: CN104267378A
• 公開（公告）日: 2015-01-07
• 發明（設計）人: 勵盼攀
• 申請（專利權）人: 勵盼攀
• 主分類號: G01S7/02
• 分類號: G01S7/02
• 代理公司: 暫無信息
• 代理人: 暫無信息

技術領域

本發明涉及干涉相位圖仿真領域，特別涉及一種干涉相位仿真方法。

背景技術

干涉相位的模擬是基于干涉相位與兩次成像的路徑差的比例關系來進行的。在已知DEM數據以及衛星軌道參數的前提下，干涉相位仿真的關鍵在于尋找目標點在雷達成像時刻的主從衛星的位置、以及地理坐標位置之間三者之間的對應點，從而可根據路徑差及雷達波長求出干涉相位。

發明內容

有鑒于此，本發明提供一種新的干涉相位仿真方法，通過將星載SAR測繪帶內成像幾何規律等效為n階多項式描述的目標點的地理坐標系中的經緯度值隨高程值的變化規律，從而簡化地理編碼的運算過程，提高了運算效率，節約了運算時間。

本發明提供了一種干涉相位仿真方法，包括：

步驟1：分別利用設定的n個假定高程h_0，1…h_0，n，生成主衛星雷達圖像坐標系中柵格點的地理坐標系中的前向地理編碼多項式系數{c_k}和{d_k}，其中，

ζ=c0+c1h0,l+...+cnh0,lnξ=d0+d1h0,l+...+dnh0,ln,]]>(ζ，ξ)為柵格點對應的經緯度值，1≤l≤n；

步驟2：根據目標點的主衛星雷達圖像坐標計算成像時刻的主衛星位置矢量S₁；

步驟3：設定成像時刻目標點高程值h₀的值；

步驟4：通過步驟1得到的前向地理編碼多項式系數{c_k}和{d_k}，得到高程值h₀對應的目標點的地理坐標系中的經緯度值(ζ₀，ξ₀)；

步驟：5：對DEM模型h＝z(ζ，ξ)進行內插，得到目標點對應的經緯度值(ζ₀，ξ₀)在DEM模型中的實際高程值h，并判斷實際高程值h和設定高程值h₀之差是否小于門限，若小于門限，則DEM模型匹配成功，進入步驟6，否則重新設定成像時刻目標點高程值h₀的值，并返回步驟4；

步驟6：根據DEM匹配成功時的實際高程值h，根據主衛星軌道參數、主衛星位置矢量S₁，通過前向地理編碼方法計算目標點的笛卡爾坐標；

步驟7：根據步驟6得到的目標點的笛卡爾坐標，以及從衛星軌道參數，采用后向地理編碼計算目標點在從圖像中的坐標，并由此獲得目標點在成像時刻所對應的從衛星的位置矢量S₂；