高级检索

留言板

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

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

基于改进的像元二分模型估测郁闭度

丹宇卓 石晶明 李心怡 罗帅 李明泽

丹宇卓, 石晶明, 李心怡, 罗帅, 李明泽. 基于改进的像元二分模型估测郁闭度[J]. 北京林业大学学报, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
引用本文: 丹宇卓, 石晶明, 李心怡, 罗帅, 李明泽. 基于改进的像元二分模型估测郁闭度[J]. 北京林业大学学报, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
Dan Yuzhuo, Shi Jingming, Li Xinyi, Luo Shuai, Li Mingze. Estimation of canopy closure based on improved dimidiate pixel model[J]. Journal of Beijing Forestry University, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
Citation: Dan Yuzhuo, Shi Jingming, Li Xinyi, Luo Shuai, Li Mingze. Estimation of canopy closure based on improved dimidiate pixel model[J]. Journal of Beijing Forestry University, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263

基于改进的像元二分模型估测郁闭度

doi: 10.13332/j.1000-1522.20180263
基金项目: 国家重点研发计划项目(2017YFD0600902)
详细信息
    作者简介:

    丹宇卓。主要研究方向:遥感与地理信息系统。Email:18846823073@163.com 地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学林学院

    通讯作者:

    李明泽,博士,教授。主要研究方向:遥感与地理信息系统。Email:mingzelee@163.com 地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学教务处

