留言板

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

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

基于Google Earth Engine与机器学习的黄土梯田动态监测

李万源 田佳 马琴 金学娟 杨泽康 杨鹏辉

李万源, 田佳, 马琴, 金学娟, 杨泽康, 杨鹏辉. 基于Google Earth Engine与机器学习的黄土梯田动态监测[J]. 浙江农林大学学报, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
引用本文: 李万源, 田佳, 马琴, 金学娟, 杨泽康, 杨鹏辉. 基于Google Earth Engine与机器学习的黄土梯田动态监测[J]. 浙江农林大学学报, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
LI Wanyuan, TIAN Jia, MA Qin, JIN Xuejuan, YANG Zekang, YANG Penghui. Dynamic monitoring of loess terraces based on Google Earth Engine and machine learning[J]. Journal of Zhejiang A&F University, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
Citation: LI Wanyuan, TIAN Jia, MA Qin, JIN Xuejuan, YANG Zekang, YANG Penghui. Dynamic monitoring of loess terraces based on Google Earth Engine and machine learning[J]. Journal of Zhejiang A&F University, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673

基于Google Earth Engine与机器学习的黄土梯田动态监测

doi: 10.11833/j.issn.2095-0756.20200673
基金项目: 国家自然科学基金资助项目(31960330)
详细信息
    作者简介: 李万源(ORCID: 0000-0002-6185-9762),从事水土保持研究。E-mail: 626251037@qq.com
    通信作者: 田佳(ORCID: 0000-0001-7732-4925),副教授,博士,从事水土保持研究。E-mail: 517153555@qq.com
  • 中图分类号: S157.2

