Python：如何打开 RAIN1 Vaisala sigmet IRIS 文件？已经尝试过 Py-ART、wradlib 和 xradar 产生以下错误

Question

我正在尝试从维萨拉 sigmet iris 雷达打开二进制雷达数据文件。我尝试打开的文件是指 IRIS 产品而不是原始数据。

对于来自命令的 Py-Art：

r = pyart.io.read('path_to_filename/filename.RN1ZPL1')
I get the error:
Traceback (most recent call last):
  File "/snap/pycharm-community/344/plugins/python-ce/helpers/pydev/pydevconsole.py", line 364, in runcode
    coro = func()
  File "<input>", line 1, in <module>
  File "/path/venv/lib/python3.10/site-packages/pyart/io/auto_read.py", line 125, in read
    return read_sigmet(filename, **kwargs)
  File "/path/venv/lib/python3.10/site-packages/pyart/io/sigmet.py", line 146, in read_sigmet
    sigmet_data, sigmet_metadata = sigmetfile.read_data(full_xhdr=full_xhdr)
  File "pyart/io/_sigmetfile.pyx", line 133, in pyart.io._sigmetfile.SigmetFile.read_data
  File "/path/venv/lib/python3.10/site-packages/numpy/ma/core.py", line 8437, in __call__
    result = self._func.__call__(*args, **params).view(MaskedArray)
ValueError: negative dimensions are not allowed

对于来自命令的 wralib：

r = wrl.io.open_iris_dataset('path_to_filename/filename.RN1ZPL1')
I get the error:
Traceback (most recent call last):
  File "/snap/pycharm-community/344/plugins/python-ce/helpers/pydev/pydevconsole.py", line 364, in runcode
    coro = func()
  File "<input>", line 1, in <module>
  File "/path/venv/lib/python3.10/site-packages/wradlib/io/iris.py", line 527, in open_iris_dataset
    return open_radar_dataset(filename_or_obj, engine=IrisBackendEntrypoint, **kwargs)
  File "/path/Projects/venv/lib/python3.10/site-packages/wradlib/io/xarray.py", line 345, in open_radar_dataset
    groups = _get_iris_group_names(filename_or_obj)
  File "/path/venv/lib/python3.10/site-packages/xradar/io/backends/iris.py", line 3950, in _get_iris_group_names
    with opener(filename, loaddata=False) as ds:
TypeError: 'bool' object is not callable

对于来自命令的 xradar 和 fsspec：

file = fsspec.open_local(path + 'filename.RN1ZPL1')
radar_data = xd.io.open_iris_datatree(file).xradar.georeference()
and i get:
Traceback (most recent call last):
  File "/snap/pycharm-community/344/plugins/python-ce/helpers/pydev/pydevconsole.py", line 364, in runcode
    coro = func()
  File "<input>", line 1, in <module>
  File "/path/venv/lib/python3.10/site-packages/xradar/io/backends/iris.py", line 4091, in open_iris_datatree
    sweeps = _get_iris_group_names(filename_or_obj)
  File "/path/venv/lib/python3.10/site-packages/xradar/io/backends/iris.py", line 3950, in _get_iris_group_names
    with opener(filename, loaddata=False) as ds:
TypeError: 'bool' object is not callable

Answer 1

我刚刚发现你可以使用 gdal (https://gdal.org/index.html) 或 QGIS (https://www.qgis.org/en/site/)。