专利名称：基于差分滤波的多波段InSAR相位解缠方法
技术领域：
本发明涉及多波段InSAR(Multi-band Interferometric Synthetic ApertureRadar，多波段合成孔径雷达干涉测量)技术，特别涉及一种基于差分滤波的多波 段InSAR相位解缠方法。
背景技术：
利用多波段InSAR技术，可将不同波段(频率)的InSAR数据进行联合处理， 提高对地形陡峭地区混叠干涉相位的解缠精度和可靠性，从而提高获取DEM(Digital Elevation Model，数字高程模型)的精度和可靠性。对于多波段InSAR技术而言，相位解缠仍然是干涉处理的一个关键环节。传统的 相位解缠方法大都基于如下假设相邻像素的相位差不超过半个周期η。而实际上，对于 地形起伏较大地区并采用较短波长、较长基线进行干涉成像的情况，很容易出现干涉相位 的欠采样和干涉条纹的混叠。对于这种干涉条纹，很难甚至不能进行高效、高精度的相位解缠。

发明内容
本发明的目的是提供一种由最长波长开始的差分滤波多波段InSAR相位解缠方 法，利用该方法可以由最长波长的干涉图解缠开始，把不同波段的干涉数据进行组合处理， 逐步精化解缠结果，从而提高相位解缠的可靠性和精度。为实现上述目的，本发明采用以下技术方案由最长波长的干涉数据出发，在对长波长、低频率干涉图进行有效解缠的基础上， 根据波长之间的比例关系，由解缠的低频率干涉图生成短波长、高频率干涉图的基准干涉 图，再对短波长、高频率干涉图与基准干涉图进行差分处理得到差分干涉图，从而降低短波 长、高频率干涉图的条纹频率；对差分干涉图进行滤波和相位解缠之后，再由基准干涉图和 差分干涉图的解缠结果得到短波长、高频率干涉图的解缠结果；并依此类推，一直到最短波 长的干涉图解缠完毕。本发明包括以下步骤①获取所有波长的干涉图Interi,其中，i = 1,2,…，N，N表示波段数，并对最长 波长、最低频率干涉图Inter1进行滤波处理；②对滤波后的最长波长、最低频率干涉图Inter1进行相位解缠，得到最长波长、最 低频率干涉图Inter1的解缠结果UnInter1 ；③根据较长波长、较低频率干涉图Interi的解缠结果UnInteri构造较短波长、较 高频率干涉图Interj的基准干涉图，其中，i = 1,2,…，N，N表示波段数；j = 2，…， N，且j > i ；当由最短波长、最低频率干涉图的解缠结果UnInter1构造基准干涉图时，i = 1 ；④对较短波长、较高频率干涉图Interj与基准干涉图进行差分处理，得到差分干涉图 DInterj DInterj = Interj-Inter]；⑤对差分干涉图DInterj进行滤波和相位解缠之后，再由基准干涉图和差分干涉 图的解缠结果得到较短波长、较高频率干涉图的解缠结果；⑥根据较短波长、较高频率干涉图的解缠结果，进而对更短波长、更高频率干涉图 进行解缠，并依此类推，直到最短波长的干涉图解缠完毕。所述的步骤①中，采用频率域低通滤波方法对最长波长、最低频率干涉图Inter1 进行滤波处理。在所述的步骤①中，采用空间域矢量滤波方法对最长波长、最低频率干涉图 Inter1进行滤波处理。在所述的步骤②中，采用质量图区域生长、直接相位积分和移动曲面相位拟合相 结合的自适应相位解缠方法对最长波长、最低频率干涉图Inter1进行相位解缠。采用上述技术方案的本发明，能够较好将不同波段干涉相位进行组合处理，充分 利用较长波长干涉图频率较低、相位解缠较容易的优势和较短波长干涉图频率较高、精度 较高的优势，有效解决较短波长情况下陡峭地区混叠干涉条纹的相位解缠难题，提高相位 解缠的可解性和解缠精度。


