PyCINRAD模块处理雷达基数据指南（发布1.7.0版本）(修复经纬度标注问题)

单纯xiao可爱

eeeee 发表于 2021-4-10 14:49
组合反射率暂时不支持拼图

这样子，谢谢老师回复。

sskkyy2008

请教老师1.62版本 ：想从基数据的REF中提取某个经纬度的dbz 程序代码这么写对吗？
print(path)
                f = CinradReader(path) #老版本数据
                station=path.split("/")[-1].split('_') [3]

                fstaionlon=f.stationlon
                f.set_code(station)
                f._update_radar_info()
                f.stationlat = 30.113
                f.stationlon = 102.859
                fdata=f.data
                r = f.get_data(0, 150, 'REF')
               
#                r2, xout, yout = cinrad.grid.grid_2d(r['REF'].values,r['longitude'].values,r['latitude'].values)
#                print(r2)
#                lon=(r['longitude'])
#                print('---------------------')
#                print(lon)
#                print(r.variables['longitude'])
               
#                tree = spatial.KDTree(np.dstack((r.lon.ravel(),r.lat.ravel()))[0])
                tree = spatial.KDTree(np.dstack((r['longitude'],r['latitude']))[0])
               
               
                treedata=tree.data
               
#                lonnum=101
#                latnum=26.5
                lonnum=102.865
                latnum=30.096
                t2=tree.query(np.array([lonnum,latnum]))
                print(t2)
                t3=np.unravel_index(t2[1],r['REF'].values.shape)
                print(t3)
                t4=r['REF'].values[t3]
                print(t4)
这里的t4 是否是提取的反射率值

单纯xiao可爱

eeeee 发表于 2020-4-8 20:41
这里是一个历史遗留问题，要素的名字没有处理好，暂时手动覆盖一下吧
data = f.get_data()
data.dtype  ...

遇到问题没有多加思考，后面注意，已解决。（ps:用python的老师是真的厉害）

单纯xiao可爱

单纯xiao可爱 发表于 2021-4-13 19:42
老师我又又来打扰了，之前您之前说的这个设置，现在好像不大管用了，呜呜，咋办呀

画图问题大概解决了，程序问题不太懂。。。我是又菜又菜菜。还是写代码的老师厉害。
之前我用的是2D PRODUCT，看见前辈之前有设置dtype，就没思考粘贴了一下，出现了问题。
然后我又试了3DREF，没想太多又用dtype，无果
最后用的CREF。。。这才是我最先在SWAN检查过的数据。。这就画出来了。
老师大概已经改过了还是咋样，用CREF不需要设置dtype。

藿稥丶誫炁

老师您好！ 感谢您开发的模块~ 现在使用时有两个问题想请教一下您。
我使用的时.cinradreadar读取CD雷达数据，可以读取。
想使用to_netcdf导出数据时遇到问题：1.是否只能按仰角get_data后才能进行输出？不能直接将所有仰角一并输出？2.当使用f.get_data(1, 50, 'REF') 读取某仰角数据后使用to_netcdf报错：TypeError: Invalid value for attr 'scan_time': datetime.datetime(2020, 5, 14, 14, 59, 50) must be a number, a string, an ndarray or a list/tuple of numbers/strings for serialization to netCDF files。

还有一个问题：cinradreadar 是否不可使用to_pyart模块功能？
谢谢~

青帘影

老师您好，我在windows系统下，用anaconda虚拟环境安装完成后pycharm可以检索到cinrad插件，但是import时出错：import命令为：
import cinrad
from cinrad.visualize import PPI
from cinrad.visualize import Section
from cinrad.visualize import RHI
报错代码为：
ImportError: cannot import name 'RHI' from 'cinrad.visualize' (C:\ProgramData\Anaconda3\envs\pyart_env\lib\site-packages\cinrad-1.6.3-py3.8-win-amd64.egg\cinrad\visualize\__init__.py)
请问这是什么原因导致的？

200862584@qq.co

楼主，现在新格式的Z_RADR_I_Z9531_20210412062633_O_DOR_SAD_CAP_FMT.bin，就这种文件有个FMT的是不是不能正常读取

蓝色随风

您好，能请问一下您用的是哪位作者什么方法做的HCA就是粒子相态识别

初见101

请问您这种情况是无法读取我们这的数据吗？
Traceback (most recent call last):
  File "E:/Ls/cinrad/main.py", line 11, in <module>
    f = cinrad.io.CinradReader(path)  # 老版本数据
  File "D:\ProgramData\Anaconda\lib\site-packages\cinrad\io\level2.py", line 96, in __init__
    raise RadarDecodeError(
cinrad.error.RadarDecodeError: Unable to determine the file type. Use `radar_type` keywordto specify the radar type.

航在路

老师您好，我用1.63版读取CC双偏振雷达数据，显示的可用变量没差分反射率等双偏振变量，请问这是什么原因导致的？
g.available_product(0)
['REF', 'VEL', 'SW', 'TREF']