3月15日近10年来最强的沙尘暴袭击了北京。关于此次沙尘暴的天气分析已经非常多了,本文不会分析相关的天气背景,主要从技术方面讲一下如何利用satpy处理卫星数据,从卫星视角看一下此次沙尘过程的演变。
本文的数据为 Himawari-8 静止卫星L1b产品:
以下是日期解析和数据处理函数:
因为所处理的数据比较多,上述 process 函数直接将处理后的结果保存为 .tiff 格式文件。如果你想单独看某一个时刻的卫星图像,可以在 save_dataset 之前使用 show 显示图片。
注意:satpy 官方暂时不支持 Himawari-8 卫星的 L1b 数据接口,需要单独添加接口。
利用satpy绘制真彩色图非常方便,给定 composite 参数即可,同时给定经纬度范围限制图片显示范围。
由于需要绘制的图形比较多,为了加快绘图速度,使用 dask 并行绘图:
3月15日0500UTC Himawari-8真彩色图
其实从真彩色图上已经能够看出沙尘的发展了。此外,satpy 还提供了 dust 合成产品,主要是通过不同的光谱通道对粒子反射特性进行计算。
3月15日0500UTC 沙尘合成产品
多说几句,除了 真彩色图 和 dust 的合成产品之外,satpy 还支持很多合成产品,比如 fog 和 convection 等,处理方法是类似,只需要更改 composite 参数即可。
THE END