1.一种基于飞行数据的无人机非量测型相机检校方法,其特征在于,包含如下步骤:
S1:获取无人机正常作业的飞行数据,其中,飞行数据包括影像和机载GPS数据;
S2:采用SURF特征提取算法计算每张影像的特征点及特征描述子;
S3:根据机载GPS数据和每张影像的特征描述子,对特征点进行粗匹配;
S4:对粗匹配后的两张影像的特征点进行精匹配;
S5:对精匹配后的特征点对应的两张影像相对定向,计算强连通图,并得到候选检校影像和验证影像;
S6:采用自检校SFM算法,进行检校参数,验证是否通过,是则结束,否则返回步骤S1。
2.根据权利要求1所述的无人机非量测型相机检校方法,其特征在于:步骤S1中,每张影像对应一个机载GPS数据,并且,机载GPS数据为摄影时无人机非量测型相机在空中的位置。
3.根据权利要求1所述的无人机非量测型相机检校方法,其特征在于:
步骤S2中:
将影像从RGB色彩空间转换到全色色彩空间,然后利用SURF特征提取算法计算每张影像的特征点及特征描述子;
步骤S3中:
S31:根据机载GPS数据构建Delaunay三角网,将此Delaunay三角网转换为无向图得到影像之间的拓扑关系,再利用拓扑排序算法计算任意两张影像间的图距离;
S32:针对图距离小于4的任意两张影像,计算其两两特征描述子之间的距离,将满足匹配距离准则的两个特征点称为一组同名特征点对;
步骤S4中:
当两张影像之间的同名特征点对的数量超过16,则利用RANSAC算法计算这两张影像的基础矩阵,再利用基础矩阵计算同名特征点对的极线距离,当极线距离大于3,则删除该同名特征点对;
步骤S5中:
S51:删除后,当两张影像之间的同名特征点对的数量超过16,则对这两张影像进行相对定向,每一组同名特征点对对应两条光束和一个交会角,当交会角满足最佳交会条件时,将此交会角对应的同名特征点对称为有效同名特征点对,有效同名特征点对的数量即为这两张影像的连接强度;
S52:删除连接强度小于128的影像间的连通关系,设置搜索的最大距离为3,利用广度优先搜索算法,查找连接强度值最大和次大的连通图,将强度值最大连通图作为候选检校影像,并且,将强度值次大连通图作为验证影像;
步骤S6中:
利用候选检校影像间的同名特征点对,采用自检校的SfM方法,计算出候选检校影像的焦距、像主点和畸变参数,即获取该影像的无人机非量测相机的检校参数,再利用验证影像检查检校参数的合理性。
4.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S2中,每张影像特征点和特征描述子的计算方法进一步包括以下步骤:
①RGB色彩空间转换到全色色彩空间
I=0.2989×R+0.5870×G+0.1140×B;
式中,R为红色波段的强度,G为绿色波段的强度,B为蓝色波段的强度,I为全色色彩空间的强度;
②对经过步骤1转换的影像,利用SURF特征提取算法计算其特征点及特征描述子,每个特征点为一个二维向量,每个特征描述子为一个64维的向量。
5.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S31中,任意两张影像间的图距离的计算方法包括以下步骤:
①飞行数据中的每张影像对应一个机载GPS数据(X,Y,Z),利用机载GPS数据的平面坐标(X,Y)构建Delaunay三角网,Delaunay三角网的每个顶点对应一张影像,此Delaunay三角网可以认为是带权无向图,Delaunay三角网中的边对应带权无向图的无向边且权为1;
②Delauany三角网中任意两个顶点的图距离定义为其在带权无向图中的最短路径的权的和,针对每个Delaunay三角网中的每个顶点,利用拓扑排序算法计算出该顶点与其他顶点的图距离。
6.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S32中,匹配距离准则包括以下步骤:
①令待判定的两张影像分别为A和B,特征点分别为{f0 … fi … fn}和{g0 … gj … gm},特征描述子分别为{a0 … ai … an}和{b0 … bj … bm},其中fi和gj为2维向量,ai和bj为64维向量,fi与ai一一对应,gj与bj一一对应;
②任取A中的一个特征描述子ai,计算其与{b0 … bj … bm}的所有特征描述子的距离dk,令距离最小的两个特征描述子为bu和bv,对应的距离为du和dv(du≤dv),如果du小于0.7且du/dv<0.8,则称gu为fi的最优匹配特征点;
dk=|ai-bk|
③将A和B互换,重复步骤②,如果fi与gj互为最优匹配特征点,则称fi与gj为同名特征点对。
7.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S51中,连接强度的计算方法进一步包括以下步骤:
①如果两个影像的同名特征点对的数量大于16,则进行相对定向,此时可获得两张影像对应的相对姿态信息;
②根据相对姿态信息和同名特征点对进行前方交会,得到对应的两条光束和一个交会角;
③如果两条光束的交会角大于20°且小于90度,则称其满足最佳交会条件,这两条光束对应的同名特征点对称为有效同名特征点对;
④统计任意两张影像上的有效同名特征点对的数量n,则此影像间的连接强度y(n)如式;
8.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S52中,查找候选检校影像进一步还包括以下步骤:
①删除步骤S3生成的Delaunay三角网里连接强度小于128的影像间的边,确保弱连接的影像不会参与到后续的检校计算里;
②设置搜索的最大距离为3,以每张影像为搜索起点,采用广度优先搜索算法进行区域增长搜索;
③统计每次搜索到的影像集的连接强度值,将最大连接强度值对应的影像集作为候选检校影像,将次大连接强度值对应的影像集作为验证影像。
9.根据权利要求3所述的无人机非量测型相机检校方法,其特征在于:步骤S6中,自检校SFM算法还包括以下步骤:
①设置畸变模型公式
Δx=(x-x0)(k1r2+k2r4+k3r6)+
p1[r2+2(x-x0)2]+2p2(x-x0)(y-y0)+α(x-x0)+β(y-y0)
Δy=(y-y0)(k1r2+k2r4+k3r6)+p2[r2+2(y-y0)2]+2p1(x-x0)(y-y0)
其中:Δx,Δy为像点改正值;x,y为像方坐标系下的像点坐标;x0,y0为像主点;
②候选检校影像的每个特征点构成一个共线条件方程组,将所有共线条件方程组合并为一个大的待求解方程组,采用Levenberg-Marquardt算法求得待求解方程组的最优解;
③验证影像的所有特征点形成一个共线条件方程组,将求解得到的最优解带入此共线条件方程组中,如果中误差小于1个像素,则接收此最优解,否则将候选检校影像丢弃,转步骤7重新选取候选检校影像和验证影像;
④上一步中得到的最优解即为获取候选检校影像的无人机非量测相机的检校参数。