一种无人机航拍图像快速拼接方法

技术领域

本发明涉及图像拼接技术领域，具体是一种无人机航拍图像快速拼接方法。

背景技术

近年来，随着无人机技术的发展，无人机被广泛应用于各种地理信息数据采集，在吹填项目的应用也有很多先例，说明其作业有着明显的优势；但受航空相机传感器尺寸的限制，存在着视场与分辨率间的矛盾。单张无人机航拍图像的覆盖范围较小，不能满足实际应用需求，因此如何快速高效的将多幅无人机航拍图像拼接成一幅过渡自然的无缝图片成为了近年来研究的热点。

现有技术中，直接对无人机获取的航拍图像进行拼接，容易产生无效拼接，为此，现提供一种无人机航拍图像快速拼接方法。

发明内容

本发明的目的在于提供一种无人机航拍图像快速拼接方法。

本发明的目的可以通过以下技术方案实现：一种无人机航拍图像快速拼接方法，包括以下步骤：

步骤S1：将航拍区域平等划分为航拍区块，对航拍区块进行遍历检测，检测航拍区块是否存在相邻的航拍区块，若不存在相邻的航拍区块则无需对航拍区块进行处理，若存在相邻的航拍区块，则随机挑选其中一个航拍区块进行扩充处理，使相邻的两个航拍区块之间存在重合区域；

步骤S2：分别对航拍区块进行航拍，获取航拍图像，生成航拍图像集；每一个航拍区块都有其对应的航拍图像集；

步骤S3：获取航拍图像集中的一个航拍图像，作为主体图像，再将航拍图像集中除主体图像以外航拍图像标记为样本图像，并生成样本图像集；

步骤S4：根据主体图像建立二维坐标系，通过二维坐标系对主体图像进行分割，生成图像块；在图像块中选取主像素块；设定辐射函数，根据辐射函数获取主像素块的辐射像素块；由主像素块和辐射像素块构成判定像素；获取多个判定像素，生成判定像素集；

步骤S5：将主体图像依次跟样本图像集中的样本图像进行比对，获取样本图像中的判定像素个数；并根据样本图像中的判定像素个数判断主体图像是否能与该样本图像进行拼接；若判断主体图像与该样本图像可以进行拼接，则将该样本图像与主体图像进行拼接，生成新的主体图像；再将新的主体图像依次与样本图像集中的样本图像进行比对，若符合条件，则进行拼接，直至样本图像集变为空集，生成区块图像；获取所有航拍区块的区块图像，生成区块图像集；

步骤S6：依照生成区块图像的拼接方法，将区块图像集拼接成区域图像，即航拍区域的整体图像。

进一步的，划分航拍区域的过程包括：

获取航拍区域的长度和宽度，将长度和宽度分别等分为N份，获得N*N个航拍区块；

设置扩充百分比；对航拍区块进行遍历检测，若检测到航拍区块存在相邻的航拍区块，则随机挑选其中一个航拍区块，作为扩充区块，将扩充区块与相邻航拍区块接触边的邻边长度增添扩充长度。

进一步的，生成航拍图像集的过程为：

将航拍区块范围制作成KML文件提前导入飞控软件，飞控软件根据相关飞行参数，自动规划出航线后，无人机执行航拍，获取航拍图像，生成航拍图像集。

进一步的，获取主体图像和样本图像集的过程包括：

获取航拍图像集中的一个航拍图像，作为主体图像，再将航拍图像集中除主体图像以外航拍图像标记为样本图像，并生成样本图像集。

进一步的，生成判定像素集的过程为：

以主体图像左下角为原点，建立二维坐标系；

设置分辨值R，设定R个像素点为一个单位；

设定切割值，根据切割值生成切割函数对二维坐标系进行切割，获取图像块；

设定间隔值，根据间隔值生成间隔集；通过Random函数随机获取间隔集中的数值依次作为x值和y值，生成坐标；

将生成坐标操作循环u次后，获得u个坐标，生成坐标集；

依次将坐标集中的坐标所在图像块标记为主像素块；

根据坐标设定辐射函数，根据辐射函数获取主像素块的辐射像素块，再由主像素块和辐射像素块构成判定像素，生成判定像素集；

将主像素块及其辐射像素块标记为判定像素，生成判定像素集。

进一步的，判断样本图像是否能与主体图像进行拼接的过程包括：

建立样本图像的二维坐标系，对样本图像进行分割，获取样本图像块；

将位于样本图像块中心的坐标标记为中心坐标；

依次将样本图像块作为样本主像素块，再根据中心坐标、辐射函数和范围值获取辐射坐标，并将辐射坐标所在的样本图像块标记为样本辐射像素块；将样本主像素块及其样本辐射像素块标记为像素组，生成像素组集；