Estimation of canopy closure based on improved dimidiate pixel model

  • 摘要: 目的郁闭度是森林资源调查的一个重要因子,它不仅可以反映森林冠层的郁闭程度和树木利用空间的程度,并且能够指示林分密度。遥感为区域和全球尺度精确估测郁闭度提供了前所未有的契机,使得大面积郁闭度监测制图成为可能。本文旨在利用资源三号卫星影像数据结合模型对区域尺度上落叶松林进行郁闭度估测。方法以黑龙江省佳木斯市桦南县孟家岗林场落叶松人工林为研究对象,首先采用像元二分模型对植被覆盖度进行估算,通过探寻植被覆盖度与郁闭度的关系,对像元二分模型进行改进,通过计算不同累积频率归一化植被指数(NDVI)的值作为模型参数,对比不同参数取值时模型拟合效果得到最优模型,并利用模型估测郁闭度。结果拟合结果表明,2%累积频率下模型拟合效果最好,模型R2为0.838,RMSE为0.045,最后利用该模型得出孟家岗林场落叶松人工林郁闭度分布图。结论利用改进的像元二分模型可以较为准确地估测郁闭度。本研究对探索我国北方落叶松人工林郁闭度遥感估测提供了更有效的途径,同时也为森林资源和参数的调查提供一定的参考依据。
  • 图  1  黑龙江省佳木斯市孟家岗林场落叶松样地分布

    Figure  1.  Distribution of Larix gmelinii sample plots in Mengjiagang, Heilongjiang Province

    图  2  资源三号卫星影像的预处理

    Figure  2.  ZY-3 data preprocessing

    图  3  孟家岗林场落叶松人工林植被覆盖度反演

    Figure  3.  Vegetation coverage inversion map of Larix gmelini in Mengjiagang,Heilongjiang Province

    图  4  累计频率的估测值与真实值的R2和RMSE

    Figure  4.  R2 and RMSE of the estimated value and the true value at the cumulative frequency

    图  5  孟家岗林场落叶松人工林郁闭度反演

    Figure  5.  Canopy closure inversion map of Larix gmelini in Mengjiagang,Heilongjiang Province

    表  1  样地郁闭度数据

    Table  1.   Canopy closure data of sample plot

    样地号 Sample plot No.郁闭度 Canopy closure样地号 Sample plot No.郁闭度 Canopy closure样地号 Sample plot No.郁闭度 Canopy closure
    10.81130.79250.73
    20.74140.86260.82
    30.76150.81270.88
    40.79160.79280.81
    50.84170.83290.84
    60.81180.84300.86
    70.79190.87310.84
    80.77200.87320.81
    90.73210.86330.81
    100.85220.85340.82
    110.80230.74
    120.83240.80
    下载: 导出CSV

    表  2  不同累计频率的拟合效果

    Table  2.   Fitting effects of different accumulative frequency

    森林类型 Forest type累计频率 Cumulative frequencyNDVIvNDVIsR2RMSE
    落叶松人工林 Larix gmelinii plantation 0.0050.840.04 0.7210.101
    0.010.830.0430.7820.051
    0.020.820.0450.8380.045
    0.050.810.0540.8030.055
    下载: 导出CSV
  • [1] Fiala A C S, Garman S L, Gray A N. Comparison of five canopy cover estimation techniques in the western Oregon Cascades[J]. Forest Ecology and Management, 2006, 232(1−3): 188−197.
    [2] Jennings S B. Assessing forest canopies and understorey illumination: canopy closure, canopy cover and other measures[J]. Forestry: an International Journal of Forest Research, 1999, 72(1): 59−74. doi:  10.1093/forestry/72.1.59
    [3] Ruiliang P U, Gong P. Wavelet transform applied to EO-1 hyperspectral data for forest LAI and crown closure mapping[J]. Remote Sensing of Environment, 2004, 91(2): 212−224. doi:  10.1016/j.rse.2004.03.006
    [4] Korhonen L, Korhonen K T, Rautiain M, et al. Estimation of forest canopy cover: a comparison of field measurement techniques[J]. Silva Fennica, 2006, 40(4): 577−588.
    [5] Stumpf K A. The estimation of forest vegetation cover descriptions using a vertical densitometer [J/OL]. Society of American Foresters Convention, 1994. [2018−06−02]. https://www.researchgate.net/publication/237510073.
    [6] North P R J. Estimation of fAPAR, LAI, and vegetation fractional cover from ATSR-2 imagery[J]. Remote Sensing of Environment, 2002, 80(1): 114−121. doi:  10.1016/S0034-4257(01)00292-9
    [7] 刘大伟, 孙国清, 庞勇, 等. 利用Landsat TM数据对森林郁闭度进行遥感分级估测[J]. 遥感信息, 2006, 13(1):41−42, 24. doi:  10.3969/j.issn.1000-3177.2006.01.011

    Liu D W, Sun G Q, Pang Y, et al. Estimation of forest closure using Landsat TM data[J]. Remote Sensing Information, 2006, 13(1): 41−42, 24. doi:  10.3969/j.issn.1000-3177.2006.01.011
    [8] 杜晓明, 蔡体久, 琚存勇. 采用偏最小二乘回归方法估测森林郁闭度[J]. 应用生态学报, 2008, 19(2):273−277.

    Du X M, Cai T J, Ju C Y. Estimation of forest canopy closure by using partial least square regression[J]. Chinese Journal of Applied Ecology, 2008, 19(2): 273−277.
    [9] 郑冬梅, 曾伟生, 智长贵, 等. 三峡库区森林郁闭度的遥感定量估测[J]. 中南林业科技大学学报, 2013, 33(9):1−4, 26. doi:  10.3969/j.issn.1673-923X.2013.09.001

    Zheng D M, Zeng W S, Zhi C G, et al. Remote sensing estimation of forest canopy closure in forests of Three Gorges Reservoir Region[J]. Journal of Central South University of Forestry & Technology, 2013, 33(9): 1−4, 26. doi:  10.3969/j.issn.1673-923X.2013.09.001
    [10] 谭炳香, 李增元, 陈尔学, 等. Hyperion高光谱数据森林郁闭度定量估测研究[J]. 北京林业大学学报, 2006, 28(3):95−101. doi:  10.3321/j.issn:1000-1522.2006.03.017

    Tan B X, Li Z Y, Chen E X, et al. Estimating forest crown closure using Hyperion hyperspectral data[J]. Journal of Beijing Forestry University, 2006, 28(3): 95−101. doi:  10.3321/j.issn:1000-1522.2006.03.017
    [11] Armando A, Compagna L . SAT-based model-checking for security protocols analysis[J]. International Journal of Information Security, 2008, 7(1): 3−32. doi:  10.1007/s10207-007-0041-y
    [12] Gupta A. Formal hardware verification methods: a survey[J]. Formal Methods in System Design, 1992,1 (2−3): 151−238.
    [13] 吴瑶, 秦凯伦. 林龄对长白落叶松胸径生长的影响[J]. 林业科技, 2016, 41(6):8−11. doi:  10.3969/j.issn.1001-9499.2016.06.003

    Wu Y, Qin K L. Effects of forest age on DBH growth of Larix olgensis[J]. Forestry Science & Technology, 2016, 41(6): 8−11. doi:  10.3969/j.issn.1001-9499.2016.06.003
    [14] 郭孝玉. 长白落叶松人工林树冠结构及生长模型研究[D]. 北京: 北京林业大学, 2013.

    Guo X Y. Crown structure and growth model for Larix olgensis plantation[D]. Beijing: Beijing Forestry University, 2013.
    [15] 祁有祥, 骆汉, 赵廷宁. 基于鱼眼镜头的林冠郁闭度简易测量方法[J]. 北京林业大学学报, 2009, 31(6):60−66.

    Qi Y X, Luo H, Zhao T N. Simplified approach to measure canopy closure based on fish lenses[J]. Journal of Beijing Forestry University, 2009, 31(6): 60−66.
    [16] 刘玉安, 黄波, 程涛, 等. 基于像元二分模型的淮河上游植被覆盖度遥感研究[J]. 水土保持通报, 2012, 32(1):93−97, 267.

    Liu Y A, Huang B, Cheng T, et al. Vegetation coverage in upper Huaihe River Basin based on binary pixel model of remote sensing[J]. Bulletion of Soil and Water Conservation, 2012, 32(1): 93−97, 267.
    [17] 陈妍, 宋豫秦, 王伟. 基于随机森林回归的草场植被盖度反演模型研究: 以新疆阿勒泰地区布尔津县为例[J]. 生态学报, 2018, 38(7):2384−2394.

    Chen Y, Song Y Q, Wang W. Grassland vegetation cover inversion model based on random forest regression: a case study in Burqin County, Altay, Xinjiang Uygur Autonomous Region[J]. Acta Ecologica Sinica, 2018, 38(7): 2384−2394.
    [18] Li X B, Chen Y H, Chen Z J, et al. Estimating vegetation fractional coverage of steppe based on multi-spatial resolution remotely sensed images[C]//IEEE International symposium on Geoscience and Remote Sensing. Denver: IEEE, 2006.
    [19] 汤旭光. 基于激光雷达与多光谱遥感数据的森林地上生物量反演研究[D]. 沈阳: 中国科学院东北地理与农业生态研究所, 2013.

    Tang X G. Estimation of forest aboveground biomass by integrating ICESat/GLAS waveform and TM date[D]. Shenyang: Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, 2013.
    [20] 陈文灿. 郁闭度模型的建立及应用[J]. 华东森林经理, 2015, 29(2):63−65. doi:  10.3969/j.issn.1004-7743.2015.02.018

    Chen W C. Establishment and application of canopy closure model[J]. East China Forest Management, 2015, 29(2): 63−65. doi:  10.3969/j.issn.1004-7743.2015.02.018
    [21] 国家林业局. 森林资源规划设计调查技术规程: GB/T 26424—2010[S]. 北京: 中国标准出版社, 2010.

    State Forestry Administration. Technical regulations for inventory for forest management planning and design: GB/T 26424—2010[S]. Beijing: China Standards Press, 2010.
    [22] 邹杰, 诸葛祥东. 森林郁闭度及其测定方法[J]. 黑龙江科技信息, 2011(35):290. doi:  10.3969/j.issn.1673-1328.2011.35.296

    Zou J, Zhuge X D. Forest canopy closure and its measuring methods[J]. Heilongjiang Science and Technology Information, 2011(35): 290. doi:  10.3969/j.issn.1673-1328.2011.35.296
  • [1] 李鹏飞, 郭小平, 顾清敏, 张昕, 冯昶栋, 郭光.  基于可见光植被指数的乌海市矿山排土场坡面植被覆盖信息提取研究 . 北京林业大学学报, 2020, 42(6): 102-112. doi: 10.12171/j.1000-1522.20190252
    [2] 吕娇, MustaqShah, 崔义, 徐程扬.  土壤紧实度和凋落物覆盖对城市森林土壤持水、渗水能力的影响 . 北京林业大学学报, 2020, 42(8): 102-111. doi: 10.12171/j.1000-1522.20190476
    [3] 孙钊, 潘磊, 孙玉军.  基于无人机影像的高郁闭度杉木纯林树冠参数提取 . 北京林业大学学报, 2020, 42(10): 20-26. doi: 10.12171/j.1000-1522.20190386
    [4] 邵艳莹, 吴秀芹, 张宇清, 秦树高, 吴斌.  内蒙古地区植被覆盖变化及其对水热条件的响应 . 北京林业大学学报, 2018, 40(4): 33-42. doi: 10.13332/j.1000-1522.20170414
    [5] 惠刚盈, 张连金, 胡艳波, 王宏翔, 张弓乔.  林分拥挤度及其应用 . 北京林业大学学报, 2016, 38(10): 1-6. doi: 10.13332/j.1000-1522.20160107
    [6] 李瀚之, 余新晓, 樊登星, 孙佳美, 常玉, 梁洪儒.  模拟降雨条件下黑麦草覆盖度对褐土坡面径流流速的影响 . 北京林业大学学报, 2014, 36(5): 94-98. doi: 10.13332/j.cnki.jbfu.2014.05.014
    [7] 尤号田, 邢艳秋, 冉慧, 王蕊, 霍达.  基于LiDAR 点云能量信息的樟子松郁闭度反演方法 . 北京林业大学学报, 2014, 36(6): 30-35. doi: 10.13332/j.cnki.jbfu.2014.06.009
    [8] 曾涛, 琚存勇, 蔡体久, 刘文彬, 姚月锋.  利用变量投影重要性准则筛选郁闭度估测参数 . 北京林业大学学报, 2010, 32(6): 37-41.
    [9] 李晓松, 李增元, 高志海, 白黎娜, 王琫瑜, 李世明.  基于Hyperion植被指数的干旱地区稀疏植被覆盖度估测 . 北京林业大学学报, 2010, 32(3): 95-100.
    [10] 武鹏飞, 王茂军, 张学霞.  北京市植被绿度与城市热岛效应关系研究 . 北京林业大学学报, 2009, 31(5): 54-60.
    [11] 信忠保, 余新晓, 甘敬, 王小平, 李金海.  黄河中游河龙区间植被覆盖变化与径流输沙关系研究 . 北京林业大学学报, 2009, 31(5): 1-8.
    [12] 祁有祥, 骆汉, 赵廷宁.  基于鱼眼镜头的林冠郁闭度简易测量方法 . 北京林业大学学报, 2009, 31(6): 60-66.
    [13] 周金星, JunYang, 董林水, 张旭东, 左力.  青藏铁路唐古拉山南段沿线植被多样性及盖度特征分析 . 北京林业大学学报, 2008, 30(3): 24-30.
    [14] 张学霞, 朱清科, 吴根梅, 孙丹珠.  数码照相法估算植被盖度 . 北京林业大学学报, 2008, 30(1): 164-169.
    [15] 蔡体久, 琚存勇.  毛乌素沙地植被盖度估测方法研究 . 北京林业大学学报, 2008, 30(增刊1): 292-295.
    [16] 胡理乐, 闫伯前, 孙一荣, 朱教君, .  从径向生长量分析红松幼树生长光需求 . 北京林业大学学报, 2008, 30(2): 147-150.
    [17] 余新晓, 王雄宾, 张晓明, .  北京山区防护林适宜林分结构研究 . 北京林业大学学报, 2008, 30(supp.2): 103-106.
    [18] 徐文铎, 李吉跃, 段爱国, 张求慧, 吴斌, 温俊宝, 刘金福, 常德龙, 王云琦, 匡文慧, 谭炳香, 郑凌凌, 赵燕东, 程占红, 齐春辉, 
    ZHAOGuang-jie, 冯夏莲, 张灿, 邹大林, 李贤军, 李雪萍, 刘常富, 王玉涛, 鲁绍伟, 李雪峰, 张建国, 翟洪波, LUOWen-sheng, 张路平, 何正权, 宋湛谦, 温俊宝, 韩士杰, 白陈祥, 张树文, 何兴元, 李增元, 朱天辉, 王玉杰, 吴庆利, 何承忠, 韩烈保, 何友均, 李吉跃, 余新晓, 吴斌, 洪伟, 赵广杰, 黄文豪, 林秦文, 姜伟, 骆有庆, 陈发菊, 匡秋明, 童书振, 刘凤芹, 何静, 陈尔学, 梁小红, 郭忠玲, FurunoTakeshi, 张志毅, 张养贞, 骆有庆, ]魏晓霞, ]陈玮, 李俊清, 郑兴波, 庞勇, RENQian, 许志春, 胡伟华, 张璧光, 张军, 赵桂玲, 张振明, 曾会明, 许志春, 李颖, 梁宏伟, 安新民, 崔国发, 杨凯, 宋国正, 曹川健, 郑杰, 李凤兰, 侯伟, 刘君, PaulWolfgang, 李福海, 赵广亮, 雷渊才, 李考学, 姚永刚, 田桂芳, 董建生, 张有慧, 张全来, 李永波, 赫万成, 李长明, 张世玺.  Hyperion高光谱数据森林郁闭度定量估测研究 . 北京林业大学学报, 2006, 28(3): 95-101.
    [19] 李全发, 詹亚光, 李梅, 马文辉, 杨海龙, 龙玲, 刘震, 宋小双, 侯亚南, 王明枝, 熊瑾, 张一平, 符韵林, 殷亚方, 黄国胜, 张秋英, 李慧, 
    王保平, 李景文, 李景文, 饶良懿, 杨晓晖, 杜华强, 窦军霞, 尹立辉, 徐峰, 朱金兆, 韩海荣, 吕建雄, 李妮亚, 秦瑶, 陆熙娴, 李吉跃, 李发东, 范文义, 刘文耀, 耿晓东, 赵敏, 朱金兆, 王雪军, 梁机, 李俊清, 张克斌, 李俊清, 陈晓阳, 王洁瑛, 齐实, 赵宪文, 毕华兴, 秦素玲, 于贵瑞, 陈晓阳, 康峰峰, 李云, 沈有信, 孙玉军, 陈素文, 刘桂丰, 李黎, 倪春, 唐黎明, 刘雪梅, 乔杰, 欧国强, 慈龙骏, 李凤兰, 朱国平, 蒋建平, 宋献方, 赵双菊, 马钦彦, 李伟, 文瑞钧, 张桂芹, 李伟, 韦广绥, 王玉成, 刘伦辉, 王雪, 魏建祥, 任海青, 黎昌琼, 丁霞, 孙涛, 周海江, 宋清海, 李慧, 杨谦, 张万军, , 刘莹, 孙志强, 孙晓敏, 李宗然, 
    木 材 表 面 粗 糙 度 的 分 析 . 北京林业大学学报, 2005, 27(1): 14-18.
    [20] 武三安, 胡建忠, 杜晓, 莫秋云, 陈玮, 贾黎明, 刁一伟, 李成茂, 刘晓丽, 张丽丽, 韩烈保, 张建军, 张厚江, 郭惠红, 石娟, 赵博光, 王昌俊, 申世杰, 梁波, 马履一, 李文彬, 姜笑梅, 徐文铎, 邢长山, 骆有庆, 李镇宇, 王安志, 清水晃, 宋菲, 张峻萍, 金昌杰, 苏德荣, 崔英颖, 壁谷直记, 王小平, 曾凡勇, 李海林, 李景锐, 赵林果, 石碧, 殷亚方, 沉昕, 延廣竜彦, 徐梅, 蒋艳灵, 韩瑞东, 苗毅, 陈卫平2, 韦艳葵, 胡青, 关德新, 王瀛坤, 裴铁璠, 赵永利, 严晓素, 高述民, 徐君, 李凤兰, 周军, 蒋平, 蒋平.  确定监测区域建立森林郁闭度估测方程最优样地的研究 . 北京林业大学学报, 2005, 27(6): 24-28.
  • 加载中
