PyCINRAD模块安装指南（利用Anaconda3安装）

小其其格

dataiyang 发表于 2019-7-22 11:55
你好，python新手想利用PyCINRAD处理雷达数据，安装这一步已经抓狂了！conda list查看了，需要的Cartopy、M ...

http://bbs.06climate.com/forum.p ... p;extra=&page=1

你到这个大神的帖子下去问问看，这些个问题他比较清楚些

dataiyang

小其其格 发表于 2019-7-22 21:43
http://bbs.06climate.com/forum.php?mod=viewthread&tid=88702&extra=&page=1

你到这个大神的帖子下 ...

楼主你好，安装问题已经解决了，都是很细节的问题！已经开始逐步画图了，在使用你提供的PPI.VCS.py脚本测试的时候，总是会提示TypeError: Cannot label gridlines on a AzimuthalEquidistant plot.  Only PlateCarree and Mercator plots are currently supported.我看了执行报错中提到，
  File "C:\Users\LR\Anaconda3\lib\site-packages\cartopy\mpl\gridliner.py", line 399, in _assert_can_draw_ticks
    'supported.'.format(prj=self.axes.projection))，请问你们是如何解决该问题的，非常感谢！

小其其格

dataiyang 发表于 2019-7-23 14:07
楼主你好，安装问题已经解决了，都是很细节的问题！已经开始逐步画图了，在使用你提供的PPI.VCS.py脚本测 ...

cartopy的这个功能还没有在最新的release当中包含，可以去GitHub更新一下，或者等cartopy发布下一个版本

小其其格

dataiyang 发表于 2019-7-23 14:07
楼主你好，安装问题已经解决了，都是很细节的问题！已经开始逐步画图了，在使用你提供的PPI.VCS.py脚本测 ...

1. # -*- coding: utf-8 -*-
   
2. """
   
3. Created on Sun Jan 13 21:51:35 2019
   
4. 
   
5. @author: Fulang WU
   
6. """
   
7. 
   
8. """定义一个函数，用于读取目标文件夹下的多个雷达基数据"""
   
9. #读取指定文件夹中的指定文件类型的文件名
   
10. import os
    
11. def get_filename(path,fileType):
    
12.     file_name =[]
    
13.     final_name = []
    
14.     for files in os.listdir(path):  #root为目录路径 #dirs为路径下的子目录 #files为路径下的所有非子目录
    
15.         if fileType in files:
    
16.            file_name.append(files.replace(fileType,''))#生成不带‘.bin’后缀的文件名组成的列表
    
17.            final_name = [path +item +fileType for item in file_name]#生成‘.bin’后缀的文件名组成的绝对路径列表
    
18.     return final_name #输出列表
    
19. """-----------------------------------------------------------------------------------------------------------"""
    
20. 
    
21. 
    
22. 
    
23. """读取目标文件夹下的多个雷达基数据"""
    
24. Dir = "D:/Python/Anaconda3/Scripts/Radar/" #目标文件夹
    
25. fileType = '.bin' #雷达基数据后缀名
    
26. #sa_radar_file = Dir+"Z_RADR_I_Z9571_20180304115800_O_DOR_SA_CAP.bin"
    
27. sa_radar_file = get_filename(Dir,fileType) #读取
    
28. print(sa_radar_file)
    
29. """----------------------------------------------------------------------"""
    
30. 
    
31. """---------------------------绘制雷达图像-------------------------------"""
    
32. import cinrad
    
33. from cinrad.io import CinradReader, StandardData
    
34. from cinrad.visualize import PPI
    
35. import matplotlib.ticker as mticker
    
36. from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
    
37. 
    
38. 
    
39. ele = 1 #第1个仰角
    
40. radius = 230 #绘制图像的范围大小
    
41. 
    
42. for nFiles in range(len(sa_radar_file)):
    
43.    
    
44.     f = CinradReader(sa_radar_file[nFiles]) #老版本数据
    
45.     rl = [f.get_data(i, radius, 'REF') for i in f.angleindex_r]
    
46.     vcs = cinrad.easycalc.VCS(rl)
    
47.    
    
48.    
    
49.     """
    
50.         剖面经纬度
    
51.     """
    
52.     ##1045
    