Dynamic monitoring of loess terraces based on Google Earth Engine and machine learning

  • 摘要:   目的   梯田是黄土高原最重要的水土保持措施和农业生产措施,高效、准确地获取长时间序列黄土梯田分布信息,对黄土高原的水土保持监测和评价十分重要。   方法   在Google Earth Engine(GEE)的支持下,以宁夏回族自治区固原市为研究区,使用遥感影像监督识别技术,对比随机森林(RF)、决策树(CART)、支持向量机(SVM)等3种机器学习算法的识别精度,探讨LandTrendr算法在长时间序列动态监测中的优化应用,最终获取固原市近30 a梯田分布信息。   结果   ①3种算法识别精度从大至小依次为随机森林、决策树、支持向量机。②使用随机森林算法识别梯田,基于样点检验总体精度达94.10%,Kappa系数达0.87,基于实地斑块检验总体精度达93.33%,Kappa系数达0.80。③ LandTrendr算法能有效校正时间序列中的错误值。④ 1988-2019年,固原市梯田面积减少了45.90%。⑤固原市西部的梯田使用时间较东部更长。   结论   采用本研究方法在GEE云平台可以高效、准确地遥感监测长时序、大尺度的黄土梯田。固原市近30 a梯田农业比例逐渐下降,促进了生态环境持续向好发展。图4表3参22
  • 图  1  黄土梯田动态监测流程图

    Figure  1  Flowchart of dynamic monitoring of loess terraces

    图  2  3个不同位置的原始识别结果及使用LandTrendr算法处理后的概率

    Figure  2  Classification probability of the original results and the results of using LandTrendr algorithm at 3 different positions

    图  3  1988−2019年研究区梯田面积与植被覆盖度变化

    Figure  3  Variations of annual terraces area and annual mean fractional vegetation cover in the research area from 1988−2019

    图  4  1988−2019年研究区梯田使用时间分布示意图

    Figure  4  Distribution of time to use terraces in the research area from 1988−2019

    表  1  光谱指数计算方法

    Table  1.   Calculation methods of spectral index

    光谱指数名称计算方法
    归一化植被指数BNDVI = (BnirBr)/(Bnir + Br)
    增强型植被指数BEVI = (BnirBr)/(Bnir + 6Br−7.5Bb + 1)
    归一化建筑指数BNDBI = (Bswir2Bnir)/(Bswir2 + Bnir)
    归一化湿度指数BNDMI = (BnirBswir1)/(Bnir + Bswir1)
    归一化水体指数BNDWI = (BgBnir)/(Bg + Bnir)
      说明:Br为红波段;Bg为绿波段;Bb为蓝波段;Bnir为近红     外;Bswir1为短波红外1;Bswir2为短波红外2
    下载: 导出CSV

    表  2  不同机器学习算法识别结果的样点验证精度

    Table  2.   Sample points verification accuracy of the results of different machine learning algorithms

    机器学习
    算法
    梯田的生产
    者精度/%
    梯田的用户
    精度/%
    总体精
    度/%
    Kappa
    系数
    随机森林 94.4689.0394.100.87
    决策树  78.8978.0784.400.66
    支持向量机70.8867.3678.020.52
    下载: 导出CSV

    表  3  原始识别结果及使用LandTrendr算法处理后的斑块验证精度(随机森林)

    Table  3.   Speckles verification accuracy of the original results and the results of using LandTrendr algorithm(RF)

    去除交界
    缓冲区/m
    验证像元数/
    (×104个)
    LandTrendr
    处理
    梯田的生产者
    精度/%
    梯田的用户
    精度/%
    总体精
    度/%
    Kappa
    系数
    020.66处理前78.2473.4784.730.65
    处理后80.8276.1886.380.68
    10014.43处理前77.5584.9392.040.76
    处理后81.7585.9793.330.80
      说明:去除交界缓冲区是指去除梯田与其他类型交界线缓冲区范围内的像元,减少有地理配准误差较大的像元输入。0 m代表不     去除
    下载: 导出CSV
  • [1] 杨艳芬, 王兵, 王国梁, 等. 黄土高原生态分区及概况[J]. 生态学报, 2019, 39(20): 7389 − 7397.

    YANG Yanfen, WANG Bing, WANG Guoliang, et al. Ecological regionalization and overview of the Loess Plateau [J]. Acta Ecol Sin, 2019, 39(20): 7389 − 7397.
    [2] 中华人民共和国水利部. 第1次全国水利普查水土保持情况公报[R/OL]. (2013-11-21)[2021-04-06]. http://www.mwr.gov.cn/zw/slbgb/201612/t20161222_775434.html.
    [3] 马荣亮, 朱立军, 杨晓珍, 等. 茹河流域水土保持工程隔坡梯田应用研究[J]. 水土保持通报, 2008, 28(6): 135 − 137.

    MA Rongliang, ZHU Lijun, YANG Xiaozhen, et al. Study of the bench terrace as soil and water conservation engineering in Ruhe watershed [J]. Bull Soil Water Conserv, 2008, 28(6): 135 − 137.
    [4] 胡望舒, 洪辉, 周侃, 等. 黄土丘陵区水平梯田与农民收入的关系及原因: 以宁夏回族自治区固原市为例[J]. 干旱区地理, 2013, 36(3): 536 − 544.

    HU Wangshu, HONG Hui, ZHOU Kan, et al. Relationship of level terrace and farmers’ income and its causes in loess hilly region: a case study in Guyuan [J]. Arid Land Geogr, 2013, 36(3): 536 − 544.
    [5] FENG Xiaoming, FU Bojie, LU Nan, et al. How ecological restoration alters ecosystem services: an analysis of carbon sequestration in China’s Loess Plateau [J]. Sci Rep, 2013, 3(1): 2846 − 2850. doi:  10.1038/srep02846
    [6] 余峰, 李月祥. 宁夏退耕还林工程研究[M]. 银川: 阳光出版社, 2012: 6 − 13.
    [7] 党恬敏, 穆兴民, 孙文义, 等. 高分辨率遥感影像梯田快速提取方法研究进展[J]. 人民黄河, 2017, 39(3): 85 − 89, 94. doi:  10.3969/j.issn.1000-1379.2017.03.019

    DANG Tianmin, MU Xingmin, SUN Wenyi, et al. Review of quickly discriminating approaches of terrace information based on high resolution remote sensing images [J]. Yellow River, 2017, 39(3): 85 − 89, 94. doi:  10.3969/j.issn.1000-1379.2017.03.019
    [8] XIONG Jun, THENKABAIL P S, TILTON J C, et al. Nominal 30-m cropland extent map of continental Africa by integrating pixel-based and object-based algorithms using Sentinel-2 and Landsat-8 data on Google Earth Engine [J]. Remote Sensing, 2017, 9(10): 1065 − 1091. doi:  10.3390/rs9101065
    [9] 刘凯, 彭力恒, 李想, 等. 基于Google Earth Engine的红树林年际变化监测研究[J]. 地球信息科学学报, 2019, 21(5): 731 − 739. doi:  10.12082/dqxxkx.2019.180354

    LIU Kai, PENG Liheng, LI Xiang, et al. Monitoring the inter-annual change of mangroves based on the Google Earth Engine [J]. J Geo-inf Sci, 2019, 21(5): 731 − 739. doi:  10.12082/dqxxkx.2019.180354
    [10] TELUGUNTLA P, THENKABAIL P S, OLIPHANT A, et al. A 30-m Landsat-derived cropland extent product of Australia and China using random forest machine learning algorithm on Google Earth Engine cloud computing platform [J]. ISPRS J Photogramm Remote Sensing, 2018, 144(1): 325 − 340.
    [11] 严欣荣, 张美曼, 郑亚雄, 等. 基于Sentinel-2的丛生竹林信息提取方法比较及分布特征[J]. 生态学杂志, 2020, 39(3): 1056 − 1066.

    YAN Xinrong, ZHANG Meiman, ZHENG Yaxiong, et al. Comparison of extraction methods and the distribution characteristics of cluster bamboo forest information based on Sentinel-2 [J]. Chin J Ecol, 2020, 39(3): 1056 − 1066.
    [12] GORELICK N, HANCHER M, DIXON M, et al. Google Earth Engine: planetary-scale geospatial analysis for everyone [J]. Remote Sensing Environ, 2017, 202(1): 18 − 27.
    [13] 谭深, 吴炳方, 张鑫. 基于Google Earth Engine与多源遥感数据的海南水稻分类研究[J]. 地球信息科学学报, 2019, 21(6): 937 − 947. doi:  10.12082/dqxxkx.2019.180423.

    TAN Shen, WU Bingfang, ZHANG Xin. Mapping paddy rice in the Hainan Province using both Google Earth Engine and remote sensing images [J]. J Geo-inf Sci, 2019, 21(6): 937 − 947. doi:  10.12082/dqxxkx.2019.180423.
    [14] CHEN Bangqian, XIAO Xiangming, LI Xiangping, et al. A mangrove forest map of China in 2015: analysis of time series Landsat 7/8 and Sentinel-1A imagery in Google Earth Engine cloud computing platform [J]. ISPRS J Photogramm Remote Sensing, 2017, 131(1): 104 − 120.
    [15] GONG Peng, LIU Han, ZHANG Meinan, et al. Stable classification with limited sample: transferring a 30-m resolution sample set collected in 2015 to mapping 10-m resolution global land cover in 2017 [J]. Sci Bull, 2019, 64(6): 370 − 373. doi:  10.1016/j.scib.2019.03.002
    [16] GONG Peng, LI Xuecao, ZHANG Wei. 40-Year (1978−2017) human settlement changes in China reflected by impervious surfaces from satellite remote sensing [J]. Sci Bull, 2019, 64(1): 756 − 763.
    [17] 王彦武, 牛莉婷, 张峰, 等. 黄土区高标准梯田生态服务功能及其价值[J]. 水土保持学报, 2019, 33(6): 190 − 196.

    WANG Yanwu, NIU Liting, ZHANG Feng, et al. Ecological service function and its value of high-standard terrace in loess region [J]. J Soil Water Conserv, 2019, 33(6): 190 − 196.
    [18] 柴旭荣, 李明, 周义, 等. 影像的土地覆被快速分类[J]. 遥感技术与应用, 2020, 35(2): 315 − 325.

    CHAI Xurong, LI Ming, ZHOU Yi, et al. Rapid land cover classification using Landsat time series based on the Google Earth Engine [J]. Remote Sensing Technol Appl, 2020, 35(2): 315 − 325.
    [19] 王塞, 王思诗, 樊风雷. 基于时间序列分割算法的雅鲁藏布江流域NDVI(1985−2018)变化模式研究[J]. 生态学报, 2020, 40(19): 6863 − 6871.

    WANG Sai, WANG Sishi, FAN Fenglei. Change patterns of NDVI (1985−2018) in the Yarlung Zangbo River Basin of China based on time series segmentation algorithm [J]. Acta Ecol Sin, 2020, 40(19): 6863 − 6871.
    [20] 刘宝元, 刘瑛娜, 张科利, 等. 中国水土保持措施分类[J]. 水土保持学报, 2013, 27(2): 80 − 84.

    LIU Baoyuan, LIU Yingna, ZHANG Keli, et al. Classification for soil conservation practices in China [J]. J Soil Water Conserv, 2013, 27(2): 80 − 84.
    [21] 裴杰, 牛铮, 王力, 等. 基于Google Earth Engine云平台的植被覆盖度变化长时间序列遥感监测[J]. 中国岩溶, 2018, 37(4): 608 − 616.

    PEI Jie, NIU Zheng, WANG Li, et al. Monitoring to variations of vegetation cover using long-term time series remote sensing data on the Google Earth Engine cloud platform [J]. Carsol Sin, 2018, 37(4): 608 − 616.
    [22] 王立祥, 李永平, 许强. 中国粮食问题: 宁夏粮食生产能力提升及战略储备[M]. 银川: 阳光出版社, 2015: 92 − 93.
  • [1] 金嘉乐, 孙剑, 杨利斌, 董程锋, 齐锋.  杭州村镇下垫面遥感指数与地表温度的相关性 . 浙江农林大学学报, 2020, 37(3): 563-570. doi: 10.11833/j.issn.2095-0756.20190379
    [2] 杜雨菲, 吴保国, 陈玉玲.  基于机器学习算法的广西桉树适宜性研究 . 浙江农林大学学报, 2020, 37(1): 122-128. doi: 10.11833/j.issn.2095-0756.2020.01.016
    [3] 谢巧雅, 余坤勇, 邓洋波, 刘健, 范华栋, 林同舟.  杉木人工林冠层高度无人机遥感估测 . 浙江农林大学学报, 2019, 36(2): 335-342. doi: 10.11833/j.issn.2095-0756.2019.02.015
    [4] 姚智, 张晓丽.  基于WebGIS平台的森林郁闭度遥感反演信息系统研建 . 浙江农林大学学报, 2015, 32(3): 392-398. doi: 10.11833/j.issn.2095-0756.2015.03.009
    [5] 蒋云姣, 韩轶群, 李明阳, 王艺晓.  紫金山乔木多样性遥感估测及空间分布 . 浙江农林大学学报, 2015, 32(4): 509-515. doi: 10.11833/j.issn.2095-0756.2015.04.003
    [6] 陈利, 林辉, 孙华.  基于WorldView-2影像的外来物种薇甘菊入侵遥感监测 . 浙江农林大学学报, 2014, 31(2): 185-189. doi: 10.11833/j.issn.2095-0756.2014.02.004
    [7] 张倩倩, 陈健, 江洪, 唐敏忠.  结合纹理信息Hyperion高光谱影像分类 . 浙江农林大学学报, 2013, 30(6): 880-886. doi: 10.11833/j.issn.2095-0756.2013.06.012
    [8] 王开利, 王小德, 章宬.  浙江下渚湖湿地面积变化与保护 . 浙江农林大学学报, 2013, 30(5): 784-788. doi: 10.11833/j.issn.2095-0756.2013.05.023
    [9] 崔瑞蕊, 杜华强, 周国模, 徐小军, 董德进, 吕玉龙.  近30 a安吉县毛竹林动态遥感监测及碳储量变化 . 浙江农林大学学报, 2011, 28(3): 422-431. doi: 10.11833/j.issn.2095-0756.2011.03.012
    [10] 吴春争, 冯益明, 舒清态, 李增元, 武红敢, 车腾腾.  基于高空间分辨率影像的林业小班遥感区划系统设计与实现 . 浙江农林大学学报, 2011, 28(1): 40-45. doi: 10.11833/j.issn.2095-0756.2011.01.007
    [11] 陈春雷, 武刚.  多源遥感影像的最优尺度选择 . 浙江农林大学学报, 2011, 28(1): 164-172. doi: 10.11833/j.issn.2095-0756.2011.01.027
    [12] 刘晓双, 黄建文, 鞠洪波.  高空间分辨率遥感的单木树冠自动提取方法与应用 . 浙江农林大学学报, 2010, 27(1): 126-133. doi: 10.11833/j.issn.2095-0756.2010.01.021
    [13] 张艮龙, 冯益明, 贾建华, 武红敢, 李增元.  基于SPOT-5遥感影像的小班区划探讨 . 浙江农林大学学报, 2010, 27(2): 299-303. doi: 10.11833/j.issn.2095-0756.2010.02.023
    [14] 吴见, 彭道黎.  多伦县土地利用遥感信息提取技术 . 浙江农林大学学报, 2010, 27(3): 417-423. doi: 10.11833/j.issn.2095-0756.2010.03.016
    [15] 郭徵, 江洪, 陈健, 程苗苗, 江子山, 余树全, 李土生3.  基于遥感的杭州余杭森林景观格局变化 . 浙江农林大学学报, 2010, 27(1): 36-43. doi: 10.11833/j.issn.2095-0756.2010.01.006
    [16] 顾蕾, 吴春骏, 王鑫.  基于遥感的临安市土地利用变化及驱动力分析 . 浙江农林大学学报, 2009, 26(6): 870-876.
    [17] 肖宇, 姜海瑞, 薛文杰, 王淯, 胡忠军, 徐宏发.  陕西凤县林麝生境破碎化及其景观指数评估 . 浙江农林大学学报, 2008, 25(3): 331-335.
    [18] 施拥军, 徐小军, 杜华强, 周国模, 金伟, 周宇峰.  基于BP神经网络的竹林遥感监测研究 . 浙江农林大学学报, 2008, 25(4): 417-421.
    [19] 葛静茹, 秦安臣, 赵雄伟, 孙秋华, 李明, 张锐, 贾哲.  叠加分类及其在冀太行山区土地利用类型遥感解译中的应用 . 浙江农林大学学报, 2007, 24(6): 681-685.
    [20] 吕晓平, 金根明, 赵仁友, 章华祥.  鞭节华扁叶蜂生物学习性研究 . 浙江农林大学学报, 1993, 10(1): 16-22.
  • 加载中
  • 链接本文:

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

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2021/4/730

