数据论文 I 区论文(评审中) 版本 ZH2
下载
1988–2017年黄山自然与文化遗产地森林扰动数据集
A dataset of forest disturbance on the site of Huangshan Natural and Cultural Heritage, 1988-2017
 >>
: 2018 - 11 - 16
: 2018 - 12 - 21
: 2018 - 12 - 21
657 4 0
摘要&关键词
摘要:基于Google Earth Engine(GEE)云计算平台,利用1988–2017年共计30期Landsat 5、Landsat 7、Landsat 8数据进行在线批量处理,结合森林扰动提取算法,提取30年的每年森林扰动变化;利用GEE云平台及ArcGIS 10.5软件进行结果统计及精细制图,采用流程化的处理方式,得到长时间序列的黄山自然与文化遗产地每年的森林扰动数据集。本数据集可用于黄山遗产地森林时空变化分析、遗产地的要素监测、遗产地的可持续发展等研究,能够清晰地反映黄山自然遗产地在长时间序列监测作用下的变化特征。
关键词:自然遗产地;森林扰动;GEE;黄山
Abstract & Keywords
Abstract: On the cloud computing platform Google Earth Engine (GEE), we performed batch processing on Landsat 5, Landsat 7, Landsat 8 data for the Huangshan Natural and Cultural Heritage site over the past 30 years in 1988-2017, and used forest disturbance extraction algorithm to extract annual variations of forest disturbance during this period. Then GEE cloud platform and ArcGIS10.5 software were used for fine mapping and result statistics. This dataset can be used for analyzing the temporal and spatial changes of forests, monitoring the heritage elements and studying the sustainable development of the Huangshan heritage site. It clearly reflects the changing characteristics of the site under long-term sequence monitoring.
Keywords: natural heritage ; forest disturbance; GEE; Huangshan
数据库(集)基本信息简介
数据库(集)名称1988–2017年黄山自然与文化遗产地森林扰动数据集
数据作者张少宇、朱岚巍、霍艳辉
数据通信作者朱岚巍(zhulw@radi.ac.cn
数据时间范围19882017年
地理区域森林扰动数据集范围为黄山自然与文化遗产地保护区范围。坐标:东经118°01'—118°17',北纬30°01'—30°18',包括遗产地的官方划定范围,以及黄山管理委员会规划范围,其中规划范围包含核心区范围。
空间分辨率30 m
数据量39.2 MB
数据格式*.shp,*.tif,*.lyr,*.jpg
数据服务系统网址http://www.sciencedb.cn/dataSet/handle/681
基金项目中国科学院A类先导专项(Grant No. XDA19030500)
数据库(集)组成本数据集主要包括1988–2017年长时间序列30年的黄山自然遗产地的森林扰动结果数据。这些数据保存为一个压缩文件(1988-2017年黄山自然与文化遗产地森林扰动数据集.zip),其中包含4个文件夹(地理要素、源数据、符号系统、制图)。
Dataset Profile
TitleA dataset of forest disturbance on the site of Huangshan Natural and Cultural Heritage, 1988-2017
Data corresponding authorZhu Lanwei (zhulw@radi.ac.cn)
Data authorsZhang Shaoyu, Zhu Lanwei, Huo Yanhui
Time range1988 -2017
Geographical scopeThe forest disturbance dataset covers the scope of the Huangshan Natural and Cultural Heritage Site (E 118°01'—118°17', N 30°01'—30°18'). It includes the officially delineated area of the site and the area under planning by the Huangshan Management Committee, the latter of which designates the core area.
Spatial resolution30 m
Data volume39.2 MB
Data format*.shp,*.tif,*.lyr,*.jpg
Data service system<http://www.sciencedb.cn/dataSet/handle/681>
Sources of fundingStrategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA19030500)
Dataset compositionThis dataset mainly includes data of the forest disturbance results across long-term sequences of Huangshan Natural Heritage Site from 1988 to 2017. The data set is saved as a compressed file (A dataset of forest disturbance on the site of Huangshan Natural and Cultural Heritage_1988-2017.zip), comprised of four folders (namely, geographical element, source data, symbol system, mapping) .
引 言
世界遗产主要包括自然遗产、文化遗产、混合遗产3种,是由联合国教科文组织和世界遗产委员会发起并确认的世界人民公认的、目前无法替代的财富,是全人类最具有代表性和最具有普遍价值的文物古迹及自然景观。其中自然遗产包含奇特的自然景观、完整的生态系统、原始的生物保护圈,突出体现了大自然的鬼斧神工、动植物的良好生活家园、生物进化的显著特征等。因此,世界遗产逐渐受到人类的广泛重视,越来越多的国家和人民加入保护自然遗产、保护世界遗产的行列当中。
黄山在1990年被评为世界自然与文化遗产,在2015年被评为世界地质公园,在2018年被列入世界生物圈保护区,因此黄山在中国乃至世界上都占据非常高的地位。本文利用遥感空间技术,在自然遗产的研究课题下针对中国的自然与文化遗产地–黄山进行长时间序列的监测研究,并利用云计算平台及森林扰动提取算法得到1988–2017年30年的森林扰动数据产品,并进行数据共享。相关的研究成果作为自然遗产的研究成果的一部分,反映了自然遗产的历史及现在的生存状况,同时也可为相关学者研究其他自然遗产地提供相关参考。
1   数据采集和处理方法
1.1   数据来源
本研究使用的云计算在线处理平台,具有高效、快速、方便等多个特点,解决了数据年份过多、研究区过大、数据下载繁琐等众多问题。在线获取的原始数据集是美国地址勘探局(http://glovis.usgs.gov/)提供的Landsat产品数据集。Landsat具有较高的分辨率(30 m)、较长的时间跨度、开放式使用等众多特点[1],对长时间序列监测研究具有显著优势。区域尺度的植被变化监测研究,主要使用的是Landsat影像数据 [2,3,4 ]。通过人工预览选取30年的较好的数据产品(表1),并进行相应的预处理得到可在线使用的1988–2017年的基础数据集。
表1   黄山自然遗产地遥感影像数据列表
序号成像日期卫星传感器轨道号
11988-06-03Landsat 5TM120-039
21989-11-29Landsat 5TM120-039
31990-10-15Landsat 5TM120-039
41991-11-03Landsat 5TM120-039
51992-11-21Landsat 5TM120-039
61993-12-26Landsat 5TM120-039
71994-09-24Landsat 5TM120-039
81995-10-13Landsat 5TM120-039
91996-04-22Landsat 5TM120-039
101997-11-03Landsat 5TM120-039
111998-04-28Landsat 5TM120-039
121999-12-27Landsat 5TM120-039
132000-10-10Landsat 5TM120-039
142001-02-15Landsat 5TM120-039
152002-09-30Landsat 5TM120-039
162003-02-05Landsat 5TM120-039
172004-04-20Landsat 7ETM120-039
182005-03-06Landsat 7ETM120-039
192006-05-20Landsat 5TM120-039
202007-01-07Landsat 7ETM120-039
212008-02-27Landsat 7ETM120-039
222009-03-09Landsat 5TM120-039
232010-03-28Landsat 5TM120-039
242011-05-18Landsat 5TM120-039
252012-04-26Landsat 7ETM120-039
262013-04-07Landsat 8OLI120-039
272014-10-25Landsat 7ETM120-039
282015-02-06Landsat 8OLI120-039
292016-03-28Landsat 8OLI120-039
302017-10-09Landsat 8OLI120-039
1.2   数据处理方法
近几十年来,许多遥感变化检测方法被提出、应用、对比分析和评价[2]。同时也出现了很多监测森林扰动的指数,如MODIS全球扰动指数(MGDI)、基于缨帽变化的扰动指数(DI)、归一化湿度指数(NDMI)、综合森林特征指数(IFZ)、归一化植被指数(NDVI)、归一化燃烧比指数(NBR)等[3]。但是每个指数具有不同的特点,针对不同的扰动类型具有不同的监测能力,同时众多指数对于不同的研究区也会出现一定的差异性。本文参考众多文献进行相关的分析,选取DI指数作为研究方法的基础,并采取因地制宜的方式进行改进,对黄山自然遗产地进行监测。
根据森林扰动监测算法的原理,基于森林扰动DI算法的基础上并加以改进,DI指数在根据对影像波段的进行缨帽变化后,提取其主要的分量进行数学运算,以获取主要的森林扰动信息,其DI算法的主要原理如下[5]
森林扰动指数DI = Bt-(Gt+Wt) (1)
标准化后的亮度Bt=(B-Bμ)/Bσ (2)
标准化后的绿度Gt =(G-Gμ)/Gσ (3)
标准化后的湿度Wt =(W-Wμ)/Wσ (4)
森林扰动指数平均值DIμ= Bμ- ( Gμ+ Wμ) (5)
扰动像元平均值DIm= (DIμ1+DIμ2 +...+ DIμn)/n 其中1…n为年份 (6)
扰动像元以及恢复像元判别:
扰动像元,连续三年DI值大于θ,即定义为扰动像元;如果DIn、DIn+1、DIn+2均大于θ,则DIn为第n年扰动值,此时DIn对应的该像元为扰动像元。恢复像元,DIn、DIn+1、DIn+2中DIn为扰动像元,DIn+1、DIn+2均小于θ,则DIn对应的该像元为恢复像元。根据扰动平均值DIm及相应结果进行阈值θ设置。θ阈值的调整首先参考DIm平均值的大小,一般往大于DIm的方向选择。本次实验中通过计算DIm值,且人工选取多个较好的持久森林样本参考其DI值,综合考虑后本实验初次θ阈值选取为4.0,以步长为0.5,分别多次进行实验,比较其检测变化的范围大小程度以及监测效果,经过多次试验后确定阈值的取值。
本文的研究路线如图1,利用Google earth engine云计算在线处理平台,先进行原始数据的预处理、去云、监督分类、NDVI滤波;然后结合landsat TOA数据产品进行缨帽变换、扰动提取、结果导出;以上工作大部分可依靠编程完成。在数据平台导出数据后,然后通过ArcGIS进行处理,同时利用哨兵数据、Google Earth高分辨率影像、实地验证进行结果分析和结果验证,最终生成1988–2017年共30年每年的森林扰动产品数据集。


图1   森林扰动研究技术路线
2   数据样本描述
2.1   数据组成
本数据集主要包括1988–2017年长时间序列30年的黄山自然遗产地的森林扰动结果数据,这些数据保存为一个压缩文件。压缩文件中包含4个文件夹(地理要素、源数据、符号系统、制图),总数据量为39.2 MB。地理要素文件夹中包含建筑、水系、边界、索道、公路、步道对应的矢量数据,文件为shp格式,可直接进行查看编辑;此边界为相关课题边界研究人员参考多方资料确定的黄山自然遗产地的边界范围;其他要素为本研究工作人员搜集数据,结合地形图,规划图矢量化而成。源数据文件夹中每个文件夹单独存储每一年的对应的TIF栅格数据文件,共有6个属性值,分别为未扰动森林、扰动森林、恢复森林、建筑裸岩和耕地,采用的UTM-WGS84坐标系;符号系统文件夹中为ArcGIS软件中属性值中的符号系统,此为模板方便赋值使用,格式为lyr可直接使用或编辑;制图文件夹中为1988–2017年30年的制图成果,格式为jpg。
2.2   数据样本
图2为黄山自然遗产地1988–2017年中的6期的森林扰动监测产品。通过GEE平台中的监督分类的方式,选取足够多的样本进行分类,主要分为森林、水体、建筑、裸岩、耕地,然后通过NDVI滤波再次筛选森林像元。森林扰动信息的提取是对森林像元做变化监测,根据连续3年的处理像元值的比较来判断是否为扰动像元或恢复像元。使用的方法是GEE平台中Join函数方法,查看当前年份前后年份的像元值是否大于θ,只要存在连续3年大于θ,均定义为为扰动像元。例如,起始年份1988年只查看当前年份及后一年数据即1989年像元值是否大于θ,以及终止年份只查看当前年份及左年份即2016年像元值是否大于θ,以此定义是否为扰动像元。森林扰动的主要点位集中在旅游路线的主干道两侧、山顶的宾馆附近、水系附近、较高海拔、以及森林与裸岩的交界处。主要的扰动类型为森林–裸岩、森林–建筑、森林–土地以及植被的自然变迁。例如,2006年黄山自然遗产地的保护面积为160.6平方公里,缓冲区面积为490平方公里,2006年的森林扰动面积接近7.6万平方米,恢复面积为6.8万平方米,扰动面积占据总覆盖面积的0.0473%,恢复面积占据总覆盖面积的0.0423%。

(a)


(b)


(c)


(d)


(e)


(f)

图2   黄山自然与文化遗产地6期森林扰动图
3   数据结果分析与精度验证
3.1   数据结果分析
根据结果统计可得1988–2017年30年中每年扰动面积和恢复面积的数量见图3。根据30年的森林扰动趋势可看出,在早期及2011年前后具有大面积的扰动,在2008年前后出现的大面积扰动,随后未出现较好的恢复。最近10年虽一直频繁出现不同程度的扰动,但整体占据比例很小,仍在自然变化的可接受范围之内。在每次的扰动之后均有一定幅度的面积恢复,可知主要以森林自然变化为主。选取每年的具体扰动点位进行比较来看,大多扰动点位出现在低等植物与高等植物、裸岩交界之处;由于数据质量的限制,不得不采用每年份最好的影像,于此就满足不了同一月份进行比较;同时由于黄山季节变化也会出现很大的差异,数据选取问题给提取精度造成相对的影响,造成一定程度的扰动信息误提取,导致提取精度一般;但鉴于此方法高效、快捷、可同时进行长时间序列分析,因此具有一定的普遍适用性和可扩展性,可优选其使用数据并用于其他地区研究。


图3   森林恢复和扰动面积列表
3.2   精度验证方式
采用辅助数据,2017年哨兵数据,Landsat原始数据,高分影像google earth数据进行目视解译。采用目视解译的方法,根据随机分层抽样原理[6],由于google 高分影像限制,选取2008–2017年10年的每年的扰动点,每层随机抽样30个样本;未扰动森林、非森林每层随机抽取样本100个,共500个像元点。根据精度评价混淆矩阵,计算制图精度,用户精度,总体精度,误检率,漏检率,Kappa系数加以评价。总体精度为80.66%,Kappa系数为0.77,见表1。
表1   精度混淆矩阵及评价表
预测值
真实值
持续森林非森林2008200920102011201220132014201520162017总计用户精度误检率%
持续森林946100.094.006.00
非森林397100.097.003.00
2008322530.0083.3316.67
200982230.0073.3326.67
201011126130.0086.6613.34
2011812130.0070.0030.00
201251122130.0073.3326.67
20136122130.0073.3326.67
201452530.0083.3316.67
2015112830.0093.336.67
201612930.0096.673.33
201711332230.0073.3326.67
总计13511026222721222428313222500.0
制图精度69.6388.1896.15100.096.29100.0100.091.6789.2990.3290.63100.0总体
%
kappa
漏检率%30.3711.823.800.003.700.000.002.901.719.709.400.0080.660.77
通过实地验证的方式,共选取4种典型样本点,见图4。图4a为早期扰动,现已经完全恢复,森林密集程度走势较好,左侧为该样本点DI值的走势图。由图4a可得当地森林覆盖良好;图4b为早期多次扰动,近几年逐渐恢复,由实地情况可知为较新树种覆盖良好;图4c为持续未扰动森林,一直处于较好的状态,由图可知此样本相比图4a、4b两样本明显有更好的覆盖效果;图4d为近几年扰动仍没有完全恢复样本,经过现场查看可得此处多为植被生长变化造成的自然扰动,由于是裸岩与植被混合出,容易出现监测扰动值持续较高,由此可见此处的覆盖度较差。



(a)




(b)




(c)




(d)

图4   4种典型样本点
4   数据使用方法和建议
原数据使用的TIF格式、SHP格式的文件,便于研究工作者查看和进行编辑。本文采用的是ArcGIS软件进行的编辑绘制,同时数据集中也包含了本文所固定的图层模板、便于更好的查看和制图。本数据集具有很强的代表性,可以为众多研究学者研究森林变化、植被变化、自然遗产地监测提供一定程度的参考。同时本研究数据集可直接用于黄山遗产地保护、规划、可持续发展等研究。
致 谢
非常感谢能够参与到自然遗产的保护研究工作当中,能够有机会致力于遗产的保护和发展;同时感谢杨瑞霞副研究员在验证工作中给予指导,感谢杜文杰、王普、陈亚亚等同学在研究过程中提供给予很多帮助。
[1]
汤冬梅, 樊辉, 张瑶. Landsat时序变化检测综述[J]. 地球信息科学学报, 2017, 19(08):1069-1079.
[2]
HEALEY S P , COHEN W B , SPIES T A , et al. The Relative Impact of Harvest and Fire upon Landscape-Level Dynamics of Older Forests: Lessons from the Northwest Forest Plan[J]. Ecosystems, 2008, 11(7): 1106-1119.
[3]
杨辰, 沈润平. 森林扰动遥感监测研究进展[J]. 国土资源遥感, 2015, 27(01):1-8.
[4]
张连华, 庞勇, 岳彩荣, 李增元. 基于缨帽变换的景洪市时间序列Landsat影像森林扰动自动识别方法研究[J]. 林业调查规划, 2013, 38(02):6-12+19.
[5]
Healey S P , Cohen W B , YANG Z Q, et al. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection[J]. Remote Sensing of Environment, 2005, 97(3): 301-310.
[6]
张宇欣, 倪静, 杨存建, 等. 基于扰动指数的岷江上游森林扰动时间序列研究——以理县为例[J]. 浙江林业科技, 2017, 37(03): 48-53.
数据引用格式
张少宇, 朱岚巍, 霍艳辉. 1988–2017年黄山自然与文化遗产地森林扰动数据集[DB/OL]. Science Data Bank, 2018. (2018-11-16). DOI: 10.11922/sciencedb.681.
稿件与作者信息
论文引用格式
张少宇, 朱岚巍, 霍艳辉. 1988–2017年黄山自然与文化遗产地森林扰动数据集[J/OL]. 中国科学数据, 2018. (2018-12-14). DOI: 10.11922/csdata.2018.0076.zh.
张少宇
ZHANG Shaoyu
主要承担工作:处理程序的编写、数据处理工作、结果分析、数据制图等工作。
(1994—),男,河南省驻马店人,硕士研究生,研究方向为“一带一路”自然遗产地时空变化遥感监测。
朱岚巍
ZHU Lanwei
主要承担工作:实地验证、数据下载、统计资料、相关资料的收集等工作。
zhulw@radi.ac.cn
(1980—),女,内蒙古省赤峰人,副研究员,研究方向为数字自然遗产研究、微波遥感研究、遥感应用研究等。
霍艳辉
HUO Yanhui
主要承担工作:验证思路设计、内业验证、实地验证、数据整理等工作。
(1994—),男,河北省邢台市人,硕士研究生,研究方向为“一带一路”自然要素遥感监测、多源数据珊瑚礁提取研究等。
出版历史
I区发布时间:2018年12月21日 ( 版本ZH2
参考文献列表中查看
中国科学数据
csdata