数据
在 JAXA 的 P-Tree System 注册一个账号,注册后你的邮箱会收到一封有 FTP 账号密码的邮件,登录 FTP 即可检索下载 Himawari-8 与 Himawari-9 的数据了。 在 FTP 中,进入 /jma/netcdf 路径下,就是可以直接使用的 L3 级数据。
文件命名规则为: NC_H{nn}_{YYYYMMDD_hhmm}_R/r{bb}_FLDK.{xxxxx}_{yyyyy}.nc
nn: Himawari 系列卫星编号,08为 Himawari 8,09为 Himawari 9YYYYMMDD_hhmm: 观测的时间bb: 数据内频段的数量xxxxx: 扫描像素数(横向分辨率),2801 为 5km 分辨率,7001 为 2km 分辨率。yyyyy: 扫描线数(纵向分辨率),2401 为 5km 分辨率,6001 为 2km分辨率。
bb前面的 R/r表示区域,R为全圆盘,r为日本区域。
绘图
读取 albedo_01、albedo_02、albedo_03三个变量,构成 RGB 数组,再调整以下即可得到真彩色图。
python
import numpy as npimport matplotlib.pyplot as pltimport xarray as xrFILE = "NC_H09_20260604_0500_R21_FLDK.02801_02401.nc"ds = xr.open_dataset(FILE)b = ds['albedo_01'].valuesg = ds['albedo_02'].valuesr = ds['albedo_03'].values# to rgb imagergb = np.stack([r, g, b], axis=-1)# normalizergb = (rgb - np.nanmin(rgb)) / (np.nanmax(rgb) - np.nanmin(rgb))# gamma correctiongamma = 2.2rgb = np.power(rgb, 1/gamma)ratio = 28/24fig = plt.figure(figsize=(10, 10/ratio))plt.imshow(rgb)plt.axis('off')plt.tight_layout(pad=0)plt.show()
7000*6000 分辨率圆盘图局部(1500*1500)