葵花(Himawari)8/9卫星数据处理

news2024/10/6 4:09:38

李国春

葵花8(含葵花9,下同)是静止气象卫星,数据范围是星下点为140°E赤道上空的大圆盘。空间分辨率有500m、1km和2km三种。虽然其空间分辨率不高,但是有极高的时间分辨率,观测数据刷新仅10分钟。这对一些生态、农业和大范围自然灾害的动态检测具有非常高的应用价值。

本文介绍如何使用RSD处理葵花8数据。

一、数据下载

下载地址(ftp.ptree.jaxa.jp),登陆后进入jma,然后有两个目录,一个是hsd,里面存放的是扩展名为dat的数据。另一个是netcdf,里面存储的是nc格式数据。

1. dat格式数据

进入hsd目录,里面一大堆的压缩文件。某一个文件的文件名可能如下:

HS_H09_20230226_0350_B04_FLDK_R10_S0210.DAT.bz2

文件名中 FLDK ,表示是大圆盘数据,是我们需要的,没有 FLDK 我们暂时不予考虑。

B04 表示该数据集是4波段数据,文件名中的这个字段可以是B01~B16,代表1~16波段。

R10 是分辨率指示。有R05、R10和R20三种,分别代表葵花8的500m、1km和2km分辨率数据。

S02 是大圆盘分割条带的指示。S01~S10表示从大圆盘分割出来的某一条。

可见,下载大圆盘各波段不同分辨率的数据需要下载一大堆的文件。

注意:下载的文件解压后都放在一个目录中,不要每个文件放在单独的一个目录

2. nc 格式的数据

在日本气象局jma目录下进入netcdf目录,直接下载nc格式的数据文件。实际上这些nc格式的数据是 hdf5格式。

二、处理500m的大圆盘数据

1. 打开葵花8 500m 大圆盘数据

启动RSD,点击图1的菜单命令

 图1 打开葵花8大圆盘数据的菜单命令                                                

点击后,从打开文件对话框选择一个文件打开,如:

HS_H09_20230226_0350_B03_FLDK_R05_S0310.DAT

注意这个文件的文件名中一定要包含 R05 ,表示500m数据,另外要包含 FLDK ,表示是大圆盘数据。点击打开后,RSD自动搜索该500m分辨率数据的全部条带和全部波段(500m只有3波段)数据。然后出现投影对话框,直接点击确定。RSD自动投影后,打开结果见图2。

图2 等经纬度网格的葵花8 500m数据

2. 向现有任务添加葵花8 500m数据

如果需要在特定的区域使用特定的投影处理葵花8数据,可以在投影对话框中更改投影和框架范围。也可以在已经处理中的任务添加葵花8数据。例如,中国区域500m的Albers投影框架。从图3的菜单命令选择添加葵花8 500m数据。

图3 向现有任务添加葵花8数据的菜单命令

添加结果见图4。

 图4 中国大陆地区添加葵花8 500m数据

3. 保存投影后的葵花8 500m数据

可以从菜单命令选择保存为RSD格式文件,或者从层列表勾选葵花8数据层,点击 任务管理->导出数据->导出GeoTIFF格式文件 ,或者右击该层,在弹出菜单选择 导出GeoTIFF格式文件 。或者也可以选择导出其它格式文件。

三、处理1km的大圆盘数据

同样使用图1的菜单命令,选择1km数据。数据文件可以选择如:

HS_H09_20230226_0350_B02_FLDK_R10_S0510.DAT

R10表示是1km数据,只选择一个文件,RSD就自动查找所有的1km大圆盘数据文件。

打开后如图5。

 图5 1km葵花8大圆盘数据打开结果

图中红框可见1km葵花8数据有3个波段,分别为1、2、4波段。

其处理它同500m数据。

四、处理2km的大圆盘数据

同样使用图1的菜单命令,选择2km数据。数据文件可以选择如:

HS_H09_20230226_0350_B05_FLDK_R20_S0210.DAT

R20表示是2km数据,只选择一个文件RSD就自动查找所有的大圆盘2km数据文件。

打开后结果图略。

2km葵花8数据有12个波段,为5~16波段。其它处理同500m数据。

五、计算葵花8数据反射率、亮温

