【WRF数据准备】基于GEE下载静态地理数据-叶面积指数LAI及绿色植被率Fpar

【WRF数据准备】基于GEE下载静态地理数据

  • 准备:WRF所需静态地理数据(Static geographical data)
    • 数据范围说明
  • 基于GEE下载叶面积指数及绿色植被率
    • GEE数据集介绍
    • 数据下载:LAI(叶面积指数)和Fpar(绿色植被率)
    • 数据处理:基于Python处理为单波段LAI数据
  • 参考

GEE的介绍可参见另一博客 【GEE学习第一期】GEE介绍、注册及基本使用,本博客主要介绍基于GEE平台下载所需静态地理数据。

准备:WRF所需静态地理数据(Static geographical data)

主要的所需静态地理数据如下:

  • Albedo
  • 叶面积指数LAI: Leaf area index
  • 光合有效辐射分数Fpar: Green vegetation fraction
  • 不透水面积ISA: Impervious Surface Area
  • 土地利用Land use

GEE的相关静态地理数据介绍可参见另一博客-【GEE数据库】WRF常用数据集总结。

数据范围说明

在确定研究区域范围(此博客以粤港澳大湾区为例)后,可下载处理所需的静态地理数据。
在这里插入图片描述
经纬度范围:

  • 网格1:经度(106 120.58) 纬度(16.626 29.748)
  • 网格2:经度(108.43 118.15) 纬度(18.813 27.561)
  • 网格3:经度(111.346 115.477) 纬度(21.486 24.402)

更新,采用Lambert投影确定的研究区范围如下:
在这里插入图片描述
网格经纬度范围:上:33 下:14 左:102 右:123,下载数据年份选择2020年。

基于GEE下载叶面积指数及绿色植被率

叶面积指数LAI(leaf area index) 为单位土地面积绿色叶片的单面面积总和,即:叶面积指数=叶片总面积/土地面积,它与植被的密度、结构(单层或复层)、树木的生物学特性(分枝角、叶着生角、耐荫性等)和环境条件(光照、水分、土壤营养状况)有关,是表示植被利用光能状况和冠层结构的一个综合指标。

叶面积指数的测定有直接方法、间接方法。由于直接方法具有一定的破坏性,耗时耗力,并且无法反映大面积、大范围内的植被LAI分布情况,间接测定方法尤其是其中的遥感法,以其具有的快速、实时、大面积、不受植被类型影响等优势,成为应用最为广泛的LAI监测手段。

随着遥感技术的飞速发展,越来越多高时间分辨率、高空间分辨率、多光谱、高光谱遥感卫星数据实现了LAI监测。

GEE数据集介绍

Google Earth Engine ——MYD15A2H V6 MODIS综合叶面积指数(LAI)和光合有效辐射分数(FPAR)产品是一个500米分辨率的8天综合数据集。

  • 叶面积指数(LAI)
  • 光合有效辐射吸收比率(fraction of absorbed photosynthetically active radiation, FAPAR/FPAR)是描述太阳光在冠层辐射传输过程中植被吸收比例的参量, 一般定义为植被对波长在400 nm—700 nm 间太阳辐射能量的吸收比率。

数据介绍官网-MYD15A2H.061: Aqua Leaf Area Index/FPAR 8-Day Global 500m
在这里插入图片描述
GEE中代码:

ee.ImageCollection('MODIS/061/MYD15A2H')

数据下载:LAI(叶面积指数)和Fpar(绿色植被率)

1、选取空间范围

//var geometry = ee.FeatureCollection('文件地址/Boundary');
var geometry = ee.Geometry.Rectangle([-102,34,-84,48],'EPSG:4326',false);
//var geometry = ee.Geometry.Rectangle([-94.9910,40.6107,-92.3703,42.57729],'EPSG:4326',false)

说明:EPSG:4326是WGS 84坐标系统,也称为全球地球坐标系统(GCS),是一种广泛使用的地理坐标系统。 这个坐标系统使用经度(longitude)和纬度(latitude)来定义地球上的位置,其中经度的范围是-180到+180,纬度的范围是-90到+90。

2、选取数据集来源:类型、时间范围等

  • ee.ImageCollection(‘MODIS/006/MCD15A3H’):加载 MODIS 版本 006 的 MCD15A3H 数据集,这个数据集提供了“叶面积指数”(LAI)的信息。
  • .select(‘Lai’):从数据集中选择“Lai”波段,提取叶面积指数。
  • .filterDate():过滤图像集合,只保留在 2012 年 1 月 1 日到 2012 年 12 月 31 日之间的图像。
