一种基于NPP数据计算全国草原产量数据的方法技术

本发明专利技术公开了一种基于NPP数据计算全国草原产量数据的方法,包括如下步骤：S1：获取全国草原年NPP栅格数据、草原类型参数X、全国草原类型矢量数据；S2：应用ArcGIS软件将S1获取的草原类型参数X转化为与全国草原类型矢量数据一致的栅格数据X

【技术实现步骤摘要】
一种基于NPP数据计算全国草原产量数据的方法

：
[0001]本专利技术涉及遥感数据处理领域，特别涉及一种基于NPP(植被净初级生产力)数据计算草原产量数据的方法。

技术介绍
：
[0002]草原产量数据是重要的草地资源属性，通常指单位面积草原在单位时间内(一般是一年)地上部分和地下部分植被累积的生物量，即重量；对草原产量数据进行计算是对草原资源进行合理利用的基础。目前草原产量计算有两种方式：1.根据草原植被调查样点产量数据和该点归一化植被指数拟合建立估产模型，此方法需要大量人工实地调查数据作为依据且估产模型受地域性及调查时间限制，在全国区域尺度上推广难度较大且需要消耗大量人力、物力以及财力；2.基于CASA等模型估算草原产量，此类模型基于植被光合有效辐射、光能转化率、植物的呼吸消耗、植被碳密度等变量来建立估产模型，此方法同样需要大量的人工实地调查数据作为依据，获取全国草原产量数据困难。

技术实现思路
：
[0003]本专利技术的目的在于提供一种避免大量人工调查，能快速、准确获得全国草原产量数据的方法。
[0004]本专利技术由如下技术方案实施：一种基于NPP数据计算全国草原产量数据的方法，其包括如下步骤：
[0005]S1：获取全国草原年NPP栅格数据、草原类型参数X、全国草原类型矢量数据；
[0006]S2：应用ArcGIS软件将S1获取的所述草原类型参数X转化为与所述全国草原类型矢量数据一致的栅格数据X
G
；
[0007]S3：应用ArcGIS软件将S1获取的所述全国草原年NPP栅格数据转化为与所述全国草原类型矢量数据的投影类型一致的栅格数据NPP
G
；
[0008]S4：应用ArcGIS进行栅格运算，利用下式，用NPP
G
除以X
G
，得到全国草原产量栅格数据Y
G
：
[0009][0010]其中，Y为全国草原产量，单位kg/m2；X为草原类型参数，无量纲；NPP为植被净初级生产力。
[0011]优选的，S1中所述全国草原年NPP数据获取步骤如下：
[0012]S11：采用MODIS卫星MOD17A3数据产品，获取全国草原每日NPP栅格数据，空间分辨率为500m，时间分辨率为1年；
[0013]S12：将S11获取到的所述全国草原每日NPP栅格数据加载到ArcGIS软件中，使用SUM工具计算得到全国草原年NPP栅格数据。
[0014]优选的，所述全国草原类型参数X数据通过计算得到：
[0015][0016]其中，BGB为草地植被地下部分(根系)生物量；LBGB为活根系生物量占总根系生物量的比例；ANPP为植被地上部分生产力；k为植物干重中C的含量。
[0017]本专利技术的优点：所述方法引入草原类型参数X，结合MODIS卫星获取的遥感NPP栅格数据进行计算，利用ArcGIS将所述数据转化为与全国草原类型矢量数据一致的栅格数据NPP
G
、X
G
，并进行栅格运算，用NPP
G
除以X
G
，直接得到全国草原产量栅格数据Y
G
：避免了大量的人工草原植被调查和计算过程，可以快速、准确得到全国草原产量栅格数据。
附图说明：
[0018]图1为本专利技术实施例1的流程图。
具体实施方式：
[0019]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0020]实施例1：基于NPP数据计算全国草原产量数据的方法，其原理如下：
[0021]草原产量与NPP数据存在下述关系：
[0022]NPP＝ANPP+BNPP
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0023]Y*≈ANPP
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0024]其中，ANPP为植被地上部分生产力，BNPP为植被地下部分生产力，Y
*
为草原产草量；又：
[0025]BNPP＝BGB
·
LBGB
·
Turnover
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0026]Turnover＝0.0009
·
ANPP+0.25
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0027]其中，BGB为草地植被地下部分(根系)生物量；LBGB为活根系生物量占总根系生物量的比例；Turnover为草地植物根系周转值。
[0028]由此得到：
[0029][0030]经过公式(5)计算得到的草原产草量Y
*
是以kg C/m2为单位的，进一步利用下式将其转化为以kg/m2的干草产量Y：
[0031][0032]其中，k为植物干重中C的含量；由此，可以得到草原产量Y与NPP及草原类型参数X的关系：
[0033][0034][0035]如图1所示，基于NPP数据计算全国草原产量数据的方法包括以下步骤：
[0036]S1：获取全国草原年NPP栅格数据、草原类型参数X、全国草原类型矢量数据；在本实施例中，全国草原年NPP栅格数据获取步骤如下：
[0037]S11：采用MODIS卫星MOD17A3数据产品MOD17A3HGF.v006版本，获取全国草原每日NPP栅格数据，该数据产品空间分辨率为500m，时间分辨率为1年，每个栅格像元都表示所在空间位置的NPP值；
[0038]S12：将S11获取到的每日NPP栅格数据加载到ArcGIS软件中，选择ArcToolbox，点击Spatial Analyst Tools(空间分析工具)并选择Local(局部分析)工具框中的Cell Statistics(像元统计数据)工具，使用SUM(求和)工具计算年NPP栅格数据。
[0039]草原类型参数X通过代入式(8)参数计算得到，其中式(8)各参数选取不唯一，本实施例中采取草原领域研究中普遍使用的系列参数，各参数选取如下：LBGB采用0.79，k采用方精云的研究，为0.45；BGB和ANPP，根据草地类型分别设定其数值，具体如下表1：
[0040]表1.实施例1的ANPP、BGB取值
[0041][0042][0043]最终本实施例的草原类型参数X如下表2；
[0044]表2全国各大草原类型参数X值
[0045]草地类型X低地盐化草甸类1.311828345改良草地0.780534006干热稀树灌草丛类0.538253329高寒草甸草原类4.222325006高寒草甸类2.230748649高寒荒漠草原类2.593587062高寒荒本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于NPP数据计算全国草原产量数据的方法，其特征在于，其包括如下步骤：S1：获取全国草原年NPP栅格数据、草原类型参数X、全国草原类型矢量数据；S2：应用ArcGIS软件将S1获取的所述草原类型参数X转化为与所述全国草原类型矢量数据一致的栅格数据X
G
；S3：应用ArcGIS软件将S1获取的所述全国草原年NPP栅格数据转化为与所述全国草原类型矢量数据的投影类型一致的栅格数据NPP
G
；S4：应用ArcGIS进行栅格运算，利用下式，用NPP
G
除以X
G
，得到全国草原产量栅格数据Y
G
：其中，Y为全国草原产量，单位kg/m2；X为草原类型参数，无量...

【专利技术属性】
技术研发人员：张全民，高俊刚，王保林，白耀华，柴洋洋，马香云，刘洋，马淑娟，
申请(专利权)人：内蒙古小草数字生态产业有限公司，
类型：发明
国别省市：