pyart中关于线性规划部分代码的疑问
我画大饼说这个周末就要改完本科论文的,结果 pyart
的线性规划部分的函数越看越迷惑,写一篇帮自己梳理一下。
欲望大过了能力,由此产生了焦虑
我画大饼说这个周末就要改完本科论文的,结果 pyart
的线性规划部分的函数越看越迷惑,写一篇帮自己梳理一下。
昨天开始改本科论文,发现了一个问题:如何给自己绘制的雷达 PPI 图添加上地图底图。
画地图底图的话肯定是用 Cartopy
或者是 cnmaps
来画,但是问题是雷达数据存放的形式是极坐标形式,而 Cartopy
或者 cnmaps
绘制地图使用的是经纬度坐标。将雷达数据由极坐标系转换为经纬度坐标系有一个小问题,就是数据量大。haversine
虽然可以根据两个地理位置的距离和角度计算经纬度,但是好像只能计算两个点之间的,一个一个点算的话太慢了。
后来我想到的解决办法是绘制两个图层,然后分别叠加到一起。matplotlib
的 Axes
类有一个 set_alpha
函数可以设置背景透明度,只需要让覆盖在上面的图层背景透明,显示出下面的图,就可以达到图层叠加的目的了。
老高让我学 WRF
,我就装 WRF
。装 WRf
的过程出现了很多问题,就开贴记录一下。
WRF 版本 V4.4.1
昨天本来打算用 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
,水几篇博客。
距离上次写日记是10天前,这期间到也没什么事,就只是忙于作业,也没遇到什么困难。昨天开始想期末的课程论文,人要麻了,又是天气学,又是气象卫星数据应用,最后还要答辩,感觉好多事情。今天一天除了上课吃饭就是窝在宿舍里,天气很晴,阳光很足,也不热,tmd,为什么我要窝在宿舍里搞作业,操。
今天下午不去办公室了,直接在宿舍一边上网课一边开摆,提前放假!
首先当然是下载好zsh
的安装包然后传到服务器上了。
在这个网站可以搜索所有的rpm
包,并下载,适用于RHEL 7的zsh
安装包下载地址如下