基于混合迭代的干涉相位仿真方法
【專利摘要】本發明提供一種基于混合迭代的干涉相位仿真方法，通過采用混合迭代的方式，保證了在進行DEM模型匹配的過程中不出現迭代發散，并能夠提高迭代的速度，從而提高了運算效率，節約了運算時間。
【專利說明】基于混合迭代的干涉相位仿真方法
【技術領域】
[0001] 本發明涉及干涉相位圖仿真領域，特別涉及一種基于混合迭代的干涉相位仿真方 法。
【背景技術】
[0002] 干涉相位的模擬是基于干涉相位與兩次成像的路徑差的比例關系來進行的。在已 知DEM數據以及衛星軌道參數的前提下，干涉相位仿真的關鍵在于尋找目標點在雷達成像 時刻的主從衛星的位置、以及地理坐標位置之間三者之間的對應點，從而可根據路徑差及 雷達波長求出干涉相位。
【發明內容】
[0003] 有鑒于此，本發明提供一種基于混合迭代的干涉相位仿真方法，通過采用混合迭 代的方式，保證了在進行DEM模型匹配的過程中不出現迭代發散，并能夠提高迭代的速度， 從而提高了運算效率，節約了運算時間。
[0004] 本發明提供了一種基于混合迭代的干涉相位仿真方法，包括：
[0005] 步驟1 :分別利用設定的η個假定高程V1-Vn，生成主衛星雷達圖像坐標系中柵 格點的地理坐標系中的前向地理編碼多項式系數{ck}和{dk}，其中，
【權利要求】
1. 一種基于混合迭代的干涉相位仿真方法，其特征在于，所述方法包括： 步驟1 :分別利用設定的η個假定高程Vi···Vn，生成主衛星雷達圖像坐標系中柵格點 的地理坐標系中的前向地理編碼多項式系數{ck}和{dk}，其中，
U〇,ξ〇)為柵格點對應的經緯度值，1彡1彡η; 步驟2:利用步驟1得到的{ck}和{dk}構造匹配目標函數fQO=ζ(ζ(hQ)，Uhtl))-Iv并劃定解區間 丨[/?)] ，其中z(ζ，ξ)為DEM模型，解區間 [4〇)，們I為DEM的最小值和最大值； 步驟3:任選E[/r，/r]42)e [#)，)] ，利用
對目標函數f0〇 進行反二次插值，其中0 = ^=M2))/s=M2))/Μ0))； 步驟4 :判斷步驟3得到的插值結果是否使得目標函數fGO的絕對值小于門限值，若 是，則DEM模型匹配成功，令實際高程值Zz= 43)，然后進入步驟9,否則判斷目標函數fQO 的絕對值是否減小且解區間減少一半，如果是，則以插值結果更新解區間邊界，然后返回步 驟3;否則進入步驟5; 步驟5 :任選/^1}e ，利用
目標函數進行割線插值； 步驟6 :判斷步驟5得到的插值結果是否使得目標函數fQO的絕對值小于門限值，若 是，則DEM模型匹配成功，令實際高程值/2 =g2),進入步驟9,否則判斷目標函數fQO的 絕對值是否減小且解區間減少一半，如果是，則以插值結果更新解區間邊界，然后返回步驟 3 ;否則進入步驟7 ; 步驟7:利用=(hT+e；)/ 2對目標函數進行均值插值； 步驟8 :判斷步驟7得到的插值結果是否使得目標函數fQO的絕對值小于門限值，若 是，則DEM模型匹配成功，令實際高程值/? = ,進入步驟9,否則以插值結果更新解區間 邊界，然后返回步驟3; 步驟9 :根據DEM匹配成功時的實際高程值h，根據主衛星軌道參數、主衛星位置矢量S1，通過前向地理編碼方法計算目標點的笛卡爾坐標； 步驟10 :根據步驟9得到的目標點的笛卡爾坐標，以及從衛星軌道參數，采用后向地理 編碼計算目標點在從圖像中的坐標，并由此獲得目標點在成像時刻所對應的從衛星的位置 矢量S2 ; 步驟11 :根據主衛星的位置S1、從衛星的位置S2以及目標點的笛卡爾坐標，計算目標點 對主從天線的波程差，得到干涉相位值。
【文檔編號】G01S7/40GK104316921SQ201410539166
【公開日】2015年1月28日 申請日期:2014年9月27日 優先權日:2014年9月27日
【發明者】勵盼攀 申請人:勵盼攀