如何用Python重命名NetCDF变量？遇只读权限报错求助

错误原因 #

使用nc.Dataset(fn)打开文件时，默认是只读模式（mode='r'），无法执行重命名变量这类修改操作，因此抛出“Write to read only”错误。

解决方案 #

方案1：以读写模式打开原文件直接修改 #

修改文件打开方式，添加mode='r+'参数获取读写权限，修改后关闭文件确保更改生效。操作前建议备份原文件，避免意外损坏。

#! /usr/bin/env python3
# coding: utf-8
import netCDF4 as nc
fn='L1AHITBR'
# 以读写模式打开文件
ds = nc.Dataset(fn, mode='r+')
print(ds.variables.keys())
ds.renameVariable('himawari_9_ahi_channel_10_brightness_temperature','himawari_8_ahi_channel_10_brightness_temperature')
# 关闭文件保存修改
ds.close()

方案2：复制文件到新文件并重命名变量（更安全） #

如果原文件是只读状态，或者NetCDF格式（如经典格式）不支持原地修改，推荐创建新文件复制所有内容，仅修改目标变量名称，避免损坏原文件。

#! /usr/bin/env python3
# coding: utf-8
import netCDF4 as nc

input_fn = 'L1AHITBR'
output_fn = 'L1AHITBR_renamed.nc'

# 上下文管理器自动处理文件关闭
with nc.Dataset(input_fn, 'r') as src, nc.Dataset(output_fn, 'w') as dst:
    # 复制全局属性
    dst.setncatts(src.__dict__)
    
    # 复制所有维度
    for name, dim in src.dimensions.items():
        dst.createDimension(name, len(dim) if not dim.isunlimited() else None)
    
    # 遍历变量，复制并重命名目标变量
    for name, var in src.variables.items():
        if name == 'himawari_9_ahi_channel_10_brightness_temperature':
            new_var_name = 'himawari_8_ahi_channel_10_brightness_temperature'
            # 创建新变量
            new_var = dst.createVariable(new_var_name, var.datatype, var.dimensions)
            # 复制变量数据
            new_var[:] = var[:]
            # 复制变量属性
            new_var.setncatts(var.__dict__)
        else:
            # 复制其他变量
            new_var = dst.createVariable(name, var.datatype, var.dimensions)
            new_var[:] = var[:]
            new_var.setncatts(var.__dict__)

额外注意事项 #

• 方案1操作前务必备份原文件，原地修改存在文件损坏风险。
• 确认当前用户对原文件有写入权限，否则即使指定r+模式也会报错。
• 若原文件为NetCDF经典格式（而非NetCDF4），原地修改可能受格式限制，优先使用方案2。

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