欢迎访问陕西师范大学学报(自然科学版)官方网站！

陕西师范大学学报(自然科学版) >

2023 , Vol. 51 >Issue 6: 25 - 36

DOI: https://doi.org/10.15983/j.cnki.jsnu.2023127

黄河流域生态保护与资源可持续利用专题 主持人：周正朝

SfM摄影测量法对黄土高原典型切沟的测量精度评价

  • 陈薪伊 1, 2 ,
  • 史扬子 1, 2 ,
  • 杨扬 , 1, 2, * ,
  • 黄婷婷 3 ,
  • 刘瑛娜 1 ,
  • 张光辉 1, 2
展开
  • 1 北京师范大学 地理科学学部 地表过程与资源生态国家重点实验室,北京 100875
  • 2 北京师范大学 地理科学学部 地理学院,北京 100875
  • 3 水利部 水土保持监测中心,北京 100053
*杨扬,女,副教授,博士生导师,主要从事土壤侵蚀与水土保持研究。E-mail:

Copy editor: 程琴娟

收稿日期: 2022-11-29

  网络出版日期: 2023-11-21

基金资助

国家自然科学基金(42130701)

国家自然科学基金(41601277)

收起

Accuracy evaluation of structure from motion (SfM) photogrammetry on the measurement of typical gullies in the Loess Plateau

  • CHEN Xinyi 1, 2 ,
  • SHI Yangzi 1, 2 ,
  • YANG Yang , 1, 2, * ,
  • HUANG Tingting 3 ,
  • LIU Yingna 1 ,
  • ZHANG Guanghui 1, 2
Expand
  • 1 State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
  • 2 School of Geography, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
  • 3 Monitoring Center of Soil and Water Conservation, Ministry of Water Resources, Beijing 100053, China

Received date: 2022-11-29

  Online published: 2023-11-21

Fold

摘要

在黄土高原六道沟小流域选择7条代表性切沟,以激光扫描法(laser scanning,LS)和实时动态差分(real-time kinematic,RTK)GPS的测量结果为基准,评估运动恢复结构(structure from motion,SfM)摄影测量法对切沟及切沟侵蚀的测量精度。与LS和RTK GPS相比,SfM摄影测量法所测各切沟DEM总体差异不大,平均高程差值为-0.162~0.436 m;SfM摄影测量法对切沟长度、表面宽度、周长和面积的测量精度较高,相对误差基本不超过±7%;所测深度和体积的相对误差较大,介于-37.4%~18.0%。与RTK GPS所测切沟侵蚀相比,SfM摄影测量法测得的切沟侵蚀空间分布和形态参数变化均十分接近。研究结果证实,SfM摄影测量法可作为一种高精度测量方法应用于黄土高原的切沟和切沟侵蚀监测。

关键词： 切沟; 切沟侵蚀; SfM摄影测量法; 激光扫描法; 实时动态差分GPS; 黄土高原

本文引用格式

陈薪伊 , 史扬子 , 杨扬 , 黄婷婷 , 刘瑛娜 , 张光辉 . SfM摄影测量法对黄土高原典型切沟的测量精度评价[J]. 陕西师范大学学报(自然科学版), 2023 , 51(6) : 25 -36 . DOI: 10.15983/j.cnki.jsnu.2023127

Abstract

The structure from motion (SfM) photogrammetry was used to measure seven typical gullies in the Liudaogou catchment located in the Loess Plateau, and the accuracy in quantifying gully and gully erosion was evaluated via comparing the results to those acquired with the laser scanning (LS) and real-time kinematic (RTK) GPS. The DEM measured with SfM photogrammetry for each gully were similar with those obtained using LS and RTK GPS, and the average elevation differences ranged from -0.162 to 0.436 m when subtracting the former from the latter. SfM photogrammetry exhibited high precisions in quantifying the length, surface width, perimeter and area of each gully, and the corresponding relative errors barely exceeded ±7%. In contrast, the measurement errors of SfM photogrammetry were relatively greater for gully depth and volume, which ranged from -37.4% to 18.0%. Regarding gully erosion monitoring, both the spatial distribution of soil erosion and gully morphological parameter changes obtained by SfM photogrammetry were similar with those by RTK GPS. These results demonstrate that SfM photogrammetry can serve as a high-precision method for gully and gully erosion monitoring in the Loess Plateau.

Key words： gully; gully erosion; structure from motion (SfM) photogrammetry; laser scanning; real-time kinematic (RTK) GPS; the Loess Plateau

切沟是指普通农具无法横跨耕作的侵蚀沟,通常切破犁底层,深入母质,宽度和深度均超过50 cm[1-2]。切沟侵蚀是指切沟发生发展过程中导致的侵蚀。作为水力侵蚀最严重的阶段,切沟侵蚀是全球河流泥沙的主要来源之一,其产沙量可占流域产沙总量的10%~94%[3]。在我国黄土高原地区,切沟侵蚀产沙量可占小流域总产沙量的50%以上,部分地区甚至可以达到80%以上[4-5]。选择准确、高效的测量方法对切沟及切沟侵蚀进行监测,对于揭示切沟侵蚀时空分布规律、阐明切沟侵蚀机理、构建和完善切沟侵蚀模型具有重要意义。
近年来,以激光扫描法(laser scanning,LS)、实时动态差分(real-time kinematic,RTK)GPS和摄影测量法为典型代表的现代土壤侵蚀测量技术开始逐渐应用于切沟及切沟侵蚀监测,均取得了良好效果[6⇓⇓-9]。LS利用激光扫描设备发射激光,根据其返回时间推算所测地物的三维坐标,进而生成DEM。该方法测量精度极高[10],结果常被用作验证其他方法的标准数据[11];然而,LS所用设备昂贵,对操作人员的专业要求较高,目前应用十分有限[12]。RTK GPS是随GPS定位精度不断突破发展而来的一项全新测量技术,兼具定位精度高(厘米级)、测量耗时短等优势[13-14];但是,该方法的实际应用受到卫星信号限制,成本较高且要求操作人员到达实地测量点,从而在一定程度上限制了其广泛应用。摄影测量法可准确测量地表形态变化,早在1980年代末便被用于切沟侵蚀研究[12]。传统的立体摄影测量法需专业人员设计专门拍摄路线并利用笨重、昂贵的量测相机对切沟进行连续拍照,总体应用有限。随着计算器视觉技术的快速发展,尤其是尺度不变特征转换等图像特征检测和匹配技术的出现,基于运动恢复结构(structure from motion,SfM)的摄影测量方法应运而生[11]。该方法利用一部普通相机从多角度对切沟进行连续拍照,然后借助若干已知地面控制点坐标便可重建所测地物的三维模型。尽管在切沟监测中的应用较晚[8],但该方法因具有便携、快速、成本低和操作简单等优点,发展势头迅猛[12,15-16]。
目前,SfM摄影测量法在国内切沟监测中的应用还较少,沟壑纵横、切沟侵蚀严重的黄土高原也不例外。本研究在黄土高原六道沟小流域选取规模不等、形态各异的7条典型切沟,以LS和RTK GPS的测量结果为基准,评估SfM摄影测量法对切沟和切沟侵蚀的测量精度,以期为SfM摄影测量法在黄土高原的实践提供参考。

1 材料与方法

1.1 研究区概况

研究区为陕西省神木市的六道沟小流域(110°21'~110°23' E,38°56'~38°51' N),总面积6.9 km2,主沟道长4.2 km。该区地处毛乌素沙地边缘,属水蚀、风蚀交错带(图1a);气候为中温带半干旱气候,多年平均降水量为437.4 mm,主要集中在夏季;植被为沙生植物与黄土区植物交错出现;区内地面组成物质主要为第四纪黄土沉积物,土壤类型包括沙黄土、黄绵土、风沙土及坝地淤土,土壤结构松散,易遭受侵蚀[17-18]。
图1 六道沟小流域(a)及所选切沟(b)的地理位置

Fig.1 Geographical locations of Liudaogou catchment (a) and gullies selected (b) for investigation

2018年,在六道沟小流域开展实地调查,选取不同规模和形态的代表性切沟共7条(图1b),于植被尚未完全萌发的初春(4月底、5月初)及植被枯萎凋落的秋末(10月)进行测量。在春季,利用SfM摄影测量法对1~6号切沟(记为1S、2S、3S、4S、5S、6S)进行测量。由于三维激光扫描及RTK GPS设备均为外借,时间所限,仅测量了其中部分切沟;前者对1~4号及6号共5条切沟进行了测量、后者只对规模最大及中等的3号和5号切沟进行了测量。在秋季,同时利用SfM摄影测量法和RTK GPS对5号和7号切沟(记为5F、7F)进行测量。在7条切沟中,仅5号切沟在春、秋季均利用RTK GPS进行了测量,通过对比可计算该切沟在当年雨季的侵蚀体积。将该计算结果作为参考基准,评估SfM摄影测量法对切沟侵蚀的测量精度。

1.2 野外切沟测量

1.2.1 SfM摄影测量法

进行SfM摄影测量之前,在切沟沟缘、沟底的局部地势平坦且无植被遮挡处布设十字形地面控制点(ground control point,GCP),各GCP均匀分布,在沟缘转折处适当加密。利用Leica D510激光测距仪测量相邻GCP之间的空间直线距离及两点连线与水平面之间的夹角,在此基础上计算两点之间的水平投影距离和相对高程,再推算各GCP的空间直角相对坐标[19]。
各切沟利用Canon EOS 70D相机进行拍摄,相机焦距设置为18 mm。为确保后期三维重建的顺利进行,相邻照片的重叠率至少为60%[20]。拍摄照片数量因切沟规模而异,介于1 320~11 600张(表1)。将各切沟照片导入软件Agisoft Metashape 1.8.2中,将清晰度级别不足的照片删除后进行三维重建,依次生成稀疏点云、稠密点云和0.01 m×0.01 m分辨率的DEM[21-22]。
表1 各切沟拍摄照片数量

Tab.1 Number of images taken for each gully

切沟编号 测量时间 照片数量/张
1 春 1 320
2 春 7 442
3 春 7 956
4 春 3 665
5 春 2 015
5 秋 1 863
6 春 11 600
7 秋 2 268

1.2.2 LS

利用激光扫描法进行野外实地测量时,先在所测区域内选择宽阔、无地物与植被遮挡的位置布设棱镜靶标,利用Leica Nova MS60三维激光扫描全站仪进行地理坐标校正,至少重复3次。校正结束后,利用该全站仪扫描切沟不同部位,并在地形剧烈变化处适当增加扫描次数。将扫描得到的点云数据导入ArcGIS 10.2,基于不规则三角网模型(triangulated irregular network,TIN)构建DEM。为与SfM摄影测量法结果进行对比,利用重采样生成0.01 m×0.01 m分辨率的DEM。

1.2.3 RTK GPS

利用中海达iRTK2经典版RTK GPS以大致50 cm的间距测量各切沟的沟缘、沟坡及沟底,并在地形剧烈变化处适当增加测量点。与LS类似,将点云数据导入ArcGIS 10.2中构建DEM,并重采样至0.01 m × 0.01 m分辨率。

1.3 切沟相对DEM计算

3种测量方法采用的坐标系不同,得到的切沟DEM无法直接进行比较。RTK GPS采用大地坐标系;LS和SfM摄影测量法的点云数据虽然均采用独立空间直角坐标系,但二者的坐标原点并不完全一致。因此,先使用ArcGIS 10.2中的Georeferencing工具,以LS或RTK GPS得到的数据为基准,对SfM摄影测量法所得DEM进行配准。初春的植被虽较稀疏,也会对切沟沟缘的识别造成一定的影响。在春季,使用3种方法对3号切沟进行了测量,计算其相对DEM之前,以RTK GPS数据为基准对LS和SfM摄影测量法所测数据进行配准;接下来,结合目视解译和局部坡度变化识别不同测量方法得到的各切沟沟缘,利用TIN生成各切沟沟缘面DEM并重采样至0.01 m × 0.01 m分辨率,将各切沟DEM减去相应的沟缘面DEM得到相对DEM后再行比较。在相对DEM中,沟底的高程为负值,且绝对值大;沟缘的高程绝对值小,在理想状况下应接近0。然而,沟缘经矢量勾绘后转化为栅格并重采样得到,在与切沟DEM做相减计算时往往会出现对应栅格高程不等于0的情况。

1.4 切沟形态参数提取与对比

切沟各形态参数通过沟缘面DEM和切沟相对DEM在ArcGIS 10.2中加以提取。基于沟缘面DEM可提取切沟长度、周长、面积和周长-面积比[23],基于切沟相对DEM可得到平均表面宽度、平均深度和体积。利用相对误差(E)评估SfM摄影测量法对各切沟形态参数的测量精度[9]:
E= m i - M i M i×100%。
式中:mi为SfM摄影测量法所测第i条切沟的形态参数;Mi为LS或RTK GPS所测第i条切沟的形态参数。

2 结果与分析

2.1 切沟测量精度

2.1.1 SfM与LS所测切沟DEM对比

图2对比了春季SfM与LS所测各切沟的相对DEM。1号(1S)、2号(2S)、3号(3S)、4号(4S)、6号(6S)切沟的高程平均差值分别为-0.077、0.120、-0.162、0.436和0.373 m。
图2 SfM(a~e)与LS(f~j)所测1S、2S、3S、4S和6S相对DEM

注:网络版为彩图。

Fig.2 Relative DEM of 1S, 2S, 3S, 4S and 6S measured by SfM (a~e) and LS (f~j)

然而,从高程差值的累积频率分布(图3)与空间分布(图4)来看,2种方法对1S和3S的测量差异总体较小,相对高程差值主要集中在-0.5~0.5 m,对应栅格占比均在85%左右。SfM较LS总体高估了4S和6S的高程,各栅格高程差值多为正值,主要集中在0~1.5 m范围内,对应栅格占比分别为70.2%和67.8%,介于-0.5~0.5 m的栅格占比分别为41.3%和43.2%。至于2S,尽管2种方法所测平均高程差值较小,但空间变异较大,相对高程差值的分布较分散(图3)。高程差值主要集中在-2~2 m,对应栅格占比为70.0%,介于-0.5~0.5 m的栅格仅占22.0%。
图3 SfM与LS所测不同切沟相对高程差值累积频率分布

Fig.3 Cumulative frequency distribution of relative elevation differences between SfM and LS for different gullies

图4 SfM与LS所测1S(a)、2S(b)、3S(c)、4S(d)和6S(e)相对高程差值空间分布

注:网络版为彩图。

Fig.4 Spatial distribution of relative elevation differences of 1S (a), 2S (b),3S (c), 4S (d) and 6S (e) between SfM and LS

2.1.2 SfM与RTK GPS所测切沟DEM对比

利用SfM和RTK GPS测量春季3号(3S)和5号(5S)切沟,以及秋季5号(5F)和7号(7F)切沟的相对DEM(图5)发现,2种方法所测各切沟相对高程的平均差值不大,依次为-0.033、-0.096、-0.128和-0.075 m。从累积频率分布(图6)与空间分布(图7)来看,SfM与RTK GPS所测5S、5F和7F的高程差值均主要集中在-0.5~0.5 m,对应栅格占比分别为95.6%、96.3%和94.7%,超过±1.0 m的栅格仅占0.1%、1.2%和0.3%。相比之下,SfM和RTK GPS对3S的测量差异略大,高程差值介于-0.5~0.5 m的栅格较少,对应占比为85.3%,甚至还有3.8%的栅格高程差值超过±1 m。
图5 SfM(a~d)与RTK GPS(e~h)所测3S、5S、5F、7F相对DEM

注:网络版为彩图。

Fig.5 Relative DEM of 3S, 5S, 5F and 7F measured by SfM (a~d) and RTK GPS (e~h)

图6 SfM与RTK GPS所测不同切沟相对高程差值累积频率分布

Fig.6 Cumulative frequency distribution of relative elevation differences between SfM and RTK GPS for different gullies

图7 SfM与RTK GPS所测3S(a)、5S(b)、5F(c)、7F(d)相对高程差值空间分布

注:网络版为彩图。

Fig.7 Spatial distribution of relative elevation differences of 3S (a),5S (b), 5F (c) and 7F (d) between SfM and RTK GPS

2.1.3 切沟形态参数对比

SfM、LS和RTK GPS所测各切沟形态参数如表2所示。配对t检验结果显示,无论与LS还是与RTK GPS相比,SfM所测各切沟形态参数均无显著差异(P>0.05)。
表2 SfM、LS和RTK GPS所测各切沟形态参数

Tab.2 Morphological parameters of each gully measured by SfM, LS and RTK GPS

切沟
编号		 测量
时间		 测量
方法		 长度/m 平均表面
宽度/m		 平均深
度/m		 周长/m 面积/m2 周长-面积
比/m-1		 体积/m3
1 春 SfM 16.57 6.10 2.77 37.82 75.98 0.50 210.53
LS 15.85 6.21 2.62 37.75 78.00 0.48 203.97
2 春 SfM 52.28 6.44 3.50 86.57 205.17 0.42 717.16
LS 50.71 6.90 3.62 81.66 217.62 0.38 786.82
3 春 SfM 73.00 8.84 1.46 177.92 600.84 0.30 875.38
LS 70.41 8.78 1.30 178.37 598.01 0.30 777.95
RTK GPS 70.75 8.81 1.42 181.11 602.09 0.30 856.32
4 春 SfM 12.53 4.54 1.58 37.33 56.22 0.66 89.01
LS 12.32 4.64 1.98 37.38 56.70 0.66 112.25
5 春 SfM 12.06 4.36 1.04 32.60 49.75 0.66 51.57
RTK GPS 11.98 4.40 0.95 31.85 49.34 0.65 46.73
5 秋 SfM 12.65 4.59 1.06 32.81 52.58 0.62 55.79
RTK GPS 12.55 4.54 0.93 34.36 52.12 0.66 48.55
6 春 SfM 31.31 2.09 0.62 46.55 35.13 1.33 21.62
LS 31.16 2.37 0.92 44.82 37.55 1.19 34.55
7 秋 SfM 15.27 2.38 0.80 31.50 25.32 1.24 20.19
RTK GPS 14.29 2.27 0.76 32.67 24.28 1.35 18.58
相比LS,SfM所测各切沟长度和周长精度最高,相对误差均不超过±6%,平均相对误差仅分别为2.7%和1.9%(表3)。SfM对切沟平均表面宽度的测量精度也较高,除6S的相对误差高达11.8%外,各切沟的相对误差均不超过±7%。相比之下,SfM对切沟深度的测量误差较大。例如,SfM较LS大面积高估了4S和6S的高程(图4d、4e),其对平均深度的低估百分比分别高达20.2%和32.6%。对于其余3条切沟,SfM对平均深度的测量误差介于-4%~13%。尽管SfM所测各切沟面积与LS结果相近,相对误差均不超过±7%;由于对深度的测量误差,除1S和2S外,SfM对3S、4S和6S的体积测量误差均超过±10%,分别为12.5%、-20.7%和—37.4%。两种方法测得的各切沟周长-面积比也十分接近(表2),但由于LS所测周长-面积比本身较小,部分切沟的相对误差超过±10%(表3)。
表3 SfM所测不同切沟形态参数的相对误差

Tab.3 Relative errors of different gully morphological parameters measured by SfM 单位:%

形态参数 与LS对比 与RTK GPS对比
最小值 最大值 平均值 最小值 最大值 平均值
长度/m 0.5 4.5 2.7 0.7 6.9 2.9
平均表面宽度/m -11.8 0.7 -4.3 -0.9 4.8 1.3
平均深度/m -32.6 12.3 -7.6 2.8 14.0 7.9
周长/m -0.3 6.0 1.9 -4.5 2.4 -1.9
面积/m2 -6.4 0.5 -3.0 -0.2 4.3 1.4
周长-面积比/m-1 -0.7 12.4 5.3 -7.5 1.5 -3.2
体积/m3 -37.4 12.5 -10.2 2.2 14.9 9.0
对同时利用RTK GPS和SfM 测量的3S、5S、5F和7F各形态参数进行统计分析发现,2种方法所测结果总体相差不大(表2)。与SfM和LS的对比类似,SfM对各切沟长度、平均表面宽度、周长、面积和周长-面积比的测量误差较小,均不超过±8%,平均相对误差依次为2.9%、1.3%、-1.9%、1.4%和-3.2%;对深度和体积的测量误差较大,但也不超过±15%,平均相对误差分别为7.9%和9.0%(表3)。

2.1.4 切沟测量精度影响因素

与LS和RTK GPS相比,SfM对切沟的测量误差主要与植被覆盖、切沟沟缘曲折程度、光照条件以及沟坡形态有关。尽管开展野外调查的初春和秋末植被较为稀疏,但其也会对LS和SfM的测量结果造成干扰,影响切沟沟缘的识别和DEM重建,使2种方法所测DEM无法完全重合、相对高程出现差异(图4)。相对而言,RTK GPS直接对地表进行测量,基本不受植被和枯落物干扰[12]。当SfM与其进行对比时,往往会高估地表高程。例如,3S沟头附近的沟底发育了1条小沟,但由于沟底植被干扰,SfM难以对其进行完整识别和测量,导致SfM在此所测高程偏大(图8)。
图8 SfM与RTK GPS所测3S相对高程差值空间分布(a)及沟底小沟实地拍摄照片(b)

注:网络版为彩图。

Fig.8 Spatial distribution of relative elevation differences of 3S between SfM and RTK GPS (a) and a photo of a small channel at the bottom of 3S (b)

总体而言,切沟的周长-面积比越大,沟缘越曲折,SfM的测量误差相应越大。例如,在SfM和LS同时开展监测的5条切沟中,6S的沟缘最为曲折,LS所测周长-面积比高达1.19 m-1,SfM所测面积和体积误差均最大(表2)。此外,6S因拍摄过多照片造成的信息冗余(表1),可能也是造成SfM测量误差的原因之一[24]。而在与RTK GPS的对比中,7F的周长-面积比最大(表2),SfM所测切沟长度、平均表面宽度和面积误差均最大。
SfM对切沟的测量精度还受到光照条件的显著影响。例如,SfM和LS所测4S的高程差值呈西低东高的态势(图4d),这主要是因为SfM测量当天下午光线充足,斜射的阳光只可到达东岸,无法进入较深的沟底。东岸的细节充足曝光,相机与所摄地物的距离较实际偏小,进行三维重建的三角测量计算时得到的高程偏大(沟深被低估),与LS结果相减为正值。相比之下,西岸光线较暗,物体远近对比较不明显,相机与地物的距离较实际偏大,SfM后期计算得到的高程偏小(沟深被高估),与LS相减为负值。可见,SfM的应用受限于天气情况,最好在阴天或多云天气进行,确保所测地物各部位的光线基本一致[19,25]。
沟坡形态主要影响SfM与RTK GPS对切沟的测量差异。对于沟坡上某些坡度较大的部位,RTK GPS无法直接进行测量,只能通过沟缘和沟底的测量结果插值获得。如果沟坡并非直形坡,而呈凸形、凹形等不规则形状,那么SfM与RTK GPS所测高程就会出现较大差异。SfM与RTK GPS共同测量的各切沟沟坡均出现了局部坡度超过60°的陡峭凸出部位(图9),SfM所测相对高程明显较RTK GPS大(图7)。可见,相比RTK GPS,SfM对地形陡峭、无法到达地物存在明显的测量优势[12]。然而,SfM对地物的三维重建高度依赖GCP的坐标信息,对于黄土高原某些沟深坡陡的切沟(如2S)来说,GCP只能布设在沟岸和沟底,可能会在一定程度上影响SfM对切沟的测量精度(图3、图4b)[19]。
图9 SfM所测3S(a)、5S(b)、5F(c)和7F(d)坡度空间分布

注:网络版为彩图。

Fig.9 Spatial distribution of slope degree of 3S (a), 5S (b), 5F (c) and 7F (d) measured by SfM

2.2 切沟侵蚀测量精度

2.2.1 SfM与RTK GPS所测切沟侵蚀对比

将SfM和RTK GPS所测5F的DEM减去5S的DEM,得到5号切沟在当年雨季的土壤侵蚀空间分布(图10)。2种方法都发现该切沟的沟头处发生侵蚀,沟尾处出现堆积,但SfM测得的切沟沟底高程变化明显较RTK GPS复杂。此外,SfM所测切沟沟缘尤其是东侧和北侧沟缘的高程变化较RTK GPS所测结果大,且主要表现为高程下降。从切沟形态参数(表4)来看,相比RTK GPS,SfM所测切沟平均表面宽度和体积变化较大,但周长变化明显偏小。对于切沟长度、平均深度和面积变化,2种方法所测结果相差不大。
图10 SfM(a)与RTK GPS(b)所测5号切沟侵蚀空间分布

注:网络版为彩图。

Fig.10 Spatial distribution of soil erosion measured by SfM (a) and RTK GPS(b) for gully No. 5

表4 SfM和RTK GPS所测5号切沟形态参数变化

Tab.4 Changes of morphological parameters of gully No.5 measured by SfM and RTK GPS

测量方法 长度/m 平均表面宽度/m 平均深度/m 周长/m 面积/m2 体积/m3
SfM 0.59 0.23 0.02 0.21 2.84 4.22
RTK GPS 0.57 0.14 -0.02 2.51 2.78 1.81

2.2.2 切沟侵蚀测量精度影响因素

与切沟监测类似,SfM对切沟侵蚀的测量精度主要与植被及其枯落物覆盖有关。春季萌发的植被会对沟缘线造成遮挡,当使用SfM测量时,解译得到的沟缘高程偏大且易向切沟内侧偏移,尤其是植被较为丰富的东侧和北侧沟缘;到了秋季,植被枯萎、叶片凋落,切沟沟缘出露。因此,相比基本不受植被干扰的RTK GPS,SfM所测切沟沟缘高程和平均表面宽度变化均较大(表4),因2种方法所测切沟面积变化相差不大,因而SfM所测切沟体积变化相应较大(表4)。然而,雨季过后,SfM所测切沟周长仅增加0.21 m,远小于RTK GPS所测结果。当使用RTK GPS进行测量时,5号切沟在雨季前后的周长-面积比基本相当,可见侵蚀并未显著改变沟缘的曲折程度。2种方法测得的切沟周长差异可能主要还是与沟缘植被有关。雨季前,由于沟缘植被干扰,SfM所测切沟沟缘较曲折,周长-面积比为0.66 m-1,周长与RTK GPS测量结果相差不大(表2)。雨季过后,尽管侵蚀使切沟有所扩张,但植被枯萎凋落后沟缘曲折程度降低,SfM所测周长-面积比降至0.62 m-1(表2),该方法测量得到的周长增加有限(表4)。因设备所限,本研究仅对1条切沟的雨季侵蚀进行了系统测量。未来有必要增加切沟测量条数,尤其是雨季期间土壤侵蚀明显的切沟,系统开展SfM对切沟侵蚀的测量精度评估。

2.2.3 切沟侵蚀影响因素

