一种宽带合成孔径雷达子带干涉数据处理方法
【专利摘要】一种宽带合成孔径雷达子带干涉数据处理方法，该方法有六大步骤：步骤一：主副图像的自适应精密配准及副图像的重采样；步骤二：子带频谱位置和带宽的优化；步骤三：主副图像的带通滤波；步骤四：子带干涉图和子带差分干涉图的生成；步骤五：参考面和地形高程(DEM)的相位补偿；步骤六：子带差分干涉图的堆栈处理。本发明在大多数情况下无需相位解缠，即可获取绝对雷达干涉相位，从而可以通过转换直接得到数字地形图或者地表形变图，消除了由于相位解缠错误引起的误差扩散；它特别适用于高山地区的数字地形测绘和大形变速率的滑坡及矿山塌陷监测。
【专利说明】一种宽带合成孔径雷达子带干涉数据处理方法
【技术领域】
[0001]本发明涉及一种宽带合成孔径雷达子带干涉数据处理方法，属于合成孔径雷达干涉【技术领域】。它可以在大多数情况下无需相位解缠，即可获取绝对雷达干涉相位，从而可以通过转换直接得到数字地形图或者地表形变图。
【背景技术】
[0002]合成孔径雷达干涉是一门根据高分辨率雷达图像的相位数据来提取地面目标三维空间信息的技术。无论是利用雷达干涉进行地形数字高程测量还是地表形变测量，都有一个相位解缠的过程。因为雷达干涉图的相位都被一个2 的角度所缠绕，所谓相位解缠，就是通过累积雷达干涉图像元间的相位梯度，从而获得像元绝对相位的过程。这一过程要能够得到正确结果，对于地形测量来说，一个必要条件是，在空间域上，相邻两像元间的相位差的绝对值不得大于;对于利用堆栈差分干涉图的时序分析方法的形变测量，除了前面那个必要条件必须得到满足外，在时间域上，同一像元相邻两时间点上的相位差的绝对值亦不得大于n。在实际工程中，尤其是在高山地区的地形测量时，或者是对具有较大形变速率梯度的滑坡和矿山的形变场进行干涉测量时，上面提到的相位差的绝对值不得大于
的条件常常是不能得到满足的，因而造成测量结果出现很大误差，有时甚至会出现无解的情况。根据雷达干涉理论，要解决这一相位解缠的问题，就是要想办法提高干涉相位测量模糊度。对于地形测量而言，其模糊度定义为:
[0003]
【权利要求】
1.一种宽带合成孔径雷达子带干涉数据处理方法，其特征在于:该方法具体步骤如下: 步骤一:主副图像的自适应精密配准及副图像的重采样 由于子带干涉技术主要是应用在高山地区的高程测量或者是具有较大形变速率的滑坡和矿山塌陷测量，在主副图像进行配准时，不仅要考虑雷达干涉系统引起的地面同一目标在主副图像上的位置的偏移，还要考虑地面目标因较大的高程或较大的形变位移引起的主副图像上的位置偏移；在这里，副图像利用Sinc函数进行重采样，Sinc函数的位移量由下式决定:
Offset(i, j) = ( A i, Aj) = (( Ai, A j)i；J)System+(( A i, A j)i;Pobject (3) 上式中，offsets j)表示了副图像中位于第i行第j列的像元相对于主图像上同一像兀位置的偏移量，它是系统偏移和目标本身因闻程或形变引起的偏移量的和；offset(i, j)由副图像s中位于(i，j)的像元的邻域与相应的主图像p中像元的邻域的协方差函数取最大值时的位移量决定，见下面的(4)和(5)式
【文档编号】G01S7/41GK103728604SQ201310582414
【公开日】2014年4月16日 申请日期:2013年11月19日 优先权日:2013年11月19日
【发明者】夏耶, 郭小方, 范景辉, 李曼 申请人:中国国土资源航空物探遥感中心