CN110160503A - 一种顾及高程的无人机景观匹配定位方法 - Google Patents
一种顾及高程的无人机景观匹配定位方法 Download PDFInfo
- Publication number
- CN110160503A CN110160503A CN201910270586.6A CN201910270586A CN110160503A CN 110160503 A CN110160503 A CN 110160503A CN 201910270586 A CN201910270586 A CN 201910270586A CN 110160503 A CN110160503 A CN 110160503A
- Authority
- CN
- China
- Prior art keywords
- triangle
- matching
- unmanned plane
- pair
- elevation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/30—Interpretation of pictures by triangulation
- G01C11/34—Aerial triangulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Multimedia (AREA)
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Abstract
本发明提供一种顾及高程的无人机景观匹配定位方法,包括分别构建基准图像和无人机拍摄的实时图像中的区域空间关系,得到相应的三角形集合,所述基准图像经配准,图中每一个点都具有三维坐标;根据相似三角形的判定定理,基于基准图像的三角形集合和实时图像的三角形集合,进行相似空间关系的粗匹配,然后根据面积比一致性和旋转角度一致性去除伪匹配三角形对,得到正确匹配的三角形对;基于正确匹配的三角形对,获取无人机拍摄的实时图像中特征点的像素坐标以及对应的平面与高程坐标,从而解算无人机的三维坐标,实现顾及高程的三维定位导航。本发明充分利用了图像中存在的空间位置关系进行匹配,顾及了无人机的高程信息,得到三维定位坐标。
Description
技术领域
本发明属于摄影测量研究技术领域,具体涉及无人机的定位方法。
背景技术
近年来,无人机由于造价低廉、起降方便、操作灵活、机动性强等优点,无人机影像具有分辨率高,纹理清晰丰富,且重叠度高,大比例尺小面积等优点,因而被广泛应用于各类测图及非测图领域,如自然资源管理与检测、土地利用动态监测、道路监测、林业草场检测、交通、军事等领域。
无人机飞行器的导航定位数据传统上主要依赖惯性导航或卫星制导设备获取。惯性制导可以实现不依赖外界信息的完全自主导航,隐蔽性好,不受外界干扰和地形影响,可以全天候工作,但定位误差也会随时间而累积,影响定位精度。卫星制导具有全球性、全天候、连续精密导航与定位能力,实时性较为出色,但抗干扰性差(电磁干扰、建筑物遮挡等)。
景象匹配导航由于测量精度高、自成体系、抗干扰性强、系统能耗小、体积小等突出特点,成为了飞行器自主定位系统的研究热点。但是现有的景象匹配定位仅关注二维坐标,无法实现无人机平面坐标和高程坐标的实时定位解算。
发明内容
本发明的目的是提供一种基于景观匹配的无人机辅助三维定位方法,该方法可以实现无人机顾及高程的三维定位导航。可以降低无人机对全球导航卫星系统(GlobalNavigation Satellite System,GNSS)的依赖,实现无人机GNSS信号受干扰或失效时的定位导航。
为达到上述目的,本发明的技术方案提供一种顾及高程的无人机景观匹配定位方法,包括以下步骤:
步骤1、分别构建基准图像和无人机拍摄的实时图像中的区域空间关系,得到相应的三角形集合,所述基准图像经与相应范围内的DEM数据进行配准,图中每一个点都具有三维坐标;
步骤2,基于空间关系相似性的匹配,包括根据相似三角形的判定定理,基于从步骤1所得基准图像的三角形集合和实时图像的三角形集合,进行相似空间关系的粗匹配,初步得到匹配的三角形对,然后根据面积比一致性和旋转角度一致性去除伪匹配三角形对,得到正确匹配的三角形对;
步骤3,解算无人机位置,包括基于步骤2所得正确匹配的三角形对,获取无人机拍摄的实时图像中特征点的像素坐标以及对应的平面与高程坐标,从而解算无人机的三维坐标,实现顾及高程的三维定位导航。
而且,步骤1中,提取区域空间关系的实现方式如下,
获取图像中地物最小外接矩形的中心点,构成特征点集,在特征点集中任选3个点构成三角形集合,并剔除2种类型的三角形,一种为等腰三角形或等边三角形,一种为最小角小于相应预设阈值的三角形。
而且,步骤2中,根据相似三角形的判定定理,进行相似空间关系的粗匹配,实现方式如下,
设步骤1所得基准图像的三角形集合GT有a个三角形,步骤1所得实时图像GS的三角形集合有b个三角形,计算相似性矩阵MTS大小为a×b,MTS(i,j)为GT中第i个三角形与GS中第j个三角形的最大角之差与最小角之差的和;根据预设的角度阈值α,若MTS(i,j)<α,则认为GT中第i个三角形与Gs中第j个三角形相似,记为三角形对所有匹配的三角形对构成集合
而且,步骤2中,根据面积比一致性去除伪匹配三角形对,实现方式如下,
计算集合Gmatching中每一组三角形对的面积比,取λ为三角形对面积比值的均值,根据预设的比例阈值λ0,认为面积比值在 (λ-λ0,λ+λ0)范围内的三角形对为正确匹配三角形对,进行剔除后得到新的匹配三角形对的集合G′matching。
而且,步骤2中,λ0建议取值范围为0.1≤λ0≤0.15。
而且,步骤2中,根据旋转角度一致性去除伪匹配三角形对,实现方式如下,
计算G′matching中每一组三角形对的最大角方向旋转角度和最小角方向旋转角度,取θmax和θmin为最大角方向旋转角度和最小角方向旋转角度的均值,根据预设的角度阈值θ0和θ′0,当一对匹配的三角形对的最大角方向旋转角度θ在(θmax-θ0,θmax+θ0)范围内,且最小角方向旋转角度θ′在(θmin-θ′0,θmin+θ′0)范围内,认为该三角形对是正确匹配的三角形对。
而且,角度阈值θ0和θ′0建议取值范围为3°~5°。
和现有技术相比,本发明具有如下优点和有益效果:
(1)本发明充分利用了图像中存在的空间位置关系,以及相应的数学理论方法,准确地找到正确匹配的三角形对;
(2)本发明在解算无人机二维平面位置的同时,顾及了无人机的高程信息,最终的解算结果是三维(包含平面坐标以及高程坐标)定位坐标。
附图说明
图1是本发明实施例的流程示意图;
图2是本发明实施例最小角方向和最大角方向的示意图。
具体实施方式
以下结合附图和实施例详细说明本发明的技术方案。
本发明实施例基于景观匹配实现无人机顾及高程的定位。首先,制备包含高程数据的底图数据,然后与底图数据进行景观匹配得到实时图像中的点所对应的三维坐标,最后通过空间后方交会实现无人机三维位置的解算。实施例的流程如图1所示,包括以下步骤:
1、分别构建基准图像和实时图像中的区域空间关系。
具体实施时,可以预先设置基准图像作为底图,例如相同区域的高分辨率遥感影像,或者是无人机影像的正射拼接影像。实施例中,采用DEM数据作为正射校正的数据源,对基准影像进行校正和配准,生成底图数据,为无人机三维坐标的解算提供基准数据。数字高程模型(Digital Elevation Model),简称DEM。
对于实时图像T,大小为w×h,先对T进行目标检测,剔除面积在设定阈值范围外的目标,并提取剩余目标的最小外接矩形的中心点,构成点集合PT。在PT中任选3个点构成三角形集合,并剔除下面2种类型的三角形:等腰三角形或等边三角形;最小角小于阈值β(即三角形的三个顶点几乎共线),最终得到实时图像的特征三角形集合GT。
对于基准图像S,大小为W×H,首先将其与相应范围内的DEM数据进行配准,使得S中的每一个点都具有三维坐标。其次对S进行目标检测,剔除面积在设定阈值范围外的目标,并提取剩余目标的最小外接矩形的中心点,构成点集合PS。在PS中,选择能被实时图像面积覆盖的任意3个点构成三角形集合,并剔除下面2种类型的三角形:等腰三角形或等边三角形;最小角小于阈值β′,得到三角形集合GS。
等腰三角形或者等边三角形有两条边或者三条边是相等的,即使两个三角形是相似的,也无法准确的判断两个三角形之间顶点的对应关系,所以要剔除等腰或等边三角形,避免出现误差。例如:任意两个等边三角形是相似的,但无法判断其中一个三角形的顶点分别与另一个三角形的哪个顶点是相对应的。
角度阈值β和β′均建议取值范围为5°~10°,当三角形的最小角小于这个范围时,三角形的三个顶点几乎位于同一条直线上。
2、区域空间关系的匹配。
本步骤通过对两个三角形集合中的三角形进行匹配,得到相匹配的三角形对,称为同名三角形对,从而获取实时图像中像素点在基准图像中所对应的同名点及其三维坐标。
具体步骤如下:
2.1根据三角形的判定定理,进行相似空间关系的粗匹配。
采用的相似三角形判定定理为:如果两个三角形的两个角分别对应相等(或三个角分别对应相等),则两个三角形相似,即两角对应相等,两个三角形相似。本发明中选择三角形的最大角和最小角作为对应的两个角。
设三角形集合GT、GS分别有a、b个三角形,计算相似性矩阵MTS大小为a×b,MTS(i,j)为GT中第i个三角形与GS中第j个三角形的最大角之差与最小角之差的和;设定角度阈值α,若MTS(i,j)<α,则认为GT中第i个三角形与GS中第j个三角形相似,记为三角形对所有匹配的三角形对构成集合
角度阈值α建议取值范围为0°~12°。
2.2根据面积比一致性去除伪匹配三角形对。
如果实时图像和基准图像中的任意三对特征点匹配,则由三对特征点分别组成的三角形一定是相似的,同时,两幅图像的分辨率比例决定了这对三角形面积之间的比例,设为λ,且这个比例对于任意匹配的三角形对都是相同的。
计算集合Gmatching中每一组三角形对的面积比,取λ为三角形对面积比值的均值,大多数的三角形对面积比值接近这个值,相差较大的三角形对则认为是伪匹配三角形对。设定比例阈值λ0,认为面积比值在(λ-λ0,λ+λ0)范围内的三角形对为正确匹配三角形对,进行剔除后得到新的匹配三角形对的集合G′matching。
比例阈值λ0建议取值范围为0.1≤λ0≤0.15。
2.3根据旋转角度一致性去除伪匹配三角形对。
实时图像与基准图像之间的旋转角度决定了匹配三角形对之间的旋转角度,但是这个旋转角度对于任意匹配的三角形对都是相同的,利用这个原理,可以对G′matching中的三角形对进行更加精细的筛选。
定义三角形中最小角所在顶点到对边中点的向量所在的方向称为三角形的最小角方向,如图2中三角形ABC的顶点B到对边中点D的向量方向三角形中最大角所在顶点到对边中点的向量所在的方向称为三角形的最大角方向,如图2中的设最大方向角的对应旋转角度为θ,最小方向角的对应旋转角度为θ′。
计算G′matching中每一组三角形对的最大角方向旋转角度θ和最小角方向旋转角度θ′,取θmax和θmin为所有三角形对最大角方向旋转角度和最小角方向旋转角度的均值,大多数三角形对的旋转角度接近这个值,相差较大的三角形对则认为是伪匹配三角形对。设定角度阈值θ0和θ0′,最大角方向和最小角方向的旋转角度均在阈值内,则认为是正确的匹配三角形对。
角度阈值θ0和θ′0建议取值范围为3°~5°。θ为最大方向角的对应旋转角度,θ′为最小方向角的对应旋转角度,从理论上来说,一组三角形对中,θ和θ′应该是相等的,但是由于图像仿射等误差,所以允许θ和θ′有一定容错范围,即一对匹配的三角形对的最大角方向旋转角度θ在(θmax-θ0,θmax+θ0)范围内,且最小角方向旋转角度θ′在(θmin-θ′0,θmin+θ′0)范围内,那么认为该三角形对是正确匹配的三角形对。
3、解算无人机位置。
通过步骤2获取了正确匹配的三角形对,进而得到相对应的匹配点对(即同名点对),则可以获取无人机实时图像中特征点的像素坐标以及对应的平面与高程坐标,再利用现有的空间后方交会方法对飞行图像进行解析,解算无人机的三维坐标,实现顾及高程的三维定位导航。
基准图像与对应区域的高程数据(DEM数据)进行了配准校正,那么基准图像上的每一个像素点都具有三维坐标,即平面地理坐标和高程坐标。通过步骤2可以获得相对应的同名点对p和P(以实时图像中特征点p和基准图像中特征点P为例),那么点P的三维坐标即点p的三维坐标,故可知特征点p的像素坐标和三维坐标。
具体实施时,以上流程可采用计算机软件技术实现自动运行。
本发明中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (7)
1.一种顾及高程的无人机景观匹配定位方法,包括以下步骤:
步骤1、分别构建基准图像和无人机拍摄的实时图像中的区域空间关系,得到相应的三角形集合,所述基准图像经与相应范围内的DEM数据进行配准,图中每一个点都具有三维坐标;
步骤2,基于空间关系相似性的匹配,包括根据相似三角形的判定定理,基于从步骤1所得基准图像的三角形集合和实时图像的三角形集合,进行相似空间关系的粗匹配,初步得到匹配的三角形对,然后根据面积比一致性和旋转角度一致性去除伪匹配三角形对,得到正确匹配的三角形对;
步骤3,解算无人机位置,包括基于步骤2所得正确匹配的三角形对,获取无人机拍摄的实时图像中特征点的像素坐标以及对应的平面与高程坐标,从而解算无人机的三维坐标,实现顾及高程的三维定位导航。
2.根据权利要求1所述顾及高程的无人机景观匹配定位方法,其特征在于:步骤1中,提取区域空间关系的实现方式如下,
获取图像中地物最小外接矩形的中心点,构成特征点集,在特征点集中任选3个点构成三角形集合,并剔除2种类型的三角形,一种为等腰三角形或等边三角形,一种为最小角小于相应预设阈值的三角形。
3.根据权利要求1所述顾及高程的无人机景观匹配定位方法,其特征在于:步骤2中,根据相似三角形的判定定理,进行相似空间关系的粗匹配,实现方式如下,
设步骤1所得基准图像的三角形集合GT有a个三角形,步骤1所得实时图像GS的三角形集合有b个三角形,计算相似性矩阵MTS大小为a×b,MTS(i,j)为GT中第i个三角形与GS中第j个三角形的最大角之差与最小角之差的和;根据预设的角度阈值α,若MTS(i,j)<α,则认为GT中第i个三角形与GS中第j个三角形相似,记为三角形对所有匹配的三角形对构成集合
4.根据权利要求3所述顾及高程的无人机景观匹配定位方法,其特征在于:步骤2中,根据面积比一致性去除伪匹配三角形对,实现方式如下,
计算集合Gmatching中每一组三角形对的面积比,取λ为三角形对面积比值的均值,根据预设的比例阈值λ0,认为面积比值在(λ-λ0,λ+λ0)范围内的三角形对为正确匹配三角形对,进行剔除后得到新的匹配三角形对的集合G′matching。
5.根据权利要求4所述顾及高程的无人机景观匹配定位方法,其特征在于:步骤2中,λ0建议取值范围为0.1≤λ0≤0.15。
6.根据权利要求3所述顾及高程的无人机景观匹配定位方法,其特征在于:步骤2中,根据旋转角度一致性去除伪匹配三角形对,实现方式如下,
计算G′matching中每一组三角形对的最大角方向旋转角度和最小角方向旋转角度,取θmax和θmin为最大角方向旋转角度和最小角方向旋转角度的均值,根据预设的角度阈值θ0和θ′0,当一对匹配的三角形对的最大角方向旋转角度θ在(θmax-θ0,θmax+θ0)范围内,且最小角方向旋转角度θ′在(θmin-θ′0,θmin+θ′0)范围内,认为该三角形对是正确匹配的三角形对。
7.根据权利要求6所述顾及高程的无人机景观匹配定位方法,其特征在于:角度阈值θ0和θ′0建议取值范围为3°~5°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910270586.6A CN110160503B (zh) | 2019-04-04 | 2019-04-04 | 一种顾及高程的无人机景观匹配定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910270586.6A CN110160503B (zh) | 2019-04-04 | 2019-04-04 | 一种顾及高程的无人机景观匹配定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110160503A true CN110160503A (zh) | 2019-08-23 |
CN110160503B CN110160503B (zh) | 2021-06-04 |
Family
ID=67639042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910270586.6A Active CN110160503B (zh) | 2019-04-04 | 2019-04-04 | 一种顾及高程的无人机景观匹配定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110160503B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111457896A (zh) * | 2020-04-20 | 2020-07-28 | 中国人民解放军空军航空大学 | 一种单像空间后方交会非迭代方法 |
CN113361552A (zh) * | 2020-03-05 | 2021-09-07 | 西安邮电大学 | 定位方法及装置 |
CN117288168A (zh) * | 2023-11-24 | 2023-12-26 | 山东中宇航空科技发展有限公司 | 一种低功耗的无人机城市建筑航拍系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021381A1 (en) * | 2001-07-25 | 2003-01-30 | Reiner Koppe | Method and device for the registration of two 3D image data sets |
CN103236083A (zh) * | 2013-05-06 | 2013-08-07 | 北京建筑工程学院 | 基于立体影像库的实时三维测量方法 |
CN107063193A (zh) * | 2017-03-17 | 2017-08-18 | 东南大学 | 基于全球卫星导航系统动态后处理技术航空摄影测量方法 |
CN107330927A (zh) * | 2017-06-01 | 2017-11-07 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 机载可见光图像定位方法 |
-
2019
- 2019-04-04 CN CN201910270586.6A patent/CN110160503B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021381A1 (en) * | 2001-07-25 | 2003-01-30 | Reiner Koppe | Method and device for the registration of two 3D image data sets |
CN103236083A (zh) * | 2013-05-06 | 2013-08-07 | 北京建筑工程学院 | 基于立体影像库的实时三维测量方法 |
CN107063193A (zh) * | 2017-03-17 | 2017-08-18 | 东南大学 | 基于全球卫星导航系统动态后处理技术航空摄影测量方法 |
CN107330927A (zh) * | 2017-06-01 | 2017-11-07 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 机载可见光图像定位方法 |
Non-Patent Citations (1)
Title |
---|
朱红等: "近景影像三角网内插点密集匹配方法", 《测绘科学》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361552A (zh) * | 2020-03-05 | 2021-09-07 | 西安邮电大学 | 定位方法及装置 |
CN111457896A (zh) * | 2020-04-20 | 2020-07-28 | 中国人民解放军空军航空大学 | 一种单像空间后方交会非迭代方法 |
CN111457896B (zh) * | 2020-04-20 | 2021-09-14 | 中国人民解放军空军航空大学 | 一种单像空间后方交会非迭代方法 |
CN117288168A (zh) * | 2023-11-24 | 2023-12-26 | 山东中宇航空科技发展有限公司 | 一种低功耗的无人机城市建筑航拍系统 |
CN117288168B (zh) * | 2023-11-24 | 2024-01-30 | 山东中宇航空科技发展有限公司 | 一种低功耗的无人机城市建筑航拍系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110160503B (zh) | 2021-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Photogrammetric processing of low‐altitude images acquired by unpiloted aerial vehicles | |
CN110675450B (zh) | 基于slam技术的正射影像实时生成方法及系统 | |
Xie et al. | Study on construction of 3D building based on UAV images | |
WO2019029099A1 (zh) | 基于图像梯度联合优化的双目视觉里程计算方法 | |
CN107808362A (zh) | 一种基于无人机pos信息与图像surf特征结合的图像拼接方法 | |
CN105809687A (zh) | 一种基于图像中边沿点信息的单目视觉测程方法 | |
CN110160503A (zh) | 一种顾及高程的无人机景观匹配定位方法 | |
CN108665499B (zh) | 一种基于视差法的近距飞机位姿测量方法 | |
CN106940181B (zh) | 一种无人机影像像控分布网构建与航片可选范围匹配方法 | |
CN113361499B (zh) | 基于二维纹理和三维姿态融合的局部对象提取方法、装置 | |
CN108917753A (zh) | 基于从运动恢复结构的飞行器位置确定方法 | |
CN109376208A (zh) | 一种基于智能终端的定位方法、系统、存储介质及设备 | |
Lee et al. | Vision-based terrain referenced navigation for unmanned aerial vehicles using homography relationship | |
Zhao et al. | RTSfM: Real-time structure from motion for mosaicing and DSM mapping of sequential aerial images with low overlap | |
CN103345765B (zh) | 基于dsp+fpga的移动平台下运动目标检测装置及其方法 | |
CN108225273A (zh) | 一种基于传感器先验知识的实时跑道检测方法 | |
CN104751451B (zh) | 基于无人机低空高分辨率影像的密集点云提取方法 | |
CN105389819A (zh) | 一种鲁棒的半标定下视图像极线校正方法及系统 | |
CN107063191B (zh) | 一种摄影测量区域网整体相对定向的方法 | |
CN104121884B (zh) | 卫星影像像元观测天顶角和方位角的计算方法 | |
CN108921896A (zh) | 一种融合点线特征的下视视觉罗盘 | |
CN110968910B (zh) | 一种双视线正交的激光雷达卫星姿态设计方法及控制系统 | |
CN110411449B (zh) | 一种航空侦察载荷目标定位方法、系统及终端设备 | |
Kim et al. | Target detection and position likelihood using an aerial image sensor | |
Wang et al. | Automated mosaicking of UAV images based on SFM method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |