一种无人机点云数据的土壤粗糙度测量方法

基于无人机点云数据土壤粗糙度测量方法

张田1,2,3, 周忠发1,2,3, 王玲玉1,2,3, 赵馨1,2,3, 张文辉1,2,3, 张淑1,2,3, 王宇1,2,3

1.贵州师范大学地理与环境科学学院/喀斯特研究院,贵阳 550001

2.贵州省喀斯特山地生态环境国家重点实验室培育基地,贵阳 550001

3.国家喀斯特石漠化防治工程技术研究中心,贵阳 550001

摘要:

耕地土壤粗糙度是影响土壤水分、微波遥感观测和植物生长等农情信息监测方面的重要因子。土壤粗糙度通常是根据野外拍摄照片进行判读,但存在判读效率低、受人为因素影响等不足。无人机低空遥感对观测地表起伏状况具有良好的敏感性,因此基于摄影测量利用无人机对地表进行拍摄,将其结果与粗糙度测量板获取的数据进行对比,探讨无人机数据测量土壤粗糙度的精确性。研究结果表明: 近景摄影测量的平均绝对误差主要集中在0.4~1.2 cm之间,平均相对误差为6.16%,均方根误差为0.40 cm。表明无人机点云摄影测量可有效运用于地表粗糙度的测量,且当分块采样的面积越小时,获得粗糙度的值越精确。

0 引言

对土壤粗糙度的定义一般从2个角度进行描述: 一个是从空气动力学概念出发,将地表风速为0的高度定义为土壤粗糙度[1]; 另一个从地形学角度出发,将地表凹凸不平的程度定义为地表粗糙度,也称地形粗糙度[2]。土壤粗糙度表征土壤表面的起伏程度,是影响雷达入射波和地面相互作用进而影响后向散射系数的参数,其测量结果的准确性直接影响土壤含水量反演精度[3⇓-5]; 在农情信息监测方面,对耕地土壤入渗参数、土壤孔隙率、地表径流与土壤侵蚀有重要影响[6⇓-8],快速、准确地测量地表粗糙度一直以来都是微波定量遥感和农情信息监测的重要内容。

目前,地表粗糙度的测量方法根据其特点可分为接触式和非接触式2种。常用的接触式测量主要有插板法、链条法和销钉高程测量[9]等。接触式的测量方法获取的是局部粗糙度特征,全域粗糙度特征的获取工作量较大,推广应用具有较大难度[10]。常用的非接触式测量方式主要包括: 超声波测量法[11]、摄影测量法[12]和激光扫描法[13-14]等。超声波测量法容易造成反射,测量精度不高; 激光扫描法主要是通过手持或三脚架固定扫描仪对目标地块进行扫描,其优点是测量精度高、效率高、无扰动[15],但其所需的设备复杂、价格昂贵,且获取的数据存在明显的阴影效应[16],在目标地块较大时,数据采集的工作量也随之增加; 相比之下,摄影测量法的精度虽不及激光扫描法高,但其耗时短、所需设备简单、影像分辨率高[17-18]等优点使其应用更具有推广性[19]。

低空遥感技术是在无人机上搭载摄像头,通过无线电远程控制无人机在低空进行拍摄[20-21]。通过无人机获取较高质量的、具有一定重叠度的密集点云,经过一定的预处理(拼接、降噪和重采样等)构建目标物体的三维模型。目前,点云数据的应用主要集中在工程应用[22]、林业调查[23-24]以及城市三维建模[25-26]等方面,现阶段将点云数据应用到土壤表面粗糙度测量的研究较少。

因此,本文结合激光扫描和摄影测量的优点,利用无人机获取研究区耕地的裸土点云,构建耕地表面三维模型,基于该三维模型模拟真实土壤表面,提取该模型表面所有点云的高程值,从而计算土壤粗糙度的各项参数。将此方法计算的土壤粗糙度参数和传统粗糙度测量板测得的结果进行对比和分析,讨论无人机点云数据对获取地表粗糙度的适用性和准确性。以期为耕地土壤粗糙度提供精准、无损、大范围、高效的测量新思路和方法,为微波定量遥感和农情信息监测等提供可靠的数据支撑。