图(5) / 表 (2)
计量
  • 文章访问数:  972
  • HTML全文浏览量:  1007
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-12
  • 修回日期:  2019-01-11
  • 网络出版日期:  2019-06-11
  • 刊出日期:  2019-06-01

基于改进的像元二分模型估测郁闭度

doi: 10.13332/j.1000-1522.20180263
    基金项目:  国家重点研发计划项目(2017YFD0600902)
    作者简介:

    丹宇卓。主要研究方向:遥感与地理信息系统。Email:18846823073@163.com 地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学林学院

    通讯作者: 李明泽,博士,教授。主要研究方向:遥感与地理信息系统。Email:mingzelee@163.com 地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学教务处

摘要: 目的郁闭度是森林资源调查的一个重要因子,它不仅可以反映森林冠层的郁闭程度和树木利用空间的程度,并且能够指示林分密度。遥感为区域和全球尺度精确估测郁闭度提供了前所未有的契机,使得大面积郁闭度监测制图成为可能。本文旨在利用资源三号卫星影像数据结合模型对区域尺度上落叶松林进行郁闭度估测。方法以黑龙江省佳木斯市桦南县孟家岗林场落叶松人工林为研究对象,首先采用像元二分模型对植被覆盖度进行估算,通过探寻植被覆盖度与郁闭度的关系,对像元二分模型进行改进,通过计算不同累积频率归一化植被指数(NDVI)的值作为模型参数,对比不同参数取值时模型拟合效果得到最优模型,并利用模型估测郁闭度。结果拟合结果表明,2%累积频率下模型拟合效果最好,模型R2为0.838,RMSE为0.045,最后利用该模型得出孟家岗林场落叶松人工林郁闭度分布图。结论利用改进的像元二分模型可以较为准确地估测郁闭度。本研究对探索我国北方落叶松人工林郁闭度遥感估测提供了更有效的途径,同时也为森林资源和参数的调查提供一定的参考依据。

