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

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

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

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

Cartopy学习使用

昨天本来打算用 matplotlib 和 Cartopy 画几张海温图的，结果意外发现我的 Arch Linux 更新软件包以后，估计是 libc 版本太新导致 matplotlib break 了，只要一调用 plt.show 或者是 plt.savefig 必崩，给我人整傻了，又是更换内核又是搞虚拟机，今天下午临时装了一个 Ubuntu 22.04 用着，白白折腾了一天的时间 🤡。

更新：我发现不是 matplotlib break 了，是 Cartopy 和新版的 Shapely 不兼容导致的。执行以下命令从源码重新安装 Shapely 即可恢复正常。

1

pip install "shapely<2" --no-binary shapely --force-reinstall

我之前就用过 Cartopy，但是并不是特别熟。今天画海温图的时候顺便探索了一下，发现了几件有趣的事情，所以分享一下。

xarray处理气象数据

最近做课程论文需要读取气象数据 (nc格式)，之前处理雷达数据时一直都是用 netCDF4 直接梭哈，但是处理全球大尺度的数据时，经常搞日平均月平均，单独取某些月做操作，这个时候 netCDF 就不够看了。正好学一学 xarray，水几篇博客。

Cython 学习和使用

写在前面

长期更新

强烈建议亲手编写代码进行操作

最近想写一些用 Python 控制 wpa_supplicant 的程序，然而找到的轮子似乎都是通过 d-Bus 实现的，没有用官方提供的接口的 (也许是因为只提供了 C 的接口吧)。于是想自己通过 Cython 在 Python 中调用这些接口，无奈 Cython 水平太低，只好再学学 Cython。

通过 Google 搜到了一本英文的手册：《Cython: A Guide for Python Programmers》，在此记录学习过程做备忘。

这篇文章会一直更新，直到我看完这本书为止，至于什么时候更新嘛，咳咳….应该不会一直拖着不更，考完了研，人就是闲

使用 pyart 绘制多雷达拼图

前言

这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制单个雷达图以及多雷达拼图的教程。

这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的，想着既然要写自己的博客了那就一块搬过来吧，当然原博客的地址还是要加的

由于中间修改过代码，并不确定代码 100% 正确

原地址：Py-ART 简易中文教程

用 pyart 绘制雷达图像并进行反射率修正

前言

这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制雷达图并进行反射率修正的教程。

这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的，想着既然要写自己的博客了那就一块搬过来吧，当然原博客的地址还是要加的

由于中间修改过代码，并不确定代码 100% 正确

原地址：Py-ART 简易中文教程

用 pyart 绘制简单的雷达图像

前言

这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制简单雷达图的教程。

这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的，想着既然要写自己的博客了那就一块搬过来吧，当然原博客的地址还是要加的

由于中间修改过代码，并不确定代码 100% 正确

原地址：Py-ART 简易中文教程