一種基于Chirp-Z變換的去平地效應方法
【技術領域】
[0001] 本發明涉及干涉相位圖的優化領域，特別涉及一種基于Chirp-z變換的去平地效 應方法。
【背景技術】
[0002] 由于InSAR本身受空間幾何關系的影響，處于方位向或距離向上高度相同的點在 干涉圖上使本應保持不變的相位差發生了變化，造成干涉紋圖的相位偏移，使干涉相位差 無法反映實際地面的高程變化，這種現象稱為平地效應。由于平地相位的存在會使干涉條 紋密集化，增加相位解纏的難度，因此，在進行相位解纏前，需要進行平地效應消除。
【發明內容】
[0003] 有鑒于此，本發明提供一種新基于Chirp-z變換的去平地效應方法，通過先對頻 譜進行FFT變換，判斷出待測頻率的主瓣位置；然后再在主瓣范圍內進行Chirp-Z變換，細 化頻譜的主瓣，以獲取精確的局部條紋頻率，從而能夠更加徹底地去平地效應。
[0004] 本發明提供了一種基于Chirp-Z變換的去平地效應方法，包括：
[0005] 步驟1 :對復干涉相位圖爐丨（/,/)的每行進行傅里葉變換得到距離向頻譜，并按 幅度求和，即：
，確定頻譜主瓣范圍并計算最大值辦= max [X! (k)]作為頻譜幅度最大值；
[0006] 步驟2 :確定所述頻譜幅度最大值隊對應的距離向條紋頻率f1;
[0007] 步驟3 :對頻譜主瓣進行Chirp-Z變換，得到距離向的估計的條紋頻率偏移量g1;
[0008] 步驟4 :獲得距離向的頻譜細化后的精確的條紋頻率fi+gl作為距離向的峰值頻 率，然后對距離向頻譜進行圓周位移，使得條紋的距離向峰值頻率移到頻譜的零頻處；
[0009] 步驟5:將圓周位移后的距離向頻譜進行傅里葉反變換，即可得到去距離向平地 的干涉相位圖朽(/,_/)，?
[0010] 步驟6:對去距離向平地的干涉相位圖的每列進行傅里葉變換得到方 位向頻譜，并按幅度求和，即：
并計算其最大值：B2 = max [X2(k)]作為頻譜幅度最大值；
[0011] 步驟7 :確定所述頻譜幅度最大值&對應的方位向條紋頻率f 2;
[0012] 步驟8 :對頻譜主瓣進行Chirp-Z變換，得到方位向的估計的條紋頻率偏移量g2;
[0013] 步驟9 :獲得方位向的頻譜細化后的精確的條紋頻率f2+g2作為方位向的峰值頻 率，然后對方位向頻譜進行圓周位移，使得條紋的方位向峰值頻率移到頻譜的零頻處；
[0014] 步驟10 :將圓周位移后的方位向頻譜進行傅里葉反變換，即可得到去距離向平地 的干涉相位圖
[0015] 其中，復干涉相位圖的行方向為方位向，列方向為距離向。
【具體實施方式】
[0016] 為使本發明實施例的目的、技術方案和優點更加清楚，下面將對本發明實施例中 的技術方案進行清楚、完整的描述，顯然，所描述的實施例是本發明一部分實施例，而不是 全部的實施例。基于本發明中的實施例，本領域普通技術人員在沒有做出創造性勞動的前 提下所獲得的所有其他實施例，都屬于本發明保護的范圍。
[0017] 本發明提供了一種基于Chirp-z變換的去平地效應方法，包括：
[0018] 步驟1:對復干涉相位圖奶〇',/)的每行進行傅里葉變換得到距離向頻譜，并按 幅度求和，即：
商定頻譜主瓣范圍并計算最大值：Bi= max [X! (k)]作為頻譜幅度最大值；
[0019] 步驟2:確定所述頻譜幅度最大值隊對應的距離向條紋頻率f1;
[0020] 步驟3 :對頻譜主瓣進行Chirp-Z變換，得到距離向的估計的條紋頻率偏移量g1;
[0021] 步驟4 :獲得距離向的頻譜細化后的精確的條紋頻率fi+gl作為距離向的峰值頻 率，然后對距離向頻譜進行圓周位移，使得條紋的距離向峰值頻率移到頻譜的零頻處；
[0022] 步驟5 :將圓周位移后的距離向頻譜進行傅里葉反變換，即可得到去距離向平地 的干涉相位圖供2 (/，y);
[0023] 步驟6:對去距離向平地的干涉相位圖的每列進行傅里葉變換得到方 位向頻譜，并按幅度求和，即：
并計算其最大值：B2 = max [X2(k)]作為頻譜幅度最大值；
[0024] 步驟7 :確定所述頻譜幅度最大值&對應的方位向條紋頻率f 2;
[0025] 步驟8 :對頻譜主瓣進行Chirp-Z變換，得到方位向的估計的條紋頻率偏移量g2;
[0026] 步驟9:獲得方位向的頻譜細化后的精確的條紋頻率f2+g2作為方位向的峰值頻 率，然后對方位向頻譜進行圓周位移，使得條紋的方位向峰值頻率移到頻譜的零頻處；
[0027] 步驟10 :將圓周位移后的方位向頻譜進行傅里葉反變換，即可得到去距離向平地 的干涉相位圖
[0028] 其中，復干涉相位圖的行方向為方位向，列方向為距離向。
[0029] 綜上所述，本發明提供的基于Chirp-Z變換的去平地效應方法所獲得的頻譜的峰 值位置更加準確，從而能夠更準確的去除平地效應。
[0030] 最后應說明的是：以上實施例僅用以說明本發明的技術方案，而非對其限制；盡 管參照前述實施例對本發明進行了詳細的說明，本領域的普通技術人員應當理解：其依然 可以對前述實施例所記載的技術方案進行修改，或者對其中部分技術特征進行等同替換， 而這些修改或者替換，并不使相應技術方案的本質脫離本發明各實施例技術方案的精神和 范圍。
【主權項】
1. 一種一種基于Chirp-Z變換的去平地效應方法，其特征在于，所述方法包括： 步驟1 :對復干涉相位圖奶的每行進行傅里葉變換得到距離向頻譜，并按幅度求 和，即，確定頻譜主瓣范圍并計算最大值辦=maxDCiGO] 作為頻譜幅度最大值； 步驟2 :確定所述頻譜幅度最大值隊對應的距離向條紋頻率f1; 步驟3 :對頻譜主瓣進行Chirp-Z變換，得到距離向的估計的條紋頻率偏移量g1; 步驟4 :獲得距離向的頻譜細化后的精確的條紋頻率匕+^作為距離向的峰值頻率，然 后對距離向頻譜進行圓周位移，使得條紋的距離向峰值頻率移到頻譜的零頻處； 步驟5 :將圓周位移后的距離向頻譜進行傅里葉反變換，即可得到去距離向平地的干 涉相位圖(^(/，刀；步驟6 :對去距離向平地的干涉相位圖灼（/,/)的每列進行傅里葉變換得到方位向頻 譜，并按幅度求和，即 f計算其最大值：B2=max[X2(k)] 作為頻譜幅度最大值； 步驟7 :確定所述頻譜幅度最大值&對應的方位向條紋頻率f2; 步驟8 :對頻譜主瓣進行Chirp-Z變換，得到方位向的估計的條紋頻率偏移量g2; 步驟9 :獲得方位向的頻譜細化后的精確的條紋頻率&+&作為方位向的峰值頻率，然 后對方位向頻譜進行圓周位移，使得條紋的方位向峰值頻率移到頻譜的零頻處； 步驟10 :將圓周位移后的方位向頻譜進行傅里葉反變換，即可得到去距離向平地的干 涉相位圖沖".)； 其中，復干涉相位圖的行方向為方位向，列方向為距離向。
【專利摘要】本發明提供一種基于Chirp-Z變換的去平地效應方法，所述方法通過先對頻譜進行FFT變換，判斷出待測頻率的主瓣位置；然后再在主瓣范圍內進行Chirp-Z變換，細化頻譜的主瓣，以獲取精確的局部條紋頻率，從而能夠更加徹底地去平地效應。
【IPC分類】G01S13/90
【公開號】CN105022062
【申請號】CN201510452273
【發明人】黃琴
【申請人】寧波高新區寧源科技服務有限公司
【公開日】2015年11月4日
【申請日】2015年7月24日