发明名称 一种基于卫星影像的无人机影像快速概略拼接方法
摘要 本发明的一种基于卫星影像的无人机影像快速概略拼接方法:使用无人机机载数码相机对地面进行近似垂直摄影,计算得到无人机影像的航摄区域的地面分辨率;对无人机影像进行降采样,使其与卫星影像具有近似的地面分辨率;对无人机和卫星影像进行双边滤波,以去除弱特征点;用降尺度空间SURF算法对无人机和卫星影像进行特征点提取;用欧氏距离进行特征点的匹配,并用模板匹配辅助优选一对正确的匹配点;利用匹配点的尺度和方向信息构建相似变换矩阵,以确定无人机影像在卫星影像中的概略位置;将多幅无人机影像配准到卫星影像中,得到无人机影像间的方位信息,完成无人机影像概略拼接。本发明无需连续的无人机影像,其计算效率高、速度快、适用广。
申请公布号 CN105550994A 申请公布日期 2016.05.04
申请号 CN201610051875.3 申请日期 2016.01.26
申请人 河海大学 发明人 李浩;孙彦玮;王莎;李珊;杨彪
分类号 G06T3/40(2006.01)I 主分类号 G06T3/40(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 一种基于卫星影像的无人机影像快速概略拼接方法,其特征在于,包括如下步骤:步骤1、使用无人机机载数码相机对地面进行近似垂直摄影,获取航摄区域的地面场景的无人机影像;通过网络获取航摄区域的卫星影像,用作无人机影像概略拼接的底图;通过获取所述相机像元尺寸、镜头焦距以及飞行航高信息,计算得到无人机影像的航摄区域的地面分辨率;步骤2、影像预处理,包括影像降采样和双边滤波:用卫星影像地面分辨率与无人机影像地面分辨率的比值,作为无人机影像的降采样比例,对无人机影像进行降采样,使其与卫星影像具有近似相同的地面分辨率;并且,用滤波器对无人机影像和卫星影像进行双边滤波;步骤3、利用降尺度空间的SURF算法对无人机影像和卫星影像进行特征点的提取:由于无人机影像和卫星影像具有近似相同的地面分辨率,所以仅需选用少量滤波器模板参与特征点提取的计算;步骤4、进行欧氏距离匹配,并用模板匹配辅助优选匹配点:利用欧氏距离对无人机影像和卫星影像的两幅图像中的特征点进行匹配,利用特征点主方向对卫星影像进行旋转,消除影像间的旋角;然后在此基础上,利用基于标准平方差的模板匹配方法加以辅助,剔除误差较大的匹配点,计算得到一对最优匹配特征点;步骤5、在所述的SURF算法中,提取的特征点具有坐标、尺度和方向等信息,利用所述的一对最优匹配特征点的像片坐标、尺度和方向信息,构建相似变换矩阵R,其矩阵计算公式如下:<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>s</mi><mo>&CenterDot;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mo>-</mo><mi>s</mi><mo>&CenterDot;</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></mtd><mtd><msub><mi>d</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>&CenterDot;</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>s</mi><mo>&CenterDot;</mo><mi>cos</mi><mi>&theta;</mi></mrow></mtd><mtd><msub><mi>d</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000914826550000011.GIF" wi="596" he="164" /></maths>其中,s表示两特征点尺度比值,θ表示两特征点主方向之间旋角,d<sub>x</sub>和d<sub>y</sub>分别表示两特征点在x方向和y方向的像片坐标之间的差;以无人机影像的左上角为坐标原点,计算出该影像的四个顶点坐标,通过相似变换矩阵求得所述四个顶点坐标在卫星影像中的坐标(x′,y′),得到无人机影像的概略位置:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "(" close = ")"><mtable><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>y</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mi>R</mi><mo>&CenterDot;</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>x</mi></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000914826550000012.GIF" wi="293" he="151" /></maths>其中,R为相似变换矩阵,(x,y)为无人机影像中的四个顶点坐标;步骤6、将航摄区域内的所有无人机影像对卫星影像进行配准,即可得到无人机影像之间的方位信息,从而完成无人机影像的概略拼接。
地址 211100 江苏省南京市江宁开发区佛城西路8号