如何编辑hdf5文件的一部分

Question

我正在尝试编辑现有 hdf5 文件中的降水率值，以便将 >= 10 的值重写为 1，将

import h5py
import numpy as np
import glob

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

for file in filenames:
    f = h5py.File(file,'r+')
    new_value = np.zeros((3600, 1800))
    new_value = new_value.astype(int)
    precip = f['Grid/precipitationCal'][0][:][:]

    for i in precip:
        for j in i:
            if j >= 10.0:
                new_value[...] = 1
            else:
                pass
    precip[...] = new_value
    f.close()

Answer 1

您似乎只是在更改数组的值，而实际上并未更新文件对象中的任何内容。另外，我会摆脱那个 for 循环 - 它很慢！试试这个：

import h5py
import numpy as np
import glob

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

for file in filenames:
    f = h5py.File(file,'r+')
    precip = f['Grid/precipitationCal'][0][:][:]

    # Replacing the for loop
    precip[precip>10.0] = 1

    # Assign values
    f['Grid/precipitationCal'][0][:][:] = precip
    f.close()

Answer 2

您似乎没有将新值写入文件，而只是将它们存储在数组中。