从SfM和RTK GPS所测5号切沟的相对DEM(图10)和各形态参数变化(表4)可以发现,该切沟在雨季前后总体变化不大。一方面,该切沟接近小流域分水岭,汇水面积小,侵蚀动力不大;另一方面,该切沟沟缘处的植被生长较好,沟底的枯落物也较丰富,在一定程度上可保护切沟,削弱降雨和径流侵蚀。

3 结论

本文以黄土高原六道沟小流域内7条规模与形态各异的典型切沟为研究对象,通过与LS和RTK GPS对比,评估了SfM摄影测量法对切沟和切沟侵蚀的测量精度。研究结果表明,SfM摄影测量法可作为一种简单、便捷的高精度方法应用于黄土高原的切沟和切沟侵蚀测量。与LS和RTK GPS相比,SfM摄影测量法所测各切沟相对DEM总体差异不大,平均高程差值介于-0.162~0.436 m;各切沟长度、表面宽度、周长和面积的相对误差基本不超过±7%;深度和体积的相对误差略大,介于-37.4%~18.0%。SfM对切沟的测量精度主要与植被覆盖、切沟沟缘曲折程度、光照条件以及沟坡形态有关。因设备用时所限,仅评估了其中1条切沟的雨季侵蚀,发现RTK GPS和SfM摄影测量法所测土壤侵蚀空间分布与切沟形态参数变化均十分接近。
致谢:诚挚感谢西华师范大学张斌教授、硕士研究生郑吉林和吴曦在野外调查工作中的协助。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
刘宝元, 杨扬, 陆绍娟. 几个常用土壤侵蚀术语辨析及其生产实践意义[J]. 中国水土保持科学, 2018, 16(1):9-16.

LIU B Y, YANG Y, LU S J. Discriminations on common soil erosion terms and their implications for soil and water conservation[J]. Science of Soil and Water Conservation, 2018, 16(1):9-16.

[2]
张光辉. 切沟侵蚀研究进展与展望[J]. 水土保持学报, 2020, 34(5):1-13.

ZHANG G H. Advances and prospects for gully erosion researches[J]. Journal of Soil and Water Conservation, 2020, 34(5):1-13.

[3]
POESEN J, NACHTERGAELE J, VERSTRAETEN G, et al. Gully erosion and environmental change:importance and research needs[J]. Catena, 2003, 50(2/3/4):91-133.

DOI

[4]
景可. 黄土高原沟谷侵蚀研究[J]. 地理科学, 1986, 6(4):340-347.

JING K. A study on gully erosion on the Loess Plateau[J]. Scientia Geographica Sinica, 1986, 6(4):340-347.

DOI

[5]
郑粉莉, 徐锡蒙, 覃超. 沟蚀过程研究进展[J]. 农业机械学报, 2016, 47(8):48-59,116.

