留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于无人机可见光影像的毛竹林郁闭度估测方法

王雨阳 王懿祥 李明哲 梁丹

王雨阳, 王懿祥, 李明哲, 梁丹. 基于无人机可见光影像的毛竹林郁闭度估测方法[J]. 浙江农林大学学报, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
引用本文: 王雨阳, 王懿祥, 李明哲, 梁丹. 基于无人机可见光影像的毛竹林郁闭度估测方法[J]. 浙江农林大学学报, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
WANG Yuyang, WANG Yixiang, LI Mingzhe, LIANG Dan. Estimation method of Phyllostachys edulis forest canopy density based on UAV visible image[J]. Journal of Zhejiang A&F University, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
Citation: WANG Yuyang, WANG Yixiang, LI Mingzhe, LIANG Dan. Estimation method of Phyllostachys edulis forest canopy density based on UAV visible image[J]. Journal of Zhejiang A&F University, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576

基于无人机可见光影像的毛竹林郁闭度估测方法

doi: 10.11833/j.issn.2095-0756.20210576
基金项目: 浙江省林业局与中国林业科学研究院合作项目(2019SY06)
详细信息
    作者简介: 王雨阳(ORCID: 0000-0002-9496-0891),从事无人机遥感等研究。E-mail: 2804649121@qq.com
    通信作者: 梁丹(ORCID: 0000-0002-6713-122X),讲师,博士,从事测量误差理论与无人机遥感研究。E-mail: liangdan812345@163.com
  • 中图分类号: S758.5

