2006-2015年蒙古国中东部地区产草量遥感反演数据集

李舸，王卷乐

1．山东理工大学，山东淄博 255049

2．中国科学院地理科学与资源研究所，资源与环境信息系统国家重点实验室，北京 100101

3．江苏省地理信息资源开发与利用协同创新中心，南京 210023

数据库（集）基本信息简介

数据库（集）名称 2006-2015年蒙古国中东部地区产草量遥感反演数据集数据作者 李舸、王卷乐数据通信作者 王卷乐（wangjl@igsnrr.ac.cn）数据时间范围 2006-2015年地理区域 蒙古国中央省、乌兰巴托市、肯特省、戈壁苏木贝尔省、中戈壁省和肯特省（42°38'-50°30'N，101°27'-113°51'E）。空间分辨率 250 m数据量 113 MB数据格式 Geotiff数据服务系统网址 http://www.sciencedb.cn/dataSet/handle/691基金项目 中国科学院A类战略性先导科技专项（XDA2003020302）；中国科学院“十三五”信息化专项科学大数据工程项目（XXH13505-07）。数据库（集）组成本数据集由2部分组成：（1）“Dataset Result.Zip”包含2006-2015逐年的蒙古国中东部 6省产草量空间分布数据。数据命名格式为YYYY.tif，YYYY 代表年份。数据量为 89.6 MB。（2）“Dataset Map.Zip”包含2006-2015逐年的蒙古国中东部6省产草量空间分布图。数据命名格式为YYYY.jpg。数据量为24 MB。

引 言

草地是一种可更新的自然资源，是陆地生态系统的主要组成部分，在气候调节、水土保持和防风固沙等方面发挥着重要作用[1-2]。产草量是衡量草地生产力的重要指标，是草原生态状况最直接的反映。蒙古国位于亚洲大陆内部，是蒙古高原的重要组成部分，其最大的生物群落为草原，约占陆地国土总面积的83.4%[3]。由于草地生态系统稳定性较差，易受外界因素影响，近年来，在气候变化和过度放牧等自然和人文因素多重作用下，蒙古国草地生态环境日益恶化，局部地区草地呈现退化趋势。表现为牧草生产力下降、水土流失等，生态屏障作用也在减弱，对蒙古高原生态环境及农牧业可持续发展构成严重威胁。及时准确地掌握产草量的时空分布及变化状况是科学合理利用、管理草地的关键，也是草原畜牧业和生态环境可持续发展的重要条件[4-6]。

目前，国内外学者对于蒙古高原产草量的研究中，多数是针对中国内蒙古地区或新疆自治区局部区域的草地监测，对蒙古国产草量的研究较少且缺乏长时间序列的连续监测[7-10]。本研究选取蒙古国中东部的6个省份为研究区，基于野外实测数据及卫星遥感数据，对该地区2006-2015年产草量进行估算，得到蒙古国中东部地区10年间产草量时空分布及变化数据集。

1 数据采集和处理方法

1.1 区域范围

本研究选取的蒙古国中东部地区6个省份分别为中央省、乌兰巴托市、肯特省、戈壁苏木贝尔省、中戈壁省和前杭爱省，研究区位置及样点分布见图1。研究区位于亚洲大陆腹地，42°38'-50°30'N，101°27'-113°51'E，是典型的干旱半干旱区，属于大陆性温带草原气候，季温差和日温差较大，降水少且主要集中在夏季。地势西高东低，平均海拔约1400 m，草地特征明显，主要草地类型包括草甸草地、典型草地和荒漠草地。区域面积为30.76万平方公里，约占蒙古国面积的19.6%，人口178万，占全国人口的一半以上。草原畜牧业发达，牲畜量1929万头只，约占国家牲畜总量的30%。

图1 研究区位置与样点分布图

1.2 基础数据

1.2.1 遥感数据

遥感数据选取MODIS传感器陆地标准产品MODIS13Q1和MODIA17A2H中的归一化植被指数（NDVI，Normalized Difference Vegetation Index）、增强型植被指数（EVI，Enhanced Vegetation Index）和净光合（PSNnet，Net Photosynthesis），并通过计算得到了修改型土壤调节植被指数（MSAVI，Modified Soil-Adjusted Vegetation Index）。NDVI、EVI和MSAVI的空间分辨率为250 m，时间分辨率为16 d；PSNnet空间分辨率为500 m，时间分辨率为8 d。数据源为美国国家航空航天局的MODIS Web 网站（https://modis.gsfc.nasa.gov/）。

1.2.2 实测与辅助数据

本研究中地面样方数据通过本课题组和蒙古国立大学联合开展的野外调查获取，采样时间为2013年和2014年牧草生长旺盛期的8月，通过齐地面刈割法对产草量进行采集。样地大小为10 m×10 m，在其内部随机选取3个0.5 m×0.5 m的样方，用GPS记录样点坐标，取样方地上部分鲜草称重并记录，3个样方平均值为该样地产草量。经后期剔除筛选处理后，最终获取产草量有效样方数据共计29个。社会经济数据从蒙古国国家统计局[11]（http://www.1212.mn/）下载获取。蒙古国行政边界数据、气象数据和草地类型等数据由蒙古国立大学提供。

