NCL(7)—台风路径图

好久不见，格点气象的逗比们恭祝大家新年猴啊！

废话不多说，先献上红包一枚，不成敬意。

---

---

     心意表达完毕就进入正题吧，今天为大家带来的是非常实用的“台风路径图的绘制技巧”，脚本的提供者是帅气的兵哥哥、气象界大咖@夏朗的芒果。推荐台风界大咖们参考参考。

---

先看一下效果图吧！观测数据(黑线，共40个时次)是采用了日本的最佳路径集数据，预报数据(彩色实线，每个起报时次分别有24h，4次预报)则是各家台风中心的路径预报结果。然后呢，不同的彩虹色对应着不同的起报时次，大方而美观有没有。由于测试数据是我们给定的，预报总时次约为观测的4倍，因此在绘图之前需要注意数据的预处理，保证在最后绘图的时候满足同一时次的数据一一对应~

---

下面我们就一起来看一下这幅图的绘制步骤。

第一步，应该不需要说了吧，做一些常规的事情：

可以看到呢，数据类型是大家已经非常熟悉的文本文件，在前几期已经有过比较详细的介绍。这里可以注意一下的是numAsciiRow这个函数，可用来获取数据的行数：

定义一系列数组来存放绘制路径所需要的经纬度信息，然后balabala对数据进行一些预处理方便后续对其进行绘制。需要注意的是，起报时次是6小时一次，每个起报时次进行24h的预报(即4个时次的预报结果)，因此存放预报经纬度的数组的设定是二维，而观测数据是一维：

先绘制底图，这里应该不需要多说了吧~底图颜色的搭配根据自己的喜好就OK：

接着是路径的绘制，观测数据共有40个时次，三家预报中心的预报数据在每个时次对应有4时次的预报，采用一个简单的循环进行绘制：

使用gsn_add_polymarker函数给观测路径加上Marker标注，使其显得更为突出：

最后再加上图右侧的标注：

最后将各部分“组件”拼接成一幅完整的图片，需要注意的是，这里并没有采用传统“一页多图”的绘制手段，而是采用了“拼接”的方式使得图片更加紧凑、美观：

好了，这期的NCL实战绘图技巧就到这里了，祝各位绘制出满意的图片！