设置阈值；分别获取判定像素集中判定像素的像素值和像素组集中像素组的像素值，进行比对，若二者的像素值相差小于或等于阈值，则将该像素组标记为判定像素，若二者的像素值相差大于阈值，则不标记；

设定合格值；获取像素组集中被标记为判定像素个数；

若像素组集中被标记为判定像素个数小于合格值，则判定主体图像与该样本图像不可进行拼接；若像素组集中被标记为判定像素个数大于或等于合格值，则判定主体图像与该样本图像可以进行拼接。

进一步的，生成区块图像的过程为：

通过SURF特征点检测算法提取主体图像和样本图像中的特征点并进行匹配，再通过SDD算法对主体图像和样本图像进行图像配准，并使用RANSAC为图像计算单应性矩阵，使主体图像和样本图像转化为同一坐标下，再通过羽化，使样本图像和主体图像融合为新的主体图像，再将新的主体图像依次与样本图像集中的样本图像进行比对，符合条件后进行拼接，直至样本图像集变为空集，生成区块图像。

进一步的，生成航拍区域的整体图像过程为：

获取各个航拍区块的区块图像，生成区块图像集，依照生成区块图像的方法，将区块图像集中的区块图像拼接成区域图像，即航拍区域的整体图像。

与现有技术相比，本发明的有益效果是：本方法通过对航拍区域的切割，获取多个航拍区块，使多台无人机同时作业，获取航拍图像，极大程度上提高了航拍效率，同时，在划分航拍区块时添加了扩充区域，防止后续图像拼接过程中发生相邻航拍区块间不存在相同特征点的情况；且在航拍图像拼接前，根据航拍图像集生成主体图像和样本图像集，获取主体图像的主像素块和辐射像素块，生成判定像素；将主体图像依次跟样本图像集中的样本图像进行比对，获取样本图像中存在的判定像素个数；并根据样本图像中的判定像素个数判断主体图像是否能与该样本图像进行拼接；当判断主体图像与该样本图像可以进行拼接时，再将该样本图像与主体图像进行拼接，生成新的主体图像；当判断主体图像与该样本图像不可以进行拼接时，则不予拼接操作；如此，在图像拼接前通过判定像素判断两个图像是否存在拼接的可能，提高了图像拼接的效率，避免图像的无效拼接。

附图说明

图1为本发明的流程图。

具体实施方式

如图1所示，一种无人机航拍图像快速拼接方法，包括以下步骤：

步骤S1：将航拍区域平等划分为航拍区块，对航拍区块进行遍历检测，检测航拍区块是否存在相邻的航拍区块，若不存在相邻的航拍区块则无需对航拍区块进行处理，若存在相邻的航拍区块，则随机挑选其中一个航拍区块进行扩充处理，使相邻的两个航拍区块之间存在重合区域；

步骤S2：分别对航拍区块进行航拍，获取航拍图像，生成航拍图像集；每一个航拍区块都有其对应的航拍图像集；

步骤S3：获取航拍图像集中的一个航拍图像，作为主体图像，再将航拍图像集中除主体图像以外航拍图像标记为样本图像，并生成样本图像集；

步骤S4：根据主体图像建立二维坐标系，通过二维坐标系对主体图像进行分割，生成图像块；在图像块中选取主像素块；设定辐射函数，根据辐射函数获取主像素块的辐射像素块；由主像素块和辐射像素块构成判定像素；获取多个判定像素，生成判定像素集；

步骤S5：将主体图像依次跟样本图像集中的样本图像进行比对，获取样本图像中的判定像素个数；并根据样本图像中的判定像素个数判断主体图像是否能与该样本图像进行拼接；若判断主体图像与该样本图像可以进行拼接，则将该样本图像与主体图像进行拼接，生成新的主体图像；再将新的主体图像依次与样本图像集中的样本图像进行比对，若符合条件，则进行拼接，直至样本图像集变为空集，生成区块图像；获取所有航拍区块的区块图像，生成区块图像集；

步骤S6：依照生成区块图像的拼接方法，将区块图像集拼接成区域图像，即航拍区域的整体图像；

需要进一步说明的是，在具体实施过程中，划分航拍区域的过程包括：

获取航拍区域的长度和宽度，将长度和宽度分别等分为N份，获得N*N个航拍区块；

设置扩充百分比K；

对航拍区块进行遍历检测，若检测到航拍区块存在相邻的航拍区块，则随机挑选其中一个航拍区块，作为扩充区块，将扩充区块与相邻航拍区块接触边的邻边长度增添扩充长度，所述扩充长度为领边原本长度*扩充百分比K；

若航拍区块A和航拍区块B相邻；

设定航拍区块A东方向的边为L

设定航拍区块B东方向的边为L

