基于机器学习、遥感和Penman-Monteith方程的农田蒸散发混合模型研究_刘燕_2022

news2024/11/19 23:35:05

基于机器学习、遥感和Penman-Monteith方程的农田蒸散发混合模型研究_刘燕_2022

  • 摘要
    • 关键词
  • 1 绪论
  • 2 数据与方法
    • 2.1 数据
    • 2.2 机器学习算法
    • 2.3 Penman-Monteith方程
    • 2.4 Medlyn公式
    • 2.5 模型性能评估
  • 3 基于机器学习算法的混合模型估算农田蒸散量的评价与比较
  • 4 利用人工神经网络算法和遥感植被指数改善PM方程估算农田蒸散量的精度
  • 5 结论

在这里插入图片描述

摘要

  本文的研究内容旨在发展能够准确获取全球或区域尺度蒸散量的混合模型,主要工作包括:
  (1)构建了六种不同的混合蒸散模型。这六种混合模型基于六种传统的机器学习算法,分别是KNN、RF、SVM、XGBoost、ANN、LSTM。然后使用17个全球分布的农田通量站点的涡度协方差通量塔的通量数据构建模型。
  (2)分析了十种不同的输入因子组合对混合模型估算农田蒸散量性能的影响。将气象数据分别和一种、两种、三种以及四种遥感数据输入到机器学习方法中,构建气孔导度模型,进而运用Penman-Monteith(PM)方程计算蒸散。
  (3)评估了六种基于机器学习的混合模型在使用相同的输入变量下估算蒸散的精度,以得到一种较优的基于机器学习的估算蒸散方法。
  (4)构建Medlyn-Penman-Monteith模型。将得到的较优的混合机器学习模型和Medlyn-PM模型比较,以得到适用于不同的环境条件的模型,提升模型模拟区域尺度农田蒸散的精度,能够为区域尺度准确估算农田蒸散提供方法支持。
  (5)研究了ANN-PM模型在干旱环境下的准确性。利用ANN-PM模型的模拟值与观测值之间的相关系数(r)来评估模型的性能。通过比较ANN-PM模型模拟蒸散值与观测值之间的时间序列图,研究模型是否能捕捉干旱站点蒸散的时间序列变化。

关键词

蒸散;机器学习;遥感;Penman-Monteith方程

1 绪论

  本文的主要研究内容如下:
  (1)基于六种不同的ML算法,使用全球17个通量站点的气象数据和中分辨率成像光谱仪(Moderate-resolution Imaging Spectroradiometer, MODIS)数据构建了估算蒸散的六种混合模型,并优化模型的参数。
  (2)分析十种不同的输入变量组合对混合模型的影响。将十种输入因子分别输入到混合模型中,比较每种模型的相关系数和均方根误差来衡量不同的变量组合。进一步比较基于机器学习算法的六种混合模型使用相同的输入变量估算蒸散的性能,得到了一种性能较优的混合模型。
  (3)将上述得到的基于ML算法的混合模型与构建的Medlyn-PM模型比较,评估其估算蒸散的性能。最后研究性能较好的混合模型在干旱环境下的准确性。

2 数据与方法

2.1 数据

2.1.1 通量站点数据
  本文使用来自FLUXNET2015数据集的5个气象因子数据和观测ET数据,日尺度上时间连续的气象数据包括气温(Temperature, Ta)、降水(Precipitation, P)、大气二氧化碳浓度(Carbon dioxide concentration, Ca)、太阳辐射(Solar radiation, SW)和饱和水汽压差(Vapor Pressure Deficit, VPD)。
2.1.2 遥感数据
  MODIS数据包括归一化差异植被指数(NDVI)、增强植被指数( EVI)、植被近红外反射率(NIRv) 和 短 波 红 外 波 段 (SWIR)。

2.2 机器学习算法

  六种机器学习算法(KNN、RF、SVM、XGBoost、ANN、LSTM)

2.3 Penman-Monteith方程

  本研究构建的六种基于ML算法的混合模型和Medlyn-PM模型都是基于PM方程的物理框架,主要差异在于Gs的计算。PM方程的计算公式如下:
在这里插入图片描述
在这里插入图片描述

2.4 Medlyn公式

  基于PM方程的ET区域尺度建模的关键是要准确地量化表面电导,这是ET建模中不确定性的最大来源。因此,Medlyn等将最优气孔行为理论与光合作用模型相结合,为表征Gs提供了一个新的生物物理框架。该Gs模型的计算方法如下:
在这里插入图片描述
在这里插入图片描述
  本文将Medlyn模型与ML方法相结合,研究哪种纯粹基于ML的Gs方法或者Medlyn模型和ML的集成可以表现得更好。

