朱文泉博士NPP反演学习交流QQ群：810280611
1 软件准备
链接：https://pan.baidu.com/s/1aF33Hq3-hi-CkThRX6OiXw
提取码：b0nf
2 数据处理
2.1 所需数据类型
植被分类数据、NDVI数据：中国科学院资源环境科学数据中心(http://www.resdc.cn)
气象数据：中国气象科学数据共享服务网（http://cdc.cma.gov.cn/）
（1）某年研究区月NDVI时间序列，包含12个月，由MODIS数据整理得到。
（2）某年研究区月太阳总辐射，包含12个月，单位为MJ/m^2，由气象数据插值得到。
（3）某年研究区中国月降水量，包含12个月，单位为mm，由气象数据插值得到。
（4）某年研究区月平均温度，包含12个月，单位为℃，由气象数据插值得到。
（5）估算植被NPP的参数文件。
（6）研究区植被分类文件。
2.2 IDM空间差值
2.2.1 数据准备
（1）降雨数据
泾河流域各站点降雨逐月数据保存为Excel 97-2003文件（*.xls）
其中气象站点代码（stationid）气象站点属性表中的站点代码“stationid”相对应（名称可以不同，但数据类型必须相同）。
（2）气象站点数据投影坐标系转换
【system toolboxes】-【data management】-【projections and transformation】-
【feature】-【project】
2.2.2 数据操作
将气象站点数据、降雨数据excel表导入ArcGIS
鼠标右击气象站点数据【Joins and Relates】-【Join】
选择用作连接依据的字段【stationsid】,然后选择要连接到图层的表【’ 2008$’ 】,选择 excel 表中要用作连接依据的字段【stationsid】，其他参数使用缺省值，单击【OK】
为将连接的数据与地理要素永久保存在一起，将数据导出到一个新要素类中
鼠标右击气象站点数据，【data】-【export data】
输出要素名，设置存储路径，单击 OK
在执行插值处理之前，查看菜单栏 【Geoprocessing】-【Geoprocessing options】-【Results are temporary by default】，是否被勾选。本操作需要保留插值结果，故不选中。
利用Excel输出IDW批量处理语句
【ArcToolBox】-【Spatial Analyst Tools】-【Interpolation】-【IDW】-【Batch】
单击右侧的【+】增加 11 行用来输入2018年12个月的数据。
在【input point features】中选择被用来进行插值的离散点数据 jiangyu；在【Z value field】中选择要参与插值的字段；在【Output cell size】中设置输出表面的栅格大小为 “500”；在 【Power】中输入 IDW 的幂值，幂值是一个正实数，其缺省值为 2；【Search radius】为缺省值，VARIABLE 12；【nput barrier polyline features】为可选项；点击【OK】。
在插值窗口中，单击 【Environments …】-【Processing Extent】，在 【Extent 】的下拉菜单中选择 【Same as layer jh_ploygen】，单击 OK，如图所示。
处理结果为
2.3 ENVI波段叠加
启动 ENVI5.3 软件，选择 【File】-【Open】，打开 12 个月的温度插值数据。
【Toolbox】工具箱中，双击【Raster Management】-【IDL】-【Layer stacking】进行波段叠加。
打开 【Layer Stacking Parameters】窗口。单击【Import File…】按钮，弹出 【Layer Stacking Input File】窗口，选中 12 个月的温度插值数据，单击【OK】。输入的文件将出现在 【Selected Files for Layers Stacking】列表中，如果不一致可通过 【Reorder Files…】按钮来调整波段顺序。数据加载进来后，会自动读出图像投影信息和像元大小。【Resampling】方法使用【Cubic Convolution】。在【output file range】，选择输出文件范围；如果选择 Inclusive，输出图像的地理范围将是所有输入文件范围的并集；如果选择Exclusive，输出图像的地理范围仅包含所有输入文件的重叠范围。此处选用缺省值 【Inclusive】。设置文件名及存储路径，单击【OK】。波段叠加后，band1 表示 1 月的气温数据，……band12 表示 12 月的气温数据。
3 插件运行
打开ENVI CLASSIC,导入数据
设置结果保存路径，成果会输出到该文件夹下面。默认情况下，会输出年度植被净初级生产力（npp_sum）和植被年平均覆盖率（veg_cov_mean），若勾选了“同时输出 NPP 及植被覆盖度时间序列文件”复选框，则还会额外输出每个月的植被 NPP及植被覆盖度。
4 学习笔记
（1）所有栅格文件的投影类型、空间分辨率及范围必须一致。
（2）静态参数文件可参考excel表格文件“MODIS_npp_para_file_植被类型_20210317”。
（3）植被分类文件可将示例分类数据进行投影类型转换、重采样、裁剪获得。（CASA模型\example\data\class_2006_albers）