1 数据源

研究区位于贵州省安顺市关岭布依族苗族自治县花江镇境内的一块裸土耕地,样区中心坐标为25°49'43.37″N,105°35'34.44″E,该样区选在平坦、开阔的平地,且经过表层翻耕,避免了因地形及表面作物导致的测量误差,选取的样地面积为20 m×16 m,将其划分为20个面积为4 m×4 m样方。

1.1 测量板实测数据

本文使用的测量板有水平和垂直2个方向的刻度线,间隔为1 cm,如图1(a)所示。将测量板垂直插入划分好的土壤样方中,拍摄角度与地面平行,每个样方中进行一次测量板拍照。内业计算利用Matlab软件对该20张测量板照片进行裁剪、校正并分别提取土壤表面轮廓线。由于粗糙度的测量结果受测量方向的影响较大,因此测量板的数据采集在a,b,c,d这4个方向随机分布,如图1(b)所示。

1.2 无人机点云摄影测量

实验使用的无人机为大疆-御2专业版,无人机的具体飞行参数如表1所示。

表1   无人机、相机及飞行参数

拍摄时间为2020年12月9日11: 00,在研究区连续拍摄了94张正射影像,坐标系为WGS-84坐标,飞行高度为10 m,对应影像的地面分辨率为0.43 cm,航向和旁向重叠度均为75%,保证了获取的每张影像之间均有足够的重叠度生成匹配点云。

2  研究方法

2.1 无人机影像处理

利用Agisoft PhotoScan摄影测量软件对无人机影像进行拼接处理,并基于内插法生成密集点云,共生成20 986 847个点,利用点云生成网络和纹理,最后构建三维实景模型。

三维重建的原理主要是实现点云的网格化处理[27],其采用的算法是运动恢复结构(structure from motion,SFM)算法,SFM算法直接将摄影测量与计算机视角算法相结合,对不同距离、不同视角的重叠影像自动提取影像特征并进行匹配[28],该算法在影像的重叠部分自动检测影像特征点(即几何相似性),并检测这些特征点在多幅影像序列中的运动特征和相机参数,在保证影像的连续性和拍摄时间间隔较短的前提下,基于每张影像的位置、方向、焦距和特征点的运动规律就可以从多幅影像中的二维坐标计算任意位置的三维坐标,从而校正影像并创建点云; 另一方面,当相机参数未知时,可以从不同影像中获取等值点的二维坐标导出其在空间中的三维坐标[29⇓⇓-32]。因此,在满足影像重叠率、地面分辨率和连续性的条件下,无需地面控制点,基于无人机在拍摄影像时写入的全球定位系统(global positioning system,GPS)方位和位置姿态系统(position and orientation system,POS)信息可生成高精度的密集点云。

基于每个样本中测量板的采样方向,对点云模型进行剖面采样,即在点云模型上截取和测量板方向相同的剖面,比较2种方法获取的土壤粗糙度差异。技术路线如图2所示。为避免采样长度和采样间隔对土壤粗糙度测量结果的影响,本文在样方内采取分块采样的方式提取土壤表面粗糙度,即将线采样模式上升为面采样模式,如图3所示。将Agisoft PhotoScan软件导出的三维表面模型导入三维点云编辑和处理软件Cloud Comparer中,在拍摄测量板的同一样方内进行点云数据的分块采样,将其剪裁为20个1 m×1 m的样块作为样本,再将裁剪后的点云数据输出为文本格式,将大地坐标系转换成影像坐标系,获取土粒之间相对位置与高程,从而得出土壤粒径的大小,计算粗糙度的各项参数。

2.2 土壤粗糙度计算

