求一份python绘制风云4b L1级数据真彩图demo

chhyu

最近领导要求绘制卫星云图，我用的是python，将14个通道数据单独取出来并绘制出来了，单个通道都是黑白的，想绘制官方那种全圆盘真彩图，一直感觉差点意思，求大佬指点，如果有demo那更好了

tqm

https://www.moexin.cn/posts/use-satpy
https://www.moexin.cn/posts/process-fy4b

chhyu

tqm 发表于 2024-11-1 19:40
https://www.moexin.cn/posts/use-satpy
https://www.moexin.cn/posts/process-fy4b

非常感谢大佬，这就去学习一下

tqm

因为还要计算绿色波段，最好直接用satpy读出来再丢进matplotlib操作了

1. import os
   
2. from glob import glob
   
3. from satpy.scene import Scene
   
4. filenames = glob('FY4B-_AGRI--_N_DISK_*')
   
5. scn = Scene(filenames, reader='agri_fy4b_l1')
   
6. composite = "true_color"
   
7. scn.load([composite])
   
8. scn = scn.resample(scn.finest_area(), resampler='native')
   
9. scn.save_dataset(composite, filename='FY4B.png')

复制代码

手动合成RGB大概可以stack计算好后的3个波段做RGB然后imshow，或者照抄satpy的GenericCompositor的RGB合成方法

chhyu

tqm 发表于 2024-11-6 21:14
因为还要计算绿色波段，最好直接用satpy读出来再丢进matplotlib操作了手动合成RGB大概可以stack计算好后的3 ...

大佬，有没有demo的，我对这里面的计算一无所知