CN109612445B - 一种基于无人机的WebGIS平台下高精度地形建立方法 - Google Patents

一种基于无人机的WebGIS平台下高精度地形建立方法 Download PDF

Info

Publication number
CN109612445B
CN109612445B CN201811541925.1A CN201811541925A CN109612445B CN 109612445 B CN109612445 B CN 109612445B CN 201811541925 A CN201811541925 A CN 201811541925A CN 109612445 B CN109612445 B CN 109612445B
Authority
CN
China
Prior art keywords
terrain
aerial vehicle
unmanned aerial
steps
point
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.)
Active
Application number
CN201811541925.1A
Other languages
English (en)
Other versions
CN109612445A (zh
Inventor
李竹有
李艳生
张社荣
王超
姜佩奇
孙洪稳
王晓莉
孔买群
张玉彬
刘山川
朱罗方
樊贵明
浦瑞
杨国华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinohydro Bureau 14 Co Ltd
Original Assignee
Sinohydro Bureau 14 Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinohydro Bureau 14 Co Ltd filed Critical Sinohydro Bureau 14 Co Ltd
Priority to CN201811541925.1A priority Critical patent/CN109612445B/zh
Publication of CN109612445A publication Critical patent/CN109612445A/zh
Application granted granted Critical
Publication of CN109612445B publication Critical patent/CN109612445B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种基于无人机的WebGIS平台下高精度地形建立方法,属于地形测绘领域,本发明的目的在于提供一种基于WebGIS的精细化且轻量化的地形建立与应用方式,其步骤为步骤1.航测基点布设与定位;步骤2.无人机图像采集;步骤3.图像识别与点云生成;步骤4.DEM模型生成;步骤5.地形瓦片数据的生成;步骤6.地形瓦片数据的Web端解析。本发明为整个施工管理可视化平台的搭建提供高精度的地形信息。

Description

一种基于无人机的WebGIS平台下高精度地形建立方法
技术领域
本发明属于无人机测绘领域,更具体的说涉及一种基于无人机的WebGIS平台下的高精度地形建立方法。
背景技术
随着我国工程建设信息化水平的不断提高,越来越多的工程项目开始采用可视化手段进行施工过程的精细化管理。通常可视化要素利用物联网平台采集,通过无线传输传至数据库,再经由后台处理传至前端,从而实现与用户之间的交互。大型工程施工管控要素众多,为了保证数据处理的效率和展示的流畅,传统的施工可视化管理平台大多采用C/S架构实现。针对当前国际通用的EPC、PMC、BOT、 PPP等管理模式,施工质量的管控人员涵盖业主、监理、总承包方、施工方,从项目部总经理到专项经理再到部门员工,均需对项目进行全面直观的把控,而传统的C/S架构系统仅针对施工现场管理人员,其本身存在时效性和地域性问题,无法实现多维度的协同管理,因此越来越无法适应当前管理模式的需要。
WebGIS是将GIS技术应用到Internet平台的产物,是GIS技术的延伸和发展。WebGIS具有运行速度快、展示效果好、服务人群多的特点,能够实现施工全范围的覆盖和全网络端的共享。地理信息数据是WebGIS的核心,传统的 3D-WebGIS平台的地理信息展示主要有两种方式,第一种方式是将处理后的卫星影像数据导入平台当中,这种方式产生的地形精度过低,往往会造成部分施工管控要素埋入地形或处于悬空状态的问题,无法与地形很好的结合;第二种方式是将地形整体导入,精度高但加载速度慢,极易造成浏览器的崩溃。
现有的基于WebGIS的施工管理可视化平台仍无一种合理有效的地形展示方式。
发明内容
本发明的目的在于提供一种基于WebGIS的精细化且轻量化的地形建立与应用方式,为整个施工管理可视化平台的搭建提供高精度的地形信息,同时尽可能降低地形加载对浏览器负载的影响。
为了实现上述目的,本发明是通过以下技术方案实现的:所述的基于无人机的WebGIS平台下高精度地形建立方法由以下几个步骤建立:
步骤1.航测基点布设与定位;
步骤2.无人机图像采集;
步骤3.图像识别与点云生成;
步骤4.DEM模型生成;
步骤5.地形瓦片数据的生成;
步骤6.地形瓦片数据的Web端解析。
优选的,所述的步骤1航测基点布设与定位包括以下两个步骤:
c.航测基点布设:航测基点在布设的过程中应尽可能将航测基点设置在施工范围的边缘区域;
航测基点坐标确定:航测基点坐标系选用当地的工程坐标系,每一个航测基点采集3次,取平均值。
优选的,所述的步骤2无人机图像采集包括以下步骤:
d.软硬件选择;
e.起飞点选择:起降点应选在平坦、空旷的地域;
f.航线规划:在最优航线计算过程中,首先需根据最大航向重叠度和最大旁向重叠度得到图像相邻边界的位移差,按照这个距离将整个区域分网,最后通过蚁群算法,将采集内容覆盖全部网格作为限定条件,计算出航线路径最优的结果。
优选的,所述的步骤3图像识别与点云生成包括以下步骤:
g.获取无人机拍照相机的标定参数,利用OpenCV进行中值滤波和直方图均匀化处理;
h.利用SIFT算法进行特征点的提取;
i.用k-d树结构辅助,利用欧氏距离作为标准,使用最邻近算法进行特征点匹配;
j.通过交叉过滤和基础矩阵过滤的方式剔除误匹配点,得到高质量的匹配点对集合;
k.对匹配出的特征点对进行稀疏重构,恢复其在三维空间中的位置;
l.对稀疏重构后的点云进行密集重构,采用PMVS算法,经过区域匹配、范围扩展、滤波,完成密集重构,得到加密点云。
优选的,所述的步骤4DEM模型生成包括以下步骤:
c.将加密后的点云导出,将点云数据转换成三角网数据并导出;
d.利用MapMatric软件转换功能将三角网数据转换成国家标准的DEM格式,并进行存储。
9.根据权利要求1所述的一种基于无人机的WebGIS平台下高精度地形建立方法,其特征是:所述的步骤5地形瓦片数据的生成包括以下步骤:
c.实现图像识别和表面模型的构建;
d.通过坐标提取、坐标转换、gdal渲染与属性映射,得到标准的地形瓦片数据。
优选的,所述的步骤1航线航测基点按照如下方式进行布置:在挖除边坡表土后开挖一直径为40cm的圆形孔或40cm×40cm的正方形孔约40cm深,用钢筋砼浇注底盘至地面高度。
优选的,所述的航线规划:在最优航线计算过程中,首先需根据最大航向重叠度和最大旁向重叠度得到图像相邻边界的位移差,其计算公式如下:
公式1:
Figure BDA0001908383760000031
公式2:
Figure BDA0001908383760000032
其中,H为相机距地面的高度,θ为相机镜头视角,pxmax、pymax为允许最大重叠影像边长,lx、ly为相机拍摄范围边长,Dxmax、Dymax为相邻边界的最大位移差。
本发明有益效果:
1.生成的地形精度相较于卫星影像生成的地形精度大大提高,完全能够满足工程建设过程中地形应用的需要。
2.生成的地形加载效率快,采用了LOD加载模式,避免了因地形整体过大导致前端负载过大的问题。
3.生成过程便捷高效,通过多次调用脚本依次对标准化的航测图片和生成的数据进行处理,可以很方便的实现地形瓦片数据的生成。
4.本发明方法步骤简单、实现方便、可操作性强,大大简化了了地形的数据过程和处理效率,而且能确保满足施工管理可视化平台的搭建要求。
附图说明
图1为本发明的监测基点埋设结构示意图;
图2为本发明的局部地形整体航拍图;
图3为本发明的无人机相邻边界最大位移差计算方法;
图4为本发明的点云数据生成效果图;
图5为本发明的DEM模型生成效果图;
图6为本发明的地形瓦片解析效果图;
具体实施方式
下面将结合本发明实施例和附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属本发明保护范围。
如图所示,本发明的实现主要包括以下步骤:
步骤1,航测基点布设与定位:
航测基点布设:航测监 测基点用于将生成的地形坐标系转换为当地工程坐标系,在布设的过程中应尽可能将航测基点设置在施工范围的边缘区域,数量不宜超过 5个。航测基点颜色需显眼、中心点需明确。
航测基点坐标确定:坐标系选用当地的工程坐标系,将移动GPS的坐标系调成当地的工程坐标系,放置在航测基点的位置进行坐标信息的采集,每一个航测基点采集3次,取平均值,当出现误差较大的采集点时,需重新采集。
监测基点按照如下方式进行布置:在挖除边坡表土后开挖一直径为40cm的圆形孔或40cm×40cm的正方形孔约40cm深,用钢筋砼浇注底盘至地面高度,并涂刷蓝色油漆作为标记,监测基点示意图如图1所示。监测基点设置在位移稳定的区域,避免在松动的表层上设点。在位移稳定的基础上,监测基点尽可能靠近监测坡体,以缩小无人机航测区域,提高监测精度。本实施例中,所埋设的监测基点为3个,整体航拍图如图2所示。
步骤2,无人机图像采集:
软硬件选择:无人机选用大疆经纬M600pro搭配X5R云台,无人机飞控软件选用DJIGS Pro。
起飞点选择:无人机的起降点应选在平坦、空旷的地域。针对面积较大的区域,从终点到着陆点的飞行时间会对图像采集的时间产生显著影响,因此需要实地考察后确定。
航线规划:运用python语言编写最优路径算法进行航线的规划,将保证在全范围覆盖的前提下使航线最短。在航线规划过程中,需保证航向重叠度最佳在 80%以上,最低不低于60%;旁向重叠度最佳在70%以上,最低不低于30%。在最优航线计算过程中,首先需根据最大航向重叠度和最大旁向重叠度得到图像相邻边界的位移差,计算方法见公式1和公式2,之后按照这个距离将整个区域分网,最后通过蚁群算法,将采集内容覆盖全部网格作为限定条件,计算出航线路径最优的结果。在图像采集过程中,云台应始终垂直向下。
公式1:
Figure BDA0001908383760000051
公式2:
Figure BDA0001908383760000052
其中,H为相机距地面的高度,θ为相机镜头视角,pxmax、pymax为允许最大重叠影像边长,lx、ly为相机拍摄范围边长,Dxmax、Dymax为相邻边界的最大位移差。
航向重叠度指无人机沿航线进行航摄时,同一航线上所拍摄的两张图像重叠区域占单张图片大小的比例;旁向重叠度指无人机沿航线进行航摄时,两航线之间所拍摄相邻的两张图像重叠区域占单张图片大小的比例。在无人机飞行前需进行航线规划,通过规定无人机飞行高度,在地图中通过无人机进行定点标出所需航摄区域,无人机会根据航向重叠度及旁向重叠度的要求按照设定的航线进行图像采集。飞行过程中应保证无人机所载相机镜头垂直朝下,与无人机的航线相垂直。根据图3,无人机的飞行高度指从无人机起飞位置计算起的上升高度;航摄区域是指包含监测边坡、监测点及监测基点在内的区域,本实施例中,无人机飞行高度为100m。
步骤3,图像识别与点云生成:运用开发的系统将所采集到的无人机航摄图像进行处理及识别,并对整体边坡进行三维重建,并通过监测基点计算出坐标转换七参数,将加密后的点云全部转换到工程坐标系,生成效果如图4所示。
获取X5R相机的标定参数,利用OpenCV开源库进行中值滤波和直方图均匀化处理。
利用SIFT算法进行特征点的提取。使用图像金字塔表达尺度空间,利用高斯差分方程,构造DoG图像金字塔。利用OpenCV进行SIFT特征点检查,提取出 SIFT特征点。
用k-d树结构辅助,利用欧氏距离作为标准,使用最邻近算法进行特征点匹配。
通过交叉过滤和基础矩阵过滤的方式剔除误匹配点,得到高质量的匹配点对集合。
对匹配的出的特征点对进行稀疏重构,恢复其在三维空间中的位置。依托 OpenCV来进行图像之间的基础矩阵和本征矩阵的计算。利用相机标定参数三角化重建匹配点坐标并对已重建的图像作集束优化,得到新增图像的投影矩阵。循环重构直到所有的图像重构完成,得到稀疏重构的点云。
对稀疏重构后的点云进行密集重构。采用PMVS算法,经过区域匹配、范围扩展、滤波,完成密集重构,得到加密点云。
步骤4,DEM模型生成:
将加密后的点云导出为.las格式,利用VTK中的Delaunay将点云数据转换成三角网数据,并导出为基于文本的ASCII格式。
打开MapMatric软件,利用软件中的DEM转换功能将ASCII格式转换成国家标准的DEM格式,并进行存储。
步骤5,地形瓦片数据的生成:
地形瓦片数据的生成采用python脚本实现,首先借助opencv-python类库实现图像识别和表面模型的构建,之后借助gdal工具包实现dem地形信息到tin 格式地形信息的转换。数字高程模型的标准数据格式为.tif格式,通过坐标提取、坐标转换、gdal渲染与属性映射,即可得到标准的流式瓦片地形数据。地形数据以xml文件为索引,统一存储在一个文件夹当中,每一级分辨率又单独存储,由hdr文件、kml文件和terrain文件构成,卫星地图的精度是米级,瓦片的精度是厘米级,精度由米级提高到了厘米级,同时整体地形导入100兆左右需要1分钟左右,用瓦片加载时只需要1-3秒,尽可能降低地形加载对浏览器负载的影响。
步骤6,地形瓦片数据的Web端解析
地形数据的Web端解析利用Cesium.js插件实现,在引入插件的时候将地形信息以“new Cesium.CesiumTerrainProvider”的形式导入,最终实现地形数据的web端解析。
步骤3、步骤4、步骤5、步骤6已通过系统开发实现,步骤3、步骤4、步骤5为解析脚本、步骤6为Web前端解析程序,过程便捷高效,通过多次调用脚本依次对标准化的航测图片和生成的数据进行处理,可以很方便的实现地形瓦片数据的生成。采用了LOD加载模式,生成的地形加载效率快,避免了因地形整体过大导致前端负载过大的问题。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。

Claims (3)

1.一种基于无人机的WebGIS平台下高精度地形建立方法,其特征是:所述的基于无人机的WebGIS平台下高精度地形建立方法包括以下步骤:
步骤1.航测基点布设与定位;
所述的步骤1包括以下步骤:
a.航测基点布设:航测基点在布设的过程中应尽可能将航测基点设置在施工范围的边缘区域;
b.航测基点坐标确定:航测基点坐标坐标系选用当地的工程坐标系每一个航测基点采集多次,取平均值;
步骤2.无人机图像采集;
所述的步骤2包括以下步骤:
a.软硬件选择;
b.起飞点选择:起降点应选在平坦、空旷的地域;
c.航线规划:在最优航线计算过程中,首先需根据最大航向重叠度和最大旁向重叠度得到图像相邻边界的位移差,按照这个距离将整个区域分网,最后通过蚁群算法,将采集内容覆盖全部网格作为限定条件,计算出航线路径最优的结果;
步骤3.图像识别与点云生成;
所述的步骤3图像识别与点云生成包括以下步骤:
a.获取无人机拍照相机的标定参数,利用OpenCV进行中值滤波和直方图均匀化处理;
b.利用SIFT算法进行特征点的提取;
c.用k-d树结构辅助,利用欧氏距离作为标准,使用最邻近算法进行特征点匹配;
d.通过交叉过滤和基础矩阵过滤的方式剔除误匹配点,得到高质量的匹配点对集合;
e.对匹配的出的特征点对进行稀疏重构,恢复其在三维空间中的位置;
对稀疏重构后的点云进行密集重构,采用PMVS算法,经过区域匹配、范围扩展、滤波,完成密集重构,得到加密点云;
步骤4.DEM模型生成;
所述的步骤4包括以下步骤:
a.将加密后的点云导出,将点云数据转换成三角网数据并导出;
b.利用MapMatric软件转换功能将三角网数据转换成国家标准的DEM格式,并进行存储;
步骤5.地形瓦片数据的生成;
所述的步骤5包括以下步骤:
a.实现图像识别和表面模型的构建;
b.通过坐标提取、坐标转换、gdal渲染与属性映射,得到标准的地形瓦片数据;
步骤6.地形瓦片数据的Web端解析。
2.根据权利要求1所述的一种基于无人机的WebGIS平台下高精度地形建立方法,其特征是:所述的步骤1航线航测基点按照如下方式进行布置:在挖除边坡表土后开挖一直径为40cm的圆形孔或40cm×40cm的正方形孔约40cm深,用钢筋砼浇注底盘至地面高度。
3.根据权利要求1所述的一种基于无人机的WebGIS平台下高精度地形建立方法,其特征是:所述的航线规划:在最优航线计算过程中,首先需根据最大航向重叠度和最大旁向重叠度得到图像相邻边界的位移差,其计算公式如下:
公式1:
Figure FDA0002930443700000021
公式2:
Figure FDA0002930443700000022
其中,H为相机距地面的高度,θ为相机镜头视角,pxmax、pymax为允许最大重叠影像边长,lx、ly为相机拍摄范围边长,Dxmax、Dymax为相邻边界的最大位移差。
CN201811541925.1A 2018-12-17 2018-12-17 一种基于无人机的WebGIS平台下高精度地形建立方法 Active CN109612445B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811541925.1A CN109612445B (zh) 2018-12-17 2018-12-17 一种基于无人机的WebGIS平台下高精度地形建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811541925.1A CN109612445B (zh) 2018-12-17 2018-12-17 一种基于无人机的WebGIS平台下高精度地形建立方法

Publications (2)

Publication Number Publication Date
CN109612445A CN109612445A (zh) 2019-04-12
CN109612445B true CN109612445B (zh) 2021-04-30

Family

ID=66009470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811541925.1A Active CN109612445B (zh) 2018-12-17 2018-12-17 一种基于无人机的WebGIS平台下高精度地形建立方法

Country Status (1)

Country Link
CN (1) CN109612445B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021035608A1 (zh) * 2019-08-29 2021-03-04 深圳市大疆创新科技有限公司 航线生成方法、地面端设备、无人机、系统和存储介质
US11709819B2 (en) 2020-09-30 2023-07-25 International Business Machines Corporation Validating test results using a blockchain network
CN114201633B (zh) * 2022-02-17 2022-05-17 四川腾盾科技有限公司 一种用于无人机视觉定位的自适应卫星图像生成方法
CN115129401B (zh) * 2022-08-29 2022-11-15 吉奥时空信息技术股份有限公司 基于Cesium实现无偏移无变形叠加多源数据的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102103202B (zh) * 2010-12-01 2013-01-09 武汉大学 一种融合影像的机载激光雷达数据半监督分类方法
KR101329460B1 (ko) * 2012-03-07 2013-11-15 인천대학교 산학협력단 Gis 정보 기반 센서 배치 시뮬레이터
CN105976426B (zh) * 2016-04-27 2019-03-05 中国电子科技集团公司第二十八研究所 一种快速的三维地物模型构建方法
CN108536829B (zh) * 2018-04-11 2019-04-26 中国中医科学院中药研究所 一种提高无人机航测数据生成瓦片地图效率的方法
CN108762300A (zh) * 2018-05-31 2018-11-06 智飞智能装备科技东台有限公司 一种用于无人机组飞行的地面控制编程系统

Also Published As

Publication number Publication date
CN109612445A (zh) 2019-04-12

Similar Documents

Publication Publication Date Title
CN109612445B (zh) 一种基于无人机的WebGIS平台下高精度地形建立方法
WO2020156273A1 (zh) 一种园林空间数字化测绘和三维可视化方法
CN111091613A (zh) 基于无人机航测的三维实景建模方法
CN114518104B (zh) 基于动态遥感监测技术的国土测绘方法、系统及存储介质
CN108109437B (zh) 一种基于地图特征的无人机自主航线提取生成方法
AU2021104249A4 (en) Pesticide application method and system based on air-ground collaboration
Laefer et al. 2015 aerial laser and photogrammetry survey of Dublin city collection record
CN111322994A (zh) 基于无人机倾斜摄影的房屋密集区大比例尺地籍测量方法
CN101335431A (zh) 基于机载激光雷达数据的架空送电线路优化选线方法
CN109670789B (zh) 一种用于生产建设项目水土保持的遥感监测系统
Lv et al. Research on the technology of LIDAR data processing
CN109238242A (zh) 基于立体大场景的输电线路路径优化选线
CN114065339A (zh) 一种基于三维可视化模型的高塔建设选址方法
CN110362895A (zh) 一种基于bim+gis技术的征地拆迁应用管理系统
CN116468869A (zh) 一种基于遥感卫星影像的实景三维建模方法、设备及介质
CN106969753B (zh) 基于电力设计应用的无人机航拍数据处理方法
Zhou et al. Application of UAV oblique photography in real scene 3d modeling
Yijing et al. Construction and analysis of 3D scene model of landscape space based on UAV oblique photography and 3D laser scanner
Williams Unmanned aerial vehicle photography: Exploring the medieval city of Merv, on the Silk Roads of Central Asia
CN108692710B (zh) 一种公路宗地测量方法及系统
CN116129064A (zh) 电子地图生成方法、装置、设备及存储介质
CN113744393B (zh) 一种多层级边坡滑坡变化监测方法
KR102587445B1 (ko) 드론을 이용하여 시계열정보가 포함된 3차원 지도의 제작 방법
CN109827503A (zh) 一种田间作物精准定位方法和装置
CN115752381A (zh) 一种基于无人机遥感技术的矿山监测方法

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
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A high-precision terrain establishment method based on drone in WebGIS platform

Effective date of registration: 20231218

Granted publication date: 20210430

Pledgee: Shenzhen Deyuan Commercial Factoring Co.,Ltd.

Pledgor: SINOHYDRO BUREAU 14 Co.,Ltd.

Registration number: Y2023980072201