一、安装
1.首先要安装Anaconda,最好安装2020年以后的3.8或更高版本(太高的版本可能后续匹配下载会比较慢)。
推荐在镜像资源下载
Index of /anaconda/archive/ | 清华大学开源软件镜像站 | Tsinghua Open Source Mirror
2.安装作为前置的cartopy库,建议使用conda命令安装,省事方便无需其他操作。
conda install cartopy -c conda-forge --yes
但如果网络环境不好,进度卡住,也可以去搜索下载cartopy及它的前置库whl文件(记得有3、4个),使用pip命令安装cartopy的几个前置库,再安装cartopy。当然版本要能对应才行,具体可以再去网上查查。安装完成后去试运行一下cartopy的代码,可能有缺的部分,会有报错提醒去下载缺失文件,解压放到指定目录即可。
3.最后安装pycwr,两种方式:git源码和pip。建议使用后者,无需其他操作,运气好10分钟安装完,运气不好的话速度会很慢(换源可能会好一点)。
使用git源码下载大概率会遇到报错,主流办法是下一个微软的编译软件解决,也可以自行再找办法。(我没有试过这个方法,是搜索报错信息得到的,好像是因为有部分c++的代码,需要用软件转换成python的。不太了解这方面,如有错误还望指正。)
二、使用
1.这个是读取文件的代码,文件后缀为.bz2,自行替换路径即可。会自动分辨雷达数据类型,当然也可以指定,具体请看官方文档。
from pycwr.io import read_auto
filename = r"./data/Z_RADR_I_Z9898_20190828181529_O_DOR_SAD_CAP_FMT.bin.bz2"
2.绘图具体用法可见官方文档,下面仅放一个案例:
from pycwr.io import read_auto
import matplotlib.pyplot as plt
from pycwr.draw.RadarPlot import Graph
filename=r"D:\zzzhqx\radarImage\RadarImage\RENY2\Z_RADR_I_Z9731_20230803055644_O_DOR_SAD_CAP_FMT.bin.bz2"
PRD = read_auto(filename)
fig, ax = plt.subplots()
graph = Graph(PRD)
graph.plot_ppi(ax, 0, "dBZ", cmap="CN_ref") ## 0代表第一层, dBZ代表反射率产品
graph.add_rings(ax, [0, 50, 100, 150, 200, 250, 300])
ax.set_title("PPI Plot", fontsize=16)
ax.set_xlabel("Distance From Radar In East (km)", fontsize=14)
ax.set_ylabel("Distance From Radar In North (km)", fontsize=14)
plt.show()
效果如图:
个人使用来说,应该matplotlib的操作都能行吧,至少像去掉横纵坐标什么的,跟matplotlib的一样。下面放官方网址:
安装方法 — pycwr 0.3.4 documentationhttps://github.com/rubyzhang1/pycwr/blob/master/README_CN.md