Estimation method of Phyllostachys edulis forest canopy density based on UAV visible image

  • 摘要:   目的  由于毛竹Phyllostachys edulis的生长特点和经营特点,使得毛竹林郁闭度在毛竹林经营中尤为重要,只有保持适宜毛竹生长的郁闭度,才能提高毛竹生产力。研究无人机可见光影像的毛竹林郁闭度估测方法,可实现实时快速获取毛竹林的郁闭度。  方法  以普通旋翼无人机可见光毛竹林影像为研究对象,基于像元的阈值分类、像元的监督分类、多尺度分割的阈值分类、多尺度分割的监督分类等4种方法,选取不同钩梢和郁闭度的样地36个,利用现有软件和MATLAB编程,对各样地的毛竹林竹冠区域进行快速提取,进而估算林分郁闭度,对比目视解译的郁闭度真值计算各方法的估算精度,利用单因素方差分析比较4种方法在不同钩梢和不同郁闭度下估算郁闭度的表现。  结果  基于像元的阈值分类、基于像元的监督分类、基于多尺度分割的阈值分类、基于多尺度分割的监督分类等4种方法的郁闭度估算总体精度依次为91.81%、92.96%、93.47%、98.86%,郁闭度估测值绝对误差依次为0.038、0.030、0.024、0.004。钩梢和郁闭度等对提取结果没有显著影响。  结论  基于多尺度分割的监督分类方法总体精度最高,估算绝对误差最小,能够满足快速、准确提取并估测毛竹林林分郁闭度的要求,且适用于不同经营类型的毛竹林。图2表6参28
  • 图  1  双峰法最优阈值示意图

    Figure  1  Schematic diagram of the double peak threshold method

    图  2  不同方法提取竹冠效果对比示意图

    Figure  2  Comparison of the bamboo crown extraction results with different methods

    表  1  毛竹林样地基本情况

    Table  1.   Basic information of plot sites

    类型郁闭度样地/
    平均胸
    径/cm
    平均树
    高/m
    坡度/
    (°)
    坡向
    钩梢 0~0.7138.919.911.6北坡 
    钩梢 0.7~0.889.039.416.1西北坡
    钩梢 0.8~1.038.659.213.3北坡 
    未钩梢128.589.617.1东南坡
    下载: 导出CSV

    表  2  不同毛竹林的最优参数组合

    Table  2.   Optimal parameter combination of different Moso bamboo stands     

    类型郁闭度分割尺度形状因子紧致度因子
    0~0.7290.30.5
    钩梢 0.7~0.8310.30.5
    0.8~1.0290.30.5
    未钩梢410.30.5
    下载: 导出CSV

    表  3  不同方法的毛竹林竹冠提取精度及郁闭度误差对比

    Table  3.   Bamboo crown extraction accuracy and canopy density error of different methods

    方法总体精度/%生产者精度/%用户精度/%郁闭度误差
    TP91.81±3.08 c93.34±3.37 c96.24±2.03 b0.038±0.026 c
    SP92.96±3.66 bc95.47±3.29 b95.64±2.64 b0.030±0.026 ab
    TM93.47±2.53 b96.10±2.24 b95.57±2.57 b0.024±0.018 b
    SM98.86±0.53 a99.15±0.40 a99.36±0.53 a0.004±0.003 a
      说明:同列不同小写字母表示差异显著(P<0.05)
    下载: 导出CSV

    表  4  4种方法下钩梢和未钩梢毛竹林竹冠区域的提取精度

    Table  4.   Extraction accuracy of the four methods with truncation and non-truncation

    类型方法总体精度/%生产者精度/%用户精度/%
    钩梢 TP92.82±2.57 Ab93.85±3.11 Ac95.96±2.32 Ab
    SP93.95±2.74 Ab96.67±1.39 Ab95.00±2.98 Ab
    TM94.05±2.25 Ab96.15±2.04 Ab95.49±2.90 Ab
    SM98.9±0.59 Aa99.14±0.45 Aa99.27±0.61 Aa
    未钩梢TP89.80±3.12 Bb92.33±3.77 Ac96.81±1.15 Ab
    SP90.97±4.52 Bb93.07±4.57 Bc96.92±0.98 Bb
    TM92.30±2.77 Bb96.01±2.70 Ab95.73±1.81 Ac
    SM98.78±0.37 Aa99.16±0.29 Aa99.55±0.22 Aa
      说明:同列不同小写字母表示相同林分类型在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同林分类型之间差异显著(P<0.05)
    下载: 导出CSV

    表  5  4种方法下钩梢和未钩梢毛竹林的郁闭度误差

    Table  5.   Canopy density error of the four methods with truncation and non-truncation

    类型方法郁闭度误差类型方法郁闭度误差
    钩梢TP0.029±0.021 Ab未钩梢TP0.052±0.029 Bc
    SP0.023±0.015 AbSP0.042±0.039 Bbc
    TM0.022±0.016 AbTM0.026±0.023 Ab
    SM0.003±0.003 AaSM0.004±0.003 Aa
      说明:同列不同小写字母表示相同林分类型在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同林分类型之间差异显著(P<0.05)
    下载: 导出CSV

    表  6  不同方法下郁闭度对毛竹林竹冠提取精度的影响

    Table  6.   Effect of canopy density on the extraction accuracy under different methods

    郁闭度方法总体精度/%生产者精度/%用户精度/%
    0~0.7TP93.70±2.44 Ab95.04±2.50 Ab95.26±2.69 Ab
    SP94.98±2.11 Ab97.02±1.37 Ab95.33±2.90 Ab
    TM93.89±2.66 Ab96.15±2.52 Ab94.70±3.32 Ab
    SM98.82±0.69 Aa99.11±0.53 Aa99.09±0.68 Aa
    0.7~0.8TP92.21±2.65 Ab92.43±3.83 Ac97.03±1.68 Ab
    SP93.62±4.31 Ab96.97±1.43 Ab95.22±4.78 Ab
    TM92.71±1.40 Ab96.11±1.43 Ab95.19±1.90 Ab
    SM98.88±0.48 Aa99.08±0.32 Aa99.42±0.48 Aa
    0.8~1.0TP90.68±1.46 Ac92.46±1.39 Ac96.14±0.74 Ab
    SP93.62±4.31 Ab96.97±1.43 Ab95.22±4.78 Ab
    TM92.71±1.40 Ab96.11±1.43 Ab95.19±1.90 Ab
    SM99.28±0.33 Aa99.48±0.22 Aa99.66±0.33 Aa
      说明:同列不同小写字母表示相同郁闭度在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同郁闭度之间差异显著(P<0.05)
    下载: 导出CSV
  • [1] HERBERT T J. Area projections of fisheye photographic lenses [J]. Agric For Meteorol, 1987, 39(2/3): 215 − 223.
    [2] 濮毅涵, 徐丹丹, 王浩斌. 基于数码相片的林冠郁闭度提取方法研究[J]. 林业资源管理, 2020(6): 153 − 160.

    PU Yihan, XU Dandan, WANG Haobin. An approach on estimating canopy closure via digital images [J]. For Resour Manage, 2020(6): 153 − 160.
    [3] 高云飞, 李智广, 杨胜天, 等. 基于SPOT 5影像的郁闭度反演方法[J]. 水土保持研究, 2012, 19(2): 267 − 270.

    GAO Yunfei, LI Zhiguang, YANG Shengtian, et al. Study on canopy density retrieval method from SPOT 5 [J]. Res Soil Water Conserv, 2012, 19(2): 267 − 270.
    [4] 李擎, 王振锡, 王雅佩, 等. 基于GF-2号遥感影像的天山云杉林郁闭度估测研究[J]. 中南林业科技大学学报, 2019, 39(8): 48 − 54.

    LI Qing, WANG Zhenxi, WANG Yapei, et al. Study on canopy density inversion of Picea schrenkiana forest based on GF-2 remote sensing image [J]. J Cent South Univ For Technol, 2019, 39(8): 48 − 54.
    [5] 李蕴雅. UAV/RS 3D像对森林信息提取方法研究[D]. 北京: 北京林业大学, 2016.

    LI Yunya. Study on the Method of Extracting Forest Information from 3D UAV/RS[D]. Beijing: Beijing Forestry University, 2016.
    [6] JONATHAN L, MARC P D, STÉPHANIE B, et al. A photogrammetric workflow for the creation of a forest canopy height model from small unmanned aerial system imagery [J]. Forests, 2013, 4(4): 922 − 944.
    [7] DANDOIS J P, ELLIS E C. High spatial resolution three-dimensional mapping of vegetation spectral dynamics using computer vision [J]. Remote Sensing Environ, 2013, 136(5): 259 − 276.
    [8] 李丹, 张俊杰, 赵梦溪. 基于FCM和分水岭算法的无人机影像中林分因子提取[J]. 林业科学, 2019, 55(5): 180 − 187.

    LI Dan, ZHANG Junjie, ZHAO Mengxi. Extraction of stand factors in UAV image based on FCM and watershed algorithm [J]. Sci Silv Sin, 2019, 55(5): 180 − 187.
    [9] 苏迪, 高心丹. 基于无人机航测数据的森林郁闭度和蓄积量估测[J]. 林业工程学报, 2020, 5(1): 156 − 163.

    SU Di, GAO Xindan. Estimation of forest canopy density and stock volume based on UAV aerial survey Data [J]. J For Eng, 2020, 5(1): 156 − 163.
    [10] 汪霖. 基于无人机高分影像的森林参数估测方法[D]. 南京: 南京林业大学, 2020.

    WANG Lin. Estimation Method of Forest Parameters based on UAV High Resolution Image[D]. Nanjing: Nanjing Forestry University, 2020.
    [11] 申景昕, 范少辉, 刘广路, 等. 毛竹林采伐林窗近地层温度时空分布特征[J]. 生态学杂志, 2020, 39(11): 3549 − 3557.

    SHEN Jingxin, FAN Shaohui, LIU Guanglu, et al. Spatiotemporal distribution characteristics of temperature on the surface layer of cutting gap of Phyllostachys edulis forest [J]. Chin J Ecol, 2020, 39(11): 3549 − 3557.
    [12] 宋庆妮, 杨清培, 刘骏, 等. 毛竹扩张对常绿阔叶林土壤氮素矿化及有效性的影响[J]. 应用生态学报, 2013, 24(2): 338 − 344.

    SONG Qingni, YANG Qingpei, LIU Jun, et al. Effects of Phyllostachys edulis expansion on soil nitrogen mineralization and its availability in evergreen broadleaf forest [J]. Chin J Appl Ecol, 2013, 24(2): 338 − 344.
    [13] 杨清培, 杨光耀, 宋庆妮, 等. 竹子扩张生态学研究: 过程、后效与机制[J]. 植物生态学报, 2015, 39(1): 110 − 124.

    YANG Qingpei, YANG Guangyao, SONG Qingni, et al. Ecological studies on bamboo expansion: process, consequence and mechanism [J]. Chin J Plant Ecol, 2015, 39(1): 110 − 124.
    [14] 孙晓艳, 杜华强, 韩凝, 等. 面向对象多尺度分割的SPOT 5影像毛竹林专题信息提取[J]. 林业科学, 2013, 49(10): 80 − 87.

    SUN Xiaoyan, DU Huaqiang, HAN Ning, et al. Multi-Scale segmentation, object-based extraction of Moso bamboo forest from SPOT 5 imagery [J]. Sci Silv Sin, 2013, 49(10): 80 − 87.
    [15] 赵敏, 赵银娣. 面向对象的多特征分级CVA遥感影像变化检测[J]. 遥感学报, 2018, 22(1): 119 − 131.

    ZHAO Min, ZHAO Yindi. Object-oriented and muti-feature hierarchical change detection based on CVA for high-resolution remote sensing imagery [J]. J Remote Sensing, 2018, 22(1): 119 − 131.
    [16] 黄慧萍, 吴炳方. 地物提取的多尺度特征遥感应用分析[J]. 遥感技术与应用, 2003(5): 276 − 281.

    HUANG Huiping, WU Bingfang. Analysis of the mutli-scale characteristiscs with objects extraction [J]. Remote Sensing Technol Appl, 2003(5): 276 − 281.
    [17] ADDINK E A, JONG S, PEBESMA E J. The importance of scale in object-based mapping of vegetation parameters with hyperspectral imagery [J]. Photogramm Eng Remote Sensing, 2007, 73(8): 905 − 912.
    [18] DRĂGUŢ L, CSILLIK O, EISANK C, et al. Automated parameterisation for multi-scale image segmentation on multiple layers [J]. ISPRS J Photogramm Remote Sensing, 2014, 88: 119 − 127.
    [19] 刘金丽, 陈钊, 高金萍, 等. 高分影像树种分类的最优分割尺度确定方法[J]. 林业科学, 2019, 55(11): 95 − 104.

    LIU Jinli, CHEN Zhao, GAO Jinping, et al. Research on the method of determining the optimal segmentation scale for tree species classification of high-resolution image [J]. Sci Silv Sin, 2019, 55(11): 95 − 104.
    [20] PUREVDORJ T, TATEISHI R, ISHIYAMA T, et al. Relationships between percent vegetation cover and vegetation indices [J]. Int J Remote Sensing, 1998, 19(18): 3519 − 3535.
    [21] 祁媛, 徐伟诚, 王林琳, 等. 基于无人机遥感影像的沙糖橘果树提取方法研究[J]. 华南农业大学学报, 2020, 41(6): 126 − 133.

    QI Yuan, XU Weicheng, WANG Linlin, et al. Study on the extraction method of sugar tangerine fruit trees based on UAV remote sensing images [J]. J South China Agric Univ, 2020, 41(6): 126 − 133.
    [22] MEYER G E, NETO J C. Verification of color vegetation indices for automated crop imaging applications [J]. Comput Electron Agric, 2008, 63(2): 282 − 293.
    [23] 汪小钦, 王苗苗, 王绍强, 等. 基于可见光波段无人机遥感的植被信息提取[J]. 农业工程学报, 2015, 31(5): 152 − 159.

    WANG Xiaoqin, WANG Miaomiao, WANG Shaoqiang, et al. Extraction of vegetation information from visible unmanned aerial vehicle images [J]. Trans Chin Soc Agric Eng, 2015, 31(5): 152 − 159.
    [24] 陈时跃. 基于无人机可见光影像的郁闭山地杉木人工林的单木与林分结构参数提取[D]. 杭州: 浙江农林大学, 2020.

    CHEN Shiyue. Extraction of Individual Tree and Stand Structure Parameters based on RGB Images UVA in Closed Canopy Mountainous Forest [D]. Hangzhou: Zhejiang A&F University, 2020.
    [25] 汪霖, 李明阳, 方子涵, 等. 基于无人机数据的人工林森林参数估测[J]. 林业资源管理, 2019(5): 61 − 67.

    WANG Lin, LI Mingyang, FANG Zihan, et al. Plantation forest parameter estimation based on UAV data [J]. For Resour Manage, 2019(5): 61 − 67.
    [26] 智献坡, 李明泽, 于颖, 等. 应用4-scale模型对人工针叶林郁闭度的反演[J]. 东北林业大学学报, 2020, 48(12): 5 − 11.

    ZHI Xianpo, LI Mingze, YU Ying, et al. Artficial coniferous forest remote sensing by using geometric optical model [J]. J Northeast For Univ, 2020, 48(12): 5 − 11.
    [27] 刘赛赛, 陈冬花, 栗旭升, 等. 基于高分一号PMS的新疆落叶松林分郁闭度遥感定量估测[J]. 西北农林科技大学学报(自然科学版), 2020, 48(7): 57 − 66.

    LIU Saisai, CHEN Donghua, LI Xusheng, et al. Quantitative estimation of stand closure density of Larix sibirica by remote sensing based on GF-1 PMS [J]. J Northwest A&F Univ Nat Sci Ed, 2020, 48(7): 57 − 66.
    [28] 严羽. 基于无人机的林分郁闭度和树高估测研究[D]. 杭州: 浙江农林大学, 2019.

    YAN Yu. Estimation of Canopy Density and Tree Height Based on UAV[D]. Hangzhou: Zhejiang A&F University, 2019.
  • [1] 王佳雨, 朱玲姣, 黄程鹏, 姜培坤, 查强威, 陈林海.  硅肥和生物质炭添加对毛竹林土壤活性硅组分的影响 . 浙江农林大学学报, 2024, 41(3): 496-505. doi: 10.11833/j.issn.2095-0756.20230366
    [2] 曹善郅, 周家树, 张少博, 姚易寒, 刘娟, 李永夫.  生物质炭基尿素和普通尿素对毛竹林土壤氧化亚氮通量的影响 . 浙江农林大学学报, 2023, 40(1): 135-144. doi: 10.11833/j.issn.2095-0756.20220254
    [3] 陈旭, 刘宗悦, 徐钧杰, 祁祥斌, 余树全.  天目山毛竹林皆伐后群落的恢复特征 . 浙江农林大学学报, 2022, 39(4): 705-716. doi: 10.11833/j.issn.2095-0756.20210595
    [4] 杨帆, 汤孟平.  毛竹林立地与结构的关系及其对生物量的影响 . 浙江农林大学学报, 2020, 37(5): 823-832. doi: 10.11833/j.issn.2095-0756.20190572
    [5] 王铮屹, 戴其林, 柏宬, 陈涵, 库伟鹏, 赵明水, 余树全.  天目山皆伐毛竹林自然更新群落类型与多样性分析 . 浙江农林大学学报, 2020, 37(4): 710-719. doi: 10.11833/j.issn.2095-0756.20190472
    [6] 谢巧雅, 余坤勇, 邓洋波, 刘健, 范华栋, 林同舟.  杉木人工林冠层高度无人机遥感估测 . 浙江农林大学学报, 2019, 36(2): 335-342. doi: 10.11833/j.issn.2095-0756.2019.02.015
    [7] 贾鹏刚, 夏凯, 董晨, 冯海林, 杨垠晖.  基于无人机影像的银杏单木胸径预估方法 . 浙江农林大学学报, 2019, 36(4): 757-763. doi: 10.11833/j.issn.2095-0756.2019.04.016
    [8] 赵赛赛, 汤孟平, 唐思嘉, 张军, 李岚, 庞春梅, 赵明水.  毛竹林分可视化研究 . 浙江农林大学学报, 2016, 33(5): 826-833. doi: 10.11833/j.issn.2095-0756.2016.05.014
    [9] 孟海月, 刘强, 吴伟光.  不同经营类型毛竹林经营效益及固碳能力分析 . 浙江农林大学学报, 2014, 31(6): 959-964. doi: 10.11833/j.issn.2095-0756.2014.06.020
    [10] 张丽景, 葛宏立.  利用MODIS数据估测毛竹林总初级生产力 . 浙江农林大学学报, 2014, 31(2): 178-184. doi: 10.11833/j.issn.2095-0756.2014.02.003
    [11] 赵晓, 吕玉龙, 王聪, 李亚丹, 杜华强.  毛竹林叶面积指数和郁闭度空间分布协同克里格估算 . 浙江农林大学学报, 2014, 31(4): 560-569. doi: 10.11833/j.issn.2095-0756.2014.04.011
    [12] 蔡彦彬, 宋照亮, 姜培坤.  岩性对毛竹林土壤硅形态的影响 . 浙江农林大学学报, 2013, 30(6): 799-804. doi: 10.11833/j.issn.2095-0756.2013.06.001
    [13] 陈孝丑, 刘广路, 范少辉, 官凤英, 苏文会, 黄金华.  连续施肥对毛竹林生长特征及生物量空间构型的影响 . 浙江农林大学学报, 2012, 29(1): 52-57. doi: 10.11833/j.issn.2095-0756.2012.01.010
    [14] 马瑞升, 杨斌, 张利辉, 刘志平.  微型无人机林火监测系统的设计与实现 . 浙江农林大学学报, 2012, 29(5): 783-789. doi: 10.11833/j.issn.2095-0756.2012.05.023
    [15] 崔瑞蕊, 杜华强, 周国模, 徐小军, 董德进, 吕玉龙.  近30 a安吉县毛竹林动态遥感监测及碳储量变化 . 浙江农林大学学报, 2011, 28(3): 422-431. doi: 10.11833/j.issn.2095-0756.2011.03.012
    [16] 叶耿平, 刘娟, 姜培坤, 周国模, 吴家森.  集约经营措施对毛竹林生长季土壤呼吸的影响 . 浙江农林大学学报, 2011, 28(1): 18-25. doi: 10.11833/j.issn.2095-0756.2011.01.004
    [17] 王冬云, 张卓文, 苏开君, 王光, 雷云飞, 林明磊, 张培, 钟庸.  广州流溪河流域毛竹林的水文生态效应 . 浙江农林大学学报, 2008, 25(1): 37-41.
    [18] 丁丽霞, 王祖良, 周国模, 杜晴洲.  天目山国家级自然保护区毛竹林扩张遥感监测 . 浙江农林大学学报, 2006, 23(3): 297-300.
    [19] 高志勤, 傅懋毅.  不同毛竹林土壤碳氮养分的季节变化特征 . 浙江农林大学学报, 2006, 23(3): 248-254.
    [20] 何福基, 吴明安, 倪荣新, 谢正成, 张建忠.  杉木种子园郁闭度对种子产量的影晌 . 浙江农林大学学报, 1995, 12(3): 311-315.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20210576

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2022/5/981

