[发明专利]一种基于卫星影像的无人机影像快速概略拼接方法

权利要求书

1.一种基于卫星影像的无人机影像快速概略拼接方法，其特征在于，包括如下步骤：

步骤1、使用无人机机载数码相机对地面进行近似垂直摄影，获取航摄区域的地面场景的无人机影像；通过网络获取航摄区域的卫星影像，用作无人机影像概略拼接的底图；通过获取所述相机像元尺寸、镜头焦距以及飞行航高信息，计算得到无人机影像的航摄区域的地面分辨率；

步骤2、影像预处理，包括影像降采样和双边滤波：用卫星影像地面分辨率与无人机影像地面分辨率的比值，作为无人机影像的降采样比例，对无人机影像进行降采样，使其与卫星影像具有近似相同的地面分辨率；并且，用滤波器对无人机影像和卫星影像进行双边滤波；

步骤3、利用降尺度空间的SURF算法对无人机影像和卫星影像进行特征点的提取：由于无人机影像和卫星影像具有近似相同的地面分辨率，所以仅需选用少量滤波器模板参与特征点提取的计算；

步骤4、进行欧氏距离匹配，并用模板匹配辅助优选匹配点：利用欧氏距离对无人机影像和卫星影像的两幅图像中的特征点进行匹配，利用特征点主方向对卫星影像进行旋转，消除影像间的旋角；然后在此基础上，利用基于标准平方差的模板匹配方法加以辅助，剔除误差较大的匹配点，计算得到一对最优匹配特征点；

步骤5、在所述的SURF算法中，提取的特征点具有像素坐标、尺度和方向信息，利用所述的一对最优匹配特征点的像素坐标、尺度和方向信息，构建相似变换矩阵R，其矩阵计算公式如下：

其中，s表示两特征点尺度比值，θ表示两特征点主方向之间旋角，d_x和d_y分别表示两特征点在x方向和y方向的像素坐标之间的差；

以无人机影像的左上角为坐标原点，计算出该影像的四个顶点的像素坐标，通过相似变换矩阵求得所述四个顶点在卫星影像中的像素坐标(x′,y′)，得到无人机影像的概略位置：

其中，R为相似变换矩阵，(x,y)为无人机影像中四个顶点的像素坐标；

步骤6、将航摄区域内的所有无人机影像对卫星影像进行配准，即可得到无人机影像之间的方位信息，从而完成无人机影像的概略拼接。

2.根据权利要求1所述的一种基于卫星影像的无人机影像快速概略拼接方法，其特征在于：在所述步骤1中，利用以下计算公式获取无人机影像的地面分辨率：

其中，GSD为地面分辨率，单位是m，a为相机像元尺寸，单位是mm，H为摄影航高，单位是m，f为摄影镜头焦距，单位是mm。

3.根据权利要求1所述的一种基于卫星影像的无人机影像快速概略拼接方法，其特征在于：

在所述步骤2中，所述的影像降采样，是指用卫星影像地面分辨率与无人机影像地面分辨率的比值，作为无人机影像的降采样比例，对无人机影像进行降采样，从而使其无人机影像与卫星影像具有近似相同的地面分辨率；

所述的双边滤波，是指对无人机影像和卫星影像进行双边滤波；双边滤波器不仅考虑像素点邻域的灰度相似关系，还需考虑像素点空间的邻近关系。

4.根据权利要求1所述的一种基于卫星影像的无人机影像快速概略拼接方法，其特征在于：

在所述步骤3中，所述的利用降尺度空间的SURF算法对无人机和卫星影像进行特征点的提取，是指在使用SURF算法进行特征点检测的过程中，为了保持特征点的尺度不变性，需要构建多组、每组多层尺度空间分别提取特征点，即使引入盒状滤波器仍然需要较大的计算量，所以仅选用少量滤波器模板构建尺度空间，进而提取特征点。

5.根据权利要求1所述的一种基于卫星影像的无人机影像快速概略拼接方法，其特征在于：所述步骤4的具体过程如下：

(4-1)所述的利用欧氏距离对无人机影像和卫星影像的两幅图像中的特征点进行匹配，其欧氏距离的计算公式如下：

其中，d为欧氏距离，m表示无人机影像的一幅图像，n表示卫星影像的一幅图像，F_m和F_n分别表示m和n两幅图像上的某个特征点的描述子，k表示特征点描述子所具有的维度，f^(m)_i和f⁽ⁿ⁾_i分别表示m和n两幅图像上的第i维向量的长度；

(4-2)如果最近距离比次近距离小于阈值，其中缺省值取0.7，则认为最近距离的两个特征点为一对匹配点，然后将所述的最近距离的两个特征点及其最近距离与次近距离的比进行记录，并按照所述的距离比从小到大的原则，对匹配点进行排序；

(4-3)从上述排序的匹配点序列中按顺序依次取匹配点，计算两匹配点主方向之间的旋角，对卫星影像进行旋转，消除两张影像之间的旋角；取特征点周围N×N像素的邻域，其中N取特征点尺度的5倍；采用标准平方差法进行模板匹配，并记录下平方差系数，该平方差系数越小，说明图像的相似度越高；标准平方差系数R的计算公式如下：

其中，T表示模板图像，I表示原始图像，x′、y′表示模板图像中的坐标，x、y表示模板图像的移动步长；重复上述(4-3)过程，从匹配点序列中选取标准平方差最小的匹配点作为最优匹配点。