2000–2017年中巴经济走廊逐年荒漠化分布数据集
作者:敏玉芳 冯克庭 康建芳 艾鸣浩
2019年8月2日
本作品收录于《中国科学数据
敏玉芳, 冯克庭, 康建芳, 艾鸣浩. 2000–2017年中巴经济走廊逐年荒漠化分布数据集[J/OL]. 中国科学数据, 2019, 4(3). (2019-08-02). DOI: 10.11922/csdata.2018.0056.zh.


    摘要&关键词

    摘要:土地荒漠化是目前世界上最为严重的生态环境问题之一。监测土地荒漠化动态变化,掌握其变化规律对防治荒漠化有很重要的意义。中巴经济走廊是“一带一路”的重要组成部分,全长约3000公里,其中1300余公里为严寒高原、干旱和荒漠区段。尤其是在南段,干旱和大面积荒漠是其主要的生态环境约束因素。本文采用荒漠化差值指数(DDI)评价中巴经济走廊荒漠化程度,以归一化植被指数(NDVI)与地表反照率(Albedo)为监测指标,通过构造Albedo-NDVI特征空间,并利用Albedo和NDVI之间负相关的关系,构建DDI公式,完成2000–2017年中巴经济走廊荒漠化分类专题数据集。并通过高分辨率Landsat数据反演的植被覆盖度验证本数据集的质量和精度。以2010年数据为例,总体评价精度达到81.67%,Kappa系数为75.42%。本数据集直观反映中巴经济走廊荒漠化程度,为定量评价此区域的荒漠化严重程度提供参考,并为中巴经济走廊沿线国家进行荒漠化防治工作以及宏观决策的制定提供基础资料。

    关键词:中巴经济走廊;荒漠化分布;Albedo-NDVI特征空间;DDI荒漠化差异指数;Albedo地表反照率

    Abstract & Keywords

    Abstract: Desertification has caused severe ecological environment problems worldwide. It is very important to monitor the dynamic change of land desertification and grasp its change rule. The China-Pakistan Economic Corridor is critical part of the One Belt And One Road. The Corridor is about 3,000 kilometers long, covering more than 1,300 kilometers of permafrost, arid and desert regions. Especially its southern section is characterized by drought and large-scale desert, which constitutes the Corridor’s major ecological constraints. In this study, desertification difference index (DDI) was used to evaluate the degree of desertification in the Corridor. Based on Normalized Difference Vegetation Index (NDVI) and Albedo, we established the Albedo-NDVI space and a DDI formula through negative correlation between Albedo and NDVI. Based on DDI, we built the dataset of desert distributions along the China-Pakistan Economic Corridor during 2000 – 2017. The quality and accuracy of the dataset were verified by vegetation coverage inversed by Landsat ETM+. Taking 2010 data as an example, the overall evaluation accuracy reached 81.67%, and the Kappa coefficient was 75.42%. The dataset reflects the degree of desertification along the Corridor, providing information for desertification prevention and control in the countries around.

    Keywords: China-Pakistan Economic Corridor; desertification; Albedo-NDVI space; desertification difference index (DDI); Albedo

    数据库(集)基本信息简介

    数据库(集)名称 2000–2017年中巴经济走廊逐年荒漠化分布数据集
    数据作者 敏玉芳,冯克庭,康建芳,艾鸣浩
    数据通信作者 敏玉芳 (myf@lzb.ac.cn)
    数据时间范围 2000–2017年
    地理区域 23°47′N–40°55′N,60°20′E–80°16′E
    空间分辨率 1 km
    数据量 569 MB
    数据格式 TIF
    数据服务系统网址 http://www.crensed.ac.cn/portal/metadata/0464e947-96a0-452d-903b-a4040d6debba
    基金项目 国家科技基础条件平台建设数据共享服务项目(Y719H71006);十三五中国科学院信息化专项(XXH13506和XXH13505)。
    数据集组成 数据集包含2000–2017年18个年份中巴经济走廊Albedo、NDVI、DDI栅格数据。每个年份一个文件夹,共18个文件夹,以年份命名。每个文件夹里包含3个文件,Albedo_年(min).tif是年最小Albedo栅格数据;NDVI_年(max).tif是年最大NDVI栅格数据;DDI_年.tif是DDI栅格数据。

    Dataset Profile

    Title A dataset of desert distributions along the China-Pakistan Economic Corridor 2000 – 2017
    Data corresponding author Min Yufang (myf@lzb.ac.cn)
    Data authors Min Yufang, Feng Keting, Kang Jianfang, Ai Minghao
    Time range From 2000 to 2017
    Geographical scope 23°47′N–40°55′N, 60°20′E–80°16′E
    Spatial resolution 1 km
    Data volume 569 MB
    Data format TIF
    Data service system http://www.crensed.ac.cn/portal/metadata/0464e947-96a0-452d-903b-a4040d6debba
    Sources of funding Data sharing fundamental program for Construction of the National Science and Technology Infrastructure Platform (Y719H71006); The 13th Five-Year Informatization Plan of Chinese Academy of Sciences (XXH13506, XXH13505).
    Dataset composition This dataset contains eighteen Albedo, NDVI and DDI raster data of the China-Pakistan Economic Corridor from the year 2000 to 2017. Each folder corresponds to a specific year, and names after the year. Each folder contains three files: “Albedo_ year (min).tif” stores the smallest Albedo raster data; “NDVI_ year (Max).tif” stores the maximum NDVI raster data; “DDI_ year.tif” stores DDI raster data.


    引 言

    荒漠化包括气候变化和人类活动在内的种种因素造成的干旱、半干旱和干旱亚湿润地区的土地退化,属于破坏性的生物地理过程,通常会造成生物多样性降低,土壤肥力下降乃至生态承载力丧失[1][2]。土地荒漠化是目前世界上最为严重的生态环境问题之一,据统计已有100多个国家正受到土地荒漠化的影响,荒漠化不仅破坏生态环境,也会削弱社会和经济的发展。监测土地荒漠化动态变化,掌握其变化规律对防治荒漠化有很重要的意义。

    中巴经济走廊起点在中国新疆喀什,终点在巴基斯坦瓜达尔港,北接“丝绸之路经济带”、南连“21世纪海上丝绸之路”,是贯通南北丝路的关键枢纽,是一条包括公路、铁路、油气和光缆通道在内的贸易走廊,也是“一带一路”的重要组成部分。中巴经济走廊全长约3000公里,其中1300余公里为严寒高原、干旱和荒漠区段。尤其是在南段,干旱和大面积荒漠是其主要的生态环境约束因素。

    遥感技术为人们提供了一种全新的荒漠化监测手段,它具有观测范围广、信息量大、数据更新快和精度高的特点,通过遥感图像解译或者定量反演都可以准确及时地获取荒漠化土地变化的信息。荒漠化在遥感影像上表现为裸地地表信息的增强和植被信息的减弱,可以采用地表反照率、地表温度、地表湿度、植被指数、植被覆盖度等指标因子表征。归一化植被指数(NDVI)是反映地表植被状态的重要生物物理参数,而地表反照率(Albedo)则是反映地表对太阳短波辐射反射特性的物理参量。随着荒漠化程度的加重,地表植被遭受严重破坏,地表植被盖度降低,生物量减少,地表粗糙度增加,在遥感图像上表现为NDVI值相应减少,地表反照率得到相应的增加。Li[3]等经过定位观测研究表明,当地表反照率达到一定数值时,会发生草地荒漠化,荒漠化发生的地表反照率阈值为0.3。Albedo 与 NDVI 这两个指标因子分别通过表征裸露地表信息、植被覆盖度信息,与荒漠化过程建立联系,具有较好的相关性。曾永年等[4] 提出不单独依靠某一个参数,而是通过NDVI和Albedo的负相关关系构造“NDVI-Albedo”特征空间,并计算出荒漠化差值指数(DDI),利用DDI来进行荒漠化分级。该方法使用简单,易于获取指标,在荒漠化的分类及分级中,较仅使用遥感光谱信息进行分类的方法精度更高,所以在近年来不同地区的荒漠化定量评价中得到了一定的应用[5][6][7][4]。本文采用DDI评价中巴经济走廊荒漠化程度,研究区如图1所示。以NDVI与Albedo为监测指标,通过构造Albedo-NDVI特征空间,并利用Albedo和NDVI之间负相关的关系,构建中巴经济走廊DDI公式,并完成2000–2017年共18期中巴经济走廊荒漠化分类专题数据集,为相关部门的治理和决策提供技术支持。


    图1 研究区范围地理位置示意图


    1 数据采集和处理方法

    1.1 数据来源与预处理

    中巴经济走廊荒漠化分布数据的制备基于2种数据:MODIS植被指数产品MOD13A3数据和地表反照率产品MCD43A3数据,数据验证采用Landsat 7数据,详细列表如表1所示。


    表1 研究采用数据列表

    序号 名称 时间 来源 类型
    1 MOD13A3 2000–2017 LPDAAC 栅格
    2 MCD43A3 2000–2017 LPDAAC 栅格
    3 Landsat 7 2010 美国地质勘探局USGS 栅格


    Albedo是遥感反演中的重要参数,指单位时间、单位面积内地表全部反射辐射通量与入射太阳总辐射通量之比。随着荒漠化程度增加,地表形态发生显著变化,集中表现在地表植被覆盖率下降,土壤有机质含量降低,土壤水分减少,地表粗糙度增加,反照率上升。MCD43A3是综合了Terra和Aqua卫星的数据产品,为天尺度L3级产品,时空分辨率是day/500 m。MCD43A3产品数据包含白空反照率和局地正午太阳角的黑空反照率,且白空反照率和黑空反照率都包含1–7波段共7个窄波段,以及3个宽波段可见光(0.3~0.7\(\mathrm{ }\mathrm{\mu }\mathrm{m}\))、近红外(0.7~5.0\(\mathrm{ }\mathrm{\mu }\mathrm{m}\))、短波(0.3~5.0\(\mathrm{ }\mathrm{\mu }\mathrm{m}\))。数据来源于美国国家航空航天局(NASA)陆地过程分布式数据档案中心(LPDAAC)。数据以HDF-EOS格式存储,采用正弦曲线投影。覆盖整个中巴走廊需要6景6个波段的数据。本研究获取2000–2017年18年MCD43A3数据共39090景。首先使用MRT工具对数据进行镶嵌和重采样,将500 m分辨率的数据重采样为1 km,并将所有数据正弦曲线投影转为WGS84。然后利用中巴走廊矢量边界和Gdal的warp工具进行数据裁剪。

    NDVI是反应植物量和植被检测的重要指标,主要应用于监测植被生长状态、植被覆盖度和消除部分辐射误差等。本文中使用1km的MODIS月合成植被指数产品MOD13A3。MOD13A3数据包括1 km月合成的NDVI、增强植被指数EVI、NDVI质量文件、EVI质量文件、可见光红光波段反射率、近红外波段反射率、可见光蓝光波段反射率、中红外波段反射率、平均观测天顶角、平均太阳天顶角、平均方位角。数据来源于LPDAAC。数据以HDF-EOS格式存储,采用Sinusoidal投影。本研究获取2000–2017年18年MOD13A3L3数据共1296景。首先使用MRT工具对数据进行镶嵌和重采样,将所有数据Sinusoidal坐标转为WGS84。然后利用中巴走廊矢量边界和Gdal的warp工具进行数据裁剪,最后得到216景研究区NDVI数据。

    本文中还使用Landsat7 ETM+数据用于数据质量评估。Landsat 7 ETM+影像包括8个波段,band 1–band 5和band 7的空间分辨率为30 m,band 6的空间分辨率为60 m,band 8的空间分辨率为15 m,时间分辨率为16天。本研究中获取2010年7–8月,云量小于5%的Landsat 7影像8景。利用ENVI软件对数据进行几何校正和大气校正,NDVI计算与统计,最后利用像元二分模型进行植被覆盖度的遥感估算,得到研究区的植被覆盖度数据。

    NDVI和Albedo原始数据有异常值和像元缺失的情况,质量不是非常可靠。需要进行无效值剔除和空间插补。NDVI数据的有效范围为−2000~10000,无效值用−3000填充。MCD43A3产品中所使用的1–6白空反照率波段的有效值范围为0~1000,无效值用32766填充。本文采用反距离加权(IDW)方法实现数据的空间差值,IDW方法参数设置如下:

    (1)距离指数:通常取值范围在0.5~3之间,本文选择值为2;

    (2)搜索半径:定义对缺失像元值进行插值的输入点,选择可变搜索半径方式,插值的最邻近输入采样点数量为12。

    1.2 数据处理步骤

    1.2.1 NDVI和Albedo数据年尺度重建

    植被覆盖度的变化是荒漠化最直观的表现,需要在一年中植被最茂盛的情况下评价这一年的荒漠化程度。所以,基于年度NDVI和年度Albedo构建的荒漠化差值指数,需要以NDVI的年度最大值,Albedo的年度最小值作为基础数据。计算过程包括以下几个方面:

    (1)Albedo数据计算

    MODIS窄波段反照率向宽波段反照率的转化使用的是Liang[8][9]研究的算法,选用短波反照率计算方法,公式如下:

    \({\alpha }_{short}=0.16{\alpha }_{1}+0.291{\alpha }_{2}+0.243{\alpha }_{3}+0.116{\alpha }_{4}+0.112{\alpha }_{5}+0.081{\alpha }_{6}-0.0015\) (1)

    其中\({\alpha }_{short}\)是短波反照率,\({\alpha }_{1}–{\alpha }_{6}\)分别代表MCD43A3中\(1–6\)波段,本研究中选用了白空反照率来计算Albedo,因为白空反照率是各个入射角的积分,更接近一般意义上的地表反照率。

    (2)NDVI年度数据计算

    采用最大值合成法(Maximum Value Compositing,MVC)。该方法在遥感数据处理方面,主要用于在一定时间段内像元的数据分析重构。计算公式:

    \({NDVI}_{i}=\underset{1\le \mathrm{j}\le \mathrm{n}}{\mathrm{max}}{NDVI}_{ij}\) (2)

    其中,\(i\)为像元名,\(j\)为\(\left[1,n\right]\)时间区间的时间点,\({NDVI}_{ij}\)指第\(i\)个像元在第\(j\)时间的NDVI值。

    (3)Albedo年度数据计算

    Albedo数据的年度重建采用最小值合成法。该方法是在某一指定的时间段内,选取像元的最小值作为新的像元值。年度Albedo的计算公式为:

    \({Albedo}_{i}=\underset{1\le \mathrm{j}\le \mathrm{n}}{\mathrm{min}}{Albedo}_{ij}\) (3)

    其中,\(i\)为像元名,\(j\)为\(\left[1,n\right]\)时间区间的时间点,\({Albedo}_{ij}\)指第\(i\)个像元在第\(j\)时间的Albedo值。

    1.2.2 基于Albedo与NDVI构建荒漠化指数

    NDVI值与植被覆盖度成正比,而Albedo与植被覆盖度成反比。根据NDVI与植被覆盖度的正相关关系,以及NDVI与Albedo的负相关关系可以得到Albedo-NDVI二维空间特征图。本文以2010年数据为例,先对NDVI和Albedo数据进行归一化处理,归一化处理公式如下:

    \(\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}=\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}/10000\) (4)

    \(\mathrm{A}\mathrm{l}\mathrm{b}\mathrm{e}\mathrm{d}\mathrm{o}=\mathrm{A}\mathrm{l}\mathrm{b}\mathrm{e}\mathrm{d}\mathrm{o}/1000\) (5)

    然后在研究区域均匀选择1500个样本点,以NDVI作为X轴,Albedo作为Y轴,构建NDVI-Albedo特征空间。特征空间散点图如图2所示,根据线性统计回归分析结果,Albedo和NDVI间的负相关关系可以用以下线性关系式表示:

    \(Albedo=-0.2303\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}+0.3258\) (6)


    图2 中巴经济走廊2010年Albedo-NDVI二维特征空间图


    如图2所示,由Albedo和NDVI构造的特征空间中,Albedo与NDVI的负相关关系式的不同位置,代表荒漠化不同阶段的状态和程度,荒漠化程度随着NDVI值的减少而增加,随着Albedo值的增加而增加,即Albedo与NDVI的负相关线性表达式可以反映荒漠化的变化趋势。根据Verstraete和Pinty[10]的研究结论,如果在代表荒漠化变化趋势的垂直方向上划分Albedo-NDVI特征空间,可以将不同的荒漠化土地有效区分开来,如图3所示,用荒漠化差值指数模型DDI来表示:

    \(DDI=\alpha ×NDVI-Albedo\) (6)


    图3 NDVI-Albedo空间中荒漠化差值指数(DDI)图形表达


    在具体应用中,常数\(\alpha \)可根据Albedo-NDVI特征空间中斜率来确定,\(k\)为特征方程的斜率。

    \(\alpha ×k=-1\) (7)

    根据上面的计算结果,Albedo-NDVI的负相关关系斜率\(k=-0.2303\),则\(\alpha =4.3422\),DDI的表达式如下:

    \(DDI=4.3422×NDVI-Albedo\) (8)

    1.2.3 荒漠化等级划分

    为荒漠化评价与制图的需要,1984年联合国粮农组织(FAO)和联合国黄金规划署(UNEP)在《荒漠化评价与制图方案》中,从植被退化、风蚀、水蚀、盐碱化等4个方面,提出了荒漠化现状、发展速率、内在危险性评价的具体定量植被,并把荒漠化按其发展程度的不同,分为轻度、中度、重度、极重度。

    [11][5]


    表2 DDI分级表

    类别 DDI值
    冰雪水体 DDI≤−0.26
    重度荒漠化 −0.26<DDI≤0.12
    中度荒漠化 0.12<DDI≤0.55
    轻度荒漠化 0.55<DDI≤1.6
    非荒漠化 1.6<DDI≤4.2


    1.2.4 荒漠化指数专题数据模拟流程

    首先重建NDVI和Albedo年数据,用最大值合成法计算年度NDVI值,用最小值合成法计算年度Albedo值;然后通过拟合Albedo和NDVI数据,构建Albedo-NDVI特征空间,得到拟合斜率k;最后构建荒漠化指数,进行荒漠化等级划分,完成基于荒漠化指数的荒漠化分级图。数据处理流程如图4所示。


    图4 荒漠化指数专题数据模拟流程


    2 数据样本描述

    中巴经济走廊逐年荒漠化分布数据集有2000–2017年18个文件,每个文件夹包括年度NDVI、Albedo、DDI栅格数据和2010年荒漠化分级图示例。全部数据的空间分辨率均是1 km,年度NDVI、年度Albedo、DDI栅格数据保存格式为tif,荒漠化分级图采用jpg文件格式保存,用户可以根据DDI数据制作荒漠化分级数据。所有数据地理坐标系为WGS1984。数据结果展示如图5。


    图5 2010年中巴经济走廊荒漠化分级图


    3 数据质量控制和评估

    本文采用高分辨率数据来评价荒漠化分级数据的质量。以2010年数据为例,选取8景Landsat 7 影像数据在小范围上进行验证。8景影像空间分辨率为30 m,平均云量小于5%,数据质量良好。DDI是通过NDVI年度最大值和Albedo年度最小值计算而来,所以本文用处于植被生长旺季的8月份Landsat数据作为验证数据,反演年度植被覆盖度最大值。运用ENVI软件对数据进行预处理,主要包括配准、几何校正、图像增强和FLAASH大气校正。首先计算8景影像的NDVI值,然后进行植被盖度的估算。植被覆盖度计算公式如下:

    \(vfc=\frac{NDVI-NDV{I}_{soil}}{NDV{I}_{veg}-NDV{I}_{soil}}\) (9)

    式中:\(vfc\) 为研究区的植被覆盖度,\(NDV{I}_{soil}\)为本研究区域纯沙漠的NDVI值,\(NDV{I}_{veg}\)为完全植被覆盖区域的NDVI值。一般以5%置信度截取NDVI的上下阈值分别近似代表\(NDV{I}_{veg}\)和\(NDV{I}_{soil}\)进行计算。

    本文参考其他研究工作,将中巴经济走廊地区的荒漠化程度按植被覆盖度大概划分为:非荒漠化(高植被覆盖,\(vfc\)>60%)、轻度荒漠化(中等植被盖度,30%~60%)、中度荒漠化(中低植被覆盖,15%~30%)、重度荒漠化(低植被覆盖,<15%)。验证时,在每个影像上选取各级别荒漠化验证点30个,由于Landsat跟Modis数据分辨率不一致,所以在选取验证点时,在Landsat上选取小范围且植被覆盖度差别不大的区域作为一个验证点,与DDI分级结果进行比较(表3),并计算Kappa系数。表4是基于此方法验证的数据评价精度表,总体评价精度达到81.67%,Kappa系数为75.42%,两种分辨率的数据评价结果基本一致,反映出本文的数据在区域尺度上评价中巴经济走廊荒漠化具有可行性。


    表3 DDI分类混淆矩阵

    DDI分类\Landsat分类 非荒漠化 轻度荒漠化 中度荒漠化 重度荒漠化 总和
    非荒漠化 61 5 5 2 73
    轻度荒漠化 2 50 2 3 57
    中度荒漠化 5 2 38 4 49
    重度荒漠化 3 7 4 47 61
    总和 71 64 49 56 240


    表4 单一类别精度及总体精度分析

    DDI分类\Landsat分类 用户精度 生产者精度 总体分类精度 Kappa系数
    非荒漠化 0.8356 0.8592 0.8167 0.7542
    轻度荒漠化 0.8772 0.7813
    中度荒漠化 0.7755 0.7755
    重度荒漠化 0.7705 0.8393


    4 数据价值

    虽然荒漠化的影响因素复杂,使用目的多样化,造成指标因子繁杂,但所有的指标最终都可以归结到最基本的影响因子——植被。植被既是土地覆盖质量的综合体现,也是影响生态环境中水土保持等方面的主体因素。本文基于易获取的遥感影像数据,使用反映植被信息的NDVI与反映土壤、水分等信息的Albedo,构建中巴经济走廊的DDI模型,直观反映中巴经济走廊荒漠化程度。为定量评价荒漠化严重程度提供参考,为中巴经济走廊沿线国家进行荒漠化防治工作以及宏观决策的制定提供基础资料。

    5 数据使用方法和建议

    2000–2017年中巴经济走廊荒漠化分布相关数据保存为栅格tif格式。ArcGIS、QGIS、ENVI、ERDAS等常用的GIS与遥感软件均支持该数据的读取和操作。

    致 谢

    感谢LPDAAC提供的MODIS数据,感谢USGS提供的Landsat 7数据。感谢国家特殊环境、特殊功能观测研究台站共享服务平台给予的观测和数据方面指导与支持。

    参考文献

    1. OH K, JEONG Y, LEE D, et al. Determining Development Density Using the Urban Carrying Capacity Assessment System[J]. Landscape & Urban Planning, 2005, 73(1): 1-15.
    2. 王新源, 杨小鹏, 陈翔舜, 等. 荒漠化监测指标体系的构建——以甘肃省为例[J]. 生态经济, 2016, 32(7): 174-177, 182.
    3. SHENG GL,HARAZONO Y,OIKAWA T, et al. Grassland Desertification By Grazing and the Resulting Micrometeorological Changes in Inner Mongolia[J]. Agricultural & Forest Meteorology, 2000, 102(2): 125-137.
    4. ^ 4.0 4.1 曾永年, 向南平, 冯兆东, 等. Albedo-NDVI特征空间及沙漠化遥感监测指数研究[J]. 地理科学, 2006, 26(1): 75-81.
    5. ^ 5.0 5.1 MA Z,XIE Y,JIAO J, et al. The Construction and Application of an Aledo-ndvi Based Desertification Monitoring Model[J]. Procedia Environmental Sciences, 2011, 10: 2029-2035.
    6. 李燕,周游游,胡宝清,等. 基于3S技术的南北方典型区荒漠化演变特征对比研究[J]. 广西师范学院学报(自然科学版), 2017, 34(1): 82-90.
    7. 官雨薇. 基于遥感影像的全球荒漠化指数构建及趋势分析[D]: 成都: 电子科技大学, 2015.
    8. LIANG S. Narrowband to Broadband Conversions of Land Surface Albedo I : Algorithms[J]. Remote Sensing of Environment, 2001, 76(2): 213-238.
    9. LIANG S,SHUEY CJ,RUSS AL, et al. Narrowband to Broadband Conversions of Land Surface Albedo: Ii. Validation[J]. Remote Sensing of Environment, 2003, 84(1): 25-41.
    10. VERSTRAETE MM,PINTY B. Designing Optimal Spectral Indexes for Remote Sensing Applications[J]. Ieee Transactions on Geoscience & Remote Sensing, 1996, 34(5): 1254-1265.
    11. 毋兆鹏, 王明霞, 赵晓. 基于荒漠化差值指数(DDI)的精河流域荒漠化研究[J]. 水土保持通报, 2014, 34(4): 188-192.

    数据引用格式

    敏玉芳, 冯克庭, 康建芳, 艾鸣浩. 2000–2017年中巴经济走廊逐年荒漠化分布数据集[DB/OL]. 国家特殊环境、特殊功能观测研究台站共享平台, 2018. DOI: 10.12072/casnw.046.2018.db.