CyanideCN/PyCINRAD

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三个都是基数据

更新了问题