一种DEM辅助的SAR图像高精度配准方法技术

本发明专利技术提供一种DEM辅助的SAR图像高精度配准方法。技术方案是：对同一场景获得的两幅SAR图像，计算主图像各像元对应的高程；然后在主图像上选取适量的像元作为控制点，通过相关系数法得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像和辅图像之间的方位向配准偏移量、距离向配准偏移量；然后通过最小二乘拟合，拟合一个与像元坐标、像元对应高程相关的多项式；最后将主图像各像元的参数代入多项式，即可得到主图像像元与辅图像对应像元之间的方位向、距离向配准偏移量，完成SAR图像配准。本发明专利技术解决了在长基线在长基线和复杂地形区域情况下，现有SAR图像二维函数配准法的局部配准精度低问题。

【技术实现步骤摘要】
一种DEM辅助的SAR图像高精度配准方法
本专利技术涉及微波遥感
，特别涉及一种DEM(DigitalElevationModel数字高程模型)辅助的SAR(SyntheticApertureRadar，合成孔径雷达)图像高精度配准方法。
技术介绍
InSAR(InterferometrySyntheticApertureRadar，合成孔径雷达干涉测量)是一种对地观测技术，它通过对雷达卫星在不同观测几何下获取的多幅SAR复图像进行干涉处理，得到地表的三维信息。SAR图像配准是InSAR处理流程的第一步，SAR图像配准的精度对于后续的处理步骤以及最终的InSAR产品都至关重要，通常要求SAR图像配准的精度达到亚像素级。SAR图像配准时，一般利用两幅SAR图像，可以任意选取其中一幅图像为主图像，另外一幅图像为辅图像。目前最常用的SAR图像配准方法是二维函数法(参见文献：《RegistrationofInSARComplexImagesandInterferogramEnhancement》YANGQING-youetc,RemoteSensingTechnologyandApplication，Vol.3,No.2,May1999，第123页至第125页)。它首先采用某种优化过程(如相关系数最大化、条纹清晰度最大化等)选取主图像、辅图像中的控制点，再以控制点坐标及控制点处的配准偏移量构造二维形变函数来描述主辅图像各像元之间的偏移，进而得到每个像元处的配准偏移量。该方法操作简单，在地形平坦区域和短基线情况下能够达到SAR图像配准的精度要求。但是，由于在长基线情况下，配准偏移量对场景高程起伏敏感，导致真实的像元之间的偏移量起伏与场景地形起伏呈高相关，显然，此时一个二维函数是不足以对真实的配准偏移量进行高精度拟合的。因此在复杂地形区域和长基线情况下，二维函数法会由于拟合精度不够而造成局部配准精度低的现象。
技术实现思路
本专利技术要解决的技术问题是，提供一种DEM辅助的SAR图像高精度配准方法，有效解决了在长基线和复杂地形区域情况下，现有SAR图像二维函数配准法的局部配准精度低现象。本专利技术的技术方案是：已知对同一场景获得的两幅SAR图像，首先选取任意一幅为主图像，则另一幅为辅图像。通过主图像的轨道参数和对应的DEM计算主图像各像元对应的高程；然后在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像和辅图像之间的方位向配准偏移量、距离向配准偏移量；然后再根据控制点的坐标、控制点像元对应的高程，以及控制点的方位向配准偏移量、距离向配准偏移量，通过最小二乘拟合，在方位向、距离向分别得到一个与像元坐标、像元对应高程相关的多项式；最后将主图像各像元的坐标以及像元对应的高程代入多项式，即可得到主图像像元与辅图像对应像元之间的方位向、距离向配准偏移量，完成SAR图像配准。本专利技术的有益效果是：通过DEM和主图像轨道参数计算主图像各像元的高程，并将像元对应高程引入到计算配准偏移量的多项式中，计算得到更高精度的配准偏移量。这样将配准偏移量与SAR图像场景对应的地形关联起来，可以更高精度的拟合真实的配准偏移量，解决了在长基线在长基线和复杂地形区域情况下，现有SAR图像二维函数配准法的局部配准精度低问题。附图说明图1是本专利技术的原理流程示意图；图2为主图像和图像对应高程图；图3(a)为利用本专利技术方法后所得相干系数图，图3(b)为利用二维函数法后所得相干系数图；图4(a)为利用本专利技术方法后所得局部干涉条纹，图4(b)为利用二维函数法后所得局部干涉条纹；具体实施方式图1是本专利技术的原理流程示意图。下面结合图1进行详细说明。步骤一：通过主图像的轨道参数和对应的DEM计算主图像各像元对应的高程。将DEM地理格网上每一点的地理坐标(B,L,H)(经度，纬度，高程)代入到主图像轨道参数和成像参数条件下的距离多普勒模型中进行处理，具体方法参照论文《UtilizationofSpaceborneSARDataforMapping》(IEEETransactionsongeoscienceandremotesensing，Vol.22,No.2,March1984，第107页至第108页)，可求得DEM地理格网的各网点对应在主图像中的位置，从而得到主图像各像元(i,j)(i为方位向坐标，j为距离向坐标)对应高程h(i,j)。步骤二：首先在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像、辅图像之间的方位向配准偏移量、距离向配准偏移量。具体方法参照论文《RegistrationofInSARComplexImagesandInterferogramEnhancement》(RemoteSensingTechnologyandApplication，Vol.3,No.2,May1999，第123页至第125页)。其中控制点的分布原则应满足场景对应地形的各高程区间区域均分布有控制点，控制点数量具体视图像大小而定，一般不少于100个。步骤三：然后根据控制点的坐标、控制点像元对应的高程，以及控制点的配准偏移量，通过最小二乘拟合，在方位向、距离向分别拟合一个与坐标、高程相关的多项式。由上述处理步骤可知各控制点的坐标(i′,j′)，控制点像元对应的高程h(i′,j′)，控制点对应的方位向配准偏移量a(i′,j′)，距离向配准偏移量r(i′,j′)。则拟合的配准偏移量计算多项式为：其中apq、bpq、cpq、dpq(p＝0,1,...n；q＝0,1,...n)为多项式的待定系数，将所有控制点的坐标(i′,j′)、高程h(i′,j′)、方位向配准偏移量a(i′,j′)、距离向配准偏移量r(i′,j′)代入公式(1)，通过最小二乘法拟合得到这些待定系数。m、n为拟合多项式阶数，通常在m＝2、n＝2阶时可以得到足够的精度和较高的运算效率。步骤四：将主图像各像元的坐标(i,j)以及对应的高程h(i,j)代入多项式(1)，即可得到主图像各像元与对应辅图像各像元之间的方位向、距离向配准偏移量，完成SAR图像配准。本专利技术的实施例所用的实测数据如图2(a)所示，图2(a)是一幅高程起伏较大的地形区域的原始的SAR图像，将该图像作为主图像。此SAR图像是从现有的PALSAR系统获得的。辅图像也是从现有的PALSAR系统获得的同一场景的原始SAR图像，只是主图像和辅图像获取时间不同。本实施例的DEM数据采用的是全球公开的SRTM-CDEM数据，DEM格网大小为30米×30米。图2(b)为通过步骤一得到的SAR图像各像元对应的高程图。图3(a)为利用本专利技术方法配准后所得的相干系数，图3(b)为利用现有的二维函数法配准后所得的相干系数图。从图3可以看出利用本专利技术方法配准后所得相干系数整体高于二维函数法配准所得相干系数。尤其在图的上部分，二维函数法所得的相干系数较低，这说明在此区域二维函数法配准误差较大，而本专利技术的方法所得相干系数在该区域明显高于二维函数法，且保持较高的相干系数，说明本专利技术方法的配准精度较高。图4(a)为通过本专利技术的方法配准本文档来自技高网...