2.5 模型性能评估

  r、RMSE、MAE

3 基于机器学习算法的混合模型估算农田蒸散量的评价与比较

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

4 利用人工神经网络算法和遥感植被指数改善PM方程估算农田蒸散量的精度

4.1 基于Medlyn公式的混合模型
4.2 干旱气候下ET估计模型的评估
4.3 Medlyn公式的参数优化
4.4 实验结果分析
4.4.1 ANN-PM模型与Medlyn-PM模型的比较
4.4.2 评估ANN-PM模型在干旱气候下估算蒸散的性能

5 结论

  主要结论如下:
  (1)不同模型的最优参数:基于KNN模型的K值为5,基于RF模型和XGBoost模型的max_depth分别为8和7,基于ANN模型和LSTM模型的隐含层数都为2、每层神经元数分别是48和40,基于Medlyn模型的待定系数g0=0.06、g1=3.94。
  (2)与使用一种RS因子的混合模型相比,使用两种RS因子的混合模型呈现了较好的结果(RMSE=18.60-26.29 W m-2,r=0.87-0.96);使用三种或四种RS因子的混合模型表现出的性能和使用两种RS因子的混合模型的性能相似。
  (3)基于不同的ML算法的混合模型中,ANN-PM模型表现出较强地估算蒸散的能力(RMSE=18.67-20.69 W m-2,r=0.90-0.94)。
  (4)ANN-PM模型和Medlyn-PM模型相比,ANN-PM模型在估算农田蒸散方面具有较好的性能(RMSE=19.23-19.71 W m-2,r=0.93),拟合效果好,精度较高且适用范围广。
  (5)ANN-PM模型能较好地捕捉干旱站点(DE-Rus,US-Tw2,US-Tw3和US- Twt)ET的时间序列变化。ANN-PM模型在干旱站点估算蒸散的性能是合理的,其在干旱站点预测的ET与观测值的平均r是0.87,所有站点的平均r是0.86。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1466404.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

信息矩阵、hessian矩阵与协方差矩阵

文章目录 协方差矩阵联合概率密度hessian矩阵marginalize 本节探讨信息矩阵、hessian矩阵与协方差矩阵的关系,阐明边缘化的原理。 一个简单的示例,如下: 来自 David Mackay. “The humble Gaussian distribution”. In: (2006). 以及手写vio第…

CSS基础(下)

