ENVI学习总结——基于改进的 CASA 模型反演NPP_envi计算gpp-程序员宅基地

技术标签: 数据库  NPP  

NPP软件下载百度云链接:NPP计算软件

基于改进的 CASA 模型反演 NPP

数据为地理空间数据云提供的 TM 影像,其具体信息如图 1 所示。该实验中所使用的 TM 数据已进行过 QUAC 快速大气校正,来消除大气和光照等因素对地物反射的影响。Landsat 主题成像仪(TM)是 Landsat4 和 Landsat5 携带的传感器,每 16 天扫描同一地区,即 16 天覆盖全球一次。 TM 影像包含 7 个波段,波段 1-5 和波段 7 的空间分辨率为 30 米,波段 6(热红外波段)的空间分辨率为 120 米

一、数据准备

在软件运行之前,需要准备以下数据月平均温度栅格文件,单位为℃,由气象数据插值得到,时间范围与 NDVI 一致。 月总降水量栅格文件,单位为 mm,由气象数据插值得到,时间范围与 NDVI 一致。

月太阳总辐射 :栅格文件:单位为 MJ/m2,由气象数据插值得到,时间范围与 NDVI 一致。

NDVI 时间序列数据:栅格文件,由遥感数据计算得到。可以是一个时间序列,如:12 个月的 NDVI 数据。

植被类型图:栅格文件,确定各植被类型的空间分布。 静态参数文件该文件配置各植被类型的 NDVI 最大值、NDVI 最小值、SR 最大值、SR 最小值以及最大光能利用率(gC/MJ)。

二、设置结果文件存放路径

设置结果存放文件夹,所有结果都会输出到该文件夹下面。默认情况下,会输出年度植被净初级生产力(npp_sum)和植被年平均覆盖率(veg_cov_mean),如果勾选了“同时输出 NPP 及植被覆盖度时间序列文件”复选框,则还会额外输出每个月的植被 NPP及植被覆盖度。

三、实验过程

1.数据来源及说明

遥感数据:来源于地理科学数据网(www.csdn.store)所提供的 2020年的 NDVI 影像。该数据集经过最大值合成法 MVC(Maximum Value Composites) 处理,空间分辨率为 1km,时间分辨率为逐月。

气象数据:来自地理科学数据网(www.csdn.store),包括 2020 年的月降水量、月平均气温和月总太阳辐射数据,共覆盖东北地区的 107 个气象站点。为保持气象和遥感数据在空间上的一致性,利用插值工具将点数据转换为空间分辨率为 1km 的面栅格数据。文中所有数据均使用以 WGS 84 为基准面的 Albers 等面积圆锥投影。

植被类型图:植被类型图为源于地理科学数据网(​​​​​​​www.csdn.store)的 2020年30m精度中国植被类型图。

地理科学数据网中富含各种地理科学数据,比如道路网数据、NPP净初级生产力数据、土地利用数据(二级分类)、乡镇边界、NDVI数据、地铁站点、气象数据(降雨量、气温、蒸散量、辐射、湿度、日照时数、风速)、径流量数据、30米土地利用二级分类数据、12.5米高精度DEM高程、2020年行政区划、10、30米NDVI数据、30m精度NPP净初级生产力、夜间灯光数据、统计年鉴、道路网、水系、景区医院等poi、三级流域矢量边界、地质灾害分布数据、30m土壤理化性质数据集、30mGDP柵格数据、30m人口柵格数据、30m精度植被类型数据、30m精度连续年份土地利用数据、GPP初级生产力数据、农田作物类型分布数据、自然保护区分布数据、高精度遥感影像数据、1:10万沙漠沙地分布矢量数据、蒸散发数据、降雨量气温蒸发数据、地表径流量、水汽压、土壤侵蚀数据等。

2.数据处理

(1)月平均温度空间插值(克里金插值法必做) 启动 ArcGIS,并加载气象(温度和降水)站点数据和东北地区边界数据,在 ArcGIS 中加载温度数据(excel 格式)。气象(温度和降水)站点数据属性表与温度数据(excel 格式)属性表建立连接。

(2)降水量插值

(3)辐射量插值

(4)启动 ENVI 软件,选择 File→Open,打开 12 个月的温度插值数据。在 Toolbox 工具箱中,双击 Raster Management→Layer stacking 进行波段叠加,打开 Layer Stacking Parameters 窗口。单击 Import File…按钮,弹出 Layer Stacking Input File 窗口,选中 12 个月的温度插值数据,单击 OK。输入的文件将出现在 Selected Files for Layers Stacking 列表中,如果不一致可通过 Reorder Files…按钮来调整波段顺序。数据加载进来后,会自动读出图像投影信息和像元大小。Resampling(重采样)方法使用缺省值,即 Nearest Neighbor。单击 Inclusive 和 Exclusive,选择输出文件范围。如果选择 Inclusive,输出图像的地理范围将是所有输入文件范围的并集;如果选择 Exclusive,输出图像的地理范围仅包含所有输入文件的重叠范围。此处选用缺省值 Inclusive。设置文件名及存储路径,单击 OK。波段叠加后,band1 表示 1 月的气温数据,band2 表示 2 月的气温数据。