【技术保护点】
1.一种DEM辅助的SAR图像高精度配准方法，DEM是指数字高程模型，SAR是指合成孔径雷达，已知对同一场景获得的两幅SAR图像，首先选取任意一幅为主图像，则另一幅为辅图像，其特征在于：通过主图像的轨道参数和对应的DEM计算主图像各像元对应的高程；然后在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像和辅图像之间的方位向配准偏移量、距离向配准偏移量；然后再根据控制点的坐标、控制点像元对应的高程，以及控制点的方位向配准偏移量、距离向配准偏移量，通过最小二乘拟合，在方位向、距离向分别得到一个与像元坐标、像元对应高程相关的多项式；最后将主图像各像元的坐标以及像元对应的高程代入多项式，即可得到主图像像元与辅图像对应像元之间的方位向、距离向配准偏移量，完成SAR图像配准。

【技术特征摘要】
1.一种DEM辅助的SAR图像高精度配准方法，DEM是指数字高程模型，SAR是指合成孔径雷达，已知对同一场景获得的两幅SAR图像，首先选取任意一幅为主图像，则另一幅为辅图像，其特征在于：通过主图像的轨道参数和对应的DEM计算主图像各像元对应的高程；然后在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像和辅图像之间的方位向配准偏移量、距离向配准偏移量；然后再根据控制点的坐标、控制点像元对应的高程，以及控制点的方位向配准偏移量、距离向配准偏移量，通过最小二乘拟合，在方位向、距离向分别得到一个与像元坐标、像元对应高程相关的多项式；最后将主图像各像元的坐标以及像元对应...

【专利技术属性】
技术研发人员：甘本对，余安喜，孙造宇，董臻，张永胜，何峰，何志华，金光虎，张启雷，李德鑫，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：湖南,43