• ISSN 1673-5722
  • CN 11-5429/P

高起伏区无人机摄影测量精度评估−以西秦岭光盖山-迭山断裂为例

张波 高泽民 王爱国 刘小丰 郑龙 袁道阳

张波,高泽民,王爱国,刘小丰,郑龙,袁道阳,2022. 高起伏区无人机摄影测量精度评估−以西秦岭光盖山-迭山断裂为例. 震灾防御技术,17(2):326−339. doi:10.11899/zzfy20220213. doi: 10.11899/zzfy20220213
引用本文: 张波,高泽民,王爱国,刘小丰,郑龙,袁道阳,2022. 高起伏区无人机摄影测量精度评估−以西秦岭光盖山-迭山断裂为例. 震灾防御技术,17(2):326−339. doi:10.11899/zzfy20220213. doi: 10.11899/zzfy20220213
Zhang Bo, Gao Zemin, Wang Aiguo, Liu Xiaofeng, Zheng Long, Yuan Daoyang. Accuracy Assessment of UAV Photogrammetry in High-Relief Area−A Case Study from Guanggaishan-Dieshan Fault in West Qinling Mountain[J]. Technology for Earthquake Disaster Prevention, 2022, 17(2): 326-339. doi: 10.11899/zzfy20220213
Citation: Zhang Bo, Gao Zemin, Wang Aiguo, Liu Xiaofeng, Zheng Long, Yuan Daoyang. Accuracy Assessment of UAV Photogrammetry in High-Relief Area−A Case Study from Guanggaishan-Dieshan Fault in West Qinling Mountain[J]. Technology for Earthquake Disaster Prevention, 2022, 17(2): 326-339. doi: 10.11899/zzfy20220213

高起伏区无人机摄影测量精度评估−以西秦岭光盖山-迭山断裂为例

doi: 10.11899/zzfy20220213
基金项目: 国家自然基金青年基金(41602225);第二次青藏高原综合科学考察研究项目(2019QZKK0901);中国地震局地震预测研究所基本科研项目(2021IESLZ06、2020IESLZ02、2019IESLZ03、2018IESLZ02)
详细信息
    作者简介:

    张波,男,生于1986年。博士,副研究员。主要从事新生代构造与活动构造研究。E-mail:kjwxn999@163.com

    通讯作者:

    王爱国,男,生于1972年。研究员。主要从事地震地质、工程地震及数值模拟方面的研究。E-mail:waguo2008@163.com