(5)NDVI 时间序列

在 ENVI 中打开 2020 年 NDVI 时间序列数据时,由于背景值为-9999,像元值范围在[0,1](NDVI 值范围在[-1,1],该数据已经去除由于云、积雪等影响而出现的负值),数值之间差距太大,图像呈现为黑色。在 Toolbox 工具箱中,双击 Raster Management→Edit ENVI Header 工具,在 Data Ignore Value 文本框中填入-9999,忽略背景值影响。

(6)植被类型图

根据实验需求,在 ENVI 中打开经处理后的东北地区植被类型图,如图所示。植被类型主要有:针叶林、阔叶林、针阔混交林、灌丛、草地、栽培植被、沼泽、荒漠和非植被。

(7)静态参数文件生成

引用朱文泉教授研究结果,配置 9 类植被类型的 NDVImax、NDVImin、SRmax、SRmin 和 Emax(理想状态下最大光能利用率)参数。 其中 NDVImax 和 SRmax 的计算需要东北地区植被类型图和 NDVI 时间序列最大值数据。

①NDVI 时间序列最大值计算方法:在 ENVI 中打开 NDVI 时间序列数据 ,在 Toolbox 工具箱中,双击 Band Algebra→Band Math 工具,打开 Band Math 对话框。在 Enter an expression(运算表达式输入框)中输入表达式:b1>b2>b3>b4>b5>b6>b7>b8>b9>b10>b11>b12 。单击 Add to List 按钮,将表达式添加到列表中。单击 OK 按钮,打开 Variables to Bands Pairings 对话框,为运算表达式中各个变量赋予图像文件或者图像波段。在 Variables used in expression 列表框中选择变量 b1,在 Available Bands List 中为 b1 指定一个波段。利用同样的方法分别为所有变量指定波段。单击 Choose 按钮,选择文件名及路径保存结果,单击 OK 按钮,执行运算。

b1:选择 northeast_ndvi_2001 文件的第 1 个波段

b2:选择 northeast_ndvi_2001 文件的第 2 个波段

b3:选择 northeast_ndvi_2001 文件的第 3 个波段

b4:选择 northeast_ndvi_2001 文件的第 4 个波段

b5:选择 northeast_ndvi_2001 文件的第 5 个波段

b6:选择 northeast_ndvi_2001 文件的第 6 个波段

b7:选择 northeast_ndvi_2001 文件的第 7 个波段

b8:选择 northeast_ndvi_2001 文件的第 8 个波段

b9:选择 northeast_ndvi_2001 文件的第 9 个波段

b10:选择 northeast_ndvi_2001 文件的第 10 个波段

b11:选择 northeast_ndvi_2001 文件的第 11 个波段

b12:选择 northeast_ndvi_2001 文件的第 12 个波段

② 启动 ENVI Classic 经典版本,单击实用函数→NPP 估算,在植被 NPP 计算设置窗口中,选择配置静态参数按钮。在 npp 及植被覆盖度计算静态参数设置窗口中,选择计算 NDVImax,SRmax 按钮,依次选择东北地区植被类型图和 NDVI 时间序列最大值。在输入分类精度窗口中,输入 70,单击确定。得到 NDVImax 和 SRmax 结果,如图所示。

引用朱文泉研究结果,手动输入 NDVImin、SRmin 和 Emax 数据,如表 1 所示。完成后单击导出数据按钮,设置文件名(命名规则:静态参数设置 2020)和存放路径,生成 cfg 格式的静态参数文件。

(8)NPP 反演

启动 ENVI Classic 经典版本,单击实用函数→NPP 估算,在植被 NPP 计算设置窗口中,依次选入准备好的数据,单击完成,如图所示。

程序运行完成后,会弹出 npp 及植被覆盖度计算完成提示框(图 30)。 NPP 计算结果 经过程序计算后,会生成 4 个文件,如下所示: npp_sum:年度植被净初级生产力npp_time_series:每个月的植被 NPP veg_cov_mean:植被年平均覆盖率 veg_cov_time_series:每个月的植被覆盖度

实验结果制图输出

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/Remote_Sensinger/article/details/122738551

智能推荐

mac终端运行mysql -uroot -p时出现-bash: mysql: command not found_在terminal内输入mysql -u root -p,提示“comman not found”.-程序员宅基地

文章浏览阅读7.5k次。原因:这是由于系统默认会查找/usr/bin下的命令,如果这个命令不在这个目录下,当然就找不到命令解决方法:我们需要做的就是映射一个链接到/usr/bin目录下,相当于建立一个链接文件。1. 知道mysql命令或mysqladmin命令的完整路径,比如mysql的路径是:/usr/local/mysql/bin/mysql,2. 执行命令:$sudo ln -s /usr/lo_在terminal内输入mysql -u root -p,提示“comman not found”.

