ENVI中不同土地利用类型NPP的提取

在一些文献中会看到不同土地利用类型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 常用的函数：

种类	操作函数
基本运算	加（+）、减（-）、乘（*）、除（/）
三角函数	正弦sin（x）、余弦cos（x）、正切tan（x）
反正弦asin（x）、反余弦acos（x）、反正切atan（x）
双曲正弦sinh（x）、双曲余弦cosh（x）、双曲正切tanh（x）
关系和逻辑运算符	小于（LT）、小于等于（LE）、等于（EQ）、不等于（NE）、大于等于（GE）、大于（GT）
AND、OR、NOT、XOR
最小值运算符（<</span>）和最大值运算符（>）
其他数学函数	指数（^）和自然指数（exp（x））
自然对数（alog（x））
以10为底的对数（alog10（x））
整型取整——round（x）、ceil（x）、和floor（x）
平方根（sqrt（x））
绝对值（abs（x））

三.不同土地利用类型提取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 -

戳原文，更有料！