好久不见,格点气象的逗比们恭祝大家新年猴啊!
废话不多说,先献上红包一枚,不成敬意。
心意表达完毕就进入正题吧,今天为大家带来的是非常实用的“台风路径图的绘制技巧”,脚本的提供者是帅气的兵哥哥、气象界大咖@夏朗的芒果。推荐台风界大咖们参考参考。
先看一下效果图吧!观测数据(黑线,共40个时次)是采用了日本的最佳路径集数据,预报数据(彩色实线,每个起报时次分别有24h,4次预报)则是各家台风中心的路径预报结果。然后呢,不同的彩虹色对应着不同的起报时次,大方而美观有没有。由于测试数据是我们给定的,预报总时次约为观测的4倍,因此在绘图之前需要注意数据的预处理,保证在最后绘图的时候满足同一时次的数据一一对应~
下面我们就一起来看一下这幅图的绘制步骤。
第一步,应该不需要说了吧,做一些常规的事情:
可以看到呢,数据类型是大家已经非常熟悉的文本文件,在前几期已经有过比较详细的介绍。这里可以注意一下的是numAsciiRow这个函数,可用来获取数据的行数:
定义一系列数组来存放绘制路径所需要的经纬度信息,然后balabala对数据进行一些预处理方便后续对其进行绘制。需要注意的是,起报时次是6小时一次,每个起报时次进行24h的预报(即4个时次的预报结果),因此存放预报经纬度的数组的设定是二维,而观测数据是一维:
先绘制底图,这里应该不需要多说了吧~底图颜色的搭配根据自己的喜好就OK:
接着是路径的绘制,观测数据共有40个时次,三家预报中心的预报数据在每个时次对应有4时次的预报,采用一个简单的循环进行绘制:
使用gsn_add_polymarker函数给观测路径加上Marker标注,使其显得更为突出:
最后将各部分“组件”拼接成一幅完整的图片,需要注意的是,这里并没有采用传统“一页多图”的绘制手段,而是采用了“拼接”的方式使得图片更加紧凑、美观:
好了,这期的NCL实战绘图技巧就到这里了,祝各位绘制出满意的图片!
联系客服