一种基于Chirp-Z变换的去平地效应方法
【专利摘要】本发明提供一种基于Chirp-Z变换的去平地效应方法，所述方法通过先对频谱进行FFT变换，判断出待测频率的主瓣位置；然后再在主瓣范围内进行Chirp-Z变换，细化频谱的主瓣，以获取精确的局部条纹频率，从而能够更加彻底地去平地效应。
【专利说明】-种基于Ch i rp-z变换的去平地效应方法
【技术领域】
[0001] 本发明涉及干涉相位图的优化领域，特别涉及一种基于Chirp-Z变换的去平地效 应方法。
【背景技术】
[0002] 由于InSAR本身受空间几何关系的影响，处于方位向或距离向上高度相同的点在 干涉图上使本应保持不变的相位差发生了变化，造成干涉纹图的相位偏移，使干涉相位差 无法反映实际地面的高程变化，这种现象称为平地效应。由于平地相位的存在会使干涉条 纹密集化，增加相位解缠的难度，因此，在进行相位解缠前，需要进行平地效应消除。
【发明内容】
[0003] 有鉴于此，本发明提供一种新基于Chirp-Z变换的去平地效应方法，通过先对频 谱进行FFT变换，判断出待测频率的主瓣位置；然后再在主瓣范围内进行Chirp-Z变换，细 化频谱的主瓣，以获取精确的局部条纹频率,从而能够更加彻底地去平地效应。
[0004] 本发明提供了一种基于Chirp-Z变换的去平地效应方法，包括：
[0005] 步骤1 :对复干涉相位图奶〇_,/)的每行进行傅里叶变换得到距离向频谱，并按 幅度求和，即：
【权利要求】
1. 一种一种基于Chirp-Z变换的去平地效应方法，其特征在于，所述方法包括： 步骤1 :对复干涉相位图奶(；/,/)的每行进行傅里叶变换得到距离向频谱，并按幅度求 和，即：
I确定频谱主瓣范围并计算最大值：B1 = max [X1GO] 作为频谱幅度最大值； 步骤2 :确定所述频谱幅度最大值B1对应的距离向条纹频率; 步骤3 :对频谱主瓣进行Chirp-Z变换，得到距离向的估计的条纹频率偏移量gl ; 步骤4 :获得距离向的频谱细化后的精确的条纹频率fi+gl作为距离向的峰值频率，然 后对距离向频谱进行圆周位移，使得条纹的距离向峰值频率移到频谱的零频处； 步骤5 :将圆周位移后的距离向频谱进行傅里叶反变换，即可得到去距离向平地的干 涉相位图炉2 (/，/); 步骤6 :对去距离向平地的干涉相位图供2〇',/)的每列进行傅里叶变换得到方位向频 谱，并按幅度求和，即：.
，并计算其最大值：B2 = max [X2(k)] 作为频谱幅度最大值； 步骤7 :确定所述频谱幅度最大值B2对应的方位向条纹频率f2 ; 步骤8 :对频谱主瓣进行Chirp-Z变换，得到方位向的估计的条纹频率偏移量g2 ; 步骤9 :获得方位向的频谱细化后的精确的条纹频率f2+g2作为方位向的峰值频率，然 后对方位向频谱进行圆周位移，使得条纹的方位向峰值频率移到频谱的零频处； 步骤10 :将圆周位移后的方位向频谱进行傅里叶反变换，即可得到去距离向平地的干 涉相位图供(U); 其中，复干涉相位图的行方向为方位向，列方向为距离向。
【文档编号】G01S13/90GK104267398SQ201410539147
【公开日】2015年1月7日 申请日期:2014年9月27日 优先权日:2014年9月27日
【发明者】励盼攀 申请人:励盼攀