土壤表面起伏度的发生和演变是一个随机分布的过程,即其统计量不随空间尺度的变化而变化。因此,一般是利用其统计特性表征其土壤粗糙度[33],包括均方根高度(root mean square height,RMSH)和自相关长度(autocorrelation length,ACL),二者分别在垂直和水平2个方向上对土壤粗糙度的大小进行描述[34],RMSH反映了土壤的平均高度偏离的程度,描述的是随机过程中在垂直方向上各个点之间相互独立的特征。而ACL描述了土壤的高度在水平方向上的变化,反映的是地表高度在水平方向上的关联性。

2.2.1 RMSH

一般当地表的粗糙度越大时,RMSH值就越大[35]。

2.2.2 ACL

ACL是估计地表2点之间相互独立的标准[36]。ACL是根据测量板上x轴方向上2点之间的间隔来定义,一般是通过利用一组数据的归一化自相关函数ρ(x)计算ACL值。通常,如果2点间的距离大于ACL值,则认为这2点的高度值在统计学意义上是相互独立的。

2.2.3 精度验证

为定量描述无人机近景摄影测量方法的精度,以测量板的计算结果为基准值,衡量无人机点云摄影测量在剖面上的测量精度,选取绝对误差(absolute error,AE)、相对误差(relative error,RE)和均方根误差(root mean square error,RMSE)表示测量结果的误差,三者的值越小表明测量的精度越高[37]。

3 结果与分析

3.1 测量结果分析

3.1.1 测量板法与点云剖面测量结果分析

由采集得到的测量板数据的计算结果表明,研究区土壤地表具有较好的空间连续性,其RMSH和ACL的变异系数分别为56%和18%,这表明样区空间异质性较强。一般情况下,光滑土壤表面的RMSH<0.6 cm,粗糙土壤表面的RMSH≥0.6 cm[38-39]。如图4所示,该方法测量出的实验样区的RMSH值大多在0.8~1.4 cm之间,表明该区域的土壤粗糙度较大,即土壤表面较粗糙。ACL值主要集中30~55 cm之间,表明单个样本中各高度点之间的相互独立性较大。基于测量板的测量方向在点云摄影测量上截取剖面并计算剖面上土壤粗糙度的RMSH和ACL,通过对比2种方法的计算结果发现,RMSH和ACL的总体变化趋势一致,表明将无人机点云数据运用于地表粗糙度的测量具有可行性。

3.1.2 精度验证

以测量板法的结果为基准,对基于无人机点云面状样本为基础截取的线状剖面上的土壤粗糙度参数进行精度验证,摄影测量法的平均AE值为0.68 cm,其范围主要集中在0.4~1.2 cm之间,平均RE值为6.16%,RMSE值为0.40 cm; 整体误差较小; 二者的结果对比如表2所示。利用测量板法和近景摄影测量法对同一地块的土壤粗糙度进行测量,结果发现摄影测量法的测量结果比测量板法小。

将无人机点云摄影测量和测量板所得到的同一样方上土壤表面的轮廓线进行对比,如图5所示,测量板和无人机点云的三维建模方法所取得的二维样线趋势一致,说明无人机点云数据可有效应用于土壤表面粗糙度测量,但二者在mm尺度上仍具有一定的差异,其原因主要在于测量板主观因素较大、无人机在飞行过程中受外界以及机身的抖动等影响造成。2种方法获得的土壤表面轮廓线在水平方向上起伏较大,这也是造成其ACL值较大的原因。

3.2 无人机面状点云摄影测量结果分析

一般情况下,采样数目和采样间隔对土壤粗糙度测量结果的影响较大。且采样长度越长,获得的粗糙度的值越精确,但另一方面,内业的工作量将成倍增加,且采样长度和间隔的影响仍然不可避免。面状采样模式可有效地避免这一误差,且可以提取土壤表面更微观的粗糙度信息。由图6可以看出,RMSH值分别集中在0.02~0.03 cm之间,总体变化较为稳定,变化幅度减小; ACL值主要集中在1.0~2.0 cm之间,RMSH和ACL与截取剖面的计算结果相差较大,原因主要为摄影测量中i为样本中点云的数量,而非测量板上刻度的序列号,二者相差较大,造成整体的变化幅度较大。