一 CSS样式重置 【面试题】:你知道浏览器的兼容性问题有哪些?你进行过样式重置吗?进行过样式标准化吗? 样式重置reset/样式标准化normalize /*******第1步:样式重置(标准化):将浏览器提供的默认样式统一化 实用化***…

UnityWebGL 设置全屏

这是Unity导出Web默认打开的页面尺寸 修改后效果 修改 index.html 文件 1.div元素的id属性值为"unity-container",宽度和高度都设置为100%,意味着该div元素将占据整个父容器的空间。canvas元素的id属性值为"unity-canvas"&#xff…

CAN转WIFI

一、 产品概述 SG-CAN-WIFI 是专为 CAN 总线网络与无线 IP 网络(WLAN 或 Wi-Fi)之 间或多个 CAN 总线网络之间通过无线 IP 网络(WLAN 或 Wi-Fi)传输 CAN 总 线数据而设计,无线 IP 网络(WLAN 或 Wi-Fi&…

计算机网络面经-从浏览器地址栏输入 url 到显示主页的过程?

大概的过程比较简单,但是有很多点可以细挖:DNS解析、TCP三次握手、HTTP报文格式、TCP四次挥手等等。 DNS 解析:将域名解析成对应的 IP 地址。TCP连接:与服务器通过三次握手,建立 TCP 连接向服务器发送 HTTP 请求服务器…

JavaSE-05笔记【面向对象02】

文章目录 1. 类之间的关系2. is-a、is-like-a、has-a2.1 is-a2.2 is-like-a2.3 has-a 3. Object类3.1 toString()3.2 finalize()(了解即可)3.3 与 equals 方法 4. package 和 import4.1 package4.2 import4.3 JDK 常用开发包 5. 访问权限控制5.1 privat…

CLion的bundled MinGW能用在VSCode上吗?

跟着前辈做一个项目,用的极海的MCU,主要用到SPI和USB功能。 官方提供的SDK中的例程有 Eclipse/ Keil/ IAR 版本。 前辈根据Eclipse版本的工程信息文件(.project 和.cproject) 看里面链接到了哪些文件,然后自己手动写…

Ansible service 模块 该模块用于服务程序的管理

目录 参数将服务设置为 自启动检查端口关闭服务再次查看端口 参数 arguments #命令行提供额外的参数 enabled #设置开机启动。 name #服务名称 runlevel #开机启动的级别,一般不用指定。 sleep #在重启服务的过程中,是否等待。如在服务关闭以后等待2秒再…

gnss尾矿库安全监测系统是什么

【TH-WY1】GNSS尾矿库安全监测系统是一种利用全球导航卫星系统(GNSS)技术对尾矿库进行安全监测的系统。尾矿库是矿山企业的重要设施之一,用于存放矿山开采过程中产生的尾矿。由于尾矿库具有高能势和复杂的地质环境,存在溃坝、滑坡…

Python:函数

目录 前言: 一、函数的定义 二、函数的调用 三、函数的分类 四、全局变量和局部变量 五、函数的参数 5.1 位置参数 5.2 默认值参数 5.3 可变参数 5.4 关键字参数 5.5 命名关键字参数 5.6 参数的组合 六、函数的递归 前言: 函数就是一个过程…

nginx+keepalived实现nginx高可用集群以及nginx实现Gateway网关服务集群

一、前言 1、简介 Nginx作为一款高性能的Web服务器和反向代理服务器,被广泛使用。且现如今很多高并发场景需要后端服务集群部署,因此nginx也需要支持集群部署从而避免单点故障的问题。 本文将详细介绍使用 KeepalivedNginx 来实现Nginx的高可用集群和N…

Leetcode2583. 二叉树中的第 K 大层和

Every day a Leetcode 题目来源:2583. 二叉树中的第 K 大层和 解法1:层序遍历 排序 先使用层序遍历计算出树的每一层的节点值的和,保存在数组 levelSum 中。然后将数组进行排序,返回第 k 大的值。需要考虑数组长度小于 k 的边…

mysql group by分组后查询无数据补0

mysql经常会用到Group By来进行分组查询,但也经常会遇到一个问题,就是不满足条件的数据就不会显示,如图总共有五个业务,业务状态为3的就不会显示: 因此,想要实现,即使没有数据,也想让count显示出0而不是空的效果&…

day40打卡

day40打卡 343. 整数拆分 状态表示 ​ dp[i] 表示将正整数i拆分成至少两个正整数的和之后,这些正整数的最大乘积 状态转移方程 ​ i > 2 时,对正整数i拆出的第一个正整数是j,则有: 将i拆分为 j 和 i-j,且 i-j…

CUDA编程 - Nsight system Nsight compute 的安装和使用 - (1) 学习记录

Nsight system & Nsight compute 安装和下载地址一、Nsight Systems1.1、主要应用1.2、比较常用的分析1.2.1、情况一1.2.2、情况二 二、Nsight Compute2.1、主要应用2.2、比较常用的分析2.2.1、情况一2.2.2、情况二 三、两者的比较3.1、Nsight systems3.2、Nsight compute …

Runaway Queries 管理:提升 TiDB 稳定性的智能引擎

在数字化系统扮演重要角色的今天,数据库稳定性成为企业关注的核心问题。对于重要计算机系统而言,突发的性能下降可能对业务造成不可估量的损失。为了稳定数据库性能,用户可以从管理流程入手规范变更的测试,或者利用产品手段减少预…

su模型在3d里渲染效果如何---模大狮模型网

SketchUp模型在其他3D应用程序中进行渲染可以获得非常好的效果,取决于您所选择的渲染引擎和技术水平。下面是一些常见的渲染效果和技巧: 一:光照和阴影 通过添加适当的光源和调整阴影设置,可以使SketchUp模型在渲染中呈现出更加真…

python递归算法

递归算法 一、嵌套调用的过程二、递归的基本原则1、递归的基本原则2、无限递归调用3、正常递归调用4、阶乘问题5、力扣:231. 2 的幂6、力扣面试题 08.05. 递归乘法7、力扣、326. 3 的幂8、力扣342. 4的幂 一、嵌套调用的过程 def show1():print("show 1 run s…

第2讲:C语言数据类型和变量

第2讲:C语言数据类型和变量 目录1.数据类型介绍1.1字符型1.2整型1.3浮点型1.4 布尔类型1.5 各种数据类型的长度1.5.1 sizeof 操作符1.5.2 数据类型长度1.5.3 sizeof 中表达式不计算 2.signed 和 unsigned3.数据类型的取值范围4. 变量4.1 变量的创建4.2 变量的分类 5…

Leetcode日记 2583. 二叉树中的第 K 大层和

Leetcode日记 2583. 二叉树中的第 K 大层和 题目:解题思路:代码实现制作不易,感谢三连,谢谢啦 题目: 给你一棵二叉树的根节点 root 和一个正整数 k 。 树中的 层和 是指 同一层 上节点值的总和。 返回树中第 k 大的层和…