葵花8数据在dat文件中已经给出了计算辐亮度、反射率和亮温的系数。也可以在这里直接计算以节省的处理时间。图1中有计算各个分辨率数据反射率、亮温的单独的菜单。虽然提供了不同的菜单,但是与打开通道计数使用的是同样的数据。操作过程也相同。

比如打开2km数据的反射率和亮温,打开一个中国大陆的2km Albers投影的框架。按照图6的菜单命令添加葵花8数据,结果见图7。

 

 图6 添加葵花8数据反射率、亮温的菜单

图7 计算葵花8数据的反射率和亮温

RSD中算得的反射率放大了10000倍,亮温放大了100倍。例如图7中间红框中,第1、2两个通道(波段5、6)是反射波段,其结果是反射率,值2394表示反射率值为0.2394。3~12通道(波段7~16)是发射波段,值29850表示亮温值为298.50°K。

500m数据与1km数据均为反射波段,结果与2km数据类似,恕不赘述。

六、处理NetCDF格式葵花8数据

在jma目录下进入netcdf目录可以下载NetCDF格式数据。有2km和5km两种。使用图1中 Himawari8/9 nc 格式 菜单命令可以直接打开扩展名为nc的数据。打开结果见图8。

图8 打开NetCDF格式的葵花8数据

NetCDF格式的葵花8数据包含全部16个波段的数据。前六个数据通道为反射率波段,后面10个为黑体辐射亮温。

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

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

相关文章

数据结构中“树”的全面讲解

