python：GEDI 波形数据提取

已于 2023-08-25 19:29:25 修改
阅读量796 收藏 3
点赞数
分类专栏： GEE生态遥感 文章标签： python 开发语言
于 2023-07-22 14:41:56 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35591253/article/details/131867400
版权
GEE生态遥感 专栏收录该内容
22 篇文章 58 订阅 ¥99.90 ¥299.90
订阅专栏
这篇博客介绍了如何使用Python处理和可视化NASA的GEDI激光雷达数据，包括配置环境、读取数据、提取波形数据并创建可视化。文中提供了代码详细解释及完整代码示例。
摘要由CSDN通过智能技术生成

作者：CSDN @ _养乐多_

在这篇博客中，我们将介绍如何使用 Python 处理和可视化 GEDI（Global Ecosystem Dynamics Investigation）激光雷达数据。GEDI 是 NASA（美国国家航空航天局）推出的激光雷达地球观测任务，用于获取全球各地生态系统的三维结构信息。本文将以提取研究区内植被波形数据为例，结果如下图所示，

在这里插入图片描述

文章目录

一、配置环境

1.1 示例文件

本教程使用的 GEDI .h5 文件：2019 年 6 月 19 日的 GEDI L1B 第 2 版观测（轨道 02932，子轨道02）。
使用下面的链接直接从LPDAAC数据库下载文件，但是需要NASA Earthdata 登录帐户才能下载本教程中使用的数据。
可以通过下面的链接创建一个帐户并下载：