Accuracy Assessment of UAV Photogrammetry in High-Relief Area−A Case Study from Guanggaishan-Dieshan Fault in West Qinling Mountain

  • 摘要: 经过近10年的迅速发展,无人机摄影测量已成为活动构造研究的常用方法之一。但对于无人机摄影测量的精度评估,尤其是高起伏地区的精度评估存在不足。为此,选择白龙江北岸光盖山-迭山断裂沿线的黑峪寺、化马村,开展无人机摄影测量,并构建正射影像(DOM)和数字地表模型(DSM),配合差分GPS测绘进行校正和精度验证。通过对比实测控制点和图像提取点分析点精度,通过对比实测剖面与提取剖面分析剖面精度。研究结果表明,未经控制点校正的图像提取点与实测点存在较大误差,水平误差为5~8 m,垂直误差为几十米至上百米,但通过少数控制点校正后,点精度可达20 cm以内;6条实测剖面与提取剖面(提取自控制点校正后的图像)平均垂直精度总体为分米级,即0.16~0.65 m,标准差为0.13~0.69 m,略低于低起伏区的精度,对于测量条件恶劣的高起伏区,该精度是可接受的;异常高的垂直误差常出现在地形突变、低矮植被密集、行走困难等测量条件不理想位置。图像控制点中心点的准确识别、提取剖面线的修正准确性等因素也会影响精度评估的可靠性。
  • 随着活动构造研究的精细定量化发展,无人机摄影测量已成为常规研究手段(毕海芸等,2017)。近10年来,随着无人机摄影测量在活动构造研究领域中的应用,提高了活动构造测量精细化程度和效率(Westoby等,2012王朋涛等,2016)。

    目前对无人机摄影测量精度评估的研究较少,艾明等(2018)利用无人机摄影测量获取茶卡盆地内部活动构造定量参数,并评估其精度,发现当控制点数量增加时,误差随之下降,但当控制点数量达12~15个后,图像精度提高不明显;当飞行高度为40~60 m时,图像精度最高,90%置信度的数据精度为10~15 cm,垂直精度为5~10 cm。杨海波等(2016)在河西走廊洪水坝河东侧阶地上开展了跨断裂无人机摄影测量,发现90%置信区间内实测数据与提取数据的高程差为10~15 cm。

    上述2个代表性实例主要基于地形相对平坦的低起伏区。对于高起伏区,无人机摄影测量的应用虽较多,但精度评估工作相对较少(张波等,2018)。高起伏区通常测量条件不理想、地形起伏度高、行走困难、植被密集,图像精度能否达到低起伏区的水平,亟待系统的精度评估。

    本文研究点位于白龙江流域(图1),地形起伏度超过3 000 m,发育了白龙江断裂和光盖山-迭山断裂,发生过多次大地震、滑坡、泥石流等地质灾害(陈洪凯等,1997袁道阳等,2007余斌等,2010梁学战等,2010俞晶星等,2012刘兴旺等,2015常直杨等,2015苏琦等,2016张波等,2018)。在该地区,传统测量手段(差分GPS、全站仪等)无法开展大面积的精细测绘工作,为此,本文对白龙江北岸光盖山-迭山断裂沿线的黑峪寺、化马村开展了无人机摄影测量工作,并结合差分GPS法进行精度评估。

    图 1  西秦岭白龙江流域地形地貌及测点分布
    Figure 1.  Topography and measure sites distribution in Bailongjiang watershed, West Qinling

    无人机摄影测量是基于移动摄影重建技术快速构建高分辨率数字地形的方法(Westoby等,2012毕海芸等,2017),测量原理如图2所示。首先获取大量一定重合度的照片,然后基于尺度不变特征变换算法(Lowe,2004Moreels等,2007)进行特征检测、跟踪与匹配,得到相机的空间位置及场景结构;再利用迭代光速平差方法(Harwin等,2012Mancini等,2013Bemis等,2014Javernick等,2014Lucieer等,2014)确定相机的空间位置及场景结构,得到稀疏点云框架;最后在影像之间根据逐像素搜索匹配和多视角立体摄影测量原理得到大量匹配点,建立更密集的点云,拼接后得到高分辨率的数字地形(Johnson等,2014毕海芸等,2017艾明,2018)。

    图 2  无人机摄影测量原理
    Figure 2.  Schematic diagram of UAV photogrammetry

    本文使用大疆Inspire 1二代专业级无人机获取照片(图3),该无人机配备SONY EXMOR相机,焦距为3.61 mm,有效像素为1 240万,传感器尺寸为1/2.3寸,视场角为94°。野外采集时间主要安排在冬季(黑峪寺测量时间为2018年12月2日,化马村测量时间为2016年11月26日),此时季节性乔木、灌木等已落叶,可最大程度降低植被影响。黑峪寺和化马村测量情况如表1所示。

    表 1  黑峪寺和化马村测量情况
    Table 1.  Measurements of Heiyusi and Huama site
    工作点飞行高度/
    m
    测量面积/
    km2
    图像重叠度/
    %
    照片数量/张控制点数目/个点云密度/
    点·m−2
    正射影像(DOM)
    分辨率/cm
    数字地表模型
    (DSM)分辨率/cm
    黑峪寺124.00.725约701 1601231.44.4617.8
    化马村98.60.385约70427747.03.6514.6
    下载: 导出CSV 
    | 显示表格
    图 3  黑峪寺无人机摄影测量构建的DOM和DSM
    Figure 3.  DOM and DSM of Heiyusi site constructed by UAV Photogrammetry

    在黑峪寺、化马村分别布设12、7个控制点,控制点为50 cm×50 cm正方形。控制点尽量保持在水平和垂直方向均匀分布。使用天宝GeoExplorer 6000系列差分GPS定位仪测量控制点坐标并校正无人机点云数据,其改正精度为厘米级。

    数据处理使用Photoscan软件(1.2.5版),处理流程完全参照艾明(2018)使用的流程,主要步骤包括照片对齐、导入控制点、加密点云、建立网格、生成纹理、构建DSM和DOM等。

    从未经控制点校正的DOM上提取控制点坐标,与实测坐标进行对比,二者的差异即为水平误差。从未经控制点校正的DSM上提取相应的高程,与实测高程进行对比,二者差异即为垂直误差。

    由于地形起伏不平,飞行高度会发生变化,高度大时采集的图像控制点不清晰。为确定图像控制点中心坐标,本文在DOM图像上大致勾出标靶点范围,然后在Arcmap里使用属性表计算标靶点范围中心点的经纬度坐标(Xcenter,Ycenter),再将经纬度坐标和差分GPS测量坐标同时添加到Arcmap中,新建线图层,用线连接实测点和计算点,在属性列表中计算线的长度,代表水平误差。同时,用经纬度坐标提取相应DSM高程,将提取到的高程和实测高程进行对比,二者差异为垂直误差。

    断层陡坎是最常见的断错地貌,对断层陡坎的精细刻画有助于活动构造的定量研究。 本文通过对比分析差分GPS实测剖面和从DSM上提取剖面的差异,研究高起伏区无人机摄影测量对地形剖面的显示特征。首先将差分GPS实测剖面投影到DOM和DSM上,由于DOM经过拼接校正,其上投影剖面和实测剖面位置可能存在一定偏差,此时需将投影到DOM上的剖面位置微调到实测位置;然后使用调整后的剖面点从DSM上提取三维坐标,形成剖面;最后与实测剖面进行对比。本研究在黑峪寺和化马村分别对3条剖面进行实测,将各自修正后的提取剖面与实测剖面进行对比,分析二者差异,从而分析高起伏区无人机摄影测量对地形的细节显示特征和偏差。

    黑峪寺位于由光盖山-迭山断裂控制的大型山间盆地内。在黑峪寺以北的盆地边缘跨断层开展无人机摄影测量,测量范围为长1.5 km、宽0.5 km的近似矩形。生成的DOM像元为4.46 cm,DSM像元为17.8 cm。由于图像边缘畸变严重,本文对图像边缘进行裁剪。黑峪寺无人机摄影测量构建的DOM和DSM如图3所示。

    2.1.1   点精度

    (1)控制点坐标对比——差分GPS实测控制点坐标与未校正DOM和DSM上提取点坐标对比

    从未校正的DOM和DSM上识别并提取控制点坐标,并与差分GPS实测控制点坐标进行对比,结果如表2所示。由表2可知,12个控制点的水平误差为5.64~7.76 m,标准差为0.53 m;垂直误差为111.38~114.79 m,标准差为1.08 m。未经控制点校正的DOM和DSM水平误差与垂直误差均较大,绝对精度低,尤其是垂直方向误差可达百米。水平方向和垂直方向的标准差较小,说明方向误差较集中。产生误差的原因可能是系统偏差,即GNSS接收机系统的定位差异。

    表 2  黑峪寺未校正图像提取控制点与差分GPS实测控制点水平误差和垂直误差
    Table 2.  Horizontal and vertical errors between uncorrected image extraction points and DGPS measured points
    控制点未校正图像提取控制点差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    1104.186 256 033.924 494 02 265.156104.186 241 133.924 436 32 152.306.54112.85
    2104.186 157 033.923 752 02 236.890104.186 139 133.923 690 02 125.357.08111.54
    3104.185 135 033.924 406 02 240.789104.185 117 633.924 343 72 128.937.1111.86
    4104.184 937 033.925 042 02 246.660104.184 919 933.924 979 42 133.697.12112.97
    5104.184 051 033.924 445 02 216.075104.184 030 333.924 377 22 104.697.76111.38
    6104.184 389 033.925 279 02 244.573104.184 371 833.925 216 12 131.537.16113.05
    7104.184 360 033.926 269 02 275.363104.184 346 333.926 211 82 160.586.47114.79
    8104.183 976 033.926 130 02 266.602104.183 961 133.926 071 92 152.416.59114.19
    9104.183 428 033.925 507 02 263.471104.183 411 833.925 446 62 150.936.86112.54
    10104.183 103 033.926 348 02 279.589104.183 087 833.926 292 72 165.806.29113.79
    11104.182 400 033.926 389 02 288.628104.182 384 833.926 333 72 175.406.29113.23
    12104.182 916 033.927 011 02 311.029104.182 903 333.926 961 32 196.415.64114.62
    下载: 导出CSV 
    | 显示表格

    (2)控制点坐标对比——差分GPS实测控制点坐标与从DOM和DSM(经6个控制点校正)上提取点坐标对比

    在12个控制点中随机选择6个点对DOM和DSM进行校正,另外6个点作为检验点。在校正后的DOM和DSM上识别6个检验点,提取其坐标信息,与差分GPS实测坐标进行对比,结果如表3所示,表中检验点编号与控制点编号一致。由表3可知,6个检验点实测坐标与提取坐标的水平误差为8~20 cm,标准差为4 cm;垂直误差为7~22 cm,标准差为5 cm。经6个控制点校正后,误差显著降低,水平方向和垂直方向误差均在22 cm以内。

    表 3  经6个控制点校正后提取的检验点坐标与实测坐标对比
    Table 3.  Comparison between six control points-corrected test points and measured points
    检验点从DOM和DSM提取坐标差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    2104.186 140 033.923 689 02 125.46104.186 139 133.923 690 02 125.350.150.11
    4104.184 920 033.924 981 02 133.82104.184 919 933.924 979 42 133.690.20.13
    6104.184 371 033.925 217 02 131.74104.184 371 833.925 216 12 131.530.140.22
    8104.183 960 033.926 072 02 152.54104.183 961 133.926 071 92 152.410.080.13
    11104.182 385 033.926 335 02 175.32104.182 384 833.926 333 72 175.400.170.07
    12104.182 903 033.926 960 02 196.22104.182 903 333.926 961 32 196.410.150.19
    下载: 导出CSV 
    | 显示表格

    综上所述,在高起伏区进行无人机摄影测量时,使用少量控制点进行图像校正,可显著提高水平和垂直方向绝对精度,总体精度可控制在20 cm以内。

    2.1.2   剖面精度

    本文在两级地貌面上利用差分GPS测量了3条地形剖面(图3(c)、图4)。剖面P1(图4(a))位于山前坡洪积台地上,剖面P2、P3(图4(b))位于冲沟阶地上。同时,在未校正DSM、6个控制点校正的DSM和12个控制点校正的DSM上分别提取相同位置的地形剖面。对比实测剖面和提取剖面(图57),得到以下结论:

    图 4  黑峪寺地形剖面分布
    Figure 4.  Distribution of topographical profiles at Heiyusi site
    图 5  黑峪寺P1实测剖面与提取剖面对比
    Figure 5.  Comparison of P1 in Heiyusi between measured profiles and extraction profiles
    图 6  黑峪寺P2实测剖面与提取剖面对比
    Figure 6.  Comparison of P2 in Heiyusi between measured profiles and extraction profiles
    图 7  黑峪寺P3实测剖面与提取剖面对比
    Figure 7.  Comparison of P3 in Heiyusi between measured profiles and extraction profiles

    (1)实测剖面和提取剖面形态一致,细节清晰完整,各级地形陡坎、梯田、缓坡等可被完整识别,肉眼可见的地形细节在各剖面均有显示。

    (2)未校正DSM上提取的剖面与实测剖面具有较大的垂向偏差,而校正后的DSM上提取的剖面与实测剖面重合度较好。通过分析3条剖面上所有点的提取高度与实测高度,发现6个控制点校正剖面与差分GPS剖面的平均垂直误差分别为0.17 m(剖面P1)、0.16 m(剖面P2)和0.33 m(剖面P3),标准差分别为0.11 m(剖面P1)、0.26 m(剖面P2)和0.17 m(剖面P3);12个控制点校正剖面与差分GPS剖面的平均垂直误差分别为0.20 m(剖面P1)、0.33 m(剖面P2)和0.36 m(剖面P3),标准差分别为0.13 m(剖面P1)、0.31 m(剖面P2)和0.16 m(剖面P3)。3条剖面的平均垂直误差和标准差均较小,说明提取剖面(校正后)可如实反映真实的地形。

    (3)当控制点数量由6个增至12个时,提取剖面的垂直误差未显著提升,说明在高地形起伏区开展摄影测量时,使用少数几个控制点校正图像可达到较高的精度,实现地形的高精度测绘。

    (4)垂直误差较大的部位往往是地形坡度陡变的部位,特别是在剖面P2的地形陡坎处,垂直误差由约0.2 m升至1.6 m(图6(c))。在地形陡变(陡坎改造成梯田)部位,低矮植被密集,行走困难,实际测量时需绕行,测量条件的不理想可能造成了异常高的垂直误差。

    化马村位于宕昌岷江左岸,多级大型洪积扇被断层断错,形成一系列左旋兼具逆冲的断错地貌,断层露头上的断层擦痕侧伏角为20°,显示断层以水平分量为主,兼具倾滑运动(张波等,2018)。化马村无人机摄影测量范围为不规则形状,面积为0.385 km2。无人机平均飞行高度为98.6 m,共拍摄427张照片。生成的DOM像元为3.65 cm,DSM像元为14.6 cm。分析前裁剪了边缘畸变严重部分。

    2.2.1   差分GPS实测控制点坐标与未校正DOM和DSM上提取点坐标对比

    由于无人机与差分GPS定位系统存在差异,在未校正DOM和DSM上提取控制点和实测控制点存在较大的系统误差,水平误差为7.65~7.97 m,标准差为0.1 m;垂直误差为35.79~37.74 m,标准差为0.69 m(表4)。方向误差的标准差均较小,说明方向误差较集中。化马村垂直误差(35.79~37.74 m)远小于黑峪寺(111.38~114.79 m),这可能与无人机定位系统的不稳定、测量时间间隔较长(约2年)等因素有关。

    表 4  化马村未校正图像提取控制点与差分GPS实测控制点水平误差和垂直误差
    Table 4.  Directional error at Huama site between extraction points from uncorrected images and measured points
    控制点未校正图像提取控制点差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    1104.540 325 033.741 176 01 556.67104.540 28833.741 114 41 594.247.6537.57
    2104.541 607 033.741 010 01 532.95104.541 56733.740 947 41 569.417.9036.46
    3104.540 270 033.740 160 01 512.82104.540 22633.740 100 61 550.567.7537.74
    4104.540 635 033.739 285 01 467.24104.540 58533.739 227 61 503.287.8636.04
    5104.540 810 033.738 124 01 469.00104.540 76133.738 065 91 506.097.8937.09
    6104.541 360 033.739 435 01 500.42104.541 31533.739 375 31 537.037.8436.60
    7104.542 265 033.740 881 01 536.25104.542 22433.740 818 01 572.047.9735.79
    下载: 导出CSV 
    | 显示表格
    2.2.2   剖面精度

    使用差分GPS在化马村测量了3条剖面(P1~P3),剖面P1和P3位于二级洪积台地上,剖面P2位于河道上(图8(b)图9)。同时,从未校正DSM、7个控制点校正DSM上的相同位置提取地形剖面,对比分析实测剖面和提取剖面,发现以下特征(图10图12):

    图 8  化马村无人机摄影测量
    Figure 8.  UAV Photogrammetry of Huama site
    图 9  化马村地形剖面分布
    Figure 9.  Distribution of Topographical profiles in Huama site
    图 10  化马村P1实测剖面与提取剖面对比
    Figure 10.  Comparison of P1 in Huama between measured profiles and extraction profiles
    图 11  化马村P2实测剖面与提取剖面对比
    Figure 11.  Comparison of P2 in Huama between measured profiles and extraction profiles
    图 12  化马村P3实测剖面与提取剖面对比
    Figure 12.  Comparison of P3 in Huama between measured profiles and extraction profiles

    (1)实测剖面与提取剖面总体形态、细节特征具有较好的一致性。剖面P1、P3(图10、图12)上不同剖面对多级梯田和地形陡坡的显示具有较好的同步性。

    (2)未校正剖面与实测剖面具有较大的垂直误差,而使用7个控制点校正后的剖面和实测剖面重合较好。校正剖面的平均垂直误差分别为0.25 m(剖面P1)、0.65 m(剖面P2)和0.30 m(剖面P3),垂直误差的标准差分别为0.22 m(剖面P1)、0.69 m(剖面P2)和0.34 m(剖面P3)。剖面P2平均垂直误差达0.65 m,标准差达0.69 m,显著优于未校正图像的绝对误差,对高地形陡坎的高度测量影响较小。

    (3)垂直误差异常高的部位往往位于地形坡折处或测量条件不理想的部位,如剖面P1上的横坐标150 m附近(图10),野外实测时绕道、植被等不利因素导致了异常高的垂直误差。

    (4)化马村摄影测量精度评估得到的结论与黑峪寺类似,即提取剖面与实测剖面具有相似的剖面形态、完整的细节展现。在高起伏区,少量控制点校正后的剖面平均垂直精度可达分米级,可较准确地反映真实地形地貌。但当垂直误差达分米级后,再增加控制点数量难以提高绝对精度。在地形坡度陡变(如梯田、地形陡坎)、植被密集、行走困难等测量条件不理想部位,易出现异常高的垂直误差。

    本文介绍了在高起伏区开展无人机摄影测量的2个实例,发现经少数控制点校正后,DOM和DSM可达较高的绝对精度和相对精度,满足高精度测量需求。无人机摄影测量在高地形起伏区仍适用,基于此,研究人员可快速获取此类地区大面积高精度数字地形,满足多行业工作需要。此外,需注意以下方面:

    (1)在黑峪寺图像校正任务里,使用12个控制点得到的剖面垂直精度未得到提高,反而略有下降。这与艾明(2018)在茶卡盆地(低起伏区)开展的摄影测量精度评估略有差异,该研究发现绝对精度随控制点的增加逐渐提高,超过12~15个控制点后精度提高不明显。

    (2)本文得到的最高单点垂直精度为0.08 m,最高的剖面平均垂直精度为0.16 m(黑峪寺剖面P2)。最高的平均剖面垂直精度与杨海波等(2016)在低起伏的河西走廊地区相当,但本文得到的6条剖面平均垂直精度(0.17、0.16、0.33、0.25、0.65、0.30 m)和对应的标准差(0.11、0.26、0.17、0.22、0.69、0.34 m)低于已有研究得到的垂直精度。剖面上异常高的误差通常与地形陡变同时存在,与魏占玉等(2015)的认识一致,即平缓地形的精度高于陡峭地形的精度。虽误差略高,但对于地形复杂、植被密集、自然和人为改造强烈的高起伏区而言,仍为可接受的误差范围。在如此复杂的多因素影响环境中,难以消除这些误差。另外,本文未发现杨海波等(2016)和Johnson等(2014)所提到的倾斜变形。

    (3)在高起伏区,影响误差的因素是多方面的。本文将部分细节应用到评估工作中,可能会有利于精度的提高。由于高起伏区无法在同一高度按设计航线飞行,高度过高时,DOM上标靶的中心点有时无法准确识别。对于不清楚的控制点图像,本文首先通过Arcmap编辑工具大致勾画标靶范围,然后计算标靶范围中心点的经纬度坐标,最后根据经纬度坐标提取对应数字地表模型高度。通过计算的中心点可能较肉眼目视估计的中心点准确。

    (4)提取剖面与实测剖面进行对比时,差分GPS剖面投影到DSM上的测线可能与实际测线有偏差,此时需结合DOM、Google Earth影像上加载的野外测量轨迹、地形环境等特征,将测线修正到实际测量位置,再根据修正测线在DSM上提取剖面。该方法将有效提高实测剖面和提取剖面的重合度,提高垂直精度分析的可靠性。同时,由于地形过于复杂,提取剖面使用的测线可能难以精确修正到实际测线位置,从而形成新的误差,如化马村剖面P2可能受到了测线修正不够精确的影响。

    (5)高植被对重建图像的影响较大。本文研究点尽量避开了高大植被密集区,测量时间选择在冬季落叶期,最大程度地减小了植被的影响。部分低植被区可能是导致剖面上异常高的垂直误差因素之一。总体来看,本文较好地评估了高起伏区无人机摄影测量的点精度和剖面精度,验证了无人机摄影测量在高起伏地区的可用性和实用性,并分析了可能影响精度的原因。

    本文通过光盖山-迭山断裂沿线的2个研究点,分析了高起伏区无人机摄影测量精度,得到以下结论。

    (1)经少数控制点校正后,DOM和DSM的点精度可达20 cm以内,符合高精度测量要求。

    (2)6条实测剖面与从DSM上提取剖面的平均垂直精度为0.16~0.65 m,标准差为0.13~0.69 m。在测量条件恶劣的高起伏地区,该平均垂直精度和标准差是可接受的。

    (3)异常高的垂直误差常出现在地形突变、低矮植被密集、行走困难等测量条件不理想的位置。同时,精度分析时图像控制点中心点的准确识别、提取剖面线的修正准确性等因素也会影响精度评估的可靠性。

    致谢 感谢中国地震局兰州地震研究所张晴毅等硕士研究生在野外工作中的协助,感谢审稿人从专业测绘角度提出的宝贵意见。

  • 图  1  西秦岭白龙江流域地形地貌及测点分布

    Figure  1.  Topography and measure sites distribution in Bailongjiang watershed, West Qinling

    图  2  无人机摄影测量原理

    Figure  2.  Schematic diagram of UAV photogrammetry

    图  3  黑峪寺无人机摄影测量构建的DOM和DSM

    Figure  3.  DOM and DSM of Heiyusi site constructed by UAV Photogrammetry

    图  4  黑峪寺地形剖面分布

    Figure  4.  Distribution of topographical profiles at Heiyusi site

    图  5  黑峪寺P1实测剖面与提取剖面对比

    Figure  5.  Comparison of P1 in Heiyusi between measured profiles and extraction profiles

    图  6  黑峪寺P2实测剖面与提取剖面对比

    Figure  6.  Comparison of P2 in Heiyusi between measured profiles and extraction profiles

    图  7  黑峪寺P3实测剖面与提取剖面对比

    Figure  7.  Comparison of P3 in Heiyusi between measured profiles and extraction profiles

    图  8  化马村无人机摄影测量

    Figure  8.  UAV Photogrammetry of Huama site

    图  9  化马村地形剖面分布

    Figure  9.  Distribution of Topographical profiles in Huama site

    图  10  化马村P1实测剖面与提取剖面对比

    Figure  10.  Comparison of P1 in Huama between measured profiles and extraction profiles

    图  11  化马村P2实测剖面与提取剖面对比

    Figure  11.  Comparison of P2 in Huama between measured profiles and extraction profiles

    图  12  化马村P3实测剖面与提取剖面对比

    Figure  12.  Comparison of P3 in Huama between measured profiles and extraction profiles

    表  1  黑峪寺和化马村测量情况

    Table  1.   Measurements of Heiyusi and Huama site

    工作点飞行高度/
    m
    测量面积/
    km2
    图像重叠度/
    %
    照片数量/张控制点数目/个点云密度/
    点·m−2
    正射影像(DOM)
    分辨率/cm
    数字地表模型
    (DSM)分辨率/cm
    黑峪寺124.00.725约701 1601231.44.4617.8
    化马村98.60.385约70427747.03.6514.6
    下载: 导出CSV

    表  2  黑峪寺未校正图像提取控制点与差分GPS实测控制点水平误差和垂直误差

    Table  2.   Horizontal and vertical errors between uncorrected image extraction points and DGPS measured points

    控制点未校正图像提取控制点差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    1104.186 256 033.924 494 02 265.156104.186 241 133.924 436 32 152.306.54112.85
    2104.186 157 033.923 752 02 236.890104.186 139 133.923 690 02 125.357.08111.54
    3104.185 135 033.924 406 02 240.789104.185 117 633.924 343 72 128.937.1111.86
    4104.184 937 033.925 042 02 246.660104.184 919 933.924 979 42 133.697.12112.97
    5104.184 051 033.924 445 02 216.075104.184 030 333.924 377 22 104.697.76111.38
    6104.184 389 033.925 279 02 244.573104.184 371 833.925 216 12 131.537.16113.05
    7104.184 360 033.926 269 02 275.363104.184 346 333.926 211 82 160.586.47114.79
    8104.183 976 033.926 130 02 266.602104.183 961 133.926 071 92 152.416.59114.19
    9104.183 428 033.925 507 02 263.471104.183 411 833.925 446 62 150.936.86112.54
    10104.183 103 033.926 348 02 279.589104.183 087 833.926 292 72 165.806.29113.79
    11104.182 400 033.926 389 02 288.628104.182 384 833.926 333 72 175.406.29113.23
    12104.182 916 033.927 011 02 311.029104.182 903 333.926 961 32 196.415.64114.62
    下载: 导出CSV

    表  3  经6个控制点校正后提取的检验点坐标与实测坐标对比

    Table  3.   Comparison between six control points-corrected test points and measured points

    检验点从DOM和DSM提取坐标差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    2104.186 140 033.923 689 02 125.46104.186 139 133.923 690 02 125.350.150.11
    4104.184 920 033.924 981 02 133.82104.184 919 933.924 979 42 133.690.20.13
    6104.184 371 033.925 217 02 131.74104.184 371 833.925 216 12 131.530.140.22
    8104.183 960 033.926 072 02 152.54104.183 961 133.926 071 92 152.410.080.13
    11104.182 385 033.926 335 02 175.32104.182 384 833.926 333 72 175.400.170.07
    12104.182 903 033.926 960 02 196.22104.182 903 333.926 961 32 196.410.150.19
    下载: 导出CSV

    表  4  化马村未校正图像提取控制点与差分GPS实测控制点水平误差和垂直误差

    Table  4.   Directional error at Huama site between extraction points from uncorrected images and measured points

    控制点未校正图像提取控制点差分GPS实测控制点水平误差/m垂直误差/m
    经度/°纬度/°高度/m经度/°纬度/°高度/m
    1104.540 325 033.741 176 01 556.67104.540 28833.741 114 41 594.247.6537.57
    2104.541 607 033.741 010 01 532.95104.541 56733.740 947 41 569.417.9036.46
    3104.540 270 033.740 160 01 512.82104.540 22633.740 100 61 550.567.7537.74
    4104.540 635 033.739 285 01 467.24104.540 58533.739 227 61 503.287.8636.04
    5104.540 810 033.738 124 01 469.00104.540 76133.738 065 91 506.097.8937.09
    6104.541 360 033.739 435 01 500.42104.541 31533.739 375 31 537.037.8436.60
    7104.542 265 033.740 881 01 536.25104.542 22433.740 818 01 572.047.9735.79
    下载: 导出CSV
  • [1] 艾明, 2018. 高精度摄影测量方法在活动构造定量参数获取中的应用——以茶卡盆地两侧断裂为例. 北京: 中国地震局地质研究所.

    Ai M., 2018. Application of high-precision photogrammetry method in obtaining quantitative parameters of active tectonics-A case study of the faults on the two sides of the Caka Basin. Beijing: Institute of Geology, China Earthquake Administrator. (in Chinese)
    [2] 艾明, 毕海芸, 郑文俊等, 2018. 利用无人机摄影测量技术提取活动构造定量参数. 地震地质, 40(6): 1276—1293

    Ai M. , Bi H. Y. , Zheng W. J. , et al. , 2018. Using unmanned aerial vehicle photogrammetry technology to obtain quantitative parameters of active tectonics. Seismology and Geology, 40(6): 1276—1293. (in Chinese)
    [3] 毕海芸, 郑文俊, 曾江源等, 2017. SfM摄影测量方法在活动构造定量研究中的应用. 地震地质, 39(4): 656—674 doi: 10.3969/j.issn.0253-4967.2017.04.003

    Bi H. Y. , Zheng W. J. , Zeng J. Y. , et al. , 2017. Application of SfM photogrammetry method to the quantitative study of active tectonics. Seismology and Geology, 39(4): 656—674. (in Chinese) doi: 10.3969/j.issn.0253-4967.2017.04.003
    [4] 常直杨, 王建, 白世彪等, 2015. 白龙江流域河流纵剖面与基岩侵蚀模型特征. 山地学报, 33(2): 183—190

    Chang Z. Y. , Wang J. , Bai S. B. , et al. , 2015. Morpho-tectonic analysis of the Bailongjiang drainage basin. Mountain Research, 33(2): 183—190. (in Chinese)
    [5] 陈洪凯, 李吉均. 1997. 白龙江流域第四纪以来地貌发育基本模式研究. 重庆交通学院学报, 16(1): 15—20

    Chen H. K., Li J. J., 1997. General approach on geomorphologic evolution in Bailongjiang basin since Quaternary. Journal of Chongqing Jiaotong Institute, 16(1): 15—20. (in Chinese)
    [6] 梁学战, 陈洪凯, 唐红梅, 2010. 青藏高原东部边缘地区第四纪泥石流发育特性研究. 重庆交通大学学报(自然科学版), 29(6): 978—983

    Liang X. Z. , Chen H. K. , Tang H. M. , 2010. Studies on debris flow development properties in Quaternary on the eastern edge of Qinghai-Tibet Plateau. Journal of Chongqing Jiaotong University (Natural Sciences), 29(6): 978—983. (in Chinese)
    [7] 刘兴旺, 袁道阳, 邵延秀等, 2015. 甘肃迭部—白龙江南支断裂中东段晚第四纪构造活动特征. 地球科学与环境学报, 37(6): 111—119 doi: 10.3969/j.issn.1672-6561.2015.06.010

    Liu X. W. , Yuan D. Y. , Shao Y. X. , et al. , 2015. Characteristics of late Quaternary tectonic activity in the middle-eastern segment of the southern branch of Diebu-Bailongjiang fault, Gansu. Journal of Earth Sciences and Environment, 37(6): 111—119. (in Chinese) doi: 10.3969/j.issn.1672-6561.2015.06.010
    [8] 苏琦, 梁明剑, 袁道阳等, 2016. 白龙江流域构造地貌特征及其对滑坡泥石流灾害的控制作用. 地球科学, 41(10): 1758—1770

    Su Q. , Liang M. J. , Yuan D. Y. , et al. , 2016. Geomorphic features of the Bailongjiang river drainage basin and its relationship with geological disaster. Earth Science, 41(10): 1758—1770. (in Chinese)
    [9] 王朋涛, 邵延秀, 张会平等, 2016. sUAV摄影技术在活动构造研究中的应用——以海原断裂骟马沟为例. 第四纪研究, 36(2): 433—442 doi: 10.11928/j.issn.1001-7410.2016.02.18

    Wang P. T. , Shao Y. X. , Zhang H. P. , et al. , 2016. The application of sUAV Photogrammetry in active tectonics: Shanmagou site of Haiyuan fault, for example. Quaternary Sciences, 36(2): 433—442. (in Chinese) doi: 10.11928/j.issn.1001-7410.2016.02.18
    [10] 魏占玉, Ramon A, 何宏林等, 2015. 基于SfM方法的高密度点云数据生成及精度分析. 地震地质, 37(2): 636—648 doi: 10.3969/j.issn.0253-4967.2015.02.024

    Wei Z. Y. , Ramon A. , He H. L. , et al. , 2015. Accuracy analysis of terrain point cloud acquired by “Structure from Motion” using aerial photos. Seismology and Geology, 37(2): 636—648. (in Chinese) doi: 10.3969/j.issn.0253-4967.2015.02.024
    [11] 杨海波, 杨晓平, 黄雄南等, 2016. 移动摄影测量数据与差分GPS数据的对比分析——以祁连山北麓洪水坝河东岸断层陡坎为例. 地震地质, 38(4): 1030—1046 doi: 10.3969/j.issn.0253-4967.2016.04.018

    Yang H. B. , Yang X. P. , Huang X. N. , et al. , 2016. Data comparative analysis between SfM data and DGPS data: a case study from fault scarp in the east bank of Hongshuiba river, northern margin of the Qilian Shan. Seismology and Geology, 38(4): 1030—1046. (in Chinese) doi: 10.3969/j.issn.0253-4967.2016.04.018
    [12] 余斌, 杨永红, 苏永超等, 2010. 甘肃省舟曲8.7特大泥石流调查研究. 工程地质学报, 18(4): 437—444 doi: 10.3969/j.issn.1004-9665.2010.04.001

    Yu B. , Yang Y. H. , Su Y. C. , et al. , 2010. Research on the giant debris flow hazards in Zhouqu county, Gansu Province on August 7, 2010. Journal of Engineering Geology, 18(4): 437—444. (in Chinese) doi: 10.3969/j.issn.1004-9665.2010.04.001
    [13] 俞晶星, 郑文俊, 袁道阳等, 2012. 西秦岭西段光盖山-迭山断裂带坪定-化马断裂的新活动性与滑动速率. 第四纪研究, 32(5): 957—967 doi: 10.3969/j.issn.1001-7410.2012.05.13

    Yu J. X. , Zheng W. J. , Yuan D. Y. , et al. , 2012. Late Quaternary active characteristics and slip-rate of Pingding-Huama fault, the eastern segment of Guanggaishan-Dieshan fault zone (West Qinling Mountain). Quaternary Sciences, 32(5): 957—967. (in Chinese) doi: 10.3969/j.issn.1001-7410.2012.05.13
    [14] 袁道阳, 雷中生, 何文贵等, 2007. 公元前186年甘肃武都地震考证与发震构造探讨. 地震学报, 29(6): 654—663

    Yuan D. Y., Lei Z. S., He W. G., et al., 2007. Textual research of Wudu earthquake in 186 B. C. in Gansu Province, China and discussion on its causative structure. Acta Seismologica Sinica, 29(6): 654—663. (in Chinese)
    [15] 张波, 王爱国, 袁道阳等, 2018. 基于多源遥感解译和野外验证的断裂几何展布——以西秦岭光盖山-迭山南麓断裂为例. 地震地质, 40(5): 1018—1039 doi: 10.3969/j.issn.0253-4967.2018.05.005

    Zhang B. , Wang A. G. , Yuan D. Y. , et al. , 2018. Fault geometry defined by multiple remote sensing images interpretation and field verification: A case study from southern Guanggaishan-Dieshan fault, Western Qinling. Seismology and Geology, 40(5): 1018—1039. (in Chinese) doi: 10.3969/j.issn.0253-4967.2018.05.005
    [16] Bemis S. P. , Micklethwaite S. , Turner D. , et al. , 2014. Ground-based and UAV-based photogrammetry: a multi-scale, high-resolution mapping tool for structural geology and paleoseismology. Journal of Structural Geology, 69: 163—178. doi: 10.1016/j.jsg.2014.10.007
    [17] Harwin S. , Lucieer A. , 2012. Assessing the accuracy of georeferenced point clouds produced via multi-view stereopsis from unmanned aerial vehicle (UAV) imagery. Remote Sensing, 4(6): 1573—1599. doi: 10.3390/rs4061573
    [18] Javernick L. , Brasington J. , Caruso B. , 2014. Modeling the topography of shallow braided rivers using structure-from-motion photogrammetry. Geomorphology, 213: 166—182. doi: 10.1016/j.geomorph.2014.01.006
    [19] Johnson K. , Nissen E. , Saripalli S. , et al. , 2014. Rapid mapping of ultrafine fault zone topography with structure from motion. Geosphere, 10(5): 969—986. doi: 10.1130/GES01017.1
    [20] Lowe D. G. , 2004. Distinctive image features from scale-invariant keypoints. International Journal of Computer Vision, 60(2): 91—110. doi: 10.1023/B:VISI.0000029664.99615.94
    [21] Lucieer A. , De Jong S. M. , Turner D. , 2014. Mapping landslide displacements using structure from motion (SfM) and image correlation of multi-temporal UAV photography. Progress in Physical Geography: Earth and Environment, 38(1): 97—116. doi: 10.1177/0309133313515293
    [22] Mancini F. , Dubbini M. , Gattelli M. , et al. , 2013. Using unmanned aerial vehicles (UAV) for high-resolution reconstruction of topography: the structure from motion approach on coastal environments. Remote Sensing, 5(12): 6880—6898. doi: 10.3390/rs5126880
    [23] Moreels P. , Perona P. , 2007. Evaluation of features detectors and descriptors based on 3D objects. International Journal of Computer Vision, 73(3): 263—284. doi: 10.1007/s11263-006-9967-1
    [24] Westoby M. J. , Brasington J. , Glasser N. F. , et al. , 2012. ‘Structure-from-Motion’ photogrammetry: a low-cost, effective tool for geoscience applications. Geomorphology, 179: 300—314. doi: 10.1016/j.geomorph.2012.08.021
  • 期刊类型引用(4)

    1. 吴海波. 基于无人机技术的河道划界及精度分析. 地下水. 2024(01): 247-249 . 百度学术
    2. 张艳婷. 无人机航测技术在水利工程河道泥沙沉积量监测中的应用. 水利技术监督. 2023(08): 31-35+99 . 百度学术
    3. 王辽,谢虹,朱俊文,薛善余,张正模,张波. 三种断层位移测量工具对比分析. 地震工程学报. 2023(05): 1133-1144 . 百度学术
    4. 彭虎. 无人机航摄在水电工程大比例尺地图测绘应用研究. 水利科学与寒区工程. 2023(11): 9-12 . 百度学术

    其他类型引用(0)

  • 加载中
图(12) / 表(4)
计量
  • 文章访问数:  166
  • HTML全文浏览量:  68
  • PDF下载量:  16
  • 被引次数: 4
出版历程
  • 收稿日期:  2019-07-19
  • 刊出日期:  2022-06-30

目录

/

返回文章
返回