1.3 数据处理

MODIS数据预处理：MODIS 13Q1和MODIA 17A2H产品为HDF格式，通过MRT（MODIS Reprojection Tool）软件批处理功能对数据进行拼接处理。本研究获取的MODIS数据产品投影方式为正弦投影（Sinusoidal），为方便数据查看与后期处理，利用MRT软件对数据进行重投影，转换为UTM（Universal Transverse Mercator Projection）投影。

植被指数参数提取：MODIS数据经格式转换处理后得到TIF格式的数据文件，利用ArcGIS中Analysis Tools模块下Clip工具叠加研究区行政边界数据图层，对重投影的栅格数据进行裁剪处理。通过MVC方法得到月合成NDVI和EVI植被指数，以及净光合PSNnet。MSAVI指数则通过RED波段和NIR波段计算得到。根据MODIS产品用户手册中各指数有效值域范围指南，将超出值域范围的无效值赋值为空值。将地面样点数据矢量化，并将经ArcGIS空间插值后形成的降水数据与各遥感指数栅格图层进行叠加，提取对应样点像元值。随机选取 80%的样点数据用于建模，20%的样点数据用于精度验证。

遥感指数计算方式如公式（1）-（4）所示：

式中，ρNIR、ρRED、ρBLUE分别为近红外波段、红光波段、蓝色波段的反射率，GPP为植物总初级生产力，Leaf_MR为叶自养呼吸消耗Froot_MR为根自养呼吸消耗。

产草量估算模型建立：基于预处理后的野外实测数据和降水数据分别建立其与 NDVI、EVI、MSAVI和PSNnet 4种遥感指数之间的线性、指数和多元拟合关系，得到12种回归模型。通过决定系数、平均相对误差（REE）、平均绝对误差（RMSE）等指标对模型进行精度评价，见公式（5）-（6）。最终确定基于MSAVI的指数模型为最优模型，利用最优模型对研究区产草量进行估算。

式中，REE和RMSE分别表示平均相对误差和平均绝对误差；N为样点数；分别表示实测产草量、估算产草量和平均实测产草量，单位为g/m2。

2 数据样本描述

经过数据处理，得到2006-2015年蒙古国中东部地区6省产草量时空分布及变化数据集，数据空间分辨率为250 m，单位为kg/hm2。利用ArcGIS中Spatial Analyst模块下Reclassify工具对数据进行重分类，划分为6个等级描述产草量分布状况。图2展示了研究区6省10年间产草量时空分布。

图2 2006-2015年研究区产草量空间分布图

3 数据质量控制和评估

3.1 质量控制

为提高数据集精度和质量，本研究在样点数据采集、数据预处理、模型构建与精度验证等过程都进行了严格的质量控制。地面数据采样路线对不同类型草地均可实现良好覆盖，采样时考虑随机性和代表性，采样方法遵循相应调查技术规范。遥感数据产品均已经过辐射校正、几何校正、去云等处理，参考MODIS Web官方发布的MODIS数据指南对各遥感指数进行处理。利用两种时空分辨率的遥感数据参与产草量估算研究，两种数据单独建模，通过分析线性、指数和多元函数3种类型共计12个模型的精度指标，选取基于250 m分辨率植被指数的最优模型进行产草量反演。

3.2 质量评估

本研究野外采样点29个，随机选取其中的20%样点用于精度验证，80%样点用于建模。通过产草量估算数据有效性验证结果（表1）可以看出，所有模型均通过极显著性检验（Sig<0.01）。R2和RMSE分别为决定系数和平均绝对误差，R2越大，RMSE越小，说明模型的拟合效果越好。各类模型中指数函数模型精度最高，为72.75%；其次为多元线性模型，为70.25%；一元线性函数模型精度最低，为69.75%。各模型相关系数显示，指数函数模型R2为0.67，明显高于多元线性模型的0.59和一元线性模型的 0.56。结果表明指数函数模型更适合该地区产草量估算研究。对各指数模型进一步分析发现，基于MSAVI的指数模型模拟效果最佳，R2为0.72，RMSE为279.09 kg/hm2，模拟精度达到78%，可用于本研究区内产草量时空分布研究。

表1 产草量遥感估算有效性验证结果

4 数据使用方法和建议

本数据集为2006-2015年蒙古国中东部地区6省产草量时空分布和变化的公开共享数据产品，数据为TIF格式，空间分辨率为250 m。数据集资料来源可靠，主要基于ArcGIS等软件平台生产完成。用户可通过ArcGIS软件或ENVI图像处理软件打开使用，如需其他格式的数据文件，可借助以上软件自行转换。通过叠加研究区行政边界数据和草地类型数据统计不同区域的栅格影像值，可得到各省及各草地类型产草量的年际变化状况。由于蒙古国境外采样困难，本数据集的样点数量有限，但通过持续积累预期能够不断增加样本数据，并继续提高反演质量，为蒙古国或蒙古高原的产草量和植被分布及变化提供相关参考资料。