了解本专栏 订阅专栏 解锁全文
_养乐多_
关注
专栏目录
订阅专栏
python：GEDI 某数据点波形数据提取并保存至csv文件
养乐多的博客
07-22 886
养乐多_本文将介绍如何提取 GEDI（Global Ecosystem Dynamics Investigation）激光雷达数据中所标识激光测高数据点的波形数据，并将波形数据提取至csv文件，为了后续数据去噪、滤波和树高特征提取做准备。
python获取声音波形的关键特征
LDC，公众号【轻松学编程】
11-12 867
python获取声音波形特征 时域分析：均值,方差,标准差,最大值,均方根,偏度,峭度,波形因子, 峰值因子,脉冲因子,裕度因子 频域分析：过零率、频谱中心、频谱滚降点、MFCC (梅尔频率倒谱系数)
python 波形包络线_一个简单的波形包络提取算法
weixin_39521068的博客
12-10 2602
最近做了一个项目来提取语音信号包络波形，因此，花一些时间来研究各种包络提取算法。所谓包络检测被称为振幅解调，在许多领域都有重要的应用。它假设载波信号被确定，所以，通常的方式同步解调，优选以这样的方式的SNR，混合信号噪声抑制最强。的叫法。在信号检測领域，这样的方式通常称为“相敏检波”。锁相放大器(Lock-inAmplifier)就是这样的方式最典型的样例。假设载波比較乱。就像我如今的应用场景。...
Python 波形处理
Why Do You Run
10-08 7034
使用Jupyter Notebook进行波形分析 首先拿到了振动传感器的数据 解释数据： 我们知道机械波分为横波、竖波、以及复杂波，现实世界中全部为复杂波——上下左右都有振动。 传感器每秒记录10240条数据（**每分钟仅记录前30秒的数据**），对应X轴、Y轴、以及Z轴的数据。我们这次仅分析一天 Z轴的数据。 由于一秒10240条数据一天的数据量就是 10240*30s*60min*24h =...
python时域波形特征分析
lizz2276的博客
07-16 7152
（一）python实现信号时域统计特征提取 https://blog.csdn.net/qq_34705900/article/details/88389319 1.实验数据需求 为了对采集的压力实验数据做特征工程，需要对信号进行时域的统计特征提取，包含了均值、均方根、偏度、峭度、波形因子、波峰因子、脉冲因子、峭度因子等，现用python对其进行实现。 2.python实现 其中的输入参数含义： ①data:实验数据的DataFrame ②p1:所截取实验信号的起始采样点位置 ③p2:所截取实验
利用python的强大函数库，实现波形的小波降噪、带通滤波、时阈分析、FFT波形转换
玩火的稻草人的博客
07-21 2064
代码】利用python的强大函数库，实现波形的小波降噪、带通滤波、时阈分析、FFT波形转换。
python：GEDI HDF5 文件导入、读取文件元数据以及按 Beam 读取 SDS 元数据和子集
养乐多的博客
07-21 598
养乐多_本文将介绍GEDI HDF5 文件导入、文件元数据读取以及按 Beam 读取 SDS 元数据和子集的python代码。以及读取和处理GEDI数据过程中遇到的环境配置问题解决方法（比如cartopy安装失败）。
GEE：GEDI数据提取值到矢量区域和点
养乐多的博客
05-23 1305
本文将介绍GEDI数据集从GEE上下载到本地，并将每一个激光点的值提取为一个矢量区域，并提取值到矢量区域的方法。
pyGEDI GEDI数据的读取分析（Python版本）
dayuhaitang1的博客
03-30 949
pyGEDI GEDI数据的读取分析（Python版本）
wireshark捕获不到东西_三种方法教你用示波器快速捕获异常
weixin_33967069的博客
12-13 786
万事开头难！当你想用示波器来分析问题时，你一定有想过，我要如何才能把问题抓下来？当然，只有抓下来之后，才能进行后面种种的分析，否则一切都是空谈。本文将带你用三种最好用的方法将异常抓下来。一、滚动模式滚动模式也许你很少用，但它却是分析问题最简单、最粗暴的方法。你仅仅要做的，就是确定异常多长时间会出现，采样率是否足够。如5秒内会出现的异常，设置滚动采集7s的数据后停止，在采样率足够的前提下，我相信问题...
GEDI学习1：H5文件所需各参数及波形的提取
Moxiaix的博客
08-01 928
提取GEDIL1B的波形和参数
《Python数字信号处理应用》学习笔记——第二章 谐波
最新发布
lstef的博客
05-21 1162
之后会了解完整的FFT，它能处理复数信号（见7.9节），将rfft的结果赋值给hs，这是一个NumPy复数数值，表示的是波形中各个频率元素的振幅和相位差。和三角波一样，方波也只含有奇数谐波，这也就是其频谱中在300Hz、500Hz和700Hz这样的值处有尖峰的原因。同样 7700Hz的信号也无法与2300Hz的信号区分，9900Hz的信号无法与100Hz的区分。以5000Hz折叠，就会出现在-2100Hz，但是它在0Hz处被再次折叠，从而回到2100Hz。2. frac序列介于0~1，其频率是给定的。
Python音频处理（一）——信号，波形与频谱
热门推荐
李元静
03-06 3万+
本文目录前言波形频谱Spectrum波形对象信号对象 前言 音频处理属于大学课程《多媒体技术》，什么采样率，频谱等理论知识，小编这里会慢慢的根据实际的代码进行讲解，不会一笼统的纯理论知识一大堆摆上来，毕竟学习音频处理是为了我们处理日常生活中的工作。 关于音频处理的方式，后面的代码都会使用python语言进行处理，至于学完后可以干啥，肯定是大家非常关系的，比如音乐的降噪，模拟某个人的声音，提取音频中...
小波异常值提取代码
qq_53011270的博客
10-13 315
小波异常值提取代码
GEDI学习笔记1：数据产品简介
晚秋的专栏
08-17 8799
GEDI 科学数据产品包括描述地球 3D 特征的足迹和网格数据集。这些数据产品被分配了不同的级别，这表明了数据在收集后所经历的处理量。所有产品都是公开可用的，其中较低级别的产品（L1 和 L2）来自 NASA 的土地过程分布式主动档案中心(LPDAC)，较高级别的产品 (L3 和 L4) 来自 ORNL DAAC。数据最初传输到戈达德太空飞行中心的GEDI任务运营中心 (MOC)，该中心每周部署采购计划，然后通过科学运营中心 (SOC) 进行处理，将科学数据产品分发给上述 DAAC。 用于创建这...
15种时域特征提取的Python代码实现（平均值、方差、标准差、偏度、峭度、最大/小值、峰峰值、均方根、峰值因子、波形因子、脉冲因子、裕度因子、峭度因子、能量）
qq_41361542的博客
03-01 4157
15种时域特征提取的Python代码实现（平均值、方差、标准差、偏度、峭度、最大/小值、峰峰值、均方根、峰值因子、波形因子、脉冲因子、裕度因子、峭度因子、能量）
python剔除异常值的方法_剔除每个cycle异常值（标准差）
weixin_39725154的博客
11-20 888
波形重定后剔除每个cycle内的异常值数据格式：Cycle_name Time Altimeter_data4.0 1013914019.8631 4581.9012074.0 1013914019.9141 4582.1156234.0 1013914019.9651 4581.9575314.0 1013914020.0161 4581.8448664.0 ...
获取波形极值与间距并显示
LDC，公众号【轻松学编程】
05-07 54
1、通过signal.find_peaks获取极大值 2、获取极大值下标 3、获取极大值对应的值 4、获取极大值的下标间距（就是隔多远有一个极大值） 5、获取极大值间距的标准差、方差、均值、最大值 6、图形展示波形图并标记极大值
python实现GEDI1B数据校正
05-04
GEDI1B是GEDI (Global Ecosystem Dynamics Investigation)卫星激光测高仪的数据格式之一。GEDI1B数据是由GEDI1A数据处理得到的，其中包含了校正后的地面高度信息。 下面是一个简单的python实现GEDI1B数据校正的示例： ```python import h5py import numpy as np # 打开GEDI1B数据文件 with h5py.File('GEDI1B_sample.h5', 'r') as f: # 获取数据集 elevation = f['/BEAM0000/elevation'] latitude = f['/BEAM0000/latitude'] longitude = f['/BEAM0000/longitude'] shot_number = f['/BEAM0000/shot_number'] # 获取校正系数 bias_correction = f['/BEAM0000/bias_correction_factor'] tx_pulse = f['/BEAM0000/tx_pulse_length'] rx_gain = f['/BEAM0000/rx_gain'] # 计算校正后的高度 corrected_elevation = elevation[:] * bias_correction[:] * tx_pulse[:] * rx_gain[:] # 保存校正后的数据为新的h5文件 with h5py.File('GEDI1B_corrected.h5', 'w') as f_out: # 创建数据集 dset_elevation = f_out.create_dataset('/BEAM0000/elevation', data=corrected_elevation) dset_latitude = f_out.create_dataset('/BEAM0000/latitude', data=latitude[:]) dset_longitude = f_out.create_dataset('/BEAM0000/longitude', data=longitude[:]) dset_shot_number = f_out.create_dataset('/BEAM0000/shot_number', data=shot_number[:]) ``` 这个示例代码打开一个GEDI1B数据文件，获取数据集和校正系数，并计算出校正后的高度数据。然后将校正后的数据保存为新的h5文件。 需要注意的是，这只是一个简单的示例，实际的GEDI1B数据校正需要考虑更多的因素，比如大气校正、地形校正等。实际的校正过程需要根据具体的数据和需求进行调整。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

_养乐多_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值