后处理及个例结果分析 ¶

后处理及个例结果分析

绘图资料准备 ¶

绘图程序是基于python产生的绘图程序。在绘图之前，需要进行资料准备。

如果准备绘制原始高分辨率的资料（2km）则需执行:
```
> cd run_B5_16
> ./zz4_tank.sh
```
所产生的绘图输入文件在 data/dumfiles 里::

> cd data/dumpfiles 201910080000_tanked.50
如果只绘制稀疏化后分辨率的有云和晴空AHI卫星资料，则只需执行:
```
> cd run_B5_16
> ./ zz5_run_dump.sh
```

所产生的绘图输入文件在 data/dumfiles 里:

> cd data/dumpfiles
201910080000_grid.50
201910080000_gridcm.50

python绘图 ¶

进入plot目录，开始绘图，具体如下:

> tree -L 1 plot
plot/
├── backup
├── fulldisk.py
├── images
├── python_functions
├── rd_dump.py
└── rd_dump.py_all

开始绘图:

> python rd_dump.py   (绘制稀疏化后有云和晴空资料分布图)
> python fulldisk.py  (绘制原始高分辨率(2km)资料分布图)

产生的图形存储在images 目录下:

> tree -L 1 images
images
├── ahi_grid_201910080000_bt07.png
├── ahi_grid_201910080000_bt08.png
├── ahi_grid_201910080000_bt09.png
├── ahi_grid_201910080000_bt10.png
├── ahi_grid_201910080000_bt11.png
├── ahi_grid_201910080000_bt12.png
├── ahi_grid_201910080000_bt13.png
├── ahi_grid_201910080000_bt14.png
├── ahi_grid_201910080000_bt15.png
├── ahi_grid_201910080000_bt16.png
├── ahi_gridcm_201910080000_bt07.png
├── ahi_gridcm_201910080000_bt08.png
├── ahi_gridcm_201910080000_bt09.png
├── ahi_gridcm_201910080000_bt10.png
├── ahi_gridcm_201910080000_bt11.png
├── ahi_gridcm_201910080000_bt12.png
├── ahi_gridcm_201910080000_bt13.png
├── ahi_gridcm_201910080000_bt14.png
├── ahi_gridcm_201910080000_bt15.png
└── ahi_gridcm_201910080000_bt16.png

个例结果分析 ¶

这里我们以2019年10月08日00点00分的一次Himiwari-8卫星观测为例进行云检测结果分析。Himawari-8卫星完成一次全景扫描所需时间是10分钟，原始资料每分钟存储一次，所以每一次全景扫描每个观测通道包含10个独立文件。

现以水汽通道CH10为例，原始资料如下:

HS_H08_20191008_0100_B10_FLDK_R20_S0110.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0210.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0310.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0410.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0510.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0610.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0710.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0810.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S0910.DAT.bz2
HS_H08_20191008_0100_B10_FLDK_R20_S1010.DAT.bz2

经过原始资料格式转换，可以得到以下包括所有观测通道的国际气象组织的标准格式BUFR文件:

HS_H08_20191008_0000_B05-16_FLDK_R20_S0110.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0210.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0310.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0410.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0510.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0610.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0710.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0810.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S0910.bufr
HS_H08_20191008_0000_B05-16_FLDK_R20_S1010.bufr

经过云分析系统处理，将有云和晴空区分出来，我们可以得到以下两个BUFR文件:

ahi.201910080000.allsky.bufr  (有云卫星资料)
ahi.201910080000.bufr         (晴空卫星资料)

卫星资料云分析前后对比:

../../_images/ahi_tanked_201910080000_bt10.png

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 原始分辨率(2km)

../../_images/ahi_grid_201910080000_bt10.png

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 有云，稀疏化后分辨率(36km)

../../_images/ahi_gridcm_201910080000_bt10.png

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 晴空，稀疏化后分辨率(36km)

注解

超强台风海贝思（英语：Typhoon Hagibis，国际编号：1919，联合台风警报中心：WP202019）是2019年太平洋台风季第19个被命名的热带气旋。超强台风海贝思(Hagibis)成为2019年第三个达到日本气象厅所评定“猛烈的”形容词级别之风暴，成为今年最强的风暴。发生时间为2019年10月5日至2019年10月20日，在10月8日左右升级为强台风。
所选个例卫星观测为2019年10月8日00点，恰好在台风最强时段，经过对比分析，可以发现：1）从原始2公里高分辨率云图中能够清晰地看到台风的位置及其周围云的分布情况。2）经过云分析和稀疏化处理以后，首先可以看到观测外围资料因为质量原因被剔除，其次台风眼区资料被剔除明显。3）特别是晴空云观测信息，在台风眼区几乎完全被剔除。

不同分辨率稀疏化结果分析 ¶

不同分辨率有云卫星资料对比:

../../_images/ahi_grid_201910080000_bt101.png

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 有云，稀疏化后分辨率(06km)

../../_images/ahi_grid_201910080000_bt102.png

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 有云，稀疏化后分辨率(18km)

葵花(Himawari-8) 静止卫星红外水汽通道(Ch10)资料分布图: 有云，稀疏化后分辨率(36km)

注解

读者还可以进行：1）不同分辨率晴空卫星资料对比和 2）其他红外通道对比。此处限于篇幅，故此省略。

后处理及个例结果分析¶

绘图资料准备¶

python绘图¶

个例结果分析¶

不同分辨率稀疏化结果分析¶