风云4A,HDF文件里面14个通道里面的值代表什么，如何转为云顶亮温？

我爱气象111 · 发表于 2021-7-28 19:38:50

本人根据网上代码画了一张13通道的云图，但是并不知道怎么把值转为云顶亮温。请求各位大神教教我

edwardli · 发表于 2021-7-28 19:38:51

本帖最后由 edwardli 于 2021-7-28 22:05 编辑

      卫星观测到的是辐射能。
      卫星文件里存储的是一个“索引”（可以这么理解，因为反照率是0-100，占6位，亮温是300左右，需要占11位，两相结合），需要根据这个“索引”去HDF、AWX前部的“定标表”查找到底对应哪个具体的值（反照率92%，TBB 270K……之流）。这个过程叫“定标”calibration
      还有个问题是“定位”，尤其是HDF多数都是圆盘图DISK，即使是REGC中国区域，也有两个角扫描的是太空。而且本质上数据是按照行列存储的，比如4km分辨率的全圆盘，就是2748行*2748列。所以需要把这个数学的矩阵/数组，下标对应到地理上（经纬度），才能画出一张合格的卫星云图，进一步与其他气象数据叠加显示分析。
FY4A数据可视化流程.png

话不多说，代码来了。
Main.py

from fy4a import AGRI_L1
file_path = 'data/Satellite/FY4A-_AGRI--_N_DISK_1047E_L1-_FDI-_MULT_NOM_20210720080000_20210720081459_4000M_V0001.HDF'
sat = AGRI_L1(file_path)
C012 = sat.extract('Channel12',calibration='brightness_temperature') #1-6通道可见光的cali为'reflectance'
geo_desc = [20, 60, 100, 140, 0.01] # 顺序为南、北、西、东、分辨率，即 lat_s, lat_n, lon_w, lon_e, resolution单位°
C012 = sat.extract('Channel12',calibration='brightness_temperature',geo_desc = geo_desc)
# C012
# C012.shape
import matplotlib.pyplot as plt
fig = plt.figure(figsize=(16,9),dpi=200)
plt.imshow(C012,cmap='gray_r') #imshow是我喜欢的“概览”方法，改成其他绘图区ax的方法比如pcolormesh……即可，当然推荐cartopy声明带投影的ax。这，不是本回帖之重点。

复制代码

以下需要的类，不是我写的，by [Modabao](https://github.com/Mo-Dabao/BiteFengyun) 华风集团数据技术开发中心

fy4a.py (4.64 KB, 下载次数: 119)

Lucian · 发表于 2021-7-28 20:11:22

由卫星观测的radiance计算亮温
http://bbs.06climate.com/forum.php?mod=viewthread&tid=67798

我爱气象111 · 发表于 2021-7-29 09:23:00

edwardli 发表于 2021-7-28 19:38
卫星观测到的是辐射能。
卫星文件里存储的是一个“索引”（可以这么理解，因为反照率是0- ...

谢谢大神。

溫不暖 · 发表于 2022-4-16 11:58:58

收藏研究！赞

menseye · 发表于 2022-4-29 11:07:07

您好，能分享一下您的代码，我最近也在处理风云4的数据，我画出来圆盘图，但我一直不会投影到平面上并选取特定经纬度范围的图。

dianzimi · 发表于 2022-8-3 12:42:56

非常感谢大神

limeiman · 发表于 2022-8-11 15:20:41

{:eb502:}{:eb502:}

likailing2008 · 发表于 2023-3-4 14:50:16

非常感谢这位经常出现的大神

mjn · 发表于 2023-9-3 09:45:36

学习到了，很

感谢

		自动登录	找回密码
密码			立即注册

[求助] 风云4A,HDF文件里面14个通道里面的值代表什么，如何转为云顶亮温？

最佳答案