专利名称：一种多源传感器的遥感图像配准方法
技术领域：
本发明涉及图像处理技术领域中关于图像配准的方法，是一种针对多源传感器、多时相、多视角遥感图像的高精度、快速配准方法。
背景技术：
图像配准是多时相、多传感器、多视角的两幅或多幅图像匹配、叠加的过程，已经广泛应用于遥感图像处理、目标识别、医学影像分析、机器视觉等领域。在遥感领域，图像配准是实现异源图像融合、异质图像分析、目标变化检测与识别、图像校正、图像镶嵌等的关键步骤。图像配准方法主要分为两类区域配准法和特征配准法。(I)基于图像区域的配准基于区域的配准方法，有时也称为类相关法或者模板匹配方法，对一定大小的图像或者整幅图像，采用图像的原始灰度信息或者频域信息来度量图像之间的相似性，通过直接搜索变换空间，来估计最优几何变换模型，主要有基于灰度相似性的配准方法、基于傅立叶变换的配准方法等。其中，基于灰度相似性的配准方法，通常选择图像的某一区域或者整幅图像作为一种特殊的图像特征，因此对于此类方法，核心就是匹配标准的度量。在以往的研究中，已经出现很多种区域匹配度量标准，相关系数是一种最常用的度量标准。互信息法是由Voila和Well于1995年首次提出的。互信息法以信息论为基础，从信息熵的角度衡量两个区域的匹配程度。基于傅立叶变换的配准方法，提供了一种从频域角度衡量配准程度的方法。该方法的优点是计算速度快，能有效地去除频域独立的噪声，而这种噪声的影响是难以用空域滤波消除的。Kuglin和Hines提出了相位相关法，利用傅立叶变换将两幅待配准的图像变换到频域，然后利用互功率谱直接计算出两幅图像间的平移矢量。Chen最早将FFT用于SPOT影像的配准。1987年，De Castro和Morandi提出了扩展相位相关法，用傅立叶变换的性质实现具有旋转和平移变换的图像的配准。(2)基于图像特征的配准方法基于特征的配准方法，首先检测图像的显著特征，然后借助与特征的属性进行特征匹配，最优通过特征点对来估计变换模型；或者，直接利用特征之间的空间关系，同时完成特征匹配和变换模型估计，主要有基于特征空间的配准和基于特征描述子的配准方法
坐寸ο基于特征空间的配准方法中，Goshtasby提出了一种基于场景均一性的点匹配算法，将Clustering技术用于图像配准，Hu首先将相似不变矩应用于模式识别；Flusser和Suk导出了仿射不变矩，并成功应用于SPOT和Landsat摄像配准。基于特征描述子的配准方法，使用描述符进行配准。Goshtasby提出一种基于特征轮廓的描述符，该方法巧妙地将二值化的特征图形放到一个网格化的极坐标系中，将图形转换为矩阵进行描述。Bentoutou等人首先沿图像边缘提取Harris角点，然后采用组合不变矩来构造描述子。David Lowe设计的SIFT描述子，在多尺度空间采用DOG算子检测特征点，采用欧式距离进行测量。发明内容
本发明的目的是提供一种多源传感器的遥感图像配准方法，针对多源传感器、多 时相、多视角遥感图像进行高精度、快速配准，以解决遥感图像的自动配准问题。
为实现上述目的，本发明的技术解决方案是
一种多源传感器的遥感图像配准方法，其包括以下步骤
第一步，分别对待配准图像和参考图像进行尺度不变特征变换(SIFT)，提取特征 点，并筛选最优匹配点对；
第二步，通过随机一致性抽样算法(RANSAC)提纯初始配准点，剔除误配准点；
第三步，对图像进行分区，通过不同区域特征点的数目，计算特征点对的分布质量 权重参数，选取分布均匀的有效控制点对；
第四步，通过互信息相似判断准则，在待配准图像的控制点一定范围内搜索最优 配准点，实现控制点对的精配准，获取最优配准点对；
第五步，通过最小二乘法曲线拟和，预测非线性模型二次多项式的参数，即待配准 图像的几何形变参数，实现待配准图像与参考图像的精确配准。
所述的遥感图像配准方法，其所述第一步，包括
a.1、利用二维高斯函数核对参考图像和待配准图像进行尺度空间变换，得到多尺 度下的尺度空间表示序列。
a. 2、特征点检测，提取尺度空间极值，并通过函数拟和实现特征点的精确定位，具 体步骤如下
a. 2.1、在图像二维平面空间和高斯差分尺度空间(Difference-of-Gaussian, DOG)同时检测极值作为特征点，其中DOG算子为两个不同尺度的高斯核的差分。
a. 2. 2、尺度空间极值检测，获取相邻三个尺度空间的局部极值，保证在尺度空间 和二维图像空间都为局部极值；
a. 2. 3、去除低对比度的关键点和不稳定的边缘点，并根据拟和三维二次函数精确 定位特征点的位置和尺度；
a. 3、特征点描述，获取特征点的主方向和描述子信息，具体步骤如下
a. 3.1、提取关键点邻域像素的梯度方向，作为特征点的主方向，确定主方向的模 值和方向；
a. 3. 2、将坐标轴旋转到特征点的主方向上，确保特征描述子的旋转不变性，以特 征点选择一定大小的区域(本方法取16 X 16的区域)，分成16个4X 4的子区域，也称为种 子点，计算每个种子点八个方向上的梯度方向直方图，得到特征点的128维数据，即特征描 述子；
a. 4、计算待配准图像与参考图像中特征点的最邻近欧式距离和次邻近距离，当其 比值小于某个阈值，则被选为最优匹配点对，获得待配准图像与参考图像的初始配准点集。
所述的遥感图像配准方法，其所述第二步，包括
b.1、从初始配准点集中，随机取4个配准点对作为样本,并初始化8自由度参数变 换模型；
b. 2、在初始配准点集中检索符合参数变换模型的误差阈值低于Td( —般取O. 01) 的点集，为与样本一致的有效点集；
b. 3、根据模型误差值Td，获取最大阈值Tmax，当有效点集个数大于Tmax时，重新估计参数变换模型；
b. 4、根据模型误差值Td，获取最小阈值Tmin，当有效点集个数小于Tmin时，重新选择新的样本；
b. 5、经过N次迭代，获取初始匹配点集中的最大一致点集，用来重新估计模型，得到最优匹配点对，剔除误配准点。
所述的遥感图像配准方法，其所述第三步，包括
c.1、将整个图像划分为等4X4个等大小区域，统计每个区域中特征点的个数；
c. 2、每个区域中特征点选中概率为l/(16XNum)，其中Num为该区域特征点数，保证图像中所有特征点的选择概率的均衡；
c. 3、每个区域选择特征点的个数为Nmax/(16XNum),其中Nmax为最大选择特征点数，得到分布均匀的控制点对。
所述的遥感图像配准方法，其所述第四步，包括
d.1、以参考图像、待配准图像的控制点为中心分别选择7X7、21X21的区域作为匹配区域和待匹配区域；
d. 2、计算待匹配区域中7X7大小的图像子块与参考图像的匹配模块的互信息；
d. 3、搜索待匹配区域中与参考模块互相关系数最大的7X7的子模块，并将其中心点作为参考模块中心点的最优匹配点，实现控制点的精配准。
所述的遥感图像配准方法，其所述第五步，包括
e.1、根据精配准点对，经过最小二乘法，实现二次多项式曲线拟和，获得待配准图像相对参考图像的几何变换模型；
e. 2、根据几何变换参数，通过双线性插值法实现待配准图像的配准。
所述的遥感图像配准方法，其所述a.1步中，
二维高斯核的函数G (X，y, σ )定义如下
权利要求
1.一种多源传感器的遥感图像配准方法，其特征在于，包括以下步骤第一步，分别对待配准图像和参考图像进行尺度不变特征变换，提取特征点，并筛选最优匹配点对；第二步，通过随机一致性抽样算法提纯初始配准点，剔除误配准点；第三步，对图像进行分区，通过不同区域特征点的数目，计算特征点对的分布质量权重参数，选取分布均匀的有效控制点对；第四步，通过互信息相似判断准则，在待配准图像的控制点一定范围内搜索最优配准点，实现控制点对的精配准，获取最优配准点对；第五步，通过最小二乘法曲线拟和，预测非线性模型二次多项式的参数，即待配准图像的几何形变参数，实现待配准图像与参考图像的精确配准。
2.如权利要求1所述的遥感图像配准方法，其特征在于，所述第一步，包括a.1、利用二维高斯函数核对参考图像和待配准图像进行尺度空间变换，得到多尺度下的尺度空间表示序列。a. 2、特征点检测，提取尺度空间极值，并通过函数拟和实现特征点的精确定位，具体步骤如下a. 2.1、在图像二维平面空间和高斯差分尺度空间同时检测极值作为特征点，其中DOG算子为两个不同尺度的高斯核的差分。a. 2. 2、尺度空间极值检测，获取相邻三个尺度空间的局部极值，保证在尺度空间和二维图像空间都为局部极值；a. 2. 3、去除低对比度的关键点和不稳定的边缘点，并根据拟和三维二次函数精确定位特征点的位置和尺度；a. 3、特征点描述，获取特征点的主方向和描述子信息，具体步骤如下a. 3.1、提取关键点邻域像素的梯度方向，作为特征点的主方向，确定主方向的模值和方向；a.3. 2、将坐标轴旋转到特征点的主方向上,确保特征描述子的旋转不变性，以特征点选择16X16的区域，分成16个4X4的子区域，也称为种子点，计算每个种子点八个方向上的梯度方向直方图，得到特征点的128维数据，即特征描述子；a.4、计算待配准图像与参考图像中特征点的最邻近欧式距离和次邻近距离，当其比值小于某个阈值，则被选为最优匹配点对，获得待配准图像与参考图像的初始配准点集。
3.如权利要求1所述的遥感图像配准方法，其特征在于，所述第二步，包括b.1、从初始配准点集中，随机取4个配准点对作为样本,并初始化8自由度参数变换模型；b. 2、在初始配准点集中检索符合参数变换模型的误差阈值低于Td的点集，为与样本一致的有效点集；b.3、根据模型误差值Td，获取最大阈值Tmax，当有效点集个数大于Tmax时，重新估计参数变换模型；b.4、根据模型误差值Td，获取最小阈值Tmin，当有效点集个数小于Tmin时，重新选择新的样本；b.5、经过N次迭代，获取初始匹配点集中的最大一致点集，用来重新估计模型，得到最优匹配点对，剔除误配准点。
4.如权利要求1所述的遥感图像配准方法，其特征在于，所述第三步，包括c.1、将整个图像划分为等4X4个等大小区域，统计每个区域中特征点的个数；c.2、每个区域中特征点选中概率为l/(16XNum)，其中Num为该区域特征点数，保证图像中所有特征点的选择概率的均衡；c.3、每个区域选择特征点的个数为Nmax/(16XNum),其中Nmax为最大选择特征点数,得到分布均匀的控制点对。
5.如权利要求1所述的遥感图像配准方法，其特征在于，所述第四步，包括d.1、以参考图像、待配准图像的控制点为中心分别选择7X7、21X21的区域作为匹配区域和待匹配区域；d. 2、计算待匹配区域中7X7大小的图像子块与参考图像的匹配模块的互信息；d.3、搜索待匹配区域中与参考模块互相关系数最大的7X7的子模块，并将其中心点作为参考模块中心点的最优匹配点，实现控制点的精配准。
6.如权利要求1所述的遥感图像配准方法，其特征在于，所述第五步，包括e.1、根据精配准点对，经过最小二乘法，实现二次多项式曲线拟和，获得待配准图像相对参考图像的几何变换模型；e. 2、根据几何变换参数，通过双线性插值法实现待配准图像的配准。
7.如权利要求2所述的遥感图像配准方法，其特征在于，所述a.1步中，二维高斯核的函数G (x，y，σ)定义如下
8.如权利要求2所述的遥感图像配准方法，其特征在于，所述a.2.1步中，DOG尺度空间算子D (x，y，σ)为两个不同尺度的高斯核的差分，其定义如下D(x, y, σ) = (G(x, y, k σ )-G(x, y, σ ))*I (χ, y)=L (x, y, k σ ) -L (χ, y, σ )......(3)其中G(x,y, σ )为二维高斯核,L(x, y, σ )为方差为σ的尺度空间表示,L (x, y, k σ )为方差为ko的尺度空间表示。
9.如权利要求2所述的遥感图像配准方法，其特征在于，所述a.2. 2步中，尺度空间极值检测，获取相邻三个尺度空间的局部极值，为同一尺度空间邻域8个点和相邻两个尺度空间的9X2个点，共26个点中的极值点。
10.如权利要求2所述的遥感图像配准方法，其特征在于，所述a.3.1步中，提取的特征点主方向的模值m(x, y)定义为
11.如权利要求2所述的遥感图像配准方法，其特征在于，所述a.4步中，最优匹配点对准则为
12.如权利要求5所述的遥感图像配准方法，其特征在于，所述d.2步中，互信息是用来衡量两幅图像对应像素灰度值之间的统计相关性，值越大，相关性越强，互信息MI (A，B)定义如下
13.如权利要求5所述的遥感图像配准方法，其特征在于，所述d.3步中，将互信息作为图像匹配的测度，寻找最优匹配点，实现控制点的精配准，就是要寻找搜索空间的最大互信息位置Wtjpt
14.如权利要求6所述的遥感图像配准方法，其特征在于，所述e.1步中，二次多项式几何变换模型定义公式如下χ' = a00+a10x+a01y+a20x2+anxy+a02y2......(10)J1 = b00+b10x+b01y+b20x2+bnxy+b02y2......(11)其中
全文摘要
一种多源传感器的遥感图像配准方法，涉及图像处理技术，包括步骤对参考图像和配准图像分别进行尺度不变特征变换(SIFT)，提取特征点，计算待配准图像与参考图像中特征点的最近欧式距离和次近欧式距离，根据比值筛选最优匹配点对；通过随机一致性抽样算法剔除误配准点，筛选初配准点对；根据特征点权重系数，计算特征点对的分布质量参数，选取分布均匀的有效控制点对；以互信息相似判断准则，在待配准图像控制点搜索最优配准点，获取控制点对最优配准点对；用多项式参数变换，获取待配准图像几何形变参数，实现待配准图像与参考图像的精确配准。本发明有较快运算速度和较高配准精度，能满足多传感器、多时相、多视角遥感图像的配准需求。
文档编号G06T7/00GK103020945SQ201110280990
公开日2013年4月3日 申请日期2011年9月21日 优先权日2011年9月21日
发明者石悦, 付琨, 孙显 申请人:中国科学院电子学研究所