ZHENG F L, XU X M, QIN C. A review of gully erosion process research[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(8):48-59,116.

[6]
张姣, 郑粉莉, 温磊磊, 等. 利用三维激光扫描技术动态监测沟蚀发育过程的方法研究[J]. 水土保持通报, 2011, 31(6):89-94.

ZHANG J, ZHENG F L, WEN L L, et al. Methodology of dynamic monitoring gully erosion process using 3D laser scanning technology[J]. Bulletin of Soil and Water Conservation, 2011, 31(6):89-94.

[7]
DONG Y F, XIONG D H, SU Z A, et al. Critical topographic threshold of gully erosion in Yuanmou Dry-hot Valley in Southwestern China[J]. Physical Geography, 2013, 34(1):50-59.

DOI

[8]
KAISER A, NEUGIRG F, ROCK G, et al. Small-scale surface reconstruction and volume calculation of soil erosion in complex Moroccan gully morphology using structure from motion[J]. Remote Sensing, 2014, 6(8):7050-7080.

DOI

[9]
宋晓鹏, 张岩, 王志强, 等. 无人机摄影测量提取黄土高原切沟参数精度分析[J]. 北京师范大学学报(自然科学版), 2021, 57(5):606-612.

SONG X P, ZHANG Y, WANG Z Q, et al. Accuracy of gully morphological parameters extracted by UAV photogrammetry in the Loess Plateau[J]. Journal of Beijing Normal University (Natural Science), 2021, 57(5):606-612.

[10]
权鑫, 史红伟, 蔡强国, 等. 室内模拟降雨情形下黄土坡面沉降程度的试验研究[J]. 陕西师范大学学报(自然科学版), 2021, 49(2):98-106.

QUAN X, SHI H W, CAI Q G, et al. Experimental study on the settlement of loess slope under indoor simulated rainfall[J]. Journal of Shaanxi Normal University (Natural Science Edition), 2021, 49(2):98-106.

[11]
JAMES M R, ROBSON S. Straightforward reconstruction of 3D surfaces and topography with a camera: accuracy and geoscience application[J]. Journal of Geophysical Research, 2012, 117:F03017.

[12]
何涛, 史扬子, 杨扬, 等. 切沟侵蚀测量方法研究进展[J]. 中国水土保持科学, 2021, 19(2):151-158.

HE T, SHI Y Z, YANG Y, et al. Research progress on gully erosion measurement methods[J]. Science of Soil and Water Conservation, 2021, 19(2):151-158.

[13]
WU Y Q, CHENG H. Monitoring of gully erosion on the Loess Plateau of China using a global positioning system[J]. Catena, 2005, 63(2/3):154-166.

DOI

[14]
ZHANG Z D, LUO J Y, CHEN B. Spatially explicit quantification of total soil erosion by RTK GPS in wind and water eroded croplands[J]. Science of the Total Environment, 2020, 702:134716.

DOI

[15]
李俊利, 李斌兵, 柳方明, 等. 利用照片重建技术生成坡面侵蚀沟三维模型[J]. 农业工程学报, 2015, 31(1):125-132.

LI J L, LI B B, LIU F M, et al. Generating 3D model of slope eroded gully based on photo reconstruction technique[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(1):125-132.

[16]
ZHENG F Z, WACKROW R, MENG F R, et al. Assessing the accuracy and feasibility of using close-range photogrammetry to measure channelized erosion with a consumer-grade camera[J]. Remote Sensing, 2020, 12(11):1706.

DOI

[17]
唐克丽, 侯庆春, 王斌科, 等. 黄土高原水蚀风蚀交错带和神木试区的环境背景及整治方向[J]. 水土保持研究, 1993(2):2-15.

TANG K L, HOU Q C, WANG B K, et al. The environment background and administration way of wind-water erosion crisscross region and Shenmu experimental area on the Loess Plateau[J]. Research of Soil and Water Conservation, 1993(2):2-15.

[18]
黄婷婷, 史扬子, 曹琦, 等. 黄土高原六道沟小流域近30年来土壤侵蚀变化评价[J]. 中国水土保持科学, 2020, 18(1):8-17.

HUANG T T, SHI Y Z, CAO Q, et al. Soil erosion evaluation of Liudaogou catchment in the Loess Plateau during the past 30 years[J]. Science of Soil and Water Conservation, 2020, 18(1):8-17.

[19]
史扬子, 黄婷婷, 罗建勇, 等. SfM摄影测量法测定切沟的精度评价[J]. 水土保持研究, 2020, 27(1):39-46.

SHI Y Z, HUANG T T, LUO J Y, et al. Evaluation on the accuracy of structure from motion (SfM) photogrammetry on gully surveying[J]. Research of Soil and Water Conservation, 2020, 27(1):39-46.

[20]
LINDER W. Digital photogrammetry:a practical course[M].2nd ed.Berlin:Springer, 2006.

[21]
CARRIVICK J L, SMITH M W, QUINCEY D J. Structure from motion in the Geosciences[M]. Oxford:Wiley-Blackwell, 2016.

[22]
YANG Y, SHI Y Z, LIANG X Z, et al. Evaluation of structure from motion (SfM) photogrammetry on the measurement of rill and interrill erosion in a typical loess[J]. Geomorphology, 2021, 385:107734.

DOI

[23]
TORRI D, POESEN J, ROSSI M, et al. Gully head modelling:a Mediterranean badland case study[J]. Earth Surface Processes and Landforms, 2018, 43(12):2547-2561.

DOI

[24]
王佳文, 李彩林, 王志勇, 等. 影像数量对SFM三维重建模型测量精度的影响分析[J]. 激光杂志, 2021, 42(3):63-69.

WANG J W, LI C L, WANG Z Y, et al. Analysis on the influence of images on measurement accuracy of three-dimensional model reconstructed by SFM[J]. Laser Journal, 2021, 42(3):63-69.

[25]
WENZEL K, ROTHERMEL M, FRITSCH D, et al. Image acquisition and model selection for multi-view stereo[EB/OL].[2022-05-08].DOI:10.5194/isprsarchives-XL-5-W1-251-2013.

Options
文章导航

/

〈 〉

版权所有 © 《陕西师范大学学报(自然科学版)》编辑部
地址：中国西安陕西师范大学长安校区文汇楼C段3层　
邮编：710119
电话：029-81530883，81530879，81530885　
E-mail：lkxuebao@snnu.edu.cn
技术支持：北京玛格泰克科技发展有限公司