图(4) / 表(3)
计量
  • 文章访问数:  223
  • HTML全文浏览量:  50
  • PDF下载量:  172
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-25
  • 修回日期:  2021-04-06
  • 网络出版日期:  2021-08-09
  • 刊出日期:  2021-08-20

基于Google Earth Engine与机器学习的黄土梯田动态监测

doi: 10.11833/j.issn.2095-0756.20200673
    基金项目:  国家自然科学基金资助项目(31960330)
    作者简介:

    李万源(ORCID: 0000-0002-6185-9762),从事水土保持研究。E-mail: 626251037@qq.com

    通信作者: 田佳(ORCID: 0000-0001-7732-4925),副教授,博士,从事水土保持研究。E-mail: 517153555@qq.com
  • 中图分类号: S157.2

摘要:    目的   梯田是黄土高原最重要的水土保持措施和农业生产措施,高效、准确地获取长时间序列黄土梯田分布信息,对黄土高原的水土保持监测和评价十分重要。   方法   在Google Earth Engine(GEE)的支持下,以宁夏回族自治区固原市为研究区,使用遥感影像监督识别技术,对比随机森林(RF)、决策树(CART)、支持向量机(SVM)等3种机器学习算法的识别精度,探讨LandTrendr算法在长时间序列动态监测中的优化应用,最终获取固原市近30 a梯田分布信息。   结果   ①3种算法识别精度从大至小依次为随机森林、决策树、支持向量机。②使用随机森林算法识别梯田,基于样点检验总体精度达94.10%,Kappa系数达0.87,基于实地斑块检验总体精度达93.33%,Kappa系数达0.80。③ LandTrendr算法能有效校正时间序列中的错误值。④ 1988-2019年,固原市梯田面积减少了45.90%。⑤固原市西部的梯田使用时间较东部更长。   结论   采用本研究方法在GEE云平台可以高效、准确地遥感监测长时序、大尺度的黄土梯田。固原市近30 a梯田农业比例逐渐下降,促进了生态环境持续向好发展。图4表3参22

