pyart中关于线性规划部分代码的疑问
我画大饼说这个周末就要改完本科论文的,结果 pyart
的线性规划部分的函数越看越迷惑,写一篇帮自己梳理一下。
欲望大过了能力,由此产生了焦虑
我画大饼说这个周末就要改完本科论文的,结果 pyart
的线性规划部分的函数越看越迷惑,写一篇帮自己梳理一下。
昨天开始改本科论文,发现了一个问题:如何给自己绘制的雷达 PPI 图添加上地图底图。
画地图底图的话肯定是用 Cartopy
或者是 cnmaps
来画,但是问题是雷达数据存放的形式是极坐标形式,而 Cartopy
或者 cnmaps
绘制地图使用的是经纬度坐标。将雷达数据由极坐标系转换为经纬度坐标系有一个小问题,就是数据量大。haversine
虽然可以根据两个地理位置的距离和角度计算经纬度,但是好像只能计算两个点之间的,一个一个点算的话太慢了。
后来我想到的解决办法是绘制两个图层,然后分别叠加到一起。matplotlib
的 Axes
类有一个 set_alpha
函数可以设置背景透明度,只需要让覆盖在上面的图层背景透明,显示出下面的图,就可以达到图层叠加的目的了。
昨天本来打算用 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
,但是并不是特别熟。今天画海温图的时候顺便探索了一下,发现了几件有趣的事情,所以分享一下。
最近做课程论文需要读取气象数据 (nc
格式),之前处理雷达数据时一直都是用 netCDF4
直接梭哈,但是处理全球大尺度的数据时,经常搞日平均月平均,单独取某些月做操作,这个时候 netCDF
就不够看了。正好学一学 xarray
,水几篇博客。
长期更新
强烈建议亲手编写代码进行操作
最近想写一些用 Python 控制 wpa_supplicant 的程序,然而找到的轮子似乎都是通过 d-Bus 实现的,没有用官方提供的接口的 (也许是因为只提供了 C 的接口吧)。于是想自己通过 Cython 在 Python 中调用这些接口,无奈 Cython 水平太低,只好再学学 Cython。
通过 Google 搜到了一本英文的手册:《Cython: A Guide for Python Programmers》,在此记录学习过程做备忘。
这篇文章会一直更新,直到我看完这本书为止,至于什么时候更新嘛,咳咳….应该不会一直拖着不更,考完了研,人就是闲
这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制单个雷达图以及多雷达拼图的教程。
这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的,想着既然要写自己的博客了那就一块搬过来吧,当然原博客的地址还是要加的
由于中间修改过代码,并不确定代码 100% 正确
原地址:Py-ART 简易中文教程
这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制雷达图并进行反射率修正的教程。
这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的,想着既然要写自己的博客了那就一块搬过来吧,当然原博客的地址还是要加的
由于中间修改过代码,并不确定代码 100% 正确
原地址:Py-ART 简易中文教程
这是一篇关于使用 Pthon 库 Py-ART 读取雷达数据文件、绘制简单雷达图的教程。
这篇文章是我在去年暑假–也就是 2020年8月–做大创项目的时候写的,想着既然要写自己的博客了那就一块搬过来吧,当然原博客的地址还是要加的
由于中间修改过代码,并不确定代码 100% 正确
原地址:Py-ART 简易中文教程