差分进化算法与多对象优化问题的研究进展-程序员宅基地

文章浏览阅读888次,点赞18次,收藏18次。1.背景介绍差分进化算法(Differential Evolution, DE)是一种基于变异和重组的全局搜索优化算法,它在过去几年中得到了广泛的关注和应用。在这篇文章中,我们将深入探讨差分进化算法与多对象优化问题的研究进展,包括背景介绍、核心概念与联系、算法原理和具体操作步骤、数学模型公式详细讲解、具体代码实例和解释、未来发展趋势与挑战以及常见问题与解答。1.1 背景介绍多对象优化问...

小程序毕设作品之微信预约订座小程序毕业设计(7)中期检查报告-程序员宅基地

文章浏览阅读456次。小程序毕设作品之微信预约订座小程序毕业设计(7)中期检查报告

C语言命令编程_c语言 运行 命令行-程序员宅基地

文章浏览阅读1.5k次,点赞5次,收藏11次。C语言命令编程按下Win+R,输入CMD,进入命令行模式Ipconfig是查看本机的ip.Ping www.qq.com是查看网络是否畅通这个时候是文件被占用,退出这个应用程序就可以CD是进入一个文件夹的路径,进入路径以后,就可以直接执行System在stdlib.h中system(“C:\CMD.exe”);参数只有一个,就是输入命令行并执行所有的C程序,只要是黑窗口,都可以在命令行下执行,C黑窗口程序可以调用其他任何C语言黑窗口程序在命令行执行,调_c语言 运行 命令行

MATLAB编写自己的BP神经网络程序_编写 1 套相关的神经网络应用系统程序matlab-程序员宅基地

文章浏览阅读1.5w次,点赞32次,收藏253次。目录前言一、什么是神经网络?二、手写代码1.主代码2.归一化函数2.反归一化函数二、工具箱代码总结前言因为项目需要,自己搞了一个BP神经网络的程序,顺便做一下笔记人工智能现在的分支有1机器学习:就是通过算法从大数据中学习到规律,而后对未知的结果进行预测。 1.21监督学习:输入和输出已知,训练集中的目标是人为注入的,常见的作用是分类数据。回归分析,统计分类。 1.22非监督学习:事先没有样本,直接对数据进行建模,非监督学习最点..._编写 1 套相关的神经网络应用系统程序matlab

安装ubuntu:-程序员宅基地

文章浏览阅读38次。ubuntu

随便推点

H5技术完美实现调用手机摄像头、相册。图片上传base64,图片压缩、预览、删除以及图片旋转90度的处理--demo。_h5实现调用手机摄像头相册。完美实现图片压缩、预览、删除的功能。---demo-程序员宅基地

文章浏览阅读1w次。看了帖子上有很多关于H5调用手机摄像头、相册的例子,也有很多可靠的栗子,综合技术点针对HTML5调用用手机相机、相册的的一次代码整合。实现提供完成的Demo提供给大家,可直接应用于项目中需要图片上传的功能。因为是H5直接调用,并没有用到像Cordova或者HTML5Plus的 JavaScript库去调用手机硬件,因此在低端手机可能会调用不到摄像图、或者相册,因此这也是在H5在这里的不足,不好解决。_h5实现调用手机摄像头相册。完美实现图片压缩、预览、删除的功能。---demo

VMware 6.0 vFlash Read Cache setup-程序员宅基地

文章浏览阅读449次。VMware vFlash Read Cache allowsto leverage host local SSDs enabling a caching layer for the VMs. No other software or extra VIBs need to be installed on the host since allneces..._vmware读取外设flash

Android笔记大全链接_安卓 笔记软件 超链接-程序员宅基地

文章浏览阅读1.7k次。android笔记链接_安卓 笔记软件 超链接

ubuntu 的硬盘安装安装_ubuntu自定义硬盘安装-程序员宅基地

文章浏览阅读487次。也许是自己太笨,搞了一天用u盘安装都是重启后没有启动项,在网上找了很多的解决方法都无效,我认为肯定是启动项被写进了u盘里,所以重启才没有启动项,最终决定放弃,改用硬盘安装;将安装步骤记下来备忘 先准备两个东西EasyBCD软件和iso镜像1.安装EasyBCD_ubuntu自定义硬盘安装

java 切割字符串split用法_list.add(ss.substring(start,i));-程序员宅基地

文章浏览阅读6.4k次。java 切割字符串split用法_list.add(ss.substring(start,i));

UVa11520-程序员宅基地

文章浏览阅读498次。Problem: Fill the Square Description: 有一个矩阵,当中的元素是大写字母,但是有个前提,相邻的位置大写字母不能相同,现在这个矩阵中的某些位置已近有字母了,要你把这个矩阵填满,并且保证字典序最小。 Solution: 简单DFSDFS。但是要注意的是,如果我这种做法超时了,那就把空的位置存到一个数组中,这样dfsdfs时就可以减少栈的空间使用。 Code(JA_uva11520