在一些文献中会看到不同土地利用类型NPP的年际变化及变化趋势分析，当前商标状态为申请中。，本期演示一下如何在ENVI中根据土地利用类型图及NPP数据，提取耕地、森林、灌木及草本等不同土地利用类型的NPP值。
一.数据准备
本期数据主要用到2020年贵州土地利用类型数据、2000-2020年贵州NPP数据，数据格式：tif，坐标系：WGS_1984。在土地利用类型图中，10代表耕地，20代表森林，30代表草地，40代表灌木，50代表湿地，60代表水体，80代表人造地表，在后续掩膜制作中会用到。
二.掩膜文件的制作
要想提取不同土地利用类型的NPP，需要制作掩膜文件。
（1）打开ENVI/File/Open Image File，把土地利用类型数据加载进来：
（2）打开Basic Tools/Band Math，输入公式(B1 EQ 20)*1+(B1 NE 20)*0，点击OK。其中EQ：等于；NE：不等于，公式的意思就是把土地利用类型图20的分即森林提取出来。
注意：公式应在全英文符号下输入，否则会报错。
（3）选择土地利用类型图层，为了方便，可以点击Memory，去输出文件及文件名这一步，下同。Memory就是不输出成文件，而是作为一个Memory存在内存中，一样可以看，可以运算，但是关掉ENVI之后如果没有存，这个文件就没有了，相当于临时文件。
（4）点击File/Save File As/ Tif /Geo TIFF，把Memory文件输出到指定文件夹即为森林掩膜文件。
（5）顺便介绍一下 Band Math 常用的函数：
三.不同土地利用类型提取NPP
（1）将掩膜文件像元小调整至与NPP一致，否则没法提取不同植被类型NPP。
打开File/Open Image File，选中上步生成的森林掩膜文件及所有NPP数据。
Basic Tools/Layer Stacking/Input File加入森林掩膜文件及所有NPP数据，先把NPP数据做成一个Memory文件，点击Recorder Files按照NPP年份顺序排序，得到的Memory文件：
得到的Memory文件再和森林掩膜文件一起做一个Memory文件：
注：如果直接把森林掩膜文件和所有NPP数据只做一个Memory，后面可能会出现统计不出森林NPP值的情况。
（2）统计森林利用类型NPP值
Basic Tools/Statistics/Compute Statistics/OK
选择Histograms，点击OK，即可统计出21年森林利用类型NPP值。
点击File/Save results to text File把结果保存为文本文档，后续可以在Excel 或者 R 里作图，分析不同植被类型NPP的年际变化。
- END -
