2000–2017年中巴經濟走廊逐年荒漠化分布數據集

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年中巴經濟走廊逐年荒漠化分布數據集
數據作者	敏玉芳，馮克庭，康建芳，艾鳴浩
數據通信作者	敏玉芳 ([email protected]）
數據時間範圍	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 ([email protected])
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.

引 言

荒漠化包括氣候變化和人類活動在內的種種因素造成的乾旱、半乾旱和乾旱亞濕潤地區的土地退化，屬於破壞性的生物地理過程，通常會造成生物多樣性降低，土壤肥力下降乃至生態承載力喪失。土地荒漠化是目前世界上最為嚴重的生態環境問題之一，據統計已有100多個國家正受到土地荒漠化的影響，荒漠化不僅破壞生態環境，也會削弱社會和經濟的發展。監測土地荒漠化動態變化，掌握其變化規律對防治荒漠化有很重要的意義。

中巴經濟走廊起點在中國新疆喀什，終點在巴基斯坦瓜達爾港，北接「絲綢之路經濟帶」、南連「21世紀海上絲綢之路」，是貫通南北絲路的關鍵樞紐，是一條包括公路、鐵路、油氣和光纜通道在內的貿易走廊，也是「一帶一路」的重要組成部分。中巴經濟走廊全長約3000公里，其中1300餘公里為嚴寒高原、乾旱和荒漠區段。尤其是在南段，乾旱和大面積荒漠是其主要的生態環境約束因素。

遙感技術為人們提供了一種全新的荒漠化監測手段，它具有觀測範圍廣、信息量大、數據更新快和精度高的特點，通過遙感圖像解譯或者定量反演都可以準確及時地獲取荒漠化土地變化的信息。荒漠化在遙感影像上表現為裸地地表信息的增強和植被信息的減弱，可以採用地表反照率、地表溫度、地表濕度、植被指數、植被覆蓋度等指標因子表徵。歸一化植被指數（NDVI）是反映地表植被狀態的重要生物物理參數，而地表反照率（Albedo）則是反映地表對太陽短波輻射反射特性的物理參量。隨着荒漠化程度的加重，地表植被遭受嚴重破壞，地表植被蓋度降低，生物量減少，地表粗糙度增加，在遙感圖像上表現為NDVI值相應減少，地表反照率得到相應的增加。Li等經過定位觀測研究表明，當地表反照率達到一定數值時，會發生草地荒漠化，荒漠化發生的地表反照率閾值為0.3。Albedo 與 NDVI 這兩個指標因子分別通過表徵裸露地表信息、植被覆蓋度信息，與荒漠化過程建立聯繫，具有較好的相關性。曾永年等 提出不單獨依靠某一個參數，而是通過NDVI和Albedo的負相關關係構造「NDVI-Albedo」特徵空間，並計算出荒漠化差值指數（DDI），利用DDI來進行荒漠化分級。該方法使用簡單，易於獲取指標，在荒漠化的分類及分級中，較僅使用遙感光譜信息進行分類的方法精度更高，所以在近年來不同地區的荒漠化定量評價中得到了一定的應用。本文採用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研究的算法，選用短波反照率計算方法，公式如下：

\({\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的研究結論，如果在代表荒漠化變化趨勢的垂直方向上劃分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個方面，提出了荒漠化現狀、發展速率、內在危險性評價的具體定量植被，並把荒漠化按其發展程度的不同，分為輕度、中度、重度、極重度。

表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數據。感謝國家特殊環境、特殊功能觀測研究台站共享服務平台給予的觀測和數據方面指導與支持。

數據引用格式

敏玉芳, 馮克庭, 康建芳, 艾鳴浩. 2000–2017年中巴經濟走廊逐年荒漠化分布數據集[DB/OL]. 國家特殊環境、特殊功能觀測研究台站共享平台, 2018. DOI: 10.12072/casnw.046.2018.db.

本作品在「知識共享-署名 4.0 國際」協議下發表。