SA格式雷达数据读取错误求助
jinhuanz opened this issue · 5 comments
强大的雷达数据处理工具,三行命令读取数据,还能方便地转成nc格式👍
我这边处理其中一个数据遇到如下问题:
IndexError Traceback (most recent call last)
in
2 for f in sorted(fs):
3 print(f)
----> 4 data = PUP(f).get_data()
5 fig = PPI(data)
6 fig('./')
~/miniconda3/envs/radar/lib/python3.7/site-packages/cinrad-1.6.3-py3.7.egg/cinrad/io/level3.py in init(self, file)
41 from metpy.io.nexrad import Level3File
42
---> 43 f = Level3File(file)
44 # Because metpy interface doesn't provide station codes,
45 # it's necessary to reopen it and read the code.
~/miniconda3/envs/radar/lib/python3.7/site-packages/MetPy-1.0-py3.7.egg/metpy/io/nexrad.py in init(self, filename)
1724 # Now that we have the header, we have everything needed to make tables
1725 # Store as class that can be called
-> 1726 self.map_data = mapper(self)
1727
1728 # Process compression if indicated. We need to fail
~/miniconda3/envs/radar/lib/python3.7/site-packages/MetPy-1.0-py3.7.egg/metpy/io/nexrad.py in init(self, prod)
962 label = ''
963 if codes >> 7:
--> 964 label = self.lut_names[val]
965 if label in ('Blank', 'TH', 'ND'):
966 val = self.MISSING
IndexError: list index out of range
数据文件:
Z_RADR_I_Z9220_20210414114800_P_DOR_SA_HI_NUL_230_NUL.220.bin.zip
你这个还是PUP产品吧
<xarray.Dataset>
Dimensions: (azimuth: 364, distance: 230)
Coordinates:
* azimuth (azimuth) float64 0.001745 0.01905 0.03636 ... 6.268 0.001745
* distance (distance) float64 0.5 1.0 1.5 2.0 ... 113.5 114.0 114.5 115.0
Data variables:
VEL (azimuth, distance) float64 nan nan nan 4.5 ... nan nan nan nan
longitude (azimuth, distance) float64 117.7 117.7 117.7 ... 117.7 117.7
latitude (azimuth, distance) float64 39.05 39.05 39.06 ... 40.07 40.08
RF (azimuth, distance) float64 nan nan nan nan ... nan nan nan nan
Attributes:
elevation: 1.5
range: 115
scan_time: 2021-04-14 11:48:00
site_code: Z9220
site_name: 塘沽
site_longitude: 117.718
site_latitude: 39.043
tangential_reso: 0.5
emmm.... 看起来是的...
雷达数据小白 :(
我以为R、V、W三个都是基数据
更新了问题