53.     #sec = vcs.get_section(start_cart=(118.8, 29.8), end_cart=(119.4, 30.2)) # 传入经纬度坐标
    
54.     ##1124
    
55.     #sec = vcs.get_section(start_cart=(119.4, 30.1), end_cart=(120.0, 30.4)) # 传入经纬度坐标
    
56.     ##1158
    
57.     sec = vcs.get_section(start_cart=(119.8, 30.2), end_cart=(120.4, 30.6)) # 传入经纬度坐标
    
58.     ##1259
    
59.     #sec = vcs.get_section(start_cart=(120.6, 30.6), end_cart=(121.4, 31.0)) # 传入经纬度坐标
    
60.    
    
61.    
    
62.     fig = cinrad.visualize.PPI(rl[ele])  #0代表第1个仰角
    
63.    
    
64.     fig.plot_range_rings(radius, color='white', linewidth=1.0) #绘制圆圈
    
65.     for i in range(0,radius-30,50):
    
66.         fig.plot_range_rings(i, color='white', linewidth=1.0) #绘制圆圈
    
67.         
    
68.     fig.plot_cross_section(sec)  #绘制垂直剖面
    
69.    
    
70.     """
    
71.     设置经纬度
    
72.     """   
    
73.     liner = fig.geoax.gridlines(draw_labels=True,linewidth=2, color='gray', alpha=0.5, linestyle='--')
    
74.     liner.xlabels_top = False
    
75.     liner.ylabels_right = False
    
76.     liner.xformatter = LONGITUDE_FORMATTER
    
77.     liner.yformatter = LATITUDE_FORMATTER
    
78.     liner.xlabel_style = {'size': 18, 'color': 'white'}
    
79.     #liner.xlabel_style = {'color': 'red', 'weight': 'bold'}
    
80.     liner.ylabel_style = {'size': 18, 'color': 'white'}
    
81.    
    
82.     #保存图像
    
83.     fig(Dir)

复制代码

你用以上代码绘图试试

dataiyang

小其其格 发表于 2019-7-23 14:18
你用以上代码绘图试试

你好，测试依然出现TypeError: Cannot label gridlines on a AzimuthalEquidistant plot.  Only PlateCarree and Mercator plots are currently supported.问题，另外可以分享Z_RADR_I_Z9571_20180304115800_O_DOR_SA_CAP.bin这个文件吗？邮箱20906279@zju.edu.cn

关于cartopy我的版本是0.17.0，github中搜不到cartopy，还是说我搜索的关键字不对？

dataiyang

dataiyang 发表于 2019-7-23 14:32
你好，测试依然出现TypeError: Cannot label gridlines on a AzimuthalEquidistant plot.  Only PlateCar ...

请问这个问题能否给一个答案，谢谢！

独孤酒见

各位大拿，我在安装Metpy时，用Python setup.py install  方式安装时提示如下：
Finished processing dependencies for MetPy==0+unknown

这个提示是啥意思，是安装成功了吗？可是conda list时看不到啊。

独孤酒见

各位大拿且@小其其格帮我看看问题。
我经过严格按照小其其格的方法来安装MetPy库，看起来是应该是安装上了，因为在pip list下看到了
H:\001.png

但是诡异的问题好像出现了：
1、我在运行Skewt-plot程序时：
    from skewt import SkewT

  File "D:\2017_Project\sound_data\skewt.py", line 7, in <module>
    import metpy.calc as mpcalc

ModuleNotFoundError: No module named 'metpy'

2、更神奇的是当我退出CMD或者退出Anaconda Prompt后，再次启动这两个命令，然后输入pip list，就发现那个metpy模块不在了。

我遇到的这都是什么鬼问题啊。

slddbs

独孤酒见 发表于 2019-7-25 15:28
各位大拿且@小其其格帮我看看问题。
我经过严格按照小其其格的方法来安装MetPy库，看起来是应该是安装上了 ...

有时候，清理电脑会导致重新加载库，你可以新建， 复制，重新运行

slddbs

独孤酒见 发表于 2019-7-25 15:28
各位大拿且@小其其格帮我看看问题。
我经过严格按照小其其格的方法来安装MetPy库，看起来是应该是安装上了 ...

有时候，清理电脑会导致重新加载库，你可以新建， 复制，重新运行