var LAI = ee.ImageCollection('MODIS/006/MCD15A3H').select('Lai').filterDate(ee.Date('2012-01-01'), ee.Date('2012-12-31'))var Fpar = ee.ImageCollection('MODIS/006/MCD15A3H').select('Fpar').filterDate(ee.Date('2012-01-01'), ee.Date('2012-12-31'))

3、创建一个每月的总数集合
这段代码用于处理 MODIS 数据集中的 LAI 和 Fpar,通过按月计算每个月的平均值

  • ee.List.sequence(1, 12):生成一个包含从 1 到 12 的整数列表,表示一年中的每个月。
var months = ee.List.sequence(1, 12);
//print(months);
//Output: [1,2,3,4,5,6,7,8,9,10,11,12]
  • months.map(function(m) {…}):对每个月(m)进行迭代。
  • LAI.filter(ee.Filter.calendarRange(m, m, ‘month’)):从 LAI 图像集合中过滤出当前月份的图像。
  • .mean():计算该月份图像的平均值。
  • .set(‘month’, m):为计算出的平均图像设置一个属性,标记其对应的月份。
  • ee.ImageCollection.fromImages(…):将所有月份的平均图像组合成一个新的图像集合。
var byMonth_LAI = ee.ImageCollection.fromImages(months.map(function(m) {return LAI.filter(ee.Filter.calendarRange(m, m, 'month')).mean().set('month', m);
}).flatten());

这部分代码与计算 LAI 平均值的部分类似,只是它处理的是 Fpar 图像集合
通过相同的逻辑,过滤出每个月的图像并计算平均值,最终生成一个包含每个月 Fpar 平均值的图像集合。

