CN111444872A - 一种丹霞地貌参数测量方法 - Google Patents
一种丹霞地貌参数测量方法 Download PDFInfo
- Publication number
- CN111444872A CN111444872A CN202010252595.5A CN202010252595A CN111444872A CN 111444872 A CN111444872 A CN 111444872A CN 202010252595 A CN202010252595 A CN 202010252595A CN 111444872 A CN111444872 A CN 111444872A
- Authority
- CN
- China
- Prior art keywords
- danxia
- data
- mountain
- landform
- geomorphic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000005259 measurement Methods 0.000 claims abstract description 37
- 238000012545 processing Methods 0.000 claims abstract description 18
- 230000000007 visual effect Effects 0.000 claims abstract description 7
- 239000000178 monomer Substances 0.000 claims description 27
- 239000011435 rock Substances 0.000 claims description 15
- 230000003628 erosive effect Effects 0.000 claims description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 claims description 10
- 238000012876 topography Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 8
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 238000005530 etching Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 6
- 238000013507 mapping Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 241000720991 Illicium Species 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 229940045109 genistein Drugs 0.000 description 1
- TZBJGXHYKVUXJN-UHFFFAOYSA-N genistein Natural products C1=CC(O)=CC=C1C1=COC2=CC(O)=CC(O)=C2C1=O TZBJGXHYKVUXJN-UHFFFAOYSA-N 0.000 description 1
- 235000006539 genistein Nutrition 0.000 description 1
- ZCOLJUOHXJRHDI-CMWLGVBASA-N genistein 7-O-beta-D-glucoside Chemical compound O[C@@H]1[C@@H](O)[C@H](O)[C@@H](CO)O[C@H]1OC1=CC(O)=C2C(=O)C(C=3C=CC(O)=CC=3)=COC2=C1 ZCOLJUOHXJRHDI-CMWLGVBASA-N 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Remote Sensing (AREA)
- Software Systems (AREA)
- Astronomy & Astrophysics (AREA)
- Multimedia (AREA)
- Computer Graphics (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种丹霞地貌参数测量方法,属于无人机测绘技术领域。包括以下步骤:基于多源卫星遥感数据解译丹霞地貌;根据地形、地貌、分辨率确定数据采集区域,划定仿地飞行范围并进行航线规划;采用井字型航线和立面航线进行智能数据采集,在复杂区域增加人工数据采集;航拍图像处理;丹霞地貌参数测量及分析。本发明可以获得典型丹霞地貌全貌,解决了人工难以测量或测量不准确的问题,得到了全覆盖、高精度、多视角的丹霞地貌数据,测量参数全面、准确度高,为丹霞地貌的刻画、保护、开发利用提供了永久的测量数据。
Description
【技术领域】
本发明涉及无人机测绘技术领域,具体涉及一种丹霞地貌参数测量方法。
【背景技术】
丹霞地貌是以陡崖坡为特征的红层地貌。丹霞地貌在我国广泛分布,相对集中分布在东南、西南和西北三个地区。裸露的红色陡崖坡是丹霞地貌最重要的识别要素。由于岩层产状不同,丹霞地貌顶面形态有“斜顶、圆顶、尖顶”等多种类型;由于直接作用于陡崖坡的动力类型差异,可以分为崩塌后壁型、流水侵蚀型、风化改造型、多动力综合型等;由于红层岩性和外力作用的差异,丹霞陡崖坡常有各种洞穴、凹槽、竖向沟槽分布;此外,受间歇性构造隆升的影响,陡崖坡还表现出多层次的坡面组合,现实中丹霞地貌的坡面形态复杂多样。有学者采用面积—高度比曲线法划分丹霞地貌发育阶段,但与现在技术结合不够,定量研究不足。丹霞地貌中丹霞山数量、谷地数量、丹霞山形态、坡度、高度、体积、岩层倾角、产状、走向、丹霞谷地宽度、丹霞地貌形成过程中被蚀去的断面面积、丹霞地貌侵蚀区长度等参数都是丹霞地貌定量描述、研究、开发利用和保护的基础数据。
目前丹霞地貌测量方式主要是地面现场测量,多采用全站仪、激光测距仪等设备进行,丹霞地貌形态则多采用素描、拍照等方式记录。在丹霞地貌测量过程中,现场测量的工作量大,由于丹霞地貌复杂,山陡谷深,往往无法精准、全面地采集丹霞地貌参数,整个测量工作周期长,精度低,由于测量人员对丹霞地貌认知不同,加上丹霞地貌测量条件的限制,往往测量结果难以重复和标准化,导致所测丹霞地貌参数偏差大,难以标准。
随着测绘无人机的发展,大面积和复杂危险地形的测绘得以实现。尤其近些年来,小型智能无人机技术得到了突飞猛进地发展,具有成本低、人员安全风险低、搭载设备多样、操控简单、工作效能高等特点,目前大量应用于数字城市建设、工程地质勘察等领域,极大地提高了作业效率。然而,丹霞山岩壁陡峭,地形复杂,现有的技术无法精准、全面地采集丹霞山的图像数据,如何利用无人机测绘技术对丹霞地貌参数进行测绘,目前还没有相关报道。
【发明内容】
本发明的发明目的在于:针对上述存在的问题,提供一种丹霞地貌参数测量方法,本发明的方法结合丹霞地貌特性,利用优化的多源卫星遥感解译、无人机测绘和图像处理等技术,能够测量丹霞全貌,实现丹霞地貌重要参数标准化测量,具有效率高、测量全面准确、操作简单的优点。
为了实现上述目的,本发明采用的技术方案如下:
一种丹霞地貌参数测量方法,包括以下步骤:
S1,丹霞地貌多源卫星遥感解译:收集待测丹霞地貌区域的资料,包括多源卫星数据、地形数据和地质资料,对不同来源的卫星遥感影像进行处理,综合解译待测丹霞地貌区域的典型丹霞地貌;
S2,无人机数据采集:采用无人机作为数据采集平台;首先室内确定丹霞地貌数据采集区域,并进行航线规划,采用井字形航线和立面航线规划,将设计好的航线导入无人机后,操控无人机按照规划的航线进行无人机倾斜摄影测量,智能采集数据;地形起伏大的区域采用仿地飞行方式获取测量区域的数据;部分地形复杂区域采用人工手动飞行补充航拍数据;
S3,航拍图像处理:采用图像处理软件对无人机拍摄的照片进行去云雾处理,保证照片清晰,采用Pix4D、Context Capture软件按照流程处理无人机采集的数据,选择2000国家大地坐标系,进行空三解算,生成丹霞地貌区域的正射影像、数字地表模型、三维模型和点云;
S4,丹霞地貌参数测量及分析:对步骤S3得到的丹霞地貌模型进行切割,识别并提取丹霞山单体,进行丹霞地貌参数测量。
进一步地,所述步骤S1中的卫星数据来自landsat卫星、高分一号、高分二号、worldview等卫星,卫星遥感影像处理和综合解译包含以下过程:
S101,遥感影像处理:对多源卫星数据进行全色与多光谱波段的配准、融合,正射校正,镶嵌;landsat卫星数据采用第7、第4和第1波段进行组合,能更好地区分丹霞地貌的岩性边界;landsat卫星数据采用第4、第3和第2波段组合,形成分辨率更高的融合影像模拟自然真彩色正射影像图,与地物真实色彩相近,能更好地进行丹霞地貌目视识别;其他类型卫星根据各自波段特征进行组合,用于丹霞地貌识别。
S102,丹霞地貌解译标志:根据landsat卫星741波段组合影像,划分丹霞地貌区域,裁切保留丹霞地貌区域影像;结合裁切的丹霞地貌区域地形数据,在ENVI软件中形成卫星正射影像的三维可视化场景,采用432波段组合形成自然真彩色影像,用于圈定丹霞山单体范围;再通过其他高分辨率卫星数据,如高分一号、高分二号、worldview等进一步确定位置和范围,解译出典型的丹霞山。
进一步地,所述步骤S2中确定丹霞地貌数据采集区域的方法为:通过多源卫星遥感影像解译得到丹霞地貌的范围和位置,识别丹霞山,根据地面分辨率条件,结合奥维地图初步划定无人机数据采集区域的矢量,以kml文件形式保存。
进一步地,所述步骤S2中航线规划的方法为:将确定的采集区域kml文件导入无人机飞控系统,根据分辨率及实地地形情况,确定飞行高度及重叠率;为保证测量的精准,比例尺设置为大于1:1000;在起伏较大的区域,即正负地形高差超过200m的区域,导入该区域地形数据,采用仿地飞行模式,确保数据采集的精准;
航线的航高设置依据以下公式:
航高=(相机镜头焦距×地面分辨率)/相机像元尺寸;
航高单位为米,地面分辨率单位为米,相机镜头焦距单位为毫米,相机像元尺寸为微米;
在划定飞行区域沿主轴方向规划,采用弓字形飞行路线,考虑到负地形采集全面,航向重叠率不低于80%,旁向重叠率不低于65%。
进一步地,所述步骤S2中操控无人机智能采集数据时,按下以下条件进行:
飞行环境:在无风或微风、晴朗、能见度高的正午进行数据采集;
飞行方式:保证飞控系统与无人机连接,按照规划的航线飞行,在有CORS信号的区域,采用无人机飞控系统与CORS站连接,发射差分改正数据给无人机,保证照片POS信息的精准;如没有CORS信号,采用PPK方式进行飞行,以确保采集数据的高准确性;
数据的完整性:在高差较大区域,即正负地形高差超过200m的地形区域,因为负地形复杂,容易造成负地形或者部分丹霞山体局部区域出现数据采集重叠率低的情况,所以在完成自动航线数据采集后,根据地形情况,由人工操作无人机补充复杂负地形区域的数据采集;
数据检查:每个架次飞行完毕,导出影像和pos数据进行质量检查,主要检查照片影像质量、POS信息完备性、航向和旁向重叠率。
进一步地,所述步骤S4丹霞地貌参数测量及分析具体方法为:
S401,模型切割:对步骤S3得到的丹霞地貌数字地表模型、三维模型、点云进行分析,将点云通过Cyclone软件进行过滤处理,除去冗余数据,在Cloud Compare点云处理软件中,将识别确定的丹霞区域切割,得到典型完整的丹霞地貌;
S402,丹霞山单体提取:在Cyclone软件中对切割的丹霞地貌进行分析,有多个丹霞山时,识别丹霞山单体,采用软件中的编辑工具提取丹霞山单体;
S403,丹霞地貌参数选取:在众多丹霞地貌参数中选择最能反映丹霞地貌特征的参数,包括:丹霞山数量、谷地数量、丹霞山形态、坡度、丹霞山高度、丹霞山体积、岩层倾角、产状、走向、丹霞谷地宽度、丹霞地貌形成过程中被蚀去的断面面积、丹霞地貌侵蚀区长度。
S403,丹霞地貌参数测量:基于Cyclone和Cloud Compare软件进行丹霞地貌参数测量,对丹霞山和谷地数量、丹霞地貌系列参数进行测量,测量方法和要求如下:
丹霞山数量:在切割的丹霞地貌中,不连续的丹霞山独立算为丹霞山单体,从而计算丹霞山单体数量;
谷地数量:在切割的丹霞地貌中,独立完整的负地形为谷地,从而计算丹霞谷地数量;
丹霞山形态:水平面截取丹霞山单体底部,形成的边缘轮廓;
丹霞山坡度、产状和走向:采用Cloud Compare软件中的compass工具进行丹霞山坡度、产状和走向测量;
丹霞山高度:计算丹霞山最高处到水平面的垂直距离;
丹霞山体积:在丹霞山顶部形成水平面,计算该水平面与底部水平面之间的山体体积;
岩层倾角:选择要测量的岩层,测量其与水平面夹角;
丹霞谷地宽度:丹霞谷地两侧崖壁上缘之间的宽度;
蚀去断面面积:丹霞地貌形成过程中被蚀去的断面面积,量取蚀去断面面积;
丹霞地貌侵蚀区长度:丹霞地貌侵蚀区最大长度,量取侵蚀区最远两个端点距离;
参数分析:基于Cyclone软件进行丹霞地貌的分析,包括丹霞地貌的分解、合并、放大或缩小,产生多个平面对丹霞地貌进行分段切割分析,丹霞山单体的对比、谷地的对比。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1、本发明首先对卫星遥感影像进行处理,通过不同波段的组合划分丹霞与非丹霞区域、圈定丹霞山和谷地,从而将典型的丹霞地貌解译出来,做好区域划定,为无人机的航线规划提供基础材料,减少航线规划的人工勘测工作量,提高工作效率。
2、本发明在利用无人机采集数据时,针对丹霞山岩壁陡峭和地貌复杂的地形,采用井字型航线规划和立面航线规划,并采用弓字形飞行路线,且在正负地形高差超过200m的区域调整无人机进行仿地飞行方式获得数据,但部分区域仍然无法通过智能航线方式获取完整数据,导致数据缺失,所以对部分难以覆盖的区域采用人工手动飞行,连续拍照方式补充航拍数据,从而保证数据的齐备性,增加了测量的准确性。
3、本发明通过Pix4D、Context Capture软件进行航拍照片空三处理,获得丹霞地貌测量区域的正射影像、地表数字模型、三维模型和点云。在将点云通过Cyclone软件进行过滤处理,除去冗余数据,采用Cyclone软件中的测量工具和Cloud Compare中compass工具,根据丹霞地貌参数定义确定具体测量方法,获得该参数精准的测量数据。解决了过去测量方法无法精准测量丹霞山形态、产状、走向、岩层倾角、体积、侵蚀断面面积、侵蚀区长度等参数的问题。
4、本发明通过无人机进行丹霞地貌数据的采集,可以获得包括丹霞点和谷地在内的丹霞全貌,解决了由于丹霞山陡峭、山与谷地高差较大、地形复杂等客观原因外业人工无法测量或测量不准确的问题,得到了全覆盖、高精度、多视角的丹霞地貌数据,操作简单高效。航拍获取的数据经空三计算后生成数字正射影像、数字表面模型,三维模型和点云。切割典型丹霞地貌模型,通过点云软件在室内可标准测量丹霞地貌系列参数,测量参数全面、准确度高,还可进行所需分析,为丹霞地貌的研究、保护、开发利用提供了永久的测量数据。
【附图说明】
图1是本发明一种丹霞地貌参数测量方法的测量流程图;
图2是本发明具体实施方式中丹霞地貌图像经处理的正射影像;
图3是本发明具体实施方式中丹霞地貌图像经处理的数字地表模型;
图4是本发明具体实施方式中丹霞地貌图像经处理的三维模型;
图5是本发明具体实施方式中丹霞地貌的多个丹霞山单体和谷地单体;
图6是本发明具体实施方式中丹霞山高度的测量示意图;
图7是本发明具体实施方式中丹霞山坡度的测量示意图;
图8是本发明具体实施方式中丹霞山产状的测量示意图;
图9是本发明具体实施方式中丹霞山走向的测量示意图;
图10是本发明具体实施方式中丹霞山岩层倾角的测量示意图;
图11是本发明具体实施方式中丹霞山谷地长度的测量示意图;
图12是本发明具体实施方式中丹霞山谷地宽度测量示意图。
【具体实施方式】
本发明提供一种丹霞地貌参数测量方法,其测量的流程图见图1所示。为了更清楚地表达本发明,以下通过以桂林八角寨丹霞景区为例进行丹霞地貌参数测量,从而对本发明作进一步详细说明,其具体实施方式及步骤如下:
一种丹霞地貌参数测量方法,包括以下步骤:
S1,丹霞地貌多源卫星遥感解译:收集待测丹霞地貌区域的资料,包括多源卫星数据、地形数据和地质资料,对不同来源的卫星遥感影像进行处理,综合解译待测丹霞地貌区域的典型丹霞地貌;卫星数据来自landsat8高分一号、高分二号和worldview卫星,卫星遥感影像处理和综合解译包含以下过程:
S101,遥感影像处理:对landsat8卫星数据进行全色与多光谱波段的配准、融合,正射校正;landsat8卫星数据采用第7、第4和第1波段进行组合,能更好反映丹霞地貌的岩性边界;landsat8卫星数据采用第4、第3和第2波段组合,形成高分辨融合影像模拟自然真彩色正射影像图,与地物真实色彩相近,能更好地进行目视识别;
S102,丹霞地貌解译标志:根据landsat8卫星741波段组合影像,划分丹霞地貌区域,裁切保留丹霞地貌区域影像;结合裁切的丹霞地貌地形数据,在ENVI软件中形成卫星正射影像的三维可视化场景,采用432波段组合,圈定典型的丹霞山和丹霞谷地负地形区域;再通过高分一号、高分二号、worldview卫星数据对比,进一步确定范围、位置,解译出典型的丹霞地貌。
S2,无人机数据采集:采用精灵phantom4 RTK无人机作为数据采集平台;首先室内确定丹霞地貌数据采集区域,并进行航线规划,采用井字形航线和立面航线规划,将设计好的航线导入无人机后,操控无人机按照规划的航线进行无人机倾斜摄影测量,智能采集数据;地形起伏大的区域采用仿地飞行方式获取测量区域的数据;部分地形复杂区域采用人工手动飞行补充航拍数据;具体为:
S201,确定丹霞地貌地貌数据采集区域:通过步骤S1卫星遥感影像解译得到的范围、位置,以及地面分辨率条件,结合奥维地图初步划定无人机数据采集区域的矢量,以kml文件形式保存;
S202:航线规划的方法为:将确定的采集区域kml文件导入无人机飞控系统,根据分辨率及实地地形情况,确定飞行高度及重叠率;为保证测量的精准,比例尺设置为大于1:1000;在起伏较大的区域,即正负地形高差超过200m的区域,导入该区域地形数据,采用仿地飞行模式,确保数据采集的精准;
航线的航高设置依据以下公式:
航高=(相机镜头焦距×地面分辨率)/相机像元尺寸;
地面分辨率单位为米,相机镜头焦距单位为毫米,相机像元尺寸为微米;航高单位为米。
在划定飞行区域沿主轴方向规划,采用弓字形飞行路线,考虑到正负地形采集全面,本实施例中的飞行参数见表1。
表1八角寨丹霞地貌航线规划参数
飞行高度(m) | 航向重叠率 | 旁向重叠率 | 航测面积(km2) | |
参数 | 250 | 80% | 70% | 15 |
S203:操控无人机智能采集数据时,按下以下条件进行:
飞行环境:在无风或微风、晴朗、能见度高的正午进行数据采集;
飞行方式:保证飞控系统与无人机连接,按照规划的航线飞行,在有CORS信号的区域,采用无人机飞控系统与CORS站连接,发射差分改正数据给无人机,保证照片POS信息的精准;如没有CORS信号,采用PPK方式进行飞行,以确保采集数据的高准确性;本次飞行采用的是CORS信号方式。
数据的完整性:在高差较大,即正负地形高差超过200m的地形,由于负地形复杂,在完成自动航线数据采集后,根据地形,由人工操作无人机补充采集难以覆盖区域的数据;
数据检查:每个架次飞行完毕,导出影像和pos数据进行质量检查,主要检查照片影像质量、POS信息完备性、航向和旁向重叠率。
S3,图像处理:采用图像处理软件对照片进行去云雾处理,保证照片清晰,采用Pix4D、Context Capture软件按照流程处理无人机采集的数据,选择2000国家大地坐标系,进行空三解算,生成岩溶坡立谷地貌的正射影像、数字地表模型、三维模型和点云;得到的正射影像、数字地表模型、三维模型分别见图2-4。
S4,丹霞地貌参数测量及分析:对步骤S3得到的丹霞地貌模型进行切割,识别并提取丹霞地貌,进行丹霞地貌参数测量。
S401,模型切割:对步骤S3得到的丹霞地貌数字地表模型、三维模型、点云进行分析,将点云通过Cyclone软件进行过滤处理,除去冗余数据,在Cloud Compare点云处理软件中,将丹霞地貌与非丹霞地貌切割开来,得到典型的丹霞完整地貌;
S402,丹霞山和谷地单体提取:对切割的丹霞地貌进行分析,在Cyclone软件中识别并提取丹霞山单体和谷地单体,八角寨景区案例中提取了6个丹霞山和4个谷地。如图5。
S403,丹霞地貌参数选取:在多种丹霞参数中选择最能反映丹霞地貌特征的参数,包括:丹霞山数量、谷地数量、丹霞山形态、坡度、丹霞山高度、丹霞山体积、岩层倾角、产状、走向、丹霞地貌形成过程中被蚀去的断面面积、丹霞地貌侵蚀区长度。
S403,丹霞地貌参数测量:基于Cyclone软件、Cloud Compare软件进行丹霞地貌参数测量,对丹霞山单体数量、谷地单体数量、相关丹霞地貌参数进行测量,方法如下:
丹霞山数量:在切割的丹霞地貌中,不连续的丹霞山独立算为丹霞山单体,从而计算丹霞山单体数量,见图5,计算得到6个丹霞山和4个谷地。
谷地数量:在切割的丹霞地貌中,独立完整的负地形为谷地,从而计算丹霞谷地数量;
丹霞山形态:水平面截取丹霞山单体底部,形成的边缘轮廓;
丹霞山高度:在丹霞山单体最高处与丹霞山底部的垂直距离,计算丹霞山底平面到山顶平面的距离;见图6,测得其中一个丹霞山单体高度为230.401m;
丹霞山体积:在丹霞山顶部形成水平面,计算该水平面与底部水平面之间的山体体积。
丹霞山坡度、产状和走向:采用Cloud Compare软件中的compass工具进行丹霞山坡度、产状和走向测量,见图7-9。
岩层倾角:选择要测量的岩层,测量其与水平面夹角,见图10。
谷地长度测量:选取谷地的起点和终点,量取两点间的直线距离。见图11,测得其中一个谷地的长度为2597.470m。
谷地宽度测量:测量谷地两个崖壁的水平垂直距离。见图12,测得其中一个谷地宽度为348.829m,。
蚀去断面面积:丹霞地貌形成过程中被蚀去的断面面积。量取蚀去断面面积。
丹霞地貌侵蚀区长度:丹霞地貌侵蚀区最大长度,量取侵蚀区最远两个端点距离。
参数分析:基于Cyclone软件进行丹霞地貌的分析,包括丹霞地貌的分解、合并、放大或缩小,产生多个平面对丹霞地貌进行分段切割分析,丹霞山单体的对比、谷地的对比。
本发明的测量参数全面、准确度高,还可进行所需分析,为丹霞地貌的研究、保护、开发利用提供了永久的测量数据。
上述说明是针对本发明较佳可行实施例的详细说明,但实施例并非用以限定本发明的专利申请范围,凡本发明所提示的技术精神下所完成的同等变化或修饰变更,均应属于本发明所涵盖专利范围。
Claims (6)
1.一种丹霞地貌参数测量方法,其特征在于,包括以下步骤:
S1,丹霞地貌多源卫星遥感解译:收集待测丹霞地貌区域的资料,包括多源卫星数据、地形数据和地质资料,对不同来源的卫星遥感影像进行处理,综合解译待测丹霞地貌区域的典型丹霞地貌;
S2,无人机数据采集:采用无人机作为数据采集平台;首先室内确定丹霞地貌数据采集区域,并进行航线规划,采用井字形航线和立面航线规划,将设计好的航线导入无人机后,操控无人机按照规划的航线进行无人机倾斜摄影测量,智能采集数据;地形起伏大的区域采用仿地飞行方式获取测量区域的数据;部分地形复杂区域采用人工手动飞行补充航拍数据;
S3,航拍图像处理:采用图像处理软件对无人机拍摄的照片进行去云雾处理,保证照片清晰,采用Pix4D、Context Capture软件按照流程处理无人机采集的数据,选择2000国家大地坐标系,进行空三解算,生成丹霞地貌区域的正射影像、数字地表模型、三维模型和点云;
S4,丹霞地貌参数测量及分析:对步骤S3得到的丹霞地貌模型进行切割,识别并提取丹霞山单体,进行丹霞地貌参数测量。
2.根据权利要求1所述的一种丹霞地貌参数测量方法,其特征在于:所述步骤S1中的卫星数据来自landsat卫星、高分一号、高分二号、worldview卫星,卫星遥感影像处理和综合解译包含以下过程:
S101,遥感影像处理:对多源卫星数据进行全色与多光谱波段的配准、融合,正射校正,镶嵌;landsat卫星数据采用第7、第4和第1波段进行组合,能更好地区分丹霞地貌的岩性边界;landsat卫星数据采用第4、第3和第2波段组合,形成分辨率更高的融合影像模拟自然真彩色正射影像图,与地物真实色彩相近,能更好地进行丹霞地貌目视识别;其他类型卫星根据各自波段特征进行组合,用于丹霞地貌识别;
S102,丹霞地貌解译标志:根据landsat卫星741波段组合影像,划分丹霞地貌区域,裁切保留丹霞地貌区域影像;结合裁切的丹霞地貌区域地形数据,在ENVI软件中形成卫星正射影像的三维可视化场景,采用432波段组合形成自然真彩色影像,用于圈定丹霞山单体范围;再通过其他高分辨率卫星数据,高分一号、高分二号、worldview进一步确定位置和范围,解译出典型的丹霞山。
3.根据权利要求1所述的一种丹霞地貌参数测量方法,其特征在于,所述步骤S2中确定丹霞地貌数据采集区域的方法为:通过多源卫星遥感影像解译得到丹霞地貌的范围和位置,识别丹霞山,根据地面分辨率条件,结合奥维地图初步划定无人机数据采集区域的矢量,以kml文件形式保存。
4.根据权利要求3所述的一种丹霞地貌参数测量方法,其特征在于,所述步骤S2中航线规划的方法为:将确定的采集区域kml文件导入无人机飞控系统,根据分辨率及实地地形情况,确定飞行高度及重叠率;为保证测量的精准,比例尺设置为大于1:1000;在起伏较大的区域,即正负地形高差超过200m的区域,导入该区域地形数据,采用仿地飞行模式,确保数据采集的精准;
航线的航高设置依据以下公式:
航高=(相机镜头焦距×地面分辨率)/相机像元尺寸;
航高单位为米,地面分辨率单位为米,相机镜头焦距单位为毫米,相机像元尺寸为微米;
在划定飞行区域沿主轴方向规划,采用弓字形飞行路线,考虑到负地形采集全面,航向重叠率不低于80%,旁向重叠率不低于65%。
5.根据权利要求1所述的一种丹霞地貌参数测量方法,其特征在于,所述步骤S2中操控无人机智能采集数据时,按下以下条件进行:
飞行环境:在无风或微风、晴朗、能见度高的正午进行数据采集;
飞行方式:保证飞控系统与无人机连接,按照规划的航线飞行,在有CORS信号的区域,采用无人机飞控系统与CORS站连接,发射差分改正数据给无人机,保证照片POS信息的精准;如没有CORS信号,采用PPK方式进行飞行,以确保采集数据的高准确性;
数据的完整性:在高差较大区域,即正负地形高差超过200m的地形区域,因为负地形复杂,容易造成负地形或者部分丹霞山体局部区域出现数据采集重叠率低的情况,所以在完成自动航线数据采集后,根据地形情况,由人工操作无人机补充复杂负地形区域的数据采集;
数据检查:每个架次飞行完毕,导出影像和pos数据进行质量检查,主要检查照片影像质量、POS信息完备性、航向和旁向重叠率。
6.根据权利要求1所述的一种丹霞地貌参数测量方法,其特征在于:所述步骤S4丹霞地貌参数测量及分析具体方法为:
S401,模型切割:对步骤S3得到的丹霞地貌数字地表模型、三维模型、点云进行分析,将点云通过Cyclone软件进行过滤处理,除去冗余数据,在Cloud Compare点云处理软件中,将识别确定的丹霞区域切割,得到典型完整的丹霞地貌;
S402,丹霞山单体提取:在Cyclone软件中对切割的丹霞地貌进行分析,有多个丹霞山时,识别丹霞山单体,不连续的丹霞山独立算为丹霞山单体,采用软件中的编辑工具提取丹霞山单体;
S403,丹霞地貌参数选取:在众多丹霞地貌参数中选择最能反映丹霞地貌特征的参数,包括:丹霞山数量、谷地数量、丹霞山形态、坡度、丹霞山高度、丹霞山体积、岩层倾角、产状、走向、丹霞谷地宽度、丹霞地貌形成过程中被蚀去的断面面积、丹霞地貌侵蚀区长度;
S403,丹霞地貌参数测量:基于Cyclone和Cloud Compare软件进行丹霞地貌参数测量,对丹霞山和谷地数量、丹霞地貌系列参数进行测量,测量方法和要求如下:
丹霞山数量:在切割的丹霞地貌中,不连续的丹霞山独立算为丹霞山单体,从而计算丹霞山单体数量;
谷地数量:在切割的丹霞地貌中,独立完整的负地形为谷地,从而计算丹霞谷地数量;
丹霞山形态:水平面截取丹霞山单体底部,形成的边缘轮廓;
丹霞山坡度、产状和走向:采用Cloud Compare软件中的compass工具进行丹霞山坡度、产状和走向测量;
丹霞山高度:计算丹霞山最高处到底部水平面的垂直距离;
丹霞山体积:在丹霞山顶部形成水平面,计算该水平面与底部水平面之间的山体体积;
岩层倾角:选择要测量的岩层,测量其与水平面夹角;
丹霞谷地宽度:丹霞谷地两侧崖壁上缘之间的宽度;
蚀去断面面积:丹霞地貌形成过程中被蚀去的断面面积,量取蚀去断面面积;
丹霞地貌侵蚀区长度:丹霞地貌侵蚀区最大长度,量取侵蚀区最远两个端点距离;
参数分析:基于Cyclone软件进行丹霞地貌的分析,包括丹霞地貌的分解、合并、放大或缩小,产生多个平面对丹霞地貌进行分段切割分析,丹霞山单体的对比、谷地的对比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010252595.5A CN111444872B (zh) | 2020-03-31 | 2020-03-31 | 一种丹霞地貌参数测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010252595.5A CN111444872B (zh) | 2020-03-31 | 2020-03-31 | 一种丹霞地貌参数测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111444872A true CN111444872A (zh) | 2020-07-24 |
CN111444872B CN111444872B (zh) | 2023-11-24 |
Family
ID=71651100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010252595.5A Active CN111444872B (zh) | 2020-03-31 | 2020-03-31 | 一种丹霞地貌参数测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111444872B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112066952A (zh) * | 2020-08-05 | 2020-12-11 | 广州誉宸信息科技有限公司 | 基于无人机的路基高边坡防护检测方法和系统 |
CN112504236A (zh) * | 2020-12-12 | 2021-03-16 | 国网辽宁省电力有限公司铁岭供电公司 | 无人机航拍配电网台区测绘方法 |
CN113418448A (zh) * | 2021-06-23 | 2021-09-21 | 西北核技术研究所 | 一种破片分布检测系统和方法 |
CN114018205A (zh) * | 2021-11-18 | 2022-02-08 | 广东电网有限责任公司江门供电局 | 一种电力线路青赔面积测量方法及系统 |
CN115147729A (zh) * | 2022-06-27 | 2022-10-04 | 浙江大学 | 基于天空地数据和智能算法的固废堆填场地风险识别方法 |
CN116228782A (zh) * | 2022-12-22 | 2023-06-06 | 中国农业科学院农业信息研究所 | 基于无人机采集的麦田穗数计数方法和装置 |
CN116449859A (zh) * | 2023-04-28 | 2023-07-18 | 中交一公局第七工程有限公司 | 一种基于无人机技术的复杂山区便道规划方法及系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101606516B1 (ko) * | 2015-03-31 | 2016-03-28 | 국민대학교산학협력단 | 무인항공영상을 이용한 수목생장조사 시스템 및 그 방법 |
CN105865427A (zh) * | 2016-05-18 | 2016-08-17 | 三峡大学 | 一种基于小型无人机遥感的单体地质灾害应急调查方法 |
CN106774431A (zh) * | 2016-12-30 | 2017-05-31 | 深圳市九天创新科技有限责任公司 | 一种测绘无人机航线规划方法及装置 |
CN106969751A (zh) * | 2017-03-13 | 2017-07-21 | 西安科技大学 | 一种基于无人机遥感的采煤地表沉陷量监测计算的方法 |
RU2633642C1 (ru) * | 2016-06-08 | 2017-10-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Сибирский государственный университет геосистем и технологий" (СГУГиТ) | Способ получения, обработки, отображения и интерпретации геопространственных данных для геодезического мониторнига операционной обстановки паводковой ситуации с применением технологии дистанционного зондирования |
CA2957065A1 (en) * | 2017-02-06 | 2018-08-06 | David R. Lewis | Synthetically recreating the geostationary satellite orbital arc array with preferred elliptical orbit parameters |
CN108409410A (zh) * | 2018-03-23 | 2018-08-17 | 山西省农业科学院农业环境与资源研究所 | 一种用于果树复壮变绿的抑菌液及其制备方法 |
CN108701373A (zh) * | 2017-11-07 | 2018-10-23 | 深圳市大疆创新科技有限公司 | 基于无人机航拍的三维重建方法、系统及装置 |
CN109238239A (zh) * | 2018-09-12 | 2019-01-18 | 成都坤舆空间科技有限公司 | 基于航空摄影的数字测量三维建模方法 |
CN109697475A (zh) * | 2019-01-17 | 2019-04-30 | 中国地质大学(北京) | 一种湿地植被信息分析方法、遥感监测组件及监测方法 |
CN110132239A (zh) * | 2019-05-09 | 2019-08-16 | 苏州嘉奕晟中小企业科技咨询有限公司 | 基于cors平台的无人机测绘方法 |
CN110398246A (zh) * | 2019-07-15 | 2019-11-01 | 西安长庆科技工程有限责任公司 | 基于沙漠地区无人机航拍影像自动生成线划图的方法 |
CN110395390A (zh) * | 2019-07-17 | 2019-11-01 | 中国电建集团西北勘测设计研究院有限公司 | 一种多旋翼无人机免像控点三维建模与测图的装置及方法 |
CN110836661A (zh) * | 2019-12-04 | 2020-02-25 | 中国地质科学院岩溶地质研究所 | 一种天坑参数测量方法 |
-
2020
- 2020-03-31 CN CN202010252595.5A patent/CN111444872B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101606516B1 (ko) * | 2015-03-31 | 2016-03-28 | 국민대학교산학협력단 | 무인항공영상을 이용한 수목생장조사 시스템 및 그 방법 |
CN105865427A (zh) * | 2016-05-18 | 2016-08-17 | 三峡大学 | 一种基于小型无人机遥感的单体地质灾害应急调查方法 |
RU2633642C1 (ru) * | 2016-06-08 | 2017-10-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Сибирский государственный университет геосистем и технологий" (СГУГиТ) | Способ получения, обработки, отображения и интерпретации геопространственных данных для геодезического мониторнига операционной обстановки паводковой ситуации с применением технологии дистанционного зондирования |
CN106774431A (zh) * | 2016-12-30 | 2017-05-31 | 深圳市九天创新科技有限责任公司 | 一种测绘无人机航线规划方法及装置 |
CA2957065A1 (en) * | 2017-02-06 | 2018-08-06 | David R. Lewis | Synthetically recreating the geostationary satellite orbital arc array with preferred elliptical orbit parameters |
CN106969751A (zh) * | 2017-03-13 | 2017-07-21 | 西安科技大学 | 一种基于无人机遥感的采煤地表沉陷量监测计算的方法 |
CN108701373A (zh) * | 2017-11-07 | 2018-10-23 | 深圳市大疆创新科技有限公司 | 基于无人机航拍的三维重建方法、系统及装置 |
CN108409410A (zh) * | 2018-03-23 | 2018-08-17 | 山西省农业科学院农业环境与资源研究所 | 一种用于果树复壮变绿的抑菌液及其制备方法 |
CN109238239A (zh) * | 2018-09-12 | 2019-01-18 | 成都坤舆空间科技有限公司 | 基于航空摄影的数字测量三维建模方法 |
CN109697475A (zh) * | 2019-01-17 | 2019-04-30 | 中国地质大学(北京) | 一种湿地植被信息分析方法、遥感监测组件及监测方法 |
CN110132239A (zh) * | 2019-05-09 | 2019-08-16 | 苏州嘉奕晟中小企业科技咨询有限公司 | 基于cors平台的无人机测绘方法 |
CN110398246A (zh) * | 2019-07-15 | 2019-11-01 | 西安长庆科技工程有限责任公司 | 基于沙漠地区无人机航拍影像自动生成线划图的方法 |
CN110395390A (zh) * | 2019-07-17 | 2019-11-01 | 中国电建集团西北勘测设计研究院有限公司 | 一种多旋翼无人机免像控点三维建模与测图的装置及方法 |
CN110836661A (zh) * | 2019-12-04 | 2020-02-25 | 中国地质科学院岩溶地质研究所 | 一种天坑参数测量方法 |
Non-Patent Citations (2)
Title |
---|
ANSHUMAN BHARDWAJ ET AL: "UAV Imaging of a Martian Brine Analogue Environment in a Fluvio-Aeolian Setting", REMOTE SENSING, pages 1 - 32 * |
周福军: "高寒高海拔复杂艰险山区无人机勘察技术应用", 川藏铁路工程建造技术研讨会论文集, pages 2 - 6 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112066952A (zh) * | 2020-08-05 | 2020-12-11 | 广州誉宸信息科技有限公司 | 基于无人机的路基高边坡防护检测方法和系统 |
CN112504236A (zh) * | 2020-12-12 | 2021-03-16 | 国网辽宁省电力有限公司铁岭供电公司 | 无人机航拍配电网台区测绘方法 |
CN113418448A (zh) * | 2021-06-23 | 2021-09-21 | 西北核技术研究所 | 一种破片分布检测系统和方法 |
CN113418448B (zh) * | 2021-06-23 | 2022-03-22 | 西北核技术研究所 | 一种破片分布检测系统和方法 |
CN114018205A (zh) * | 2021-11-18 | 2022-02-08 | 广东电网有限责任公司江门供电局 | 一种电力线路青赔面积测量方法及系统 |
CN114018205B (zh) * | 2021-11-18 | 2023-06-30 | 广东电网有限责任公司江门供电局 | 一种电力线路青赔面积测量方法及系统 |
CN115147729A (zh) * | 2022-06-27 | 2022-10-04 | 浙江大学 | 基于天空地数据和智能算法的固废堆填场地风险识别方法 |
CN116228782A (zh) * | 2022-12-22 | 2023-06-06 | 中国农业科学院农业信息研究所 | 基于无人机采集的麦田穗数计数方法和装置 |
CN116228782B (zh) * | 2022-12-22 | 2024-01-12 | 中国农业科学院农业信息研究所 | 基于无人机采集的麦田穗数计数方法和装置 |
CN116449859A (zh) * | 2023-04-28 | 2023-07-18 | 中交一公局第七工程有限公司 | 一种基于无人机技术的复杂山区便道规划方法及系统 |
CN116449859B (zh) * | 2023-04-28 | 2024-03-26 | 中交一公局第七工程有限公司 | 一种基于无人机技术的复杂山区便道规划方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111444872B (zh) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111444872B (zh) | 一种丹霞地貌参数测量方法 | |
CN113034689B (zh) | 基于激光点云的地形三维模型及地形图构建方法和系统、存储介质 | |
Battulwar et al. | A state-of-the-art review of automated extraction of rock mass discontinuity characteristics using three-dimensional surface models | |
CN102074047B (zh) | 一种高精细城市三维建模方法 | |
CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
EP2356584B1 (en) | Method of generating a geodetic reference database product | |
CN111426303A (zh) | 一种岩溶坡立谷参数测量方法 | |
CN111540052B (zh) | 一种铁路沿线危岩落石的快速定位及三维重建方法 | |
CN112461205B (zh) | 基于无人机倾斜摄影测量的铁路既有线横断面制作方法 | |
CN112729130A (zh) | 卫星遥感测量树木冠层高度的方法 | |
JP2003344048A (ja) | 森林情報処理システム | |
CN114859374B (zh) | 基于无人机激光点云和影像融合的新建铁路交叉测量方法 | |
CN111189433A (zh) | 一种基于无人机航拍的岩溶峰林地貌参数测量方法 | |
CN110207676A (zh) | 一种田沟塘参数的获取方法及装置 | |
CN110836661A (zh) | 一种天坑参数测量方法 | |
Rebelo et al. | Building 3D city models: Testing and comparing Laser scanning and low-cost UAV data using FOSS technologies | |
CN110148218B (zh) | 一种大批量机载LiDAR点云数据整体优化的方法 | |
Çömert et al. | since 1961 | |
CN115601517A (zh) | 岩体结构面信息采集方法、装置、电子设备及存储介质 | |
CN112308907B (zh) | 用飞行器对坡体进行近景摄影测量的航线规划方法 | |
Bagnolo et al. | UAV-based photogrammetry for archaeological heritage site survey and 3D modeling of the sardus pater temple (Italy) | |
CN113514037A (zh) | 基于便携无人机摄影筛查的岩体露头测量方法 | |
CN117876612B (zh) | 一种基于无人机倾斜摄影的路线拟合设计方法 | |
Mora-Félix et al. | The use of RPAS for the development of land surface models for natural resources management: a review | |
CN118135137B (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 |