循环向NetCDF变量写入字符串列表时仅保留最后一个值的问题求助

我一眼就看出问题所在了——你在循环里每次都直接覆盖了v和station的全部元素，而不是只更新当前站点对应的位置，所以最后所有位置都会被最后一次循环的值覆盖。咱们来一步步修正这个问题：

问题根源分析 #

看你的循环代码：

for i in range(len(S)):
    # ... 读取文件的代码
    v[:] = nc.stringtochar(datain1)
    station[:] = nc.chartostring(v[:])

每次循环都用v[:]把整个字符变量填满当前站点的字符串，然后又用station[:]把整个站点变量全部替换，自然最后所有元素都是最后一个站点的值。

修正后的代码 #

import netCDF4 as nc
import numpy as np
import xarray as xr

# 创建站点名称列表
S = []
for i in listObsFiles:
    if i.endswith('0h2021.nc'):
        statID = i[8:13]
        S.append(statID)

# 初始化输出NetCDF文件
waveObsNC = nc.Dataset(file, 'w')

# 定义维度：nstrings对应站点数量，nchars对应每个站点ID的长度（5个字符）
nstrings = waveObsNC.createDimension('nstrings', len(S))
nchars = waveObsNC.createDimension('nchars', 5)

# 创建变量：station是字符串类型，维度为nstrings；v是字符数组（若需要中转）
station = waveObsNC.createVariable('station', 'U5', ('nstrings',))  # 用Unicode字符串更便捷，也可保留S5
v = waveObsNC.createVariable('v', 'S1', ('nstrings', 'nchars'))  # 中转用字符数组，维度对应每个站点的字符

# 循环处理每个站点，仅更新对应索引的位置
for i in range(len(S)):
    File = xr.open_dataset(f'{filepath}saveWave{S[i]}h2021.nc')
    st = File.attrs['station']  # 或直接用st = S[i]，两者值一致
    
    # 方式1：直接写入station变量（推荐，最简洁）
    station[i] = st
    
    # 方式2：若必须用字符数组中转，更新v的第i行即可
    # datain = np.array([st], dtype='S5')
    # v[i, :] = nc.stringtochar(datain)[0]

# 若使用中转方式，最后一次性转换字符串
# station[:] = nc.chartostring(v[:])

waveObsNC.close()

# 测试输出
print(station[:])

关键修改点 #

1. 维度修正：把v的维度改成('nstrings', 'nchars')，让每个站点对应一行字符，避免每次循环覆盖整个数组。
2. 索引赋值：循环中用station[i]或v[i, :]更新单个站点的位置，而非station[:]或v[:]覆盖全部元素。
3. 简化逻辑：其实不需要v这个中转变量，直接把字符串赋值给station变量即可（NetCDF4支持直接写入字符串数组），如果有版本兼容问题再用中转方式。

额外优化建议 #

你可以完全去掉循环中转的逻辑，直接把整个S列表转成数组写入，效率更高：

# 直接写入整个站点列表，无需循环读取每个文件
station[:] = np.array(S, dtype='U5')

这样就不用循环处理每个文件来获取站点名称了——毕竟你已经从文件名提取到完整的S列表了。

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