如何编辑HDF5文件部分内容？修改降水率数据集不生效求助

原代码存在3个核心错误导致修改未生效：

• 逐元素判断逻辑完全错误：代码中new_value[...] = 1是对整个new_value数组全局赋值，只要遍历过程中遇到任意一个大于等于10的数值，整个数组都会被全部设为1，完全没有实现逐元素判断的需求。
• 修改未写入HDF5文件：precip = f['Grid/precipitationCal'][0][:][:]是将HDF5数据集的内容读取到内存的numpy数组中，后续修改内存里的precip变量不会自动同步到HDF5文件，需要直接对HDF5数据集对象执行写入操作。
• 嵌套循环效率极低：numpy原生支持向量化运算，无需嵌套循环遍历每一个元素，处理速度可提升数万倍。

import h5py
import numpy as np
import glob

filenames = glob.glob("/IMERG/Exceedance/2014_E/3B-HHR.MS.MRG.3IMERG.201401*")

for file in filenames:
    # 用with上下文管理器自动管理文件句柄，避免异常时文件未关闭
    with h5py.File(file,'r+') as f:
        # 先获取HDF5数据集对象
        precip_dset = f['Grid/precipitationCal']
        # 读取第一个时间步的降水数据到内存
        precip_data = precip_dset[0, :, :]
        # 向量化二值化：>=10赋值为1，否则为0
        new_value = (precip_data >= 10.0).astype(int)
        # 直接将结果写回HDF5数据集对应位置
        precip_dset[0, :, :] = new_value

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