var byMonth_Fpar = ee.ImageCollection.fromImages(months.map(function(m) {return Fpar.filter

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

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

VantUI

官网:Vant 4 - A lightweight, customizable Vue UI library for mobile web apps. Vant组件库: 基础组件 按钮、图标、布局、提示信息等 表单组件 日历、复选框、时间选择、输入框、评分等 反馈组件 弹出框、加载、下拉菜单、消息提示、下拉刷新、滚动…

面试阿里、字节全都一面挂,被面试官说我的水平还不如应届生

测试员可以先在大厂镀金,以后去中小厂毫无压力,基本不会被卡,事实果真如此吗?但是在我身上却是给了我很大一巴掌... 所谓大厂镀金只是不卡简历而已,如果面试答得稀烂,人家根本不会要你。况且要不是大厂出来…

C#入坑JAVA MyBatis入门 CURD 批量 联表分页查询

本文,分享 MyBatis 各种常用操作,不限于链表查询、分页查询等等。 1. 分页查询 在 下文的 的「3.4 selectPage」小节,我们使用 MyBatis Plus 实现了分页查询。除了这种方式,我们也可以使用 XML 实现分页查询。 这里&#xff0c…

1-petalinux2018.3 摸索记录 -petalinux-config

一、petalinux-config的具体配置-ZYNQMP Configuration 1、Linux Compoment Selection Linux Compoment Selection,Linux组件选择. First Stage Bootloader和Auto update ps_init勾选会自动生成fsbl.elf,自动更新ps_init。 PMU Firmware平台管理单元固…

熵与信息论

经典信息论的核心概念是香农熵。假设我们得到了一个变量X的值,X的香农熵量化了我们在获悉 X的值时所能得到的平均信息量;另一种观点是将X的看作在我们获悉的值前对其不确定程度的度量。这两种观点是互补的;我们既可以将看作在我们获悉X的值前…

Ubuntu 22.04系统启动时自动运行ROS2节点

在 Ubuntu 启动时自动运行 ROS2 节点的方法 环境:Ubuntu 系统,ROS2 Humble,使用系统自带的 启动应用程序 目标:在系统启动时自动运行指定的 ROS2 节点 效果展示 系统启动后,自动运行小乌龟节点和键盘控制节点。 实践…

龙蟠科技业绩压力显著:资产负债率持续攀升,产能利用率也不乐观

《港湾商业观察》施子夫 黄懿 去年十月至今两度递表后,10月17日,江苏龙蟠科技股份有限公司(以下简称,龙蟠科技;603906.SH,02465.HK)通过港交所主板上市聆讯。 很快,龙蟠科技发布公告称,公司全…

OceanBase 安全体系解析之身份鉴别

本文作者:金长龙爱可生测试工程师,负责 DMP 产品的测试工作。 本文以MySQL为参照,详细阐述了OceanBase 在MySQL模式下的安全体系中,身份鉴别的能力,涵盖了身份鉴别机制、用户名的构成规则、密码的复杂度,以…

在Java中的动态绑定和静态绑定

动态绑定和静态绑定是两种方法调用的绑定机制静态绑定 静态绑定也称为早期绑定,是在编译时确定调用的方法。动态绑定 动态绑定也称为晚期绑定,是在运行时确定调用的方法。静态绑定用于编译时确定的方法调用,动态绑定是Java实现运行时多态的…

CISE|暴雨受邀出席第二十六届中国国际软件博览会

10月24日至26日,备受瞩目的第二十六届中国国际软件博览会(简称CISE)在国家会展中心(天津)圆满举办。CISE不仅汇聚了来自全国各地的顶尖软件企业和机构,还吸引了众多专家学者和行业精英共襄盛举,…

Cesium基础-(Entity)-(Box)

** 里边包含Vue、React框架代码详细步骤、以及代码详细解释 ** 3、Box 盒子 以下是 BoxGeometry 类的属性、方法和静态方法,以表格形式展示: 属性 属性名类型默认值描述minimumCartesian3盒子的最小 x, y, 和 z 坐标。maximumCartesian3盒子的最大 x, y, 和 z 坐标。vertex…

【PHP】PHP使用Modbus-Rut协议与RS485串口通信,向设备发送和接收数据

目录 一、前言 二、开发前说明 三、效果图 四、安装PHP扩展 五、安装phpModbus类库 六、通信逻辑 七、完整实例 一、前言 使用PHP语言与硬件设备通信交互,并向COM串口发送和接收数据。 前面写了三篇关于PHP与RS235和USB端口通信的文章,可以作为参…

现代数字信号处理I--最佳线性无偏估计 BLUE 学习笔记

目录 1. 最佳线性无偏估计的由来 2. 简单线性模型下一维参数的BLUE 3. 一般线性模型下一维参数的BLUE 4. 一般线性模型下多维参数的BLUE 4.1 以一维情况说明Rao论文中的结论 4.2 矢量参数是MVUE的本质是矢量参数中的每个一维参数都是MVUE 4.3 一般线性模型多维参数BLUE的…

视频剪辑哪个软件好用?推荐四款热门工具!!

在这个Vlog和短视频当道的互联网时代,掌握一款好用的视频剪辑软件就像拥有了打开创作世界的魔法钥匙。今天我们来聊聊视频剪辑软件,帮你成为剪辑达人哦!接下来,给大家详细介绍四款常用且各具特色的视频剪辑软件,助你轻…

算法:利用前序序列和中序序列构造二叉树

题目 链接:leetcode链接 思路分析 前序遍历的顺序是:根 左子树 右子树 中序遍历的顺序是: 左子树 根 右子树 所以,我们可以通过前序遍历获得二叉树的根 可以通过中序遍历去分割二叉树,将二叉树分割成 左子树 根…

偷懒总结篇|贪心算法|动态规划|单调栈|图论

由于这周来不及了,先过一遍后面的思路,具体实现等下周再开始详细写。 贪心算法 这个图非常好 122.买卖股票的最佳时机 II(妙,拆分利润) 把利润分解为每天为单位的维度,需要收集每天的正利润就可以,收集正利润的区间…

HarmonyOS ArkTS与C++数据类型转换

1. HarmonyOS ArkTS与C数据类型转换 本文介绍了C与TS各自数据类型与互相之间的数据类型转换,在需要使用C模块时可以快速上手对各种数据类型进行转换。 1.1. 概述 HarmonyOS的主力开发语言是ArkTS,也提供了C语言的支持,对于一些能力&#xff…

1.3 面向对象 C++面试问题

1.3.1 简述一下什么是面向对象,面向对象与面向过程的区别 什么是面向对象 面向对象(Object-Oriented Programming,OOP)是一种编程范式,它通过将现实世界中的实体抽象为“对象”来组织代码。面向对象编程关注对象及其交互&#x…

D51【python 接口自动化学习】- python基础之模块与标准库

day51 模块的导入 学习日期:20241027 学习目标:模块与标准库 -- 66 模块的导入:如何使用其他人编写好的代码功能? 学习笔记 模块的作用 导入模块的方法 # 导入模块 # 方式一 import os # 获取当前的位置 print(os.getcwd())# …

arduino uno R3更换328pb-au芯片,烧录bootloader

使用usbasp烧录器进行烧录,解压 【免费】usbsap驱动以及软件资源-CSDN文库 安装驱动 然后打开软件 界面如下 1按步骤选中芯片, ATmega328P(由于没有328PB,直接选这个也行) 2查看spi接线, 3读取芯片id&a…