图1为本发明的整体流程图；图2为本发明中采用频域低通滤波方法对最长波长、最低频率干涉图Inter1进行 滤波处理的流程图；图3为本发明中采用空间域矢量滤波方法对最长波长、最低频率干涉图Inter1进 行滤波处理的流程图；图4为本发明中采用质量图区域生长、直接相位积分和移动曲面相位拟合相结合 的自适应相位解缠方法对最长波长、最低频率干涉图Inter1和差分干涉图DInterj(j = 2，…，N)进行相位解缠的处理流程图；图5为SRTM90米间距的DEM图；图6为实验中波长=0.18m的干涉图；图7为实验中波长=0.09m的干涉图；图8为实验中波长=0.06m的干涉图；图9为波长=0. 18m干涉图的滤波结果；图10为波长=0. 18m干涉图的解缠结果；图11为根据波长=0. 18m干涉图解缠结果而构建的波长=0. 09m的基准干涉图；图12为波长=0. 09m的差分干涉图；图13为波长=0. 09m差分干涉图的滤波结果；图14为波长=0. 09m干涉图的解缠结果；图15为根据波长=0. 09m干涉图解缠结果构建的波长=0. 06m的基准干涉图；图16为波长=0. 06m的差分干涉图；图17为波长=0. 06m差分干涉图的滤波结果；图18为波长=0. 06m干涉图的解缠结果。
具体实施例方式实施例1对于不同波段的干涉数据而言，其干涉处理的难度和精度不同。当波长较长时，干 涉图频率较低，干涉相位解缠容易，但获取的高程精度较低；反之，当波长较短时，理论上能 够获取较高的高程精度，但此时容易出现干涉相位欠采样，其干涉质量较低，很难对干涉图 进行有效滤波和相位解缠处理，无法进行有效地干涉测量和高程信息获取。为了充分利用不同波段干涉数据的各自优势，本发明将不同波段的干涉数据进行 组合处理，用较长波长干涉图的解缠结果指导较短波长干涉图的相位解缠，提高了干涉相 位的可解性和解缠精度，提高获取DEM的能力和精度。具体地说，本发明由最长波长的干涉数据出发，在最长波长、最低频率干涉图进行 有效解缠的基础上，根据波长之间的比例关系，由最低频率干涉图的解缠结果构建较短波 长、较高频率干涉图的基准干涉图，再对较短波长、较高频率干涉图与基准干涉图进行差分 处理得到差分干涉图，从而降低较短波长、较高频率干涉图的条纹频率，并解决因地形起伏 较大造成的相位欠采样问题；对差分干涉图进行滤波和相位解缠之后，再由基准干涉图和 差分干涉图的解缠结果得到更短波长、更高频干涉图的解缠结果；并依此类推，一直到最短 波长的干涉图解缠完毕，其具体的实施步骤如下，如图1所示①获取所有干涉图Interi,其中，i = 1,2,…，N，N表示波段数，并采用频率域低 通滤波方法对最长波长、最低频率干涉图Inter1进行滤波处理。其中，频率域低通滤波方 法的流程图如图2所示，在滤波时首先对最长波长、最低频率干涉图Inter1进行二维离散 傅里叶变换，变换到频域后采用低通滤波器进行频域滤波，再进行二维离散傅里叶逆变换 得到干涉图的滤波结果。②对滤波后的最长波长、最低频率干涉图Inter1，可采用质量图区域生长、直接相 位积分和移动曲面相位拟合相结合的自适应相位解缠方法，先由质量图区域生长方法对高 质量的干涉相位进行解缠，再采用移动曲面相位拟合方法完成整幅干涉图的相位解缠。该 相位解缠流程如图4所示。具体地说，为了对滤波后的最长波长、最低频率干涉图Inter1进行高质量的相位
解缠，需要先计算干涉质量图，干涉质量图的计算可采用伪相干值计算公式 其中，ρ ‘ (χ, y)表示对应于干涉图像元(X，y)处的伪相干值； Φ,力表示相应 处的干涉相位值；2k+l表示计算窗口的宽度。对干涉质量图进行区域生长时，需首先给定干涉质量阈值，然后由人工给定的干 涉质量较好的种子像元出发，根据质量图生长出干涉质量较好的区域，并对该生长区域采 用直接相位积分方法进行相位解缠。直接相位积分方法的步骤如下⑴计算干涉图相邻像元间的相位差分-.I^p = φΜ ；(2)
(ii)计算差分修正值 (iii)直接相位积分，即-.φΜ = φ丨+ Δ炉。在进行移动曲面相位拟合时，可采用二次曲面进行解缠相位的拟合 移动曲面相位拟合过程如下(a)给定待拟合相位点P (x0, y0)；(b)确定点P(XQ，yQ)的邻域范围。为了简化计算过程，可以点P(XQ，yQ)为中心，给 定一个7X7的矩形窗口作为邻域范围；(c)将选取的邻域范围内的点坐标(Xi，yi)规划到以P(XQ，yQ)为原点的坐标系，即
按下式计算 此时有 (d)计算邻域内属于非屏蔽区的点数n，如果点数η > 6，列误差方程，求解移动曲 面多项式的系数。移动曲面系数的计算公式为
/-N 6 /Iv
) 7
.S-、

