1.一种卫星遥感图像配准方法,其特征在于,包括以下步骤:步骤一、拍摄当前帧卫星遥感图像,进行图像二值化处理,获得第一幅图像;
步骤二、对第一幅图像采用canny算子提取边缘,获得第二幅图像;
步骤三、对第二幅图像的每个闭合连通区域分别提取质心;
步骤四、分别以第二幅图像的每个闭合连通区域的质心为中心,提取出对应的闭合连通区域,比较所提取的闭合连通区域的面积大小,选择面积最大的闭合连通区域作为第三幅图像,选择面积最小的闭合连通区域作为第四幅图像;
步骤五、在第三幅图像中,以其质心为圆心作最大直径R1的外接圆,在外接圆内部作最大面积正方形;在第四幅图像中,以其质心为圆心作最大直径R2的外接圆,在外接圆内部作最大面积正方形;
步骤六、对第三幅图像和第四幅图像分别进行Hough变换,分别记录第三幅图像的旋转角度θ1和第四幅图像的旋转角度θ2;
步骤七、以第三幅图像的质心为基点将第三幅图像的质心位置映射到第二幅图像中,以第四幅图像的质心为基点将第四幅图像的质心位置映射到第二幅图像中;
步骤八、再对第二幅图像每个闭合连通区域分别提取质心,分别以每个质心为中心,排列组合所有正方形,将所有正方形按照面积大小进行排列,记录每个正方形中心位置、每个正方形边长大小、每个正方形所对应的外接圆面积大小、每个正方形所对应的闭合连通区域位置,最终记录该正方形中心在第二幅图像中的位置,收集上述数据得到当前帧卫星遥感图像数据库;
步骤九、根据拍摄的当前帧卫星遥感图像,给出对应的当前帧卫星遥感标准图像;对当前帧卫星遥感标准图像重复步骤二至步骤八,获得当前帧卫星遥感标准图像数据库;
步骤十、设定正方形最大边长L1、正方形所对应的外接圆最大面积S1,在当前帧卫星遥感标准图像数据库中按照已经排列好的所有正方形搜索同时满足正方形最大边长L1、正方形所对应的外接圆最大面积S1的正方形,当搜索到该正方形时,则将这两个正方形以四个顶点为基点进行匹配,计算出所对应的两幅图像的旋转角度θ3;当|θ1‑θ3|>θ4时停止搜索,其中θ4为旋转角度最大值,搜索成功,否则搜索失败;
步骤十一、设定正方形最小边长L2、正方形所对应的外接圆最小面积S2,在当前帧卫星遥感标准图像数据库中按照已经排列好的所有正方形搜索同时满足正方形最小边长L2、正方形所对应的外接圆最小面积S2的正方形,当搜索到该正方形时,则将这两个正方形以四个顶点为基点进行匹配,计算出所对应的两幅图像的旋转角度θ4;当|θ1‑θ4|<θ5时停止搜索,其中θ5为旋转角度最小值,搜索成功,否则搜索失败;
步骤十一、当搜索成功时,将当前帧卫星遥感图像与当前帧卫星遥感标准图像按照旋转角度θ3和旋转角度θ4对齐;当搜索失败时,将当前帧卫星遥感图像与当前帧卫星遥感标准图像按照旋转角度θ1和旋转角度θ2对齐;
步骤十二、利用归一化互相关方法计算对齐后的当前帧卫星遥感图像与当前帧卫星遥感标准图像的平移量,同时将这两幅图像完全对齐。
2.根据权利要求1所述的一种卫星遥感图像配准方法,其特征在于,步骤二中,对第一幅图像采用canny算子提取边缘,获得第二幅图像的具体操作步骤如下:(1)滤波:采用高斯滤波函数对原始第一幅图像进行滤波及降噪处理,二维高斯函数为:
式(1)中,G(x,y,σ)为二维高斯函数,π和e均为常数,σ为高斯函数的分布参数,(x,y)为图像坐标点,对图像进行平滑处理后得到的平滑图像为:g(x,y)=G(x,y,σ)*f(x,y) (2)式(2)中,f(x,y)为原始第一幅图像;
(2)计算梯度幅值和梯度方向
对平滑图像的像素点在水平方向和竖直方向进行求解偏导数:式中,Gx为平滑图像的像素点在水平方向的偏导数,Gy为平滑图像的像素点在水平方向的偏导数;
利用一阶有限差分来计算梯度幅值和梯度方向:式中,M(x,y)为梯度幅值,θ(x,y)为梯度方向;
(3)对梯度幅值进行非极大值抑制;
(4)采用双阈值算法检测边缘和连接边缘:选取高阈值Th和低阈值Tl,然后开始扫描图像f(x,y)对候选边缘图像中标记的候选边缘点的每一个像素点(i,j)进行,若像素点(i,j)的梯度幅值M(x,y)大于Th,则认为该点为边缘点,若像素点(i,j)的梯度幅值M(x,y)小于Tl,则认为该点不是边缘点;当梯度幅值M(x,y)处于Th与Tl之间时,则需要根据边缘连通性进行判断:若该点的相邻像素点中有边缘点,则认为该点为边缘点进行连接,否则该点为非边缘点。