图(2) / 表(6)
计量
  • 文章访问数:  594
  • HTML全文浏览量:  268
  • PDF下载量:  91
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-08-19
  • 修回日期:  2022-05-20
  • 录用日期:  2022-06-15
  • 网络出版日期:  2022-09-22
  • 刊出日期:  2022-10-20

基于无人机可见光影像的毛竹林郁闭度估测方法

doi: 10.11833/j.issn.2095-0756.20210576
    基金项目:  浙江省林业局与中国林业科学研究院合作项目(2019SY06)
    作者简介:

    王雨阳(ORCID: 0000-0002-9496-0891),从事无人机遥感等研究。E-mail: 2804649121@qq.com

    通信作者: 梁丹(ORCID: 0000-0002-6713-122X),讲师,博士,从事测量误差理论与无人机遥感研究。E-mail: liangdan812345@163.com
  • 中图分类号: S758.5

摘要:   目的  由于毛竹Phyllostachys edulis的生长特点和经营特点,使得毛竹林郁闭度在毛竹林经营中尤为重要,只有保持适宜毛竹生长的郁闭度,才能提高毛竹生产力。研究无人机可见光影像的毛竹林郁闭度估测方法,可实现实时快速获取毛竹林的郁闭度。  方法  以普通旋翼无人机可见光毛竹林影像为研究对象,基于像元的阈值分类、像元的监督分类、多尺度分割的阈值分类、多尺度分割的监督分类等4种方法,选取不同钩梢和郁闭度的样地36个,利用现有软件和MATLAB编程,对各样地的毛竹林竹冠区域进行快速提取,进而估算林分郁闭度,对比目视解译的郁闭度真值计算各方法的估算精度,利用单因素方差分析比较4种方法在不同钩梢和不同郁闭度下估算郁闭度的表现。  结果  基于像元的阈值分类、基于像元的监督分类、基于多尺度分割的阈值分类、基于多尺度分割的监督分类等4种方法的郁闭度估算总体精度依次为91.81%、92.96%、93.47%、98.86%,郁闭度估测值绝对误差依次为0.038、0.030、0.024、0.004。钩梢和郁闭度等对提取结果没有显著影响。  结论  基于多尺度分割的监督分类方法总体精度最高,估算绝对误差最小,能够满足快速、准确提取并估测毛竹林林分郁闭度的要求,且适用于不同经营类型的毛竹林。图2表6参28

