使用MetPy Level2File读取2002年.Z格式雷达数据出现多类报错

问题原因 #

你遇到的报错核心是两方面问题：

1. 旧NEXRAD Level II格式兼容问题：2002年的雷达数据属于NEXRAD Build 7之前的旧格式，采用Message 1存储反射率等数据，MetPy 1.1.0的Level2File接口原生仅支持2007年雷达系统升级后采用Message 31的新标准格式，对2008年之前的旧格式兼容性极差，你见到的Unknown message、TypeError、IndexError都是因为接口无法正确识别旧消息类型、读取到空的属性值导致的。
2. .Z格式解压与路径问题：.Z是Unix下compress工具压缩的文件格式，Windows环境下MetPy默认不支持直接读取该压缩格式，同时你代码中写的Windows路径没有加原始字符串前缀，反斜杠会被识别为转义字符，也会触发OSError: [Errno 22] Invalid argument报错。

解决方案 #

• 首先修正路径写法：给文件路径加r前缀转义反斜杠，示例：

  raddir = r"C:\Users\rconn\Documents\WORK\Cases\20020203_Lake_Effect\"
  

• 提前解压.Z文件：用7-Zip等工具把所有.Z后缀的文件提前解压为无后缀的原始二进制文件，再传入Level2File接口测试。
• 替换读取接口为Py-ART：Py-ART库对新旧NEXRAD Level II格式的兼容性远好于MetPy，你可以安装Py-ART后用pyart.io.read_nexrad_archive接口读取旧数据，提取所需的反射率、仰角、经纬度等信息后再用你现有逻辑绘图即可。
• 升级MetPy版本：如果要继续使用MetPy，可升级到最新正式版，新版已补充部分旧Level II格式的适配逻辑，可兼容部分2005年后的旧数据。

内容的提问来源于stack exchange，提问作者Ryan Connelly