-
林业遥感真实性检验与地面调查中[1-3],需要高精度的单木位置信息。样地单木位置的精确观测已有大量研究及方法探索。传统林业样地单木调查中,一般首先确立样地中心或角点从而通过拉皮尺计算每株树与中心点或角点相对距离计算出单木相对位置,此种方法通常较为费时,且在计算坐标时由于皮尺的震荡及人为观测的随机性,单木位置的坐标值误差分布具有随机性。使用全球导航卫星系统(Global Navigation Satellite System,GNSS)定位时,有使用手持GNSS接收机直接定位单木位置,也有探索实时动态载波相位差分技术(Real-time Kinematic,RTK)定位的方法,但由于每个GNSS接收机单独定位误差会在2~5 m,导致单木位置不能得到精准确立[4-7]。使用高精度全站仪(Total Station,TS)观测的坐标值系统误差在3 mm之内,但通常观测为树皮某处,也不能精确描述单木位置。
地基激光雷达(Terrestrial Laser Scanning,TLS)高精度的参数反演能力[8-13],与机载激光雷达(Airborne Laser Scanning,ALS)及人工测量数据结合可以提供精确的地面验证数据[14-16]。基于TLS数据使用Hough变换圆检测的方法在拟合单木胸径(Diameter at Breast Height,DBH)时展示出较高效率和更高精度[17],同时基于提取的直径可以获得单木DBH处圆心,此圆心可以更好地描述冠幅覆盖下的单木位置。而基于仿射变换空间校正的思想在全局校正空间点时具有很高的指导意义,通常在点校正过程中,只需要选择好校正点,使用一般仿射变换空间校正方法就可以将全局点的坐标校正。在匹配相对坐标与地理坐标时,基于影像配准的方法也有大量研究,如王书涵等[18]提出一种基于TS和快鸟影像协同的山地定位匹配方法,Sarah等[19]提出一种基于高分辨率影像的高精度冠幅制图方法。
本研究基于TLS提取的DBH圆心点,经过外业调查(In-situ Measurement,IM)的DBH数据验证筛选,将单木IM坐标位置以空间校正方法匹配至DBH中心点处。同时,基于影像林隙特征,以ALS提取的数字冠层高度模型(Canopy Height Model,CHM)产品作为基准,将匹配的IM相对坐标单木位置再匹配至空间地理坐标系下,为大面积遥感精确验证提供可靠基础单木位置数据。
-
研究区位于黑龙江佳木斯孟家岗林场(130°32′~130°52′ E,46°20′~46°30′ N),平均海拔250 m,以低山丘陵为主。地势西南低,东北高,坡度在10°~20°之间。林场主要森林类型为人工林,主要树种长白落叶松(Larix olgensis Henry),占人工林面积的60.3%。
实验中激光雷达数据分别选用2017年6月采集的TLS数据及同期飞行的ALS数据。其中,TLS采用的是Trimble TX8三维激光扫描仪,此设备测距噪声小于2 mm,最大测量距离120 m,激光波长为1 500 nm,样地扫描为一个中心站加4个边缘站的设置,并均匀布设10个圆形靶球及7个黑白靶标;ALS数据采集于中国林业科学研究院资源信息研究所研制的CAF-LiCHy机载遥感系统[20],此系统搭载了Riegl公司的全波形激光雷达扫描仪LMS-Q680i,能够提供高精度高密度离散点云数据,孟家岗实验飞行绝对高度1 000 m左右,航速50 m·s−1,脉冲重复频率300 KHz,航带旁向重叠率约62%,研究样地范围点云数据密度10 pts·m−2以上。同时为研究需要,分别对TLS及ALS数据生产了0.25 m空间分辨率的CHM产品。
针对研究内容,选择了高、中、低3个不同株数密度人工落叶松林样地。不同株数密度样地TLS扫描场景图如图1所示。
TLS及ALS数据采集的同时对4 cm以上径阶单木做了每木检尺,包括同期实测的每木DBH、树高,并记录包括枯死及弯曲等的单木生长状态。原始单木位置沿用2012年皮尺及罗盘等测定IM结果。使用索佳免棱镜CX全站仪在样地中心均匀测定约1/3单木DBH处树皮位置并记录树号,同时测定TLS的7个黑白靶标点以将TLS坐标系平面坐标匹配至TS坐标系下,此TS的测距精度在200 m测距范围内可达到2.2 mm (2 mm固定误差加0.2 mm比例误差)。如图2给出了其中一个样地TS测点及靶标点分布图。
从图2可以看出,样地中靶标设置基本均匀分布在样地中心围绕TLS中心站。TS测点基本均匀分布在样地中的设置有助于降低后续空间位置校正中的全局单木匹配误差。
表1给出了3个不同株数密度样地的概况,表中3行自上而下分别对应高中低3种样地密度。
株数密度
Stem density/(株·hm−2)树高(最小~最大/平均)
Tree height (Min~Max/Mean)/m样地大小
Sample plot size/(m×m)胸径(最小~最大/平均)
Diameter(Min~Max/Mean)/cm1 445 11.3~25.5/21.31 30×30 8.9~28.3/18.05 1 112 5.4~26.3/20.36 30×30 6.8~28.9/15.60 778 7.7~18.7/13.55 30×30 4.2~24.4/11.10 Table 1. The profile of sample plots
Accurate Measurement of Individual Tree Position Based on DBH Extraction of Terrestrial Laser Scanning
- Received Date: 2019-09-25
- Accepted Date: 2020-05-05
- Available Online: 2020-08-01
Abstract: