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