基于混合迭代的干涉相位仿真方法
【專利摘要】本發(fā)明提供一種基于混合迭代的干涉相位仿真方法，通過(guò)采用混合迭代的方式，保證了在進(jìn)行DEM模型匹配的過(guò)程中不出現(xiàn)迭代發(fā)散，并能夠提高迭代的速度，從而提高了運(yùn)算效率，節(jié)約了運(yùn)算時(shí)間。
【專利說(shuō)明】基于混合迭代的干涉相位仿真方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及干涉相位圖仿真領(lǐng)域，特別涉及一種基于混合迭代的干涉相位仿真方 法。

【背景技術(shù)】
[0002] 干涉相位的模擬是基于干涉相位與兩次成像的路徑差的比例關(guān)系來(lái)進(jìn)行的。在已 知DEM數(shù)據(jù)以及衛(wèi)星軌道參數(shù)的前提下，干涉相位仿真的關(guān)鍵在于尋找目標(biāo)點(diǎn)在雷達(dá)成像 時(shí)刻的主從衛(wèi)星的位置、以及地理坐標(biāo)位置之間三者之間的對(duì)應(yīng)點(diǎn)，從而可根據(jù)路徑差及 雷達(dá)波長(zhǎng)求出干涉相位。


【發(fā)明內(nèi)容】

[0003] 有鑒于此，本發(fā)明提供一種基于混合迭代的干涉相位仿真方法，通過(guò)采用混合迭 代的方式，保證了在進(jìn)行DEM模型匹配的過(guò)程中不出現(xiàn)迭代發(fā)散，并能夠提高迭代的速度， 從而提高了運(yùn)算效率，節(jié)約了運(yùn)算時(shí)間。
[0004] 本發(fā)明提供了一種基于混合迭代的干涉相位仿真方法，包括：
[0005] 步驟1 :分別利用設(shè)定的η個(gè)假定高程V1-Vn，生成主衛(wèi)星雷達(dá)圖像坐標(biāo)系中柵 格點(diǎn)的地理坐標(biāo)系中的前向地理編碼多項(xiàng)式系數(shù){ck}和{dk}，其中，

【權(quán)利要求】
1. 一種基于混合迭代的干涉相位仿真方法，其特征在于，所述方法包括： 步驟1 :分別利用設(shè)定的η個(gè)假定高程Vi···Vn，生成主衛(wèi)星雷達(dá)圖像坐標(biāo)系中柵格點(diǎn) 的地理坐標(biāo)系中的前向地理編碼多項(xiàng)式系數(shù){ck}和{dk}，其中，
U〇,ξ〇)為柵格點(diǎn)對(duì)應(yīng)的經(jīng)緯度值，1彡1彡η; 步驟2:利用步驟1得到的{ck}和{dk}構(gòu)造匹配目標(biāo)函數(shù)fQO=ζ(ζ(hQ)，Uhtl))-Iv并劃定解區(qū)間 丨[/?)] ，其中z(ζ，ξ)為DEM模型，解區(qū)間 [4〇)，們I為DEM的最小值和最大值； 步驟3:任選E[/r，/r]42)e [#)，)] ，利用
對(duì)目標(biāo)函數(shù)f0〇 進(jìn)行反二次插值，其中0 = ^=M2))/s=M2))/Μ0))； 步驟4 :判斷步驟3得到的插值結(jié)果是否使得目標(biāo)函數(shù)fGO的絕對(duì)值小于門限值，若 是，則DEM模型匹配成功，令實(shí)際高程值Zz= 43)，然后進(jìn)入步驟9,否則判斷目標(biāo)函數(shù)fQO 的絕對(duì)值是否減小且解區(qū)間減少一半，如果是，則以插值結(jié)果更新解區(qū)間邊界，然后返回步 驟3;否則進(jìn)入步驟5; 步驟5 :任選/^1}e ，利用
目標(biāo)函數(shù)進(jìn)行割線插值； 步驟6 :判斷步驟5得到的插值結(jié)果是否使得目標(biāo)函數(shù)fQO的絕對(duì)值小于門限值，若 是，則DEM模型匹配成功，令實(shí)際高程值/2 =g2),進(jìn)入步驟9,否則判斷目標(biāo)函數(shù)fQO的 絕對(duì)值是否減小且解區(qū)間減少一半，如果是，則以插值結(jié)果更新解區(qū)間邊界，然后返回步驟 3 ;否則進(jìn)入步驟7 ; 步驟7:利用=(hT+e；)/ 2對(duì)目標(biāo)函數(shù)進(jìn)行均值插值； 步驟8 :判斷步驟7得到的插值結(jié)果是否使得目標(biāo)函數(shù)fQO的絕對(duì)值小于門限值，若 是，則DEM模型匹配成功，令實(shí)際高程值/? = ,進(jìn)入步驟9,否則以插值結(jié)果更新解區(qū)間 邊界，然后返回步驟3; 步驟9 :根據(jù)DEM匹配成功時(shí)的實(shí)際高程值h，根據(jù)主衛(wèi)星軌道參數(shù)、主衛(wèi)星位置矢量S1，通過(guò)前向地理編碼方法計(jì)算目標(biāo)點(diǎn)的笛卡爾坐標(biāo)； 步驟10 :根據(jù)步驟9得到的目標(biāo)點(diǎn)的笛卡爾坐標(biāo)，以及從衛(wèi)星軌道參數(shù)，采用后向地理 編碼計(jì)算目標(biāo)點(diǎn)在從圖像中的坐標(biāo)，并由此獲得目標(biāo)點(diǎn)在成像時(shí)刻所對(duì)應(yīng)的從衛(wèi)星的位置 矢量S2 ; 步驟11 :根據(jù)主衛(wèi)星的位置S1、從衛(wèi)星的位置S2以及目標(biāo)點(diǎn)的笛卡爾坐標(biāo)，計(jì)算目標(biāo)點(diǎn) 對(duì)主從天線的波程差，得到干涉相位值。
【文檔編號(hào)】G01S7/40GK104316921SQ201410539166
【公開日】2015年1月28日 申請(qǐng)日期:2014年9月27日 優(yōu)先權(quán)日:2014年9月27日
【發(fā)明者】勵(lì)盼攀 申請(qǐng)人:勵(lì)盼攀