使用MetPy读取NEXRAD DPR L3产品的单位确认及数据异常排查

你犯的核心错误是没有使用Level3数据自带的缩放和偏移参数，而是自行假设了转换系数。NEXRAD的Level3产品（包括DPR降水产品）的原始16位整数，需要通过数据块里的scale和offset计算才能得到正确的物理量，不是直接除以1000。

修正后的代码 #

from metpy.io import Level3File

radar = Level3File(filename)
datadict = radar.sym_block[0][0]
data = radar.map_data(datadict['data'])
# 用数据自带的scale和offset转换为真实物理量
precip_in_per_hour = data * datadict['scale'] + datadict['offset']
# 转换为毫米每小时
precip_mm_per_hour = precip_in_per_hour * 25.4

为什么会差32倍？ #

你遇到的32倍差异，大概率是因为该DPR产品的scale参数是1/32.0（即原始整数代表1/32英寸每小时）。你之前用除以1000来转英寸，相当于把系数当成了1/1000，两者的比值是1000/32≈31.25，正好和你看到的32倍差异吻合。

关键提醒 #

不同的NEXRAD Level3产品，scale和offset参数都可能不同，绝对不能硬编码转换逻辑，必须从datadict中读取这两个值来计算真实物理量。

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