satpy 处理卫星 FY4A 数据
读取数据并画图
import os
import glob
from datetime import datetime, timedelta
from satpy.scene import Scene
from pyresample import get_area_def
import warnings warnings.filterwarnings('ignore') def fileNameSplit(fileNs):
fSplit = []
for fpn in fileNs:
fp = os.path.dirname(fpn)
fn = os.path.basename(fpn).split('_P_')[-1]
nfpn = os.path.join(fp, fn)
if not os.path.exists(nfpn):
os.rename(fpn, nfpn)
fSplit.append(nfpn)
return fSplit def time_format(info, format="%Y%m%d%H%M%S", format2="%Y%m%d%H%M00"):
dt = datetime.strptime(info, format)
strdt = dt + timedelta(hours=8)
return datetime.strftime(strdt, format2) def get_fy4a(file_paths, channel_list):
fileTime = time_format(os.path.basename(file_paths[0]).split('_')[-3])
# 创建scene对象
scn = Scene(file_paths, reader='agri_l1')
# 查看可使用的通道
if 'C03' in channel_list:
print(scn.available_dataset_names())
if 'true_color' in channel_list:
# 查看可用的合成选项
print(scn.available_composite_names())
# 读取指定通道数据
scn.load(channel_list)
# 打印加载通道后 xarray.Dataset 的数据集
# print(scn) """
>python coord2area_def.py china merc 5 54 60 139 3
### +proj=merc +lat_0=29.5 +lon_0=99.5 +ellps=WGS84 china:
description: china
projection:
proj: merc
ellps: WGS84
lat_0: 29.5
lon_0: 99.5
shape:
height: 2194
width: 2931
area_extent:
lower_left_xy: [-4397119.886334, 553583.846816]
upper_right_xy: [4397119.886334, 7135562.567523]
"""
# 画自定义区域图片
area_id = 'china'
area_name = "myarea"
proj_id = 'china_99.5_29.5'
proj4_args = '+proj=merc +lat_0=29.5 +lon_0=99.5 +ellps=WGS84'
height = 2194
width = 2931
area_extent = [-4397119.886334, 553583.846816, 4397119.886334, 7135562.567523]
# 定义地图投影和区域,然后将数据投影到该区域上
areadef = get_area_def(area_id, area_name, proj_id, proj4_args, width, height, area_extent)
china_scene = scn.resample(areadef) channel_ele = {'C03': 'vis', 'C09': 'vap', 'C12': 'ir', 'true_color': 'thr'}
for cele in channel_list:
pic_dir = r'./'
pic_name = f'fy4a_{channel_ele[cele]}_{fileTime}.png'
pic_path = os.path.join(pic_dir, pic_name)
if not os.path.exists(pic_path):
print(f'生成图片 {pic_name}')
# 保存图片
china_scene.save_dataset(cele, filename=pic_path) if __name__ == '__main__':
t1 = datetime.now() grayFilePathList = ['./FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20211020022336_20211020022753_4000M_V0001.HDF',
'./FY4A-_AGRI--_N_REGC_1047E_L1-_GEO-_MULT_NOM_20211020022336_20211020022753_4000M_V0001.HDF']
print(grayFilePathList)
channelList = ['C03', 'C09', 'C12']
get_fy4a(grayFilePathList[:1], channelList)
t2 = datetime.now()
print(t2 - t1)
print('------')
thrFilePathList = ['./FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20211020022336_20211020022753_4000M_V0001.HDF',
'./FY4A-_AGRI--_N_REGC_1047E_L1-_GEO-_MULT_NOM_20211020022336_20211020022753_4000M_V0001.HDF'] print(thrFilePathList)
channelList1 = ['true_color']
get_fy4a(thrFilePathList, channelList1)
t3 = datetime.now()
print(t3 - t2)reader='agri_l1'
版本<=0.36.0 'agri_l1' >=0.37.0 'agri_fy4a_l1' 其他请看网址
https://satpy.readthedocs.io/en/stable/index.html#reader-table
HTTPSConnectionPool(host='zenodo.org', port=443): Max retries exceeded with url: /record/1288441/files/pyspectral_atm_correction_luts_no_aerosol.tgz
# 原因下载依赖文件失败
# 手动下载地址https://zenodo.org/record/3824535/files/pyspectral_rsr_data.tgz
https://zenodo.org/record/1288441/files/pyspectral_atm_correction_luts_no_aerosol.tgzwin
pyspectral_rsr_data.tgz 解压放到下面地址
C:\Users\admin\AppData\Local\pytroll\pyspectral pyspectral_atm_correction_luts_no_aerosol.tgz 解压放到下面地址
C:\Users\admin\AppData\Local\pytroll\pyspectral\rayleigh_only linux
pyspectral_rsr_data.tgz 解压放到下面地址
/root/.local/share/pyspectral pyspectral_atm_correction_luts_no_aerosol.tgz 解压放到下面地址
/root/.local/share/pyspectral/rayleigh_onlyDon't know how to open the following files: {'./FY4A.......HDF'}
# 按我猜测,应该是satpy 和其依赖包版本的问题
h5py==3.7.0
numpy==1.21.6
satpy==0.36.0
pyresample==1.22.0
scikit-image==0.19.3
pyspectral==0.12.0参考地址
官网:https://satpy.readthedocs.io/en/stable/readers.html#available-readers
博客:https://cloud.tencent.com/developer/article/1584969牛人:https://github.com/aiwei169/FY4A/blob/main/codes/fy4a.py
satpy 处理卫星 FY4A 数据的更多相关文章
- MeteoInfo-Java解析与绘图教程(八)_java解析卫星FY-4A一级产品文件(HDF举例)
MeteoInfo-Java解析与绘图教程(八)_java解析卫星一级产品文件(HDF举例) 最近解析卫星数据遇到了一级产品,它的解析方式与之前文章说的有些不同,特此补充一下 卫星的一级产品,里面是没 ...
- 使用GDAL工具对FY3系列卫星数据进行校正
本文档主要对如何使用GDAL提供的工具对FY3系列卫星数据进行校正处理.FY3系列卫星提供的数据一般是以HDF5格式下发,一个典型的FY3A和FY3B的数据文件名如下: FY3A_MERSI_GBAL ...
- 利用TLE数据确定卫星轨道(1)-卫星轨道和TLE
1.卫星轨道 太空中卫星和天体在各种引力的作用下都在一定的轨道中周期转动着,但实际的轨道是很复杂的,一般的项目也达不到那么精确的需求(其实精确的卫星参数数据也不可能随便公开的),所以采用一阶近似的开普 ...
- PIE SDK打开静止卫星数据
1. 功能简介 静止卫星是位于地球赤道上空约3.58万km处,与地面始终保持相对静止的卫星,静止卫星的特点是覆盖区域广,具有很强的机动灵活性,能够对特定区域进行分钟级高重复观测,可快速监测灾害目标的动 ...
- 使用Geolocation校正GDAL不支持的数据
对于低分数据来说,常用的校正方式就是给定数据的经纬度查找表来进行校正.在GDAL中,这种校正方式叫Geolocation array.常用的数据有国外的MODIS数据,国内的如风云系列(FY)和海洋系 ...
- Cesium专栏-卫星轨迹
Cesium Cesium 是一款面向三维地球和地图的,世界级的JavaScript开源产品.它提供了基于JavaScript语言的开发包,方便用户快速搭建一款零插件的虚拟地球Web应用,并在性能,精 ...
- 全球DEM高程数据下载
本文主要介绍如何使用“迈高图-地图数据下载器”(以下简称:迈高图)下载全球DEM高程数据,任意下载范围自动拼接.裁剪.DEM高程数据支持下载导出为:GeoTIFF.XYZ.南方CASS等常用数据格式. ...
- MeteoInfoLab脚本示例:OMI Swath HDF数据
这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图.脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_L ...
- 蒸腾量与蒸散量(ET)数据、潜在蒸散量、实际蒸散量数据、气温数据、降雨量数据
数据下载链接:数据下载链接 引言 多种卫星遥感数据反演地表蒸腾与蒸散率(ET)产品是地理遥感生态网推出的生态环境类数据产品之一,产品包括2000-2009年三个波段RGB数据,值域0-252之 ...
- 高精度10m/30米NPP净初级生产力分布数据
数据下载链接:百度云下载链接 引言 第一性生产力是绿色植物呼吸后所剩下的单位面积单位时间内所固定的能量或所生产的有机物质,即是总第一性生产量减去植物呼吸作用所剩下的能量或有机物质.多种卫星遥感数据反 ...
随机推荐
- Vue学习:实现用户没有登陆时,访问后自动跳转登录页面
设计思路 定义路由的时候配置属性,这里使用needLogin标记访问页面是否需要登录 设置路由守卫,每个页面在跳转之前都要经过验证,校验用户信息是否存在,不存在跳转到登录页 用户登录后将用户信息存储在 ...
- 【SDOI2015】星际战争
#include<cstdio> #include<queue> using namespace std; const int M = 10000; const double ...
- Jenkins安装(Docker)版
一.jenkins安装 1.查找,下载jenkins镜像文件 启动docker,查找Jenkins镜像文件 docker search jenkins 下载Jenkins镜像文件 docker pul ...
- 【调试】ftrace(一)基本使用方法
简介 Ftrace是Linux Kernel的官方tracing系统,支持Function trace.静态tracepoint.动态Tracepoint的跟踪,还提供各种Tracer,用于统计最大i ...
- Kali配置gmssl密码算法库
Kali配置gmssl密码算法库 一.密码算法库的下载 https://github.com/guanzhi/GmSSL/releases 二.安装配置 1 解压 把刚刚下载的GmSSL 3.0.0. ...
- 原创如何给MDK5.37添加Arm Compiler 5
最新发布的MDK5.37已经不再安装Arm Compiler 5(ARMCC)编译器了,因为点击魔术棒后,在Target选项卡中选择编译器时,会看到missing:compiler version 5 ...
- Solidity8.0-02
对应崔棉大师 26-40课程https://www.bilibili.com/video/BV1yS4y1N7yu/?spm_id_from=333.788&vd_source=c81b130 ...
- 【剑指Offer】【树】二叉树中和为某一值的路径
题目:输入一颗二叉树的根节点和一个整数,打印出二叉树中结点值的和为输入整数的所有路径.路径定义为从树的根结点开始往下一直到叶结点所经过的结点形成一条路径.(注意: 在返回值的list中,数组长度大的数 ...
- 使用clipboard.js复制文字+图片到微信后图片不显示问题处理
使用clipboard.js复制文字 +图片,粘贴到微信不显示图片,而QQ可以. 解决方案:图片链接使用http,不要使用https. 使用clipboard.js实现复制功能 文字 +图片到微信客户 ...
- react修改打包后css,js,src引入的静态地址路径修改
在package.json 里边增加该配置 (当然,配置的地址可以由部署环境决定,可以是 './' , 也可以是 '../',根据需要配置即可) 查找到到"private": tr ...