一种基于插值的去平地效应方法
【技术领域】
[0001] 本发明涉及干涉相位图的优化领域，特别涉及一种基于插值的去平地效应方法。
【背景技术】
[0002] 由于InSAR本身受空间几何关系的影响，处于方位向或距离向上高度相同的点在 干涉图上使本应保持不变的相位差发生了变化，造成干涉纹图的相位偏移，使干涉相位差 无法反映实际地面的高程变化，这种现象称为平地效应。由于平地相位的存在会使干涉条 纹密集化，增加相位解缠的难度，因此，在进行相位解缠前，需要进行平地效应消除。
【发明内容】
[0003] 有鉴于此，本发明提供一种新的基于插值的去平地效应方法，通过对频谱最大值 段进行三次样条插值使得图像的频率间隔更小，从而获得更加准确的峰值位置，并获得更 加准确的条纹频率，从而能够更加彻底地去平地效应。
[0004] 本发明提供了基于插值的去平地效应方法，包括：
[0005] 步骤1 :对复干涉相位图奶〇'，/)的每行进行傅里叶变换，求取频谱，并按幅度求 和，即：
并计算其最大值：Bi=max[X1(k)];
[0006] 步骤2 :取所述最大值队处左右各两个象素点，共计五个点作为距离向待插值向 量；
[0007] 步骤3:采用三次样条插值算法对所述距离向待插值向量中每两个点之间进行 100个样点的插值计算，得到由400个样点组成的距离向插值向量￥ 1〇1);
[0008] 步骤4 :求取插值向量Yi(n)中的最大值Ci=max[Yi(n)]作为频谱幅度最大值；
[0009] 步骤5 :确定所述频谱幅度最大值Q对应的条纹频率作为距离向的峰值频率，然 后对距离向频谱进行圆周位移，使得条纹的距离向峰值频率移到频谱的零频处；
[0010] 步骤6 :将圆周位移后的距离向频谱进行傅里叶反变换，即可得到去距离向平地 的干涉相位图(/，_/);
[0011] 步骤7:对去距离向平地的干涉相位图死(/,_/)的每列进行傅里叶变换，求 取频谱，并按幅度求和，即：
，并计算其最大值：B2 = max[X2 (k)]；
[0012] 步骤8 :取所述最大值B2处左右各两个象素点，共计五个点作为方位向待插值向 量；
[0013] 步骤9:采用三次样条插值算法对所述方位向待插值向量中每两个点之间进行 100个样点的插值计算，得到由400个样点组成的方位向插值向量￥2〇1);
[0014] 步骤10:求取插值向量Y2 (n)中的最大值C2=max[Y2 (n)]作为频谱幅度最大值；
[0015] 步骤11:确定所述频谱幅度最大值C2对应的条纹频率作为方位向的峰值频率，然 后对方位向频谱进行圆周位移，使得条纹的方位向峰值频率移到频谱的零频处；
[0016] 步骤12 :将圆周位移后的方位向频谱进行傅里叶反变换，即可得到去方位向及去 距离向平地的干涉相位图
[0017] 其中，复干涉相位图的行方向为方位向，列方向为距离向。
【具体实施方式】
[0018] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将对本发明实施例中 的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是 全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前 提下所获得的所有其他实施例，都属于本发明保护的范围。
[0019] 本发明提供一种基于插值的去平地效应方法，包括：
[0020] 步骤1 :对复干涉相位图灼（/,_/)的每行进行傅里叶变换，求取频谱，并按幅度求 和，即：
，并计算其最大值：=max [X丨(k)];
[0021] 步骤2 :取所述最大值队处左右各两个象素点，共计五个点作为距离向待插值向 量；
[0022] 步骤3 :采用三次样条插值算法对所述距离向待插值向量中每两个点之间进行 100个样点的插值计算，得到由400个样点组成的距离向插值向量￥ 1〇1);
[0023] 步骤4:求取插值向量Yi(n)中的最大值C1=max[Yi(n)]作为频谱幅度最大值；
[0024] 步骤5 :确定所述频谱幅度最大值Q对应的条纹频率作为距离向的峰值频率，然 后对距离向频谱进行圆周位移，使得条纹的距离向峰值频率移到频谱的零频处；
[0025] 步骤6 :将圆周位移后的距离向频谱进行傅里叶反变换，即可得到去距离向平地 的干涉相位图^(/，刀；
[0026] 步骤7:对去距离向平地的干涉相位图炉2〇，/)的每列进行傅里叶变换，求 取频谱，并按幅度求和，即：
并计算其最大值：B2 = max [X2 (k)]；
[0027] 步骤8 :取所述最大值B2处左右各两个象素点，共计五个点作为方位向待插值向 量；
[0028] 步骤9:采用三次样条插值算法对所述方位向待插值向量中每两个点之间进行 100个样点的插值计算，得到由400个样点组成的方位向插值向量￥ 2〇1);
[0029] 步骤10:求取插值向量Y2 (n)中的最大值C2=max[Y2 (n)]作为频谱幅度最大值；
[0030] 步骤11 :确定所述频谱幅度最大值C2对应的条纹频率作为方位向的峰值频率，然 后对方位向频谱进行圆周位移，使得条纹的方位向峰值频率移到频谱的零频处；
[0031] 步骤12 :将圆周位移后的方位向频谱进行傅里叶反变换，即可得到去方位向及去 距离向平地的干涉相位图;
[0032] 其中，复干涉相位图的行方向为方位向，列方向为距离向。
[0033] 综上所述，本发明提供的基于插值的去平地效应方法所获得的频谱的峰值位置更 加准确，从而能够更准确的去除平地效应。
[0034] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽 管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然 可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换， 而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和 范围。
【主权项】
1. 一种一种基于插值的去平地效应方法，其特征在于，所述方法包括： 步骤1:对复干涉相位图奶的每行进行傅里叶变换，求取频谱，并按幅度求和， BP:卜计算其最大值办=max[XJk)]; 步骤2 :取所述最大值队处左右各两个象素点，共计五个点作为距离向待插值向量； 步骤3 :采用三次样条插值算法对所述距离向待插值向量中每两个点之间进行100个 样点的插值计算，得到由400个样点组成的距离向插值向量Yi(n); 步骤4 :求取插值向量Yi(n)中的最大值Ci=max[Yi(n)]作为频谱幅度最大值； 步骤5 :确定所述频谱幅度最大值q对应的条纹频率作为距离向的峰值频率，然后对 距离向频谱进行圆周位移，使得条纹的距离向峰值频率移到频谱的零频处； 步骤6 :将圆周位移后的距离向频谱进行傅里叶反变换，即可得到去距离向平地的干 涉相位图朽仏乃； 步骤7 :对去距离向平地的干涉相位图灼（/,_/)的每列进行傅里叶变换，求取频谱，并 按幅度求和，即并计算其最大值：B2=max[X2(k)]; 步骤8 :取所述最大值B2处左右各两个象素点，共计五个点作为方位向待插值向量； 步骤9 :采用三次样条插值算法对所述方位向待插值向量中每两个点之间进行100个 样点的插值计算，得到由400个样点组成的方位向插值向量Y2(n); 步骤10 :求取插值向量Y2 (n)中的最大值C2=max[Y2 (n)]作为频谱幅度最大值； 步骤11 :确定所述频谱幅度最大值C2对应的条纹频率作为方位向的峰值频率，然后对 方位向频谱进行圆周位移，使得条纹的方位向峰值频率移到频谱的零频处； 步骤12 :将圆周位移后的方位向频谱进行傅里叶反变换，即可得到去方位向及去距离 向平地的干涉相位图; 其中，复干涉相位图的行方向为方位向，列方向为距离向。
【专利摘要】本发明提供一种基于插值的去平地效应方法，所述方法通过对频谱的幅度最大值段进行三次样条插值使得图像的频率间隔更小，从而获得更加准确的峰值位置，并获得更加准确的条纹频率，从而能够更加彻底地去平地效应。
【IPC分类】G01S13/90
【公开号】CN105022061
【申请号】CN201510452272
【发明人】黄琴
【申请人】宁波高新区宁源科技服务有限公司
【公开日】2015年11月4日
【申请日】2015年7月24日