文章目录 一、树结构的定义与对比二、树的存储结构是什么,都有哪些存储结构?双亲表示法:孩子表示法:孩子兄弟表示法: 一、树结构的定义与对比 ​ 树结构是一种一对多的非线性结构,它是由n(n>…

第1关:ODBC程序设计

第1关:ODBC程序设计 任务描述相关知识ODBC主要功能ODBC接口主要函数ODBC应用程序开发实例DM ODBC应用程序开发总体流程DM ODBC代码编写流程DM ODBC代码编写实例 编程要求测试说明代码参考: 任务描述 本关任务:使用 ODBC 查询表中数据。 相关…

Seata分布式事务实战

目录 1 Seata1.1 Seata术语1.2 Seata AT模式1.2.1 AT模式及工作流程1.2.2 Seata-Server安装1.2.3 集成springcloud-alibaba 1.3 Seata TCC模式 2 Seata注册中心2.1 服务端注册中心配置2.2 客户端注册中心配置 3 Seata高可用 1 Seata Seata实现分布式事务基础知识:https://blog.…

“敏捷”又“精益”的新企架思想,如何指导项目实施?

就像“一千个读者眼中有一千个哈姆雷特”,企业也有一千种被看待的方式,例如消费者会用品牌和产品来描述企业;投资者会用盈利模式和盈利空间来评估企业,而“企业架构”则是企业管理者、架构师等用来精确描述企业的方式。 通过业务…

【C语言】语言篇——程序设计入门

C站的小伙伴们大家好呀!我最近在学习刘汝佳老师的《算法竞赛入门经典》,跟着这本书来学习和做习题,在这里和大家一起分享进步。下面是本书的第一部分的语言篇。 语言篇——程序设计入门 算术表达式变量及其输入顺序结构程序设计分支结构程序设…

1.1 什么是eBPF?(上)

大多数介绍eBPF的文章都是用“eBPF是一种革命性的内核技术”来描绘的。这样讲一点也不夸张。因为它允许在Linux的内核中执行沙盒程序,在不改变内核源码或加载内核模块的前提下直接地,安全地,快捷地扩展内核,并改变内核的行为。可以想像在运行时,将用户空间的eBPF程序加载到…

C4D R26 渲染学习笔记 建模篇(2):手动建模

介绍篇 C4D R26 渲染学习笔记(1):C4D版本选择和初始UI框介绍 C4D R26 渲染学习笔记(2):渲染流程介绍 C4D R26 渲染学习笔记(3):物体基本操作快捷键 建模篇 C4D R26 渲…

第一个servlet的程序

文章目录 一.Hello World的程序1.创建项目2.引入依赖3.创建目录4.编写代码5.打包程序6.部署程序7.验证程序 二.简化部署方式1.下载插件2.配置smart Tomcat插件3.测试插件 三.常见的servelt问题出现 404出现 405出现 500出现 "空白页面"出现 "无法访问此网站&quo…

图像风格迁移

一、简介 图像风格迁移是指,将一副内容图的内容,和一幅或多幅风格图的风格融合在一起,从而生成一些有意思的图片。 我们使用 TensorFlow 和 Keras 分别来实现图像风格迁移,主要用到深度学习中的卷积神经网络,即CNN。…

Transformer应用之构建聊天机器人(二)

四、模型训练解析 在PyTorch提供的“Chatbot Tutorial”中,关于训练提到了2个小技巧: 使用”teacher forcing”模式,通过设置参数“teacher_forcing_ratio”来决定是否需要使用当前标签词汇来作为decoder的下一个输入,而不是把d…

Linux:查看进程。

Linux:查看进程。 windows linux TTY如果是?说明是不是终端(控制台)启动的,而是系统内部自己启动的。 TIME是启动Linux后,这个进程一共占用了cpu多少时间00…

《Spring Guides系列学习》guide46 - guide50

要想全面快速学习Spring的内容,最好的方法肯定是先去Spring官网去查阅文档,在Spring官网中找到了适合新手了解的官网Guides,一共68篇,打算全部过一遍,能尽量全面的了解Spring框架的每个特性和功能。 接着上篇看过的gui…

《Python安全攻防:渗透测试实战指南》极致经典,学完即可包吃包住

前言 网络江湖,风起云涌,攻防博弈,从未间断,且愈演愈烈。从架构安全到被动纵深防御,再到主动防御、安全智能,直至进攻反制,皆直指安全的本质——攻防。未知攻,焉知防! 每一位网络安…

【Python】循环语句 ② ( while 嵌套循环 | 代码示例 - while 嵌套循环 )

文章目录 一、while 嵌套循环1、while 嵌套循环语法2、代码示例 - while 嵌套循环 一、while 嵌套循环 1、while 嵌套循环语法 while 嵌套循环 就是 在 外层循环 中 , 嵌套 内层循环 ; while 嵌套循环 语法格式 : while 外层循环条件:外层循环操作1外层循环操作2while 内存循…

VuePress + GitHub Actions 自动部署

文章目录 前言背景GitHub Actions简介基本概念引用 Actionworkflow 文件 自动部署创建 Action权限问题 小结参考文献 前言 我的第二本开源书籍《后台开发命令 365》上线啦。 为了方便阅读,使用 VuePress 将之前记录的后台常用 Linux 命令博文整理成一个系统的开源…

路径规划算法:基于阴阳对优化的路径规划算法- 附代码

路径规划算法:基于阴阳对优化的路径规划算法- 附代码 文章目录 路径规划算法:基于阴阳对优化的路径规划算法- 附代码1.算法原理1.1 环境设定1.2 约束条件1.3 适应度函数 2.算法结果3.MATLAB代码4.参考文献 摘要:本文主要介绍利用智能优化算法…

Compose 没有 inputType 怎么过滤(限制)输入内容?这题我会!

前言 闲话 在我之前的文章 《Compose For Desktop 实践:使用 Compose-jb 做一个时间水印助手》 中,我埋了一个坑,关于在 Compose 中如何过滤 TextField 的输入内容。时隔好几个月了,今天这篇文章就是来填这个坑的。 为什么需要…

Doris

Aggregate 模型 是相同key的数据进行自动聚合的表模型。表中的列按照是否设置了 AggregationType,分为 Key(维度列)和 Value(指标列),没有设置 AggregationType 的称为 Key,设置了 Aggregation…

散列表(哈希表)

目录 散列表 散列函数 散列表常用函数 1. 直接定址法 2. 除留余数法 2.1. exmple 3. 数字分析法 4. 平方取中法 5. 折叠法 处理冲突的方法 1. 开放定址法---线性探测 2. 二次探测法 3. 再Hash法 4. 拉链法(链地址法) 散列表(Hash table,也…

Redis缓存击穿及解决问题

缓存击穿的意思是对于设置了过期时间的key,缓存在某个时间点过期的时候,恰好这时间点对这个 Key有大量的并发请求过来,这些请求发现缓存过期- -般都会从后端DB加载数据并回设到缓存,这个时候大并发的请求可能会瞬间把DB压垮。 解决方案有两种…