English Abstract

丹宇卓, 石晶明, 李心怡, 罗帅, 李明泽. 基于改进的像元二分模型估测郁闭度[J]. 北京林业大学学报, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
引用本文: 丹宇卓, 石晶明, 李心怡, 罗帅, 李明泽. 基于改进的像元二分模型估测郁闭度[J]. 北京林业大学学报, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
Dan Yuzhuo, Shi Jingming, Li Xinyi, Luo Shuai, Li Mingze. Estimation of canopy closure based on improved dimidiate pixel model[J]. Journal of Beijing Forestry University, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
Citation: Dan Yuzhuo, Shi Jingming, Li Xinyi, Luo Shuai, Li Mingze. Estimation of canopy closure based on improved dimidiate pixel model[J]. Journal of Beijing Forestry University, 2019, 41(6): 35-43. doi: 10.13332/j.1000-1522.20180263
  • 森林郁闭度是指林冠垂直投影面积与林地面积之比,它的测定方法有很多,由于受到技术方法和条件的限定,过去的测定多采用目测法、样线法、树冠投影法和样点法进行估测[1]。这些方法精度较低,且很难对大区域尺度的郁闭度进行连续估测。随着遥感技术的发展,国际上在该领域开展了广为深入的研究[2]

    由于遥感数据波段设置的特殊性,利用不同波段遥感数据计算出来的植被指数可以有效地表征地表的森林覆盖情况,进而估算郁闭度。因而遥感数据在提取郁闭度方面的应用引起了许多学者的重视,就不同遥感变量与郁闭度的关系进行了大量的研究,结果表明郁闭度与遥感变量之间有着很强的相关性,相关的形式可能是线性或者非线性的[34]。此外,由于在特定情况下植被覆盖度与郁闭度有类似的概念,因而植被覆盖度的遥感估算方法对于郁闭度的遥感估算有一定的借鉴意义[56]

    我国学者对于郁闭度的遥感估测研究主要集中在线性和非线性回归模型两个方面。刘大伟等[7]利用Landsat TM数据的波段亮度值,归一化植被指数(NDVI),缨帽变换后的亮度、绿度和湿度等因子和郁闭度进行相关性分析,选取了相关性最高的湿度因子与训练数据建立线性回归方程反演森林郁闭度并生成等级图。杜晓明等[8]结合Landsat ETM+数据及国家森林一类清查资料,用Bootstrap方法筛选郁闭度估测数据及森林资源一类清查数据为基础[9],定量分析了遥感因子、地形因子与郁闭度的关系,利用逐步回归、最小二乘回归、主分量回归和偏最小二乘回归构建模型反演郁闭度。谭炳香等[10]利用星载EO−1Hyperion高光谱遥感数据,采用波段选择法和主成分分析法定量估测了森林郁闭度。

    然而,地表实测郁闭度的方法仅是能获得一些指定区域上的数据,不利于大范围郁闭度空间分布及变化,而利用模型监测森林郁闭度较为复杂,适合于小区域的研究,所需要的参数较多且大多数参数难以获取,从而会影响森林郁闭度的提取精度[1112]。如何用简便方法获得高精度的森林郁闭度成为亟待解决的问题。改进的像元二分模型估测森林郁闭度主要通过计算各种基本地物类型的植被指数对混合像元植被指数的贡献来估测郁闭度。本研究通过计算不同累积频率归一化植被指数(NDVI)的值作为模型参数,对比不同参数取值时模型拟合效果得到最优模型。同时,本研究利用资源三号卫星影像,空间分辨率较大,实验精度较高,可扩展应用于大尺度上的森林郁闭度反演,这为探索我国北方落叶松(Larix gmelinii)人工林郁闭度遥感估测提供了一条更有效的途径,进而为森林资源和参数的调查提供一定的参考依据,且此方法简单,精度符合要求,适于在林场等经营单位推广。

    • 研究区主要包括佳木斯市桦南县孟家岗林场,地理坐标为130°32′42″ ~ 130°52′36″E、46°20′16″ ~ 46°30′50″N。林场西与桦南县林业局的七峰、青背林场毗邻,东与桦南林业局相邻。孟家岗林场地处完达山西麓,海拔168 ~ 575 m之间,全场总面积16 733 hm2。人工造林占林地面积的76.7%,落叶松、樟子松(Pinus sylvestris)、红松(Pinus koraiensis)约占80%。孟家岗林场是一个以针叶树为主的人工林用材基地。人工造林成活率达95%以上,最高达到99%,育苗技术远远超过原苏联和日本同纬度地区的育苗水平。全场林木蓄积量64万m3,其中人工林木蓄积量占2/3,孟家岗林场落叶松样地分布图如图1所示。

      图  1  黑龙江省佳木斯市孟家岗林场落叶松样地分布

      Figure 1.  Distribution of Larix gmelinii sample plots in Mengjiagang, Heilongjiang Province

    • 本研究所用的遥感影像来源于2012年8月30日的资源三号卫星多光谱影像和2013年5月18日的资源三号卫星的多光谱影像。该数据无云覆盖,成像质量较好。冠幅是研究树冠水平伸展空间结构的指标。胸径和林分密度是影响冠幅的重要因子[13]。经调查该林分在5年间未经过抚育、采伐及其他人为干扰,林分密度未发生显著变化。不同林龄的胸径平均生长量呈现先大后小的状态,落叶松随着林龄的增加,胸径的增长幅度趋于平缓,成熟林胸径的连年生长量逐渐减小[14],由此可推断出,5年内落叶松人工林的胸径增长也并不显著。因此,林分的冠幅也未发生显著变化,故在相同的生长季节,遥感数据与样地实测数据差别5年该林分的郁闭度差异不会很大,可以将两类数据进行比较。

    • 本研究采用的样地数据采集于2017年6月,共选取30块面积为0.06 hm2的落叶松人工林样地。每块样地利用亚米级GPS定位(定位的点见图1),同时利用鱼眼相机拍摄照片,便于郁闭度真实值的求算。

    • 本研究参考了2016年一类调查所得矢量数据作为辅助数据,确定孟家岗林场的边界与落叶松人工林的范围。

    • 基于祁有祥等[15]的研究结果,本研究对鱼眼相机拍摄的照片进行处理,利用Photoshop对鱼眼相机拍摄的照片进行二值化处理,将照片亮度降低,对比度调高,勾选天空部分算出所占比例,求得东西南北中的郁闭度值,并计算出平均郁闭度值作为该样地郁闭度的真实值(表1),用于研究过程中模型的对比和检验。

      表 1  样地郁闭度数据

      Table 1.  Canopy closure data of sample plot

      样地号 Sample plot No.郁闭度 Canopy closure样地号 Sample plot No.郁闭度 Canopy closure样地号 Sample plot No.郁闭度 Canopy closure
      10.81130.79250.73
      20.74140.86260.82
      30.76150.81270.88
      40.79160.79280.81
      50.84170.83290.84
      60.81180.84300.86
      70.79190.87310.84
      80.77200.87320.81
      90.73210.86330.81
      100.85220.85340.82
      110.80230.74
      120.83240.80
    • 本文采用的是资源三号卫星拍摄的两景数据覆盖研究区,空间分辨率为5.8 m。为得到准确的光谱信息,本研究利用美国ESRI公司开发的ENVI5.3对遥感图像进行预处理。

      由于地形会产生图像几何畸变,先对图像进行无控制点的正射校正,ENVI中会默认使用自带的全球900 m空间分辨率DEM数据。生成多中心投影平面正射图像,纠正一般系统因素产生的几何畸变。再将资源三号卫星数据进行辐射定标。利用该地区的DEM影像计算出地面平均高度。此外,对影像中心进行大气校正,减少或消除大气对遥感图像的影响,得到真实的地表反射信息。最后,对两景资源三号卫星图像进行无缝镶嵌,利用孟家岗边界矢量数据对拼接后的图像进行剪裁(图2)。

      图  2  资源三号卫星影像的预处理

      Figure 2.  ZY-3 data preprocessing

    • 该研究将孟家岗林场作为研究区域。像元二分模型的原理如下:假定遥感图像中一个像元的光谱信息是由两个组分因子的光谱信息线性组合而来,而每个组分对传感器所获得的光谱值都有贡献,以其所占比例即为该组分因子所贡献信息的权重。在计算原始的植被覆盖度时,像元由植被覆盖地表与无植被覆盖地表两部分组成,所得的光谱值也只由这两个组分线性合成。对于混合像元,其遥感信息可表示为植被覆盖部分(Iv) 信息和无植被覆盖(Is)贡献的信息之和[16]

      $$I = {I_{\rm{v}}} + {I_{\rm{s}}}$$ (1)

      式中:I为混合像元的遥感信息量,Iv为全由植被覆盖的像元的遥感信息量,Is为无植被覆盖的像元的遥感信息量。

      该研究利用归一化植被指数(NDVI)作为信息量进行计算。设有植被覆盖的比例为fc,无植被覆盖的比例为1 − fc。NDVI成为了常用的指数进行植被覆盖度的估算,有利于对区域森林总体状况进行定量评估,而且方法简单易行。反映生物量、植被生长状况等信息都用该指数。陈研等[17]、Li 等[18]的研究表明,基于NDVI和不同模型对不同地区进行植被覆盖度反演估算,估算结果都很理想。陈云浩等[19]采用ASTER和MODIS两种尺度的影像,分别用NDVI、NDGI、RVI 3种指数估算了中国北方的草原植被覆盖度,通过对比分析,验证NDVI的估算精度最高。本研究利用NDVI作为参数,将NDVI代入式(1) 转化得:

      $${\rm{NDVI}} = {f_{\rm{c}}} \cdot {\rm{NDV{I_v}}} + \left( {1 - {f_{\rm{c}}}} \right) \cdot {\rm{NDV{I_s}}}$$ (2)
      $$\text{求解得}\;\;\;\;{f_{\rm{c}}}} = \frac{{{\rm{NDVI}} - {\rm{NDV{I_s}}}}}{{{\rm{NDV{I_v}}} - {\rm{NDV{I_s}}}}}$$ (3)

      式中:fc 为植被覆盖度,NDVIv为植被覆盖部分的NDVI值,NDVIs为非植被覆盖的地区的NDVI值。

    • 郁闭度与植被覆盖度的测定原理比较相似,故本研究通过修改像元二分模型,并利用改进后的模型估算林分郁闭度[20]

      郁闭度指森林中乔木树冠遮蔽地面的程度,是反映林分密度的指标。植被覆盖度指植被地上部分垂直投影面积占地面面积的比率,又称为投影盖度。对比发现,郁闭度是将植被覆盖度的分割方式做了变化,将有植被覆盖部分的非林冠分割至无植被覆盖部分,故若想利用植被盖度模型估算郁闭度,应对有植被覆盖的非林冠部分,即灌草部分的信息加以研究。

      在郁闭度模型中将林冠部分的NDVI值记为NDVI′ v,将非林冠部分的NDVI值记为NDVI′ s,由于植被类型和立地条件的不同,以及植被覆盖随季节变化和叶冠背景的干扰等因素,NDVI′ v、NDVI′ s值会随着时间和空间而改变。因此,厘清林冠部分和非林冠部分的差异,并选择合理的NDVI取值将是本研究的重点之一。

      对于人工纯林而言,不同小区域内林冠部分差异不大,而非林冠部分由于郁闭度引起的光照条件差异使得非林冠部分差异较大。根据《森林资源规划设计调查技术规程》(GB/T 26424—2010)[21],我国森林郁闭度0.20 ~ 0.39为低郁闭度,一般以0.40 ~ 0.69为中郁闭度,0.70以上为高郁闭度。本研究所涉及区域为人工林林场,且由样地实测数据发现郁闭度均大于0.7,属于密郁闭。本研究探索密郁闭状态下落叶松人工林郁闭度估测时变量的选择方法。

    • 利用改进的像元二分模型进行森林郁闭度估算的难点也在于如何确定NDVI的两个阈值的问题,利用获取到的孟家岗林场矢量图,可将落叶松人工林部分影像提取出来。通过对影像的解译发现,影像信息全部由植被覆盖与裸露土壤部分构成。根据NDVI物理意义,应选取最接近于1的值作为植被覆盖部分的信息,而最接近0的值作为裸露土壤部分的信息,在估算盖度时由此应用。

      本研究估算的是森林郁闭度,并考虑到环境噪声及部分异常区域的影响,不能直接选取0 ~ 1范围内的最大与最小值作为参数。由于郁闭度与盖度的差别为NDVIv和NDVIs贡献主体不一样,分别为有无植被覆盖部分和林冠及非林冠部分。植被覆盖部分包括乔木、灌木和草本,而林冠部分只为乔木部分,即NDVIs中加上了灌草的NDVI值,即NDVIs阈值会相应变大为NDVI′ s,而林冠部分NDVIv 的值只剩下林冠部分即乔木贡献的NDVI,即NDVIv阈值会相应变小为NDVI′ v。因此,本研究采用计算NDVI数据的频率累积值作为参数,将研究区域孟家岗林场落叶松人工林的像元NDVI值进行统计,并从小到大进行排序,根据NDVI数据频率统计表,计算NDVI数据的频率累积值。考虑到区域的差异,所对应的NDVI阈值不同,在尝试了多个累积频率不同的NDVI值代入参数后,最终设置距离累积频率最大/最小值为2%、1%、0.5%、0.2%共4个位置处的NDVI设为NDVI′ v和NDVI′ s的值作为参数构建模型,利用残差和拟合精度进行比较选取最优。

    • 本研究首先利用传统像元二分模型的方法求出研究区域的植被盖度。计算整幅图像的NDVI值,考虑到区域差异,结合研究区域的具体情况选取参数。利用感兴趣区域的研究法求出有植被覆盖地区NDVI的最大值和无植被覆盖地区即裸地区域NDVI的最小值,分别作为植被覆盖地区的NDVI和裸地的NDVI,代入公式中求出林场区域植被覆盖度的真实值。在此遥感影像中,NDVI的最小值与最大值分别为0.012和0.996,代入公式${f_{\rm{c}}} = \dfrac{{{\rm{NDVI}} - 0.012}}{{0.996 - 0.012}}$中,即

      $$\begin{aligned}{f_{\rm{c}}} =& \frac{{{\rm{NDVI}} - 0.012}}{{0.996 - 0.012}} = \frac{{{\rm{NDVI}} - 0.012}}{{0.984}} = \\ &1.016 {\rm{NDVI}} - 0.012\;2\end{aligned}$$ (4)

      利用算出的植被覆盖度绘制植被覆盖度图像(图3)。

      图  3  孟家岗林场落叶松人工林植被覆盖度反演

      Figure 3.  Vegetation coverage inversion map of Larix gmelini in Mengjiagang,Heilongjiang Province

    • 研究选取累积频率距离最大值和最小值2%和98%、5%和95%、0.5%和99.5%、1%和99%处的值,设为林冠部分NDVI值(NDVIv)和非林冠覆盖地区的NDVI值(NDVIs)后,代入公式(4),求出样地点处郁闭度的估测值。利用估测所得郁闭度值与实测的郁闭度值进行对比,得到R2、RMSE(图4)和拟合效果(表2)。

      表 2  不同累计频率的拟合效果

      Table 2.  Fitting effects of different accumulative frequency

      森林类型 Forest type累计频率 Cumulative frequencyNDVIvNDVIsR2RMSE
      落叶松人工林 Larix gmelinii plantation 0.0050.840.04 0.7210.101
      0.010.830.0430.7820.051
      0.020.820.0450.8380.045
      0.050.810.0540.8030.055

      图  4  累计频率的估测值与真实值的R2和RMSE

      Figure 4.  R2 and RMSE of the estimated value and the true value at the cumulative frequency

      由于2%所对应的R2和RMSE分别为0.838和0.045,拟合度较好而误差较小,因此我们利用2%提供的参数进行回归建模,代入参数NDVIv=0.819 459,NDVIs=0.044 601,代入公式中得到最后的模型公式:

      $$ \begin{aligned} &\;\\ &{f_{\rm{c}}} = \frac{{{\rm{NDVI}} - 0.044\;601}}{{0.819\;459 - 0.044\;601}} \end{aligned} $$

      $$\begin{aligned} &{f_{\rm{c}}} = \dfrac{{{\rm{NDVI}} - 0.044\;601}}{{0.774\;858}}\\ &{f_{\rm{c}}} = 1.29{\rm{NDVI}} - 0.057\;5 \end{aligned}$$ (5)

      根据此模型反演整个研究区域郁闭度并绘制郁闭度分布图,见图5

      图  5  孟家岗林场落叶松人工林郁闭度反演

      Figure 5.  Canopy closure inversion map of Larix gmelini in Mengjiagang,Heilongjiang Province

      对研究区内落叶松人工林郁闭度估测结果进行分析可知:该区落叶松人工林属于密林林,郁闭度大多都在0.7以上。研究区域东边郁闭度高,西边郁闭度低,北边郁闭度较高,南边郁闭度较低。

    • 相比于原始像元二分模型针对植被覆盖度进行研究,本研究通过探寻植被覆盖度与郁闭度的关系对像元二分模型进行改进。在模型参数选择时,采用计算不同累积频率NDVI的值作为土壤与灌草的混合NDVI值代入像元二分模型进行计算是本文的亮点与优势所在。此外,本研究通过对比分析,找到精度最高的累积频率参数NDVI的值,代入公式求得最佳模型。

      郁闭度本身就是一个较为粗略的数值[22],本研究提供的方法可以通过较少的实测样地点,在大区域上方便地计算出较为精确的郁闭度的值,对森林的管理与经营及森林资源的调查提供了一种更为简便的方式。基于像元二分模型是一种植被覆盖度的估算方法,对于不同分辨率、不同地区、不同数据源下都有很多研究对此进行了分析和考证。本文首次利用该方法对空间分辨率较高的资源三号卫星数据进行郁闭度的估算,提高了模型的精度,使模型具有推广性,可减少前人数据横跨多个像元导致的误差。在NDVI阈值选取时,通过大气校正,辐射校正等方法,消除大气、太阳和地形等的影响,得到较为准确的NDVI值,但植被在不同的生长阶段、不同的月份光合能力不同,含水量不同,影响NDVI的值也不同,因此是否能够推广模型运用到不同的林分以及不同的生长阶段中进行郁闭度估算,这也是我们今后研究的方向。是否可以利用高光谱影像端提取算法,计算每个像元各类地物的贡献度建立像元三分模型等也是一个值得继续深入研究的方向之一。除此之外,本研究区域的落叶松人工林长势较为相同,样地的郁闭度差别并不大,在以后的研究中需要选取不同情况的样地数据使模型更具有代表性。

      研究采用了改进的NDVI像元二分模型,利用处理后的鱼眼相机照片结果作为真实值,进行郁闭度的求算与估计,并利用少量的样地点进行模型的检验。此方法的技术路线简单,有科学依据,可操作性强,精度能满足需要,针对落叶松人工林具有一定的推广价值。本研究结论如下:

      (1)在利用像元二分模型估测落叶松人工林郁闭度时,用距离NDVImax及NDVImin为2%的位置处的NDVI值设为NDVI′ v和NDVI′ s的值作为参数构建模型精度最佳,R2和RMSE分别为0.838和0.045。

      (2)佳木斯市孟家岗林场落叶松人工林郁闭度大多都在0.7以上,属于密林地。郁闭度呈现东高西低,北高南低的分布状态。

      (3)根据不同地点的郁闭度和植被覆盖度的差异可以判断此地的林下情况,植被覆盖度值普遍比郁闭度的值高,因为植被覆盖度是植被的覆盖率而郁闭度是树冠的覆盖率。例如,研究区域的东北部的小长龙,如果植被覆盖度与郁闭度都偏大,在0.7 ~ 0.8之间,可推测该地为密郁闭度林,且林下植被较少,几乎都为森林,郁闭度接近植被覆盖度,而研究区域的西部,部分区域郁闭度较小,而植被覆盖度较大,可推测此地为草地、灌木等分布较多。

参考文献 (22)

目录

    /

    返回文章
    返回