English Abstract

王雨阳, 王懿祥, 李明哲, 梁丹. 基于无人机可见光影像的毛竹林郁闭度估测方法[J]. 浙江农林大学学报, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
引用本文: 王雨阳, 王懿祥, 李明哲, 梁丹. 基于无人机可见光影像的毛竹林郁闭度估测方法[J]. 浙江农林大学学报, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
WANG Yuyang, WANG Yixiang, LI Mingzhe, LIANG Dan. Estimation method of Phyllostachys edulis forest canopy density based on UAV visible image[J]. Journal of Zhejiang A&F University, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
Citation: WANG Yuyang, WANG Yixiang, LI Mingzhe, LIANG Dan. Estimation method of Phyllostachys edulis forest canopy density based on UAV visible image[J]. Journal of Zhejiang A&F University, 2022, 39(5): 981-988. doi: 10.11833/j.issn.2095-0756.20210576
  • 林分郁闭度是指林木树冠垂直投影面积占林地总面积的比例,是反映林分密度的重要指标,在森林资源调查中尤为重要。在传统的森林资源调查中,获取林分郁闭度的主要方法为抬头望法,该方法耗时耗力,获取的数据精度较低,且只适用于面积小、地势缓的林地,难以在地理环境复杂的大区域内应用。随着数字成像技术的发展,HERBERT[1]提出利用鱼眼相机拍摄全天空相片进行郁闭度估测,扩大了拍摄的林冠范围;濮毅涵等[2]基于普通可见光照片中的树叶、树干和天空进行分类,将估测结果与抬头望法结果进行模型比较,得到决定系数(R2)为0.77,但该方法适用于地形平坦的林区,在地形复杂的山地难以应用。为克服地形因子的限制,不少学者基于卫星遥感影像对乔木林林分郁闭度的提取进行了研究。高云飞等[3]基于SPOT 5影像,对各波段的遥感影像像元亮度值进行组合,建立反演模型,最终得到最佳组合模型的R2为0.66;李擎等[4]为提高模型精度,基于高分二号(GF-2)遥感影像,结合光谱信息、纹理特征和地形因子构建郁闭度估测模型,模型精度达89.82%。然而卫星遥感成本高、时效低、灵活性低[5],相比之下无人机可见光遥感在中小型遥感区域作业中能发挥更大的优势,满足动态森林资源监测的条件和需求,并提供更多可能性[6-8]。苏迪等[9]基于无人机可见光影像,利用主成分分析确定模型的主成分变量,建立郁闭度回归方程,估算精度为83.18%;汪霖[10]结合无人机可见光影像得到的数字正射影像(DOM)和冠层高度模型(CHM),使用阈值法对树冠区域进行提取并计算郁闭度,平均精度为92.93%。以上关于乔木林分郁闭度遥感估测的研究都取得了一定的成果,能够满足林业调查需求,然而,目前基于无人机可见光影像的毛竹林林分郁闭度估测研究为数不多。

    毛竹作为中国南方重要的笋、竹两用竹种,是最重要的森林资源之一。第九次全国森林清查结果显示:毛竹林面积为467 hm2,占竹林总面积的73%[11]。由于毛竹林的生长特点(扩鞭)[12-13]和经营特点(择伐),使得林分郁闭度在毛竹林经营中尤为重要,出笋率和采伐结果都将反映在郁闭度上,只有保持适宜毛竹生长的郁闭度,才能提高毛竹生产力。毛竹林多分布在山地丘陵地区,传统的实地调查工作量大、效率低、工作周期长、精确度低。随着无人机技术的发展,无人机林业遥感为实现低成本、高时效、高分辨率、高灵活性的动态森林资源监测提供了技术基础。基于无人机和图像识别技术的单木树冠提取已有一定进展。然而相比其他乔木林,毛竹林竹冠具有特殊性,分为钩梢与未钩梢2种形态,且竹冠鲜少呈单个状态,具有“重叠性”,目前还没有利用现有方法直接获取单株毛竹竹冠的研究。本研究对无人机可见光影像的毛竹林郁闭度估测方法进行了分析,以期为快速获取毛竹林的郁闭度提供参考。

    • 研究区位于浙江省湖州市安吉县(30°23′~30°53′N,119°14′~119°53′E),该区地处长江三角洲腹地,属于亚热带季风气候,植被类型为亚热带常绿阔叶林,年均气温为12.2~15.6 ℃,年均降水量为1 100~1 900 mm。研究区内光照充足,雨量充沛,竹资源丰富,竹林总面积达757 km2,占森林面积56.47%,其中毛竹占79.30%,被誉为“中国竹乡”。

    • 根据是否钩梢和郁闭度2个因素选设样地,共调查36个10 m×10 m的样地,其中已钩梢毛竹林样地24个,未钩梢毛竹林样地12个。利用研究区数字正射影像目视解译获得毛竹林的郁闭度,未钩梢样地为抛荒毛竹林,处于无人经营状态,郁闭度均大于0.9,故在本研究中不再细分郁闭度等级。样地基本情况见表1

      表 1  毛竹林样地基本情况

      Table 1.  Basic information of plot sites

      类型郁闭度样地/
      平均胸
      径/cm
      平均树
      高/m
      坡度/
      (°)
      坡向
      钩梢 0~0.7138.919.911.6北坡 
      钩梢 0.7~0.889.039.416.1西北坡
      钩梢 0.8~1.038.659.213.3北坡 
      未钩梢128.589.617.1东南坡
    • 数据采集所用无人机遥感系统为大疆Phantom 4 RTK,飞行系统包括3部分:飞行器、云台相机和遥控器。飞行器提供多方位视觉系统,可稳定飞行和精准悬停,进行航点飞行作业;机身装备机载D-RTK,可提供厘米级高精度准确定位,实现更精准的测绘作业。Phantom 4 RTK配备24 mm广角相机、高精度防抖云台及图像传感器、机械快门,确保成像效果。遥控器的主要功能包括作业设计、远距离信号传输和航拍高清画面实时监测等。

      2020年11月18日14:00—15:00,天气晴朗、无风、无云,采用大疆Phantom 4 RTK搭载光学相机对选定研究区域进行拍摄。无人机飞行高度为160 m,飞行范围东西走向为222.82 m,南北走向为106.28 m,航向重叠率为90%,旁向重叠率为80%,飞行路线设计为“S”型路线。拍摄共得原始照片142张,空间分辨率为3.36 cm,照片共3波段:红波段(R)、绿波段(G)、蓝波段(B)。

    • 使用大疆智图软件对无人机所获取照片进行图像处理,通过三维重建技术,生成带有空间参考信息的数字表面模型(DSM)和色调自然、无明显接痕的正射影像图。图像处理的主要流程包括:①数据预处理,对照片进行校验、筛选、对齐,剔除一些无法进行后续合成的照片;②空三加密,通过空三计算得到密集点云和纹理;③生成网格,得到带空间参考信息的数字表面模型和数字正射影像。

    • 在毛竹林高分辨率遥感影像中,由于影像信息太多、林间空隙太杂,经常会产生“同物异谱”和“同谱异物”现象,而在分类前执行分割能有效解决该问题[14-16]

      多尺度分割是从任意一像元开始,自下而上进行合并的一种分割手段,其参数设置包括分割尺度、形状因子和紧致度因子。其中形状因子的取值为[0,1],它包含了形状和光谱2个方面的意义,取值越大,表示在同质性标准中形状因子所占比例越高,光谱因子所占比例就越低,反之亦然。紧致度因子包含紧致度和光滑度,两者权重之和也为“1”。分割尺度作为各参数中最重要的一项,决定了分割所得的多边形对象内部的异质性,分割尺度越大,所生成的对象面积就越大,数目就越少,反之亦然[17]。本研究采用eCognition软件中的Estimation of Scale Paramater 2(ESP2)插件来确定分割尺度,通过设置一系列的参数进行迭代,计算局部方差以及局部方差变化率(ROC)[18-19]ROC=[VL−(VL−1)]/(VL−1)×100。其中:VL表示尺度为L时分割结果中所有对象局部方差值的方差;VL−1表示尺度为L−1时分割结果中所有对象局部方差值的方差。

      本研究对不同类型的毛竹林进行不同参数组合实验,得到不同毛竹林的最优参数组合(表2)。

      表 2  不同毛竹林的最优参数组合

      Table 2.  Optimal parameter combination of different Moso bamboo stands     

      类型郁闭度分割尺度形状因子紧致度因子
      0~0.7290.30.5
      钩梢 0.7~0.8310.30.5
      0.8~1.0290.30.5
      未钩梢410.30.5
    • ①基于像元的阈值分类(TP)。通过设定不同的特征阈值,将图像分为若干类。常用的特征值包括2类:直接来自于原始图像的灰度或彩色特征和由原始灰度或彩色值变换得到的特征值。基于无人机影像的毛竹林竹冠区域阈值分类的关键在于选取能有效区分竹冠与非竹冠的最优阈值。通过选取竹冠与非竹冠的样点,得到各样点的阈值,统计不同阈值所包含的像元数目,绘成具有明显峰值与谷值的曲线图[20-21]。2个峰值对应位置分别为竹冠与非竹冠的典型特征值,谷值为竹冠与非竹冠交界处对应的特征值,也就是最优阈值(图1)。

      图  1  双峰法最优阈值示意图

      Figure 1.  Schematic diagram of the double peak threshold method

      特征值的选取对提取精度至关重要,本研究使用ENVI软件对数字正射影像中的红、绿、蓝波段进行主成分分析,得到红光和绿光波段包含98%以上的信息量,故选取由红、绿波段变换得到的归一化绿红差异指数(INGRD)为特征值[22-23],在MATLAB中对竹冠区域进行提取,表达式为INGRD=(GR)/(G+R)。其中R为红光波段,G为绿光波段。②基于像元的监督分类(SP)。根据已知训练区提供的样本,选取并求出特征参数作为决策规则,建立判别函数对各类图像进行分类的一种方法。本研究建立竹冠区与非竹冠区2种样本,并在ENVI软件中选择最大似然法作为分类算法进行竹冠区域提取,该方法根据训练样本的均值和方差,评价待分类像元和训练样本之间的相似性进而分类,可同时考虑2个以上的波段和类别。③基于多尺度分割的阈值分类(TM)。与传统基于像元的分类不同,基于多尺度分割的阈值分类是一种面向对象分类,其分类的基本单元为影像对象,基本内容分为影像分割与影像分类2个独立模块。结合了多尺度分割的分类方法,在精确分类的同时,很好地兼顾了地物的宏观尺度和微观特征,极大地消除了传统分类带来的“椒盐效应”。在分割的基础上,本研究使用INGRD对竹冠区域进行阈值提取。④基于多尺度分割的监督分类(SM)。该分类也是在多尺度分割的基础上,基于对象进行监督分类。本研究通过eCognition软件进行多尺度分割,并选择最邻近值分类法进行分类提取,该方法通过均匀地选择样本,建立训练集,统计样本的特征信息,通过计算未分类对象与样本之间的距离进行分类,并提取竹冠区域。

    • 郁闭度是指树冠的总垂直投影面积(m2)与该样方的总面积(m2)之比,反映林分的密度和林分光能利用程度,是抚育间伐的重要指标,郁闭度=总冠幅/样方总面积。

    • 在无人机影像的基础上,根据竹冠部分的亮度、纹理和阴影等特征,在ArcGIS软件中对各个小样地进行目视解译,并以目视解译的结果作为真值对4种方法进行精度评价[24]。分别用正确识别的竹冠总面积(AC,m2)、真实的竹冠总面积(AR,m2)、识别的竹冠面积(AD,m2)、正确识别的非竹冠面积(AN,m2)、样地面积(AS,m2)、真实郁闭度(RD)、分类结果得到的郁闭度(DC)等7个指标进行精度和误差的计算。总体精度(OA)=(AC+AN)/AS×100%;用户精度(UA)=AC/AD×100%;生产者精度(PA)=AC/AR×100%;郁闭度误差(DE)=|RDDC|。

    • 使用SPSS 25.0对各样地的总体精度、用户精度、生产者精度、郁闭度误差进行统计分析,采用单因素方差分析法(one-way ANOVA)进行显著性检验,采用LSD法进行多重比较分析(P=0.05),文中数据为平均值±标准差。

    • 图2可见:无论是在钩梢毛竹林还是未钩梢毛竹林中,基于像元的阈值分类法的结果存在较多错分和漏分的情况,且提取结果存在“碎片化”情况。基于像元的监督分类法的漏分情况较少,但存在较多的错分且提取结果呈“碎片化”;基于多尺度分割的阈值分类法提取的树冠存在较严重的错分和漏分情况;基于多尺度分割的监督分类法基本解决了树冠提取结果“碎片化”的问题,提高了所提取树冠的整体性,减少了错分和漏分的情况。从表3可以看出:基于多尺度分割的监督分类法的总体精度、生产者精度和用户精度都显著高于其他3种方法(P<0.05);基于多尺度分割的监督分类法郁闭度误差最小,为0.004,显著低于其他3种方法(P<0.05),表明使用基于多尺度分割的监督分类法提取毛竹林郁闭度的结果明显优于其他3种方法。

      图  2  不同方法提取竹冠效果对比示意图

      Figure 2.  Comparison of the bamboo crown extraction results with different methods

      表 3  不同方法的毛竹林竹冠提取精度及郁闭度误差对比

      Table 3.  Bamboo crown extraction accuracy and canopy density error of different methods

      方法总体精度/%生产者精度/%用户精度/%郁闭度误差
      TP91.81±3.08 c93.34±3.37 c96.24±2.03 b0.038±0.026 c
      SP92.96±3.66 bc95.47±3.29 b95.64±2.64 b0.030±0.026 ab
      TM93.47±2.53 b96.10±2.24 b95.57±2.57 b0.024±0.018 b
      SM98.86±0.53 a99.15±0.40 a99.36±0.53 a0.004±0.003 a
        说明:同列不同小写字母表示差异显著(P<0.05)
    • 表4所示:在基于多尺度分割的监督分类法中,是否钩梢对总体精度均没有显著影响,说明基于多尺度分割的监督分类法对钩梢和未钩梢毛竹林均适用;在其他3种方法中,钩梢毛竹林总体精度显著高于未钩梢毛竹林,说明这3种方法更适用于钩梢毛竹林。无论是钩梢林分,还是未钩梢林分,基于多尺度分割的监督分类法的总体精度均显著高于其他3种方法(P<0.05)。

      表 4  4种方法下钩梢和未钩梢毛竹林竹冠区域的提取精度

      Table 4.  Extraction accuracy of the four methods with truncation and non-truncation

      类型方法总体精度/%生产者精度/%用户精度/%
      钩梢 TP92.82±2.57 Ab93.85±3.11 Ac95.96±2.32 Ab
      SP93.95±2.74 Ab96.67±1.39 Ab95.00±2.98 Ab
      TM94.05±2.25 Ab96.15±2.04 Ab95.49±2.90 Ab
      SM98.9±0.59 Aa99.14±0.45 Aa99.27±0.61 Aa
      未钩梢TP89.80±3.12 Bb92.33±3.77 Ac96.81±1.15 Ab
      SP90.97±4.52 Bb93.07±4.57 Bc96.92±0.98 Bb
      TM92.30±2.77 Bb96.01±2.70 Ab95.73±1.81 Ac
      SM98.78±0.37 Aa99.16±0.29 Aa99.55±0.22 Aa
        说明:同列不同小写字母表示相同林分类型在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同林分类型之间差异显著(P<0.05)

      表5可见:基于多尺度分割的监督分类法所得到的郁闭度误差最小,在钩梢林分中仅为0.003,在未钩梢林分中仅为0.004,显著低于其他3种方法(P<0.05)。在基于多尺度分割的监督分类法和基于多尺度分割的阈值分类法中,是否钩梢对郁闭度误差没有显著影响;在基于像元的监督分类法和基于像元的阈值分类法中,钩梢毛竹林的郁闭度误差显著低于未钩梢毛竹林(P<0.05)。

      表 5  4种方法下钩梢和未钩梢毛竹林的郁闭度误差

      Table 5.  Canopy density error of the four methods with truncation and non-truncation

      类型方法郁闭度误差类型方法郁闭度误差
      钩梢TP0.029±0.021 Ab未钩梢TP0.052±0.029 Bc
      SP0.023±0.015 AbSP0.042±0.039 Bbc
      TM0.022±0.016 AbTM0.026±0.023 Ab
      SM0.003±0.003 AaSM0.004±0.003 Aa
        说明:同列不同小写字母表示相同林分类型在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同林分类型之间差异显著(P<0.05)
    • 表6可见:4种方法对不同郁闭度毛竹林的总体精度均没有显著影响,说明这4种方法在不同郁闭度的毛竹林中均适用。

      表 6  不同方法下郁闭度对毛竹林竹冠提取精度的影响

      Table 6.  Effect of canopy density on the extraction accuracy under different methods

      郁闭度方法总体精度/%生产者精度/%用户精度/%
      0~0.7TP93.70±2.44 Ab95.04±2.50 Ab95.26±2.69 Ab
      SP94.98±2.11 Ab97.02±1.37 Ab95.33±2.90 Ab
      TM93.89±2.66 Ab96.15±2.52 Ab94.70±3.32 Ab
      SM98.82±0.69 Aa99.11±0.53 Aa99.09±0.68 Aa
      0.7~0.8TP92.21±2.65 Ab92.43±3.83 Ac97.03±1.68 Ab
      SP93.62±4.31 Ab96.97±1.43 Ab95.22±4.78 Ab
      TM92.71±1.40 Ab96.11±1.43 Ab95.19±1.90 Ab
      SM98.88±0.48 Aa99.08±0.32 Aa99.42±0.48 Aa
      0.8~1.0TP90.68±1.46 Ac92.46±1.39 Ac96.14±0.74 Ab
      SP93.62±4.31 Ab96.97±1.43 Ab95.22±4.78 Ab
      TM92.71±1.40 Ab96.11±1.43 Ab95.19±1.90 Ab
      SM99.28±0.33 Aa99.48±0.22 Aa99.66±0.33 Aa
        说明:同列不同小写字母表示相同郁闭度在不同方法之     间差异显著(P<0.05);同列不同大写字母表示相同     方法在不同郁闭度之间差异显著(P<0.05)
    • 本研究4种方法提取的毛竹林竹冠区域精度均在90%以上,而汪霖等[25]通过提取单个树冠面积,得到提取树冠区域的总体精度为93.09%,高于本研究基于像元的阈值分类法和基于像元的监督分类法,但低于基于多尺度分割的阈值分类法和基于多尺度分割的监督分类法,证明在分类前执行分割能有效提高提取精度。4种方法在竹冠区域提取精度上,基于像元的提取精度小于基于多尺度的提取精度。可能是相比于乔木林,毛竹林由于其扩鞭生长的特殊性,竹冠部分难以单个呈现,具有很强的整体性,使用传统基于像元的提取方法难以达到对整体性的要求。基于多尺度分割的监督分类法表现最好的原因是在保证对象内部同质性最大的基础上进行的分类,且结合了光谱、纹理等信息进行分类,分类依据更全面;而相较于基于多尺度分割的监督分类法,基于多尺度分割的阈值分类法由于分类规则只基于光谱信息,较为单一,故精度较低。

      在基于像元的阈值分类法、基于像元的监督分类法和基于多尺度分割的阈值分类法中,是否钩梢对竹冠区域提取的影响存在较大影响。其一是由于未钩梢毛竹林郁闭度过高,林隙较小且颜色与竹冠接近,导致光学影像的饱和度增高,对特征量的敏感度随之下降,造成误差[26-27];其二是由于未钩梢毛竹林一般处于无人经营状态,存在毛竹倒伏现象和新老竹之间的光谱差异,造成误差。但是,基于多尺度分割的监督分类法则不受毛竹林是否钩梢的影响,其原因是多尺度分割得到内部同质性最大、外部同质性最小的斑块,很好保留了分类必需的有效信息,且避免了高分辨率影像中过于丰富的光谱与纹理信息所带来的干扰,减少错分与漏分的情况,基于此进行的监督分类能在保证竹冠整体性的前提下,通过训练样本分别提取钩梢与未钩梢毛竹林的特征参数,从而提取毛竹林竹冠区域。

      在4种方法中郁闭度对竹冠区域的提取均不存在显著影响,说明4种方法对不同郁闭度的毛竹林竹冠提取都有较好的适用性。原因是不同郁闭度的毛竹林遥感影像的光谱与纹理特征差别不大,这4种方法都足以从影像中提取稳定的特征量对竹冠区域进行提取。

      就郁闭度估测结果而言,4种方法提取的郁闭度误差都小于0.04,其中基于多尺度分割的监督分类法的误差最低,郁闭度误差仅为0.004,说明基于无人机的可见光影像可以用于毛竹林的郁闭度提取,具有很高的应用价值。苏迪等[9]基于冠层高度模型数据进行主成分确定自变量建立郁闭度模型,通过检验得到模型精度为83.18%,低于本研究的估测精度,说明基于影像提取的郁闭度优于建模所得的郁闭度。严羽[28]使用标记控制分水岭算法对树冠区域与非树冠区域进行分割,基于分割得到的树冠区域面积与实测面积对比,得到郁闭度估测误差为2.33%,高于本研究基于像元的阈值分类法、基于像元的监督分类法和基于多尺度分割的阈值分类法,但低于基于多尺度分割的监督分类法。

    • 在提取毛竹林竹冠区域和郁闭度的中,本研究基于像元的阈值分类、基于像元的监督分类、基于多尺度分割的阈值分类、基于多尺度分割的监督分类等4种方法都达到了较高精度,其中基于多尺度分割的监督分类法的提取结果精度最高,竹冠区域总体精度为98.86%,郁闭度误差为0.004,显著优于其他3种方法,而且是否钩梢对于基于多尺度分割的监督分类法提取竹冠区域没有影响。在4种方法中,郁闭度对竹冠区域的提取均不存在显著影响。总体上,低廉轻便的无人机搭载可见光相机节省了大量的调查时间与精力,提高了调查效率与精度,可以用于大面积毛竹林竹冠区域和郁闭度的提取。

参考文献 (28)

目录

    /

    返回文章
    返回