[发明专利]一种由粗到精的星载合成孔径雷达图像自动配准方法

说明书

技术领域

本发明涉及遥感图像处理技术领域，是一种星载合成孔径雷达(SAR)图像的配准方法。

背景技术

星载SAR遥感数据具有全天候、全天时的独特优势，同时具有固定的重访周期，所以它比光学遥感更适合于地表监视和动态监测，随着我国合成孔径雷达(SAR)的发展，应用多时相的SAR数据进行地表和地物变化的动态监视成为十分重要的遥感应用领域。而对不同时间获取的同一地区的SAR图像进行精确配准是变化检测和动态监视的前提和基础。

目前国内外星载SAR图像的配准方法主要有两类：(1)基于区域的方法，将图像上的小窗口内像素与另一幅图像上同样尺寸窗口做统计比较。通常选用合适的相关测度，然后将满足条件的窗口中心作为控制点用于求解两幅图像之间的变换参数。这类方法匹配准确率高，可以人为控制匹配点的分布，但对图像的尺度和角度变化敏感，在低信噪比的情况下难以满足SAR影像的匹配要求；(2)基于特征的方法，主要是匹配两幅图像中相对应的特征，比如边缘特征、点特征、互信息和hausdorff距离等。这类方法自动化程度高，一定程度上对旋转和明暗变化等不敏感，但在实际应用中存在着误匹配点多，匹配点聚集，以及速度慢等缺点。

发明内容

本发明的目的是提供一种由粗到精的星载合成孔径雷达图像自动配准方法，对不同时间获取的星载合成孔径雷达(SAR)图像之间自动进行配准，以解决星载SAR图像的配准问题，提高SAR图像配准的效率和准确度。

为达到上述目的，本发明的技术方案是：

一种由粗到精实现星载合成孔径雷达图像自动配准方法，其包括：

步骤1：根据卫星的轨道数据分别计算两幅星载SAR图像四个角点的位置，并根据这些角点位置计算两幅SAR图像之间的仿射变换参数；

步骤2：对第一幅星载SAR图像进行均匀分块，并在每个图像块内提取一个正方形的特征区域；

步骤3：根据第一幅星载SAR图像的每个特征区域中心的地理位置，在第二幅SAR图像内寻找对应位置点，并以该点位置为中心提取2公里正方形的窗口区域，再按照仿射变换参数计算正方形区域在第二幅SAR图像的对应窗口区域的四点位置；

步骤4：利用得到的第二幅SAR图像上区域的四个角点位置，按照第一幅SAR图像的分辨率对第二幅SAR图像对应四个角点位置的图像采用三次卷积内插方法重采样，得到一幅临时图像区域；

步骤5：利用归一化互相关模板匹配方法对第一幅SAR图像的正方形特征区域和第二幅SAR图像的临时图像区域进行精确匹配，得到精确匹配的相关位置；

步骤6：对第一幅SAR图像的每个正方形特征区域，按照步骤3到步骤5的处理方法进行同样处理，就得到一系列精确匹配的位置点；

步骤7：利用这些精确匹配的系列位置点对第二幅SAR图像进行校正，实现两幅图像的精确配准。

所述的方法，其所述步骤1的仿射变换模型为：

x′＝a+bx+cy

y′＝d+ex+fy

其中，x、y分别表示第一幅SAR图像的行坐标、列坐标，x’、y’分别表示第二幅SAR图像的行坐标、列坐标，a、b、c、d、e、f为仿射变换的六个参数，即a表示x方向的平移量，b表示x方向的比例缩放因子，c表示y方向的比例缩放因子，d表示y方向的平移量，e表示x方向的比例缩放因子，f表示y方向的比例缩放因子。

所述的方法，其所述步骤3的处理方法是，根据第一幅SAR图像特征区域的中心位置，考虑卫星轨道的精度，以推算出特征区域对应于第一幅SAR图像中对应的区域窗口四点位置；由于第二幅SAR图像的四个角点位置已知，根据仿射变换，就可以计算出第二幅SAR图像上对应特征窗口区域的四点位置。

所述的方法，其所述步骤5中的归一化互相关模板匹配方法，是以下式表示：