海洋气象知识互动

李知恩

目前只是一名研0小白，研究方向：海洋气象
热衷于解答坛友关于WRF数值模式、Python绘图和风云卫星资料方面的问题，如果有以上疑问或者经验，欢迎和我讨论分享！留言也可以！评论也可以
因为见识浅薄，有很多疏忽和不足，欢迎大家添加我的好友批评指正
顺便安利我的摸鱼卫星公主号：木子知恩
谢谢！！（绷不住了，等级不够，还不能自己留言，那就发个帖吧）

李知恩

Python | 基于Linux系统下conda环境的部署及其迁移：https://mp.weixin.qq.com/s/SfYOaafgyUmQvCvytMrIWQ

看云淡风清

楼主厉害，想请教下怎么将FY-4 AGRI全圆盘数据转换成等经纬度格点数据，类似于ERA5数据那样的，主要是想方便后续的数据处理。网上有一些介绍，不过感觉有些复杂，不知道楼主可否介绍一下，谢谢！

李知恩

看云淡风清 发表于 2024-8-14 20:28
楼主厉害，想请教下怎么将FY-4 AGRI全圆盘数据转换成等经纬度格点数据，类似于ERA5数据那样的，主要是想方 ...

其实还是很好实现的，我可以给你我的方法，作为一个例子，你可以参考一下：
filename_cfr_l2 = "F:\computer\data\A202403250370789027\FY4B-_AGRI--_N_DISK_1330E_L2-_CLM-_MULT_NOM_20230123210000_20230123211459_4000M_V0001.NC"

rawfile = r'F:\computer\data\FullMask_Grid_4000.raw'

dim = 2748
data = np.fromfile(rawfile,dtype=float,count=dim*dim*2)
latlon = np.reshape(data,(dim, dim, 2))
lat = latlon[:,:,0]
lon = latlon[:,:,1]
data =  xr.open_dataset(filename_cfr_l2)
lat = np.ma.masked_equal(lat, 999999.9999) #space value
lon = np.ma.masked_equal(lon, 999999.9999) #space value

李知恩

李知恩 发表于 2024-8-14 23:28
其实还是很好实现的，我可以给你我的方法，作为一个例子，你可以参考一下：
filename_cfr_l2 = "F:\comp ...

有几点需要注意：
1、FullMask文件可以在风云卫星的资料下载界面获取；
2、得到的lat和lon就是正常经纬度的格点数据
3、如果要对某个变量读取要参考卫星二级产品的变量说明，用imshow画图或许会更方便；
4、如果后续还有问题可以回复，但是我下周要回老家了，下周不能频繁水贴了，难过，但是可以关注我的公主号呀！

看云淡风清

李知恩 发表于 2024-8-14 23:32
有几点需要注意：
1、FullMask文件可以在风云卫星的资料下载界面获取；
2、得到的lat和lon就是正常经纬 ...

感谢楼主的回复，棒棒哒

李知恩

一周了，终于拿到电脑了，我来水贴了

李知恩

Git学习 | 最先进的分布式版本控制系统 01
：https://mp.weixin.qq.com/s?__biz ... 7448e4702cf0286b#rd

李知恩

李知恩 发表于 2024-8-21 22:54
Git学习 | 最先进的分布式版本控制系统 01
：https://mp.weixin.qq.com/s?__biz=MzkzNzI3MzI3MA==&mid=224 ...

新推文！快关注快关注快关注我一起学习

李知恩

Python | 气象数据后处理和可视化 | 数据插值和平均：https://mp.weixin.qq.com/s?__biz ... 0&lang=zh_CN#rd