一种基于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日