为雷达 PPI 图像添加地图底图

昨天开始改本科论文，发现了一个问题：如何给自己绘制的雷达 PPI 图添加上地图底图。

画地图底图的话肯定是用 Cartopy 或者是 cnmaps 来画，但是问题是雷达数据存放的形式是极坐标形式，而 Cartopy 或者 cnmaps 绘制地图使用的是经纬度坐标。将雷达数据由极坐标系转换为经纬度坐标系有一个小问题，就是数据量大。haversine 虽然可以根据两个地理位置的距离和角度计算经纬度，但是好像只能计算两个点之间的，一个一个点算的话太慢了。

后来我想到的解决办法是绘制两个图层，然后分别叠加到一起。matplotlib 的 Axes 类有一个 set_alpha 函数可以设置背景透明度，只需要让覆盖在上面的图层背景透明，显示出下面的图，就可以达到图层叠加的目的了。