使用os.stat处理netCDF4数组时出现ValueError：路径含嵌入空字符

我有一个大型netCDF4数组数据集，想要获取其类型、属性和维度大小等基础信息。原NCAR Command Level（NCL）代码使用printVarSummary()函数实现该功能，我认为Python的os.stat()可完成相同操作，但调用时出现错误：ValueError: stat: embedded null character in path。

Python代码 #

T2_input = np.ma.core.MaskedArray()
with Dataset(dir + 'filename.nc') as file_T2:
    T2_raw = file_T2.variables['T2MEAN'][strt * 8 : (strt + doy) * 8 - 1, :, :] - 273.15
    dimsizes_T2 = T2_raw.shape
    T2_input = T2_raw[4 : (dimsizes_T2[0]-4) : 8, :, :]  
print(os.stat(T2_input))

NCL代码 #

f_T2=addfile(dir + "filename.nc","r")
T2_raw = f_T2->T2MEAN(strt*8:(strt+doy)*8-1,:,:)-273.15
dsizes_T2 = dimsizes(T2_raw)
T2_input = T2_raw(4:(dsizes_T2(0)-4):8,:,:)
delete(T2_raw)
printVarSummary(T2_input)

错误原因 #

os.stat()的作用是读取文件系统中文件的元数据（比如文件大小、修改时间），它的参数必须是文件路径字符串，但你传入的T2_input是一个numpy掩码数组（np.ma.MaskedArray），属于内存中的数据对象，不是文件路径。

当os.stat()尝试将这个数组对象转换为路径字符串时，数组内部的二进制数据包含了空字符（\0），这是文件路径中不允许的字符，因此触发了embedded null character in path错误。

而且os.stat()完全无法实现NCL中printVarSummary()的功能——后者是用来查看内存中数组/变量的维度、类型、属性等信息，和文件系统元数据无关。

正确实现方法 #

要实现类似printVarSummary()的功能，直接用numpy或netCDF4的内置方法查看数组/变量信息即可：

1. 查看数组基础信息（维度、类型等） #

# 维度大小
print("维度大小：", T2_input.shape)
# 数据类型
print("数据类型：", T2_input.dtype)
# 掩码数组的掩码元素数量
if isinstance(T2_input, np.ma.MaskedArray):
    print("被掩码的元素数量：", T2_input.mask.sum())
# 打印数组概览
print(T2_input)
# 用numpy内置函数查看详细信息
np.info(T2_input)

2. 保留并查看原netCDF变量的属性 #

如果需要查看原netCDF变量的元数据（比如单位、描述等），不要直接将变量切片为数组，而是保留netCDF变量对象：

with Dataset(dir + 'filename.nc') as file_T2:
    # 保留原netCDF变量对象
    T2_var = file_T2.variables['T2MEAN']
    # 遍历打印变量属性
    print("原变量属性：")
    for attr_name in T2_var.ncattrs():
        print(f"  {attr_name}: {getattr(T2_var, attr_name)}")
    
    # 后续数组切片操作不变
    T2_raw = T2_var[strt * 8 : (strt + doy) * 8 - 1, :, :] - 273.15
    dimsizes_T2 = T2_raw.shape
    T2_input = T2_raw[4 : (dimsizes_T2[0]-4) : 8, :, :]

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