为更直观地描述出在不同取样面积下土壤粗糙度的变化程度,在同一个样方上分别裁剪了面积为100 cm×100 cm,80 cm×80 cm,60 cm×60 cm,40 cm×40 cm,20 cm×20 cm,10 cm×10 cm的6块样本点。其测量结果如表3所示,当采样的面积逐渐减小时,RMSH值也在减小,对土壤粗糙度描述的精细度越高。ACL的比值随着采样面积的减小而出现较大波动,说明采样面积的大小对测量的结果影响较大。由图7中看出,采样面积减小时,土壤表面的起伏程度表现的更为详细,说明本文的技术手段可以获取微地貌的精准尺度信息,而所用的硬件属于常规的无人机搭载普通摄像头,所用的软件也是常见的图像处理软件,具有可推广性。

表3   不同采样面积下的土壤粗糙度参数测量结果

3.3 综合评价

测量板法成本与技术难度低、操作性强,但会对土壤表面造成一定的破坏; 基于无人机影像获取的点云数据测量地表土壤粗糙度,该方法为非接触的测量方式,不会对地表的原始形态造成破坏。除此之外,在密集点云的基础上采取了面状采样的方法进行测量,避免了因测量板的长度有限、采样间隔和采样方向(方向异质性)的不同引起的误差,提高了测量的精度。

2种方法外业数据采集和内业计算单个样本所需的具体时间如表4所示。本文采取的方法影像的后期处理所需的时间较长,尤其在拼接、建模的环节中,由于影像精度较高,数据量较大,计算机处理时间较长。分块采样的方式也使得工作量成倍增加,因而在今后的研究工作中,需要进一步的改进和优化。

表4   2种方法所需时间及精度对比

4 结论和讨论

近景摄影测量法利用无人机搭载摄像头获取待测区域的地表影像,生成密集点云数据后构建三维模型,在分块采样的基础上进行测量。与传统测量板法相比,具有如下优点:

1)稳定性及效率较高。由实验结果可知,近景摄影测量的平均绝对误差主要集中在0.4~1.2 cm之间,平均相对误差为6.16%,均方根误差为0.40 cm; 近景摄影测量面状采样的方法避免了采样长度限制、采样间隔的不同与方向异质性而造成的误差问题,且当分块采样的面积越小时,测量结果越准确。

2)使用非接触地表的测量方式。非接触式的近景摄影测量法避免了在野外测量过程中测量板对地表微结构的破坏,在地表最真实的表面进行测量,降低了测量的误差。

3)在内业计算中可替代大量野外数据采集。在获取了密集点云数据后,可根据需要自由地调节分块采样的面积和数量,有效减少野外工作的时间和成本。

4)综合测量精度、时间和效率等评价指标,基于无人机影像生成点云数据的近景摄影测量方法更高效,更适用于地表土壤粗糙度的测量,且可根据目标单位面积或影像像元大小测量相应面积的土壤粗糙度。 （原文有删减）

【作者简介】20010090327@gznu.edu.cn。

【通讯作者】周忠发(1969-),男,教授,博士生导师,研究方向为喀斯特生态环境、GIS与遥感。Email: fa6897@163.com。

【】国家自然科学基金项目“喀斯特石漠化地区生态资产与区域贫困耦合机制研究”(41661088);贵州省高层次创新型人才培养计划项目“贵州省高层次创新型人才培养计划——‘百’层次人才”(黔科合平台人才〔2016〕5674);贵州省科学技术基金资助项目“高原山区特色农作物超低空遥感特征构建与识别方法研究”(黔科合基础-ZK[2021]一般194);国家自然科学基金面上项目“地块作物生长的光学与SAR遥感同步观测响应机制研究”(42071316)

【】张田, 周忠发, 王玲玉, 赵馨, 张文辉, 张淑, 王宇. 基于无人机点云数据土壤粗糙度测量方法[J]. 自然资源遥感, 2023, 35(1): 115-122.

转自：“测绘学术资讯”微信公众号

如有侵权，请联系本站删除！