航拍区块A与航拍区块B的接触边为L

L

设定L

若选择航拍区块A作为扩充区块，则将C

若选择航拍区块B作为扩充区块，则将C

通过在原有长度的基础上增添扩充长度，使航拍区块的原有区域增添扩充区域，从而使相邻的航拍区块间存在重叠区域，以便后续无人机获取航拍区域图像的拼接；

需要进一步说明的是，在具体实施过程中，生成航拍图像集的过程为：

将航拍区块范围制作成KML文件提前导入飞控软件，同时按照地面分辨率GSD＝2.74cm/pix进行技术设计，影像数据满足1：1000比例尺的线画图、正射影像图的成图精度要求，航线按照平行于测区长边一侧布设，航摄相对航高100米，航向重叠度80％、旁向重叠度70％，保障相片有足够的重叠度；飞控软件根据相关飞行参数，自动规划出航线，飞行前再检查无人机设备正常，起飞区域安全后，无人机执行航拍，获取航拍图像p，生成航拍图像集P。

需要进一步说明的是，在具体实施过程中，获取主体图像和样本图像集的过程为：

获取航拍图像集中的一个航拍图像，作为主体图像，再将航拍图像集中除主体图像以外航拍图像标记为样本图像，并生成样本图像集；

设定航拍图像集P＝{p

通过Random函数获取1至n之间的随机数i，设定i＝1，则p

需要进一步说明的是，在具体实施过程中，生成判定像素集的过程为：

以主体图像左下角为原点，建立二维坐标系；

设置分辨值R，设定R个像素点为一个单位，即一个单位平方里存在R*R个像素点；

设定切割值Q，根据切割值Q生成切割函数x＝mQ和y＝mQ对二维坐标系进行切割，获取图像块，其中，m＝1，2，……，n；

设定间隔值j，根据间隔值生成间隔集J＝{num

通过Random函数随机获取间隔集中的数值依次作为x值和y值，生成坐标(x，y)；

将生成坐标操作循环u次后，获得u个坐标，生成坐标集D＝{d

依次将坐标集D中的坐标d

根据坐标设定辐射函数，根据辐射函数获取主像素块的辐射像素块，再由主像素块和辐射像素块构成判定像素；

当坐标为d

设定范围值f，分别获取x＝x

若发生多个辐射坐标位于同一个图像块中，则重新设定范围值f，直至辐射坐标分别位于不同的图像块中；

将主像素块及其辐射像素块标记为判定像素，生成判定像素集；

需要进一步说明的是，在具体实施过程中，判断样本图像是否能与主体图像进行拼接的过程为：

依照对主体图像建立二维坐标系和生成图像块的原理，对样本图像进行相同操作，建立样本图像的二维坐标系，对样本图像进行分割，获取样本图像块；

将位于样本图像块中心的坐标标记为中心坐标，作为代表该样本图像块的坐标；

依次将样本图像块作为样本主像素块，再根据中心坐标、辐射函数和范围值获取辐射坐标，并将辐射坐标所在的样本图像块标记为样本辐射像素块；

将样本主像素块及其样本辐射像素块标记为像素组，生成像素组集；

设置阈值z；

分别获取判定像素集中判定像素的像素值和像素组集中像素组的像素值，进行比对，若二者的像素值相差小于或等于z，则将该像素组标记为判定像素，若二者的像素值相差大于z，则不标记；

设定合格值H；

获取像素组集中被标记为判定像素个数h；

若h小于H，则判定主体图像与该样本图像不可进行拼接；

若h大于或等于H，则判定主体图像与该样本图像可以进行拼接；

需要进一步说明的是，在具体实施过程中，样本图像和主体图像进行拼接，生成区块图像的过程为：

当判定主体图像和样本图像可以进行拼接后，通过SURF特征点检测算法提取主体图像和样本图像中的特征点并进行匹配，再通过SDD算法对主体图像和样本图像进行图像配准，并使用RANSAC为图像计算单应性矩阵，使主体图像和样本图像转化为同一坐标下，再通过羽化，使样本图像和主体图像融合为新的主体图像，再将新的主体图像依次与样本图像集中的样本图像进行比对，符合条件后进行拼接，直至样本图像集变为空集，生成区块图像；

需要进一步说明的是，在具体实施过程中，生成航拍区域的整体图像过程为：

获取各个航拍区块的区块图像，生成区块图像集，依照生成区块图像的方法，将区块图像集中的区块图像拼接成区域图像，即航拍区域的整体图像；

在区块图像集中选取一张区块图像作为主体区块图像，剩余的区块图像则标记为样本区块图像，生成样本区块图像集，将主体区块图像与样本区块图像集中的样本区块图像进行一一比对，符合条件后，再进行拼接。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。