X = (ATPA) ^1ATPL 其中
(e)根据移动曲面系数求解P(X(1，y0)的解缠相位值少( ，凡)=厂；
(f)循环(a)-(e)步，直至屏蔽区内所有点都采用移动曲面拟合出相位解缠值。 需要说明的是，除此之外，也可采用枝切法、最小二乘法等其它方法进行相位解
缠，这些均为本领域普通技术人员所熟知的技术。 ③由于在基线长度、基线水平角均相同的条件下，较长波长、较低频率干涉图 Interi, i = 1，2，…，N，和较短波长、较高频率干涉图Inter」j = 2，…，N，且j > i之间对应像素的解缠干涉相位微分存在如下关系 即任意两幅干涉图上对应的解缠干涉相位微分值与波长成反比。由此，可利用较 长波长、较低频率干涉图Interi的解缠结果UnInteri按其相位微分关系构造较短波长、较 高频率干涉图Interj的基准干涉图，其中j = 2，…，N ； j > i，逐像素的干涉相位值 UnInter^ 为: 为了显示基准干涉图的构建效果，可对其进行模处理。需要说明的是，在公式 (9)和(10)中，i = 1,2,…，N，N为波段数，且当由最短波长、最低频率干涉图的解缠结果 UnInter1构造基准干涉图时，i = 1。④在步骤③的基础上，对较短波长、较高频率干涉图Interj与基准干涉图/nier；进 行差分处理，得到差分干涉图DInterj. DInterj = Interj -Inter0j。这样，通过差分处理，可将
较短波长、较高频率干涉图Intei^由较高频率转化为较低频率，从而降低干涉相位滤波的 难度和干涉相位解缠的难度。⑤对差分干涉图DIntei^进行滤波和相位解缠，得到解缠结果UnDInteiv最后再 由基准干涉图和差分干涉图的解缠结果UnDInterj计算得到较短波长、较高频干涉图的解
缠结果
UnInterj = UnDInterj + UnInter](11)⑥根据较短波长、较高频率干涉图的解缠结果，进而对更短波长、更高频率干涉图 进行解缠，并依此类推，到最短波长的干涉图解缠完毕。实施例2本实施例与实施例1不同的是，在本实施例的步骤①中，采用空间域矢量滤波方 法对最长波长、最低频率干涉图Inter1进行滤波处理。具体地说，如图3所示，它包括以下步骤(a)将干涉相位值映射为矢量空间中的单位矢量。设表示干涉图中
某像元处含噪声的干涉相位值，将其映射到矢量空间中的单位矢量后为
= cos(^(Xjy) + y)) + . η((ρ(χ^y) + y)) ； (12)(b)把矢量空间的干涉图分解为正弦分量Sin(树X，y) + n(x, y))和余弦分量 cosOO，力+；(c)分别对正弦分量 sin ( Φ (x,y)+n(x,y))和余弦分量 cos (Φ (χ, y)+η (χ, y))进 行均值滤波处理；
(d)由滤波后的正弦分量Sil^(U)和余弦分量COS^(U)计算滤波后的干涉图，令 整个矢量均值滤波可用上式表述，其中AngleO表示计算角度值 其他技术特征与实施例1相同。为验证本发明所达到的技术效果，特作以下实验图5为SRTM 90米格网间距的 DEM图，图6、图7、图8分别为由DEM仿真的不同波段的含噪声干涉图，表1为所采用DEM的 相关参数及多波段InSAR系统参数。
表权利要求
一种基于差分滤波的多波段InSAR相位解缠方法，其特征在于该方法由最长波长的干涉数据出发，在对长波长、低频率干涉图进行有效解缠的基础上，根据波长之间的比例关系，由解缠的低频率干涉图生成短波长、高频率干涉图的基准干涉图，再对短波长、高频率干涉图与基准干涉图进行差分处理得到差分干涉图，从而降低短波长、高频率干涉图的条纹频率；对差分干涉图进行滤波和相位解缠之后，再由基准干涉图和差分干涉图的解缠结果得到短波长、高频率干涉图的解缠结果；并依此类推，一直到最短波长的干涉图解缠完毕。
2.根据权利要求1所述的基于差分滤波的多波段InSAR相位解缠方法，其特征在于，它 包括以下步骤①获取所有波长的干涉图Interi,其中，i= 1,2,…，N，N表示波段数，并对最长波长、 最低频率干涉图Inter1进行滤波处理；②对滤波后的最长波长、最低频率干涉图Inter1进行相位解缠，得到最长波长、最低频 率干涉图Inter1的解缠结果UnInter1 ；③根据较长波长、较低频率干涉图Interi的解缠结果UnInteri构造较短波长、较高频 率干涉图Interj的基准干涉图/nier;，其中，i = 1，2，…，N，N表示波段数；j = 2，…，N， 且j > i ；当由最短波长、最低频率干涉图的解缠结果UnInter1构造基准干涉图时，i = 1 ；④对较短波长、较高频率干涉图Interj与基准干涉图进行差分处理，得到差分干 涉图 DInterj -.DInterj = Inter厂Inter;；⑤对差分干涉图DInterj进行滤波和相位解缠之后，再由基准干涉图和差分干涉图的 解缠结果得到较短波长、较高频率干涉图的解缠结果；⑥根据较短波长、较高频率干涉图的解缠结果，进而对更短波长、更高频率干涉图进行 解缠，并依此类推，直到最短波长的干涉图解缠完毕。
3.根据权利要求2所述的基于差分滤波的多波段InSAR相位解缠方法，其特征在于 在所述的步骤①中，采用频率域低通滤波方法对最长波长、最低频率干涉图Inter1进行滤 波处理。
4.根据权利要求2所述的基于差分滤波的多波段InSAR相位解缠方法，其特征在于 在所述的步骤①中，采用空间域矢量滤波方法对最长波长、最低频率干涉图Inter1进行滤 波处理。
5.根据权利要求2所述的基于差分滤波的多波段InSAR相位解缠方法，其特征在于 在所述的步骤②中，采用质量图区域生长、直接相位积分和移动曲面相位拟合相结合的自 适应相位解缠方法对滤波后的最长波长、最低频率干涉图Inter1进行相位解缠。
全文摘要
一种基于差分滤波的多波段InSAR相位解缠方法，该方法由最长波长的干涉数据出发，在对长波长、低频率干涉图进行有效解缠的基础上，根据波长之间的比例关系，由解缠的低频率干涉图生成短波长、高频率干涉图的基准干涉图，再对短波长、高频率干涉图与基准干涉图进行差分处理得到差分干涉图，从而降低短波长、高频率干涉图的条纹频率；对差分干涉图进行滤波和相位解缠之后，再由基准干涉图和差分干涉图的解缠结果得到短波长、高频率干涉图的解缠结果；并依此类推，一直到最短波长的干涉图解缠完毕。本发明有效解决了较短波长情况下陡峭地区混叠干涉条纹的相位解缠难题，提高相位解缠的可解性和解缠精度。
文档编号G01S13/90GK101881831SQ201010208448
公开日2010年11月10日 申请日期2010年6月24日 优先权日2010年6月24日
发明者周杨, 张红敏, 徐青, 秦志远, 靳国旺, 龚志辉 申请人:中国人民解放军信息工程大学