English Abstract

李万源, 田佳, 马琴, 金学娟, 杨泽康, 杨鹏辉. 基于Google Earth Engine与机器学习的黄土梯田动态监测[J]. 浙江农林大学学报, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
引用本文: 李万源, 田佳, 马琴, 金学娟, 杨泽康, 杨鹏辉. 基于Google Earth Engine与机器学习的黄土梯田动态监测[J]. 浙江农林大学学报, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
LI Wanyuan, TIAN Jia, MA Qin, JIN Xuejuan, YANG Zekang, YANG Penghui. Dynamic monitoring of loess terraces based on Google Earth Engine and machine learning[J]. Journal of Zhejiang A&F University, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
Citation: LI Wanyuan, TIAN Jia, MA Qin, JIN Xuejuan, YANG Zekang, YANG Penghui. Dynamic monitoring of loess terraces based on Google Earth Engine and machine learning[J]. Journal of Zhejiang A&F University, 2021, 38(4): 730-736. doi: 10.11833/j.issn.2095-0756.20200673
  • 宁夏回族自治区固原市位于黄土高塬沟壑区[1]。长期的过度放牧、不合理耕作,导致该地区植被稀疏、水土流失加剧[2],严重影响了当地社会经济发展和生态安全。梯田有效缓解了农业生产带来的水土流失问题[3],从20世纪80年代开始,固原市实施了大面积的坡改梯工程[4]。加之2000年开始实施的国家退耕还林还草工程[5],该地区的水土流失问题有所缓解,生态环境持续向好[6]。随着遥感技术的快速发展,如何从遥感影像中高效、准确、大尺度地获取梯田时空分布信息,对于指导农业生产、水土保持监测和防治水土流失具有重要的意义。传统的梯田遥感识别主要采用目视解译[7],该方法精度较高,但存在耗时耗力、成本高、方法复用性差等问题,目前更多用来采集机器学习(machine learning)的样本[8]。近年来,大部分学者采用面向对象或基于像元的监督识别技术,利用决策树(CART)、随机森林(RF)、支持向量机(SVM)、深度学习(DL)等[9-11]机器学习算法,先学习采集的样本,然后利用学习好的模型识别新的样本。面向对象技术较基于像元识别技术,不仅依靠地物的光谱特征,还利用像元和像元之间的关系提高识别精度,识别过程更加复杂,影像分辨率要求更高[7]。但是,无论采用哪种方法进行梯田遥感识别,基本上都是基于单机处理,普遍存在遥感数据获取困难、预处理复杂、性能限制等问题[9],难以开展大尺度的遥感识别研究。为了解决这些问题,Google公司借助其强大的计算资源与海量数据存储,推出了遥感云平台Google Earth Engine(GEE)[12]。借助该平台,研究人员可以极大扩展自身原有研究的覆盖范围,提供国家乃至全球尺度的研究成果[13]。目前,GEE在大尺度森林变化监测、土地利用类型分类、人类居住地动态监测等[14-16]方面应用广泛,但大尺度梯田遥感识别未见相关报道。为此,本研究在GEE平台支持下,利用Landsat时间序列数据和SRTM数字高程模型(digital elevation model,DEM),建立每年时间序列影像的百分位数特征。对比3种机器学习算法的分类精度大小,选择分类精度最高的识别结果,应用LandTrendr时序算法逐像元拟合修正时间序列,实现固原市1988−2019年度梯田动态监测的目的。研究结果可为黄土丘陵地区梯田的高效、准确识别和水土保持监测、评价提供参考。

    • 固原市(35°14′~36°31′N,105°19′~106°57′E)位于宁夏回族自治区南部的六盘山地区,辖原州区、西吉县、隆德县、泾源县、彭阳县,国土面积1.05 万km2。属大陆暖温带半干旱气候,年均气温6.3 ℃,年均降水量493.5 mm,降水量由东南向西北递减,年均蒸发量1 472.9 mm,年均无霜期152.0 d。域内地形南高北低,沟壑纵横,黄土丘陵面积达67.9%。地带性土壤以黑垆土为主,但严重的土壤侵蚀导致土壤母质层出露,黄绵土广布。植被总体上由东南的半湿润森林草原区向西北的干旱半干旱草原区过渡[4]

    • 黄土梯田动态监测的流程可分为4个主要功能模块:遥感数据加载、数据预处理、分类算法优选、序列优化。各模块从上到下,层层递进,最终实现黄土梯田动态监测(图1)。

      图  1  黄土梯田动态监测流程图

      Figure 1.  Flowchart of dynamic monitoring of loess terraces

    • 使用T1级别(质量最高)的Landsat地表反射率数据(surface reflectance, SR)。该数据产品已经过几何校正、辐射校正和大气校正,空间分辨率30 m,时间分辨率16 d。由于Landsat 5/7/8卫星的服务年限不同,1988−2011年使用Landsat 5影像,2012年使用Landsat 7影像,2013−2019年使用Landsat 8影像,共使用1 690景影像。

    • 采用30 m空间分辨率的数字高程模型,具体编号为SRTMGL1_003。

    • 地类仅分为梯田和其他2类。通过Google Earth Pro提供的高清历史影像,利用目视解译法采集样本数据。样本数据包括样点数据和斑块数据。样点数据按时间分为2010−2014年地类属性相同和2000年的样点,以满足Landsat 5/7/8不同卫星分别进行机器学习样本训练的需求。样点采集遵循以下原则:①在研究区生成5 km方形格网,以使样点分布均匀;②保持样点100 m以内属性相同。样点数据共2 673个,梯田样点1 040个,其他样点1 633个。斑块数据为6个随机分布的5 km×5 km正方形区域,参考Google Earth Pro中2019年厘米级高清遥感影像人工勾绘以及实地验证。

    • 选择Landsat对应卫星影像的红波段(Br)、绿波段(Bg)、蓝波段(Bb)、近红外(Bnir)、短波红外1(Bswir1)、短波红外2(Bswir2)6个光谱波段;再经裁边(坏像元)、光谱指数计算(计算方法如表1)、去云后,针对黄土梯田全年季相变化特点[17],统计每年度内时间序列影像百分位数特征融合影像[18],即逐像元对某一波段1 a内所有观测值取其10%、25%、50%、75%、90%百分位数,获得该像元位置该波段对应的5个指标波段;再与6个地形特征波段组合,即由数字高程计算得到的海拔、坡度、坡向,以及3个3×3、7×7、11×11像元窗口内地形起伏度波段。共计61个特征波段。

      表 1  光谱指数计算方法

      Table 1.  Calculation methods of spectral index

      光谱指数名称计算方法
      归一化植被指数BNDVI = (BnirBr)/(Bnir + Br)
      增强型植被指数BEVI = (BnirBr)/(Bnir + 6Br−7.5Bb + 1)
      归一化建筑指数BNDBI = (Bswir2Bnir)/(Bswir2 + Bnir)
      归一化湿度指数BNDMI = (BnirBswir1)/(Bnir + Bswir1)
      归一化水体指数BNDWI = (BgBnir)/(Bg + Bnir)
        说明:Br为红波段;Bg为绿波段;Bb为蓝波段;Bnir为近红     外;Bswir1为短波红外1;Bswir2为短波红外2
    • 3种机器学习算法为随机森林、决策树、支持向量机,GEE均有内建,可直接调用。另外,针对不同卫星分别进行机器学习,把样点数据分年度映射到对应合成影像并汇总(如Landsat 5包括2000、2010和2011年的样本),再按9∶1划分样本,90%的样本用于分类器训练,10%的样本用于精度验证。

    • LandTrendr算法将以年时间序列的值进行分割、逐段拟合、平滑[19],获取单个像元在整个研究时间段内的整体变化特征。具体介绍参考文献[19]。

    • 应用前文分类精度最高的机器学习算法,对研究区1988−2019年逐年进行梯田遥感识别。为减少极端天气和人类活动导致识别错误,利用地类在时间序列上连续、稳定的特征,使用LandTrendr算法[19]对识别结果的时间序列(概率为0~1的浮点)拟合平滑处理。参考中国水土保持措施分类[20],提取坡度>2°和坡度<25°区域的梯田,以减少沟壑地及塬地的误分。

    • 采用混淆矩阵的方法,以总体精度、Kappa系数、生产者精度和用户精度等指标作为识别精度评价依据。具体计算方法参考文献[18]。

    • 植被覆盖度(fractional vegetation cover, FVC)采用归一化植被指数和像元二分模型计算。具体计算方法参考文献[21]。

    • 表2为随机抽取的1 051个样点的验证结果。4种精度指标均为随机森林算法最高,决策树算法次之,支持向量机算法最小。随机森林算法基于样点检验的精度分别为:梯田的生产者精度94.46%、梯田的用户精度89.03%、总体精度94.10%、Kappa系数为0.87,都远大于另外2种算法。因此,后文采用随机森林机器学习算法进行梯田遥感识别。

      表 2  不同机器学习算法识别结果的样点验证精度

      Table 2.  Sample points verification accuracy of the results of different machine learning algorithms

      机器学习
      算法
      梯田的生产
      者精度/%
      梯田的用户
      精度/%
      总体精
      度/%
      Kappa
      系数
      随机森林 94.4689.0394.100.87
      决策树  78.8978.0784.400.66
      支持向量机70.8867.3678.020.52
    • 表3显示:去除交界100 m缓冲区后的验证精度高于未去除时(0 m)的验证精度。另外,经LandTrendr处理后梯田的生产者精度、梯田的用户精度、总体精度和Kappa系数分别为:81.75%、85.97%、93.33%、0.80,均大于LandTrendr处理前的验证精度。

      表 3  原始识别结果及使用LandTrendr算法处理后的斑块验证精度(随机森林)

      Table 3.  Speckles verification accuracy of the original results and the results of using LandTrendr algorithm(RF)

      去除交界
      缓冲区/m
      验证像元数/
      (×104个)
      LandTrendr
      处理
      梯田的生产者
      精度/%
      梯田的用户
      精度/%
      总体精
      度/%
      Kappa
      系数
      020.66处理前78.2473.4784.730.65
      处理后80.8276.1886.380.68
      10014.43处理前77.5584.9392.040.76
      处理后81.7585.9793.330.80
        说明:去除交界缓冲区是指去除梯田与其他类型交界线缓冲区范围内的像元,减少有地理配准误差较大的像元输入。0 m代表不     去除
    • 选择3个不同位置来展示LandTrendr算法拟合效果(图2),位置A原始识别结果在1994、2002、2004年被错误识别为其他类型,位置B原始识别结果在1997年被错误识别为其他类型,在2015年被错误识别为梯田类型。经LandTrendr算法处理后,这些错误类型均被校正。位置C原始识别结果与经LandTrendr算法处理后的结果均为其他类型,识别类型没有变化。

      图  2  3个不同位置的原始识别结果及使用LandTrendr算法处理后的概率

      Figure 2.  Classification probability of the original results and the results of using LandTrendr algorithm at 3 different positions

    • 经LandTrendr算法处理后的研究区梯田面积(图3)变化趋势更稳定,从1988年5 816.59 km2减少到2019年3 146.72 km2,年均减少90.85 km2·a−1。1988−2019年,研究区植被覆盖度则呈现不断增加的趋势,与梯田面积变化趋势相反。另外,处理前、处理后的梯田面积与植被覆盖度极显著(P<0.001)相关,其相关系数分别为−0.50和−0.75。

      图  3  1988−2019年研究区梯田面积与植被覆盖度变化

      Figure 3.  Variations of annual terraces area and annual mean fractional vegetation cover in the research area from 1988−2019

    • 图4显示了研究区1988−2019年梯田使用时间长短的分布。从整体上来看,梯田主要分布在六盘山山脉两侧,且西部的梯田使用时间较东部更长。从局部来看,南部的泾源县区域,梯田零星分布,使用时间相对较短;西部西吉县的沟谷条带、中部的六盘山山脉、北部原州区清水河的河谷冲积平原(红色部分)能明显区分出来。

      图  4  1988−2019年研究区梯田使用时间分布示意图

      Figure 4.  Distribution of time to use terraces in the research area from 1988−2019

    • 已有的梯田遥感监测研究[3-4]受限于单机处理性能和准确的历史样本采集,其研究内容往往时间短、区域小,限制了长时间序列、大尺度遥感监测的应用与发展。本研究使用模型迁移法,针对每一个传感器独立训练机器学习分类模型,减少了样本采集的难度,得以实现黄土梯田动态监测。然而,机器学习的识别精度主要受样本量、特征、机器学习算法的影响[7]。本研究利用多年采样法增加样本量,选取最优机器学习算法,得到较高的识别精度。另外,关于特征选取,我们前期使用了最大值、最小值、众数、中位数、平均数等多种特征融合方法,但识别精度均低于本研究的百分位数特征融合。而对于深度学习,我们在本地电脑使用相同样本集,多次构建深度学习模型并训练,然而识别精度也低于本研究的随机森林。最后引入LandTrendr算法逐像元拟合时间序列轨迹,有效校正了如图2中的异常值,提高了识别精度。而且,在斑块验证数据与样点采集时同样保留100 m空间误差时,消除2种利用类型相邻区域地理配准误差带来的系统错误后,基于2019年斑块检验总体精度93.33%,与样点验证总体精度94.10%相当,说明训练好的模型随时间迁移应用,识别性能不会降低。

    • 整体来看,研究区1988−2019年梯田面积呈减少趋势,植被覆盖度则呈现逐步增长趋势,梯田面积与植被覆盖度极显著相关(P<0.01),说明梯田面积减少有助于生态环境向好发展。局部来看,研究区在1988−1996年梯田面积年均减少69.02 km2·a−1,与该时期宁南山区逐步退耕还林还草时间一致;1997−2000年梯田面积年均增长91.60 km2·a−1,与该时期耕地面积增长趋势相同[22];在2001−2005年梯田面积下降较快,梯田面积年均减少250.51 km2·a−1,远高于1988−2019年年均减少速率90.85 km2·a−1,且2001−2005年植被覆盖度年均增长速率是1988−2019年植被覆盖度年均增长速度的4倍,与宁夏“退耕还林工程”生态政策大力实施的时间节点相符;从2007年开始,为巩固退耕还林工作,持续推进生态文明建设,研究区梯田面积下降减缓[6]。另外,研究区西部的梯田使用时间较东部更长,这可能与东部年降水量达650 mm,而西部年降雨量不到450 mm[4],在东部进行梯田退耕后有利于提高植被成活率有关。

    • 基于GEE云平台,使用随机森林机器学习算法与LandTrendr算法,可以高效、准确地实现长时间序列、大尺度的黄土梯田动态监测。相比1988年,研究区2019年梯田面积减少45.90%,植被覆盖度增长52.44%,说明近30 a梯田农业比例逐渐降低,生态环境持续向好发展。

参考文献 (22)

目录

    /

    返回文章
    返回