代码之家 › 专栏 › 技术社区 › ComputerScientist

将PyTables/HDF5文件中的所有数组从float64转换为float32

h5py pytables hdf5 numpy arrays

ComputerScientist · 技术社区 · 9 年前

我有一个包含大量子目录的PyTables文件。我有一种遍历表中所有数组数据类型的方法。它们是浮动的64;我想转换文件 在正确的位置

import h5py
import numpy as np

# filehead is a string for a file
with h5py.File(filehead, 'r+') as f:
    # Lots of stuff here ... e.g. `head` is a string

    print("/obsnorm/Standardizer/count {}".format(f[head+'/obsnorm/Standardizer/count']))
    print("count value: {}".format(f[head+'/obsnorm/Standardizer/count'].value))
    f[head+'/obsnorm/Standardizer/count'][...] = (f[head+'/obsnorm/Standardizer/count'].value).astype('float32')
    print("/obsnorm/Standardizer/count {}".format(f[head+'/obsnorm/Standardizer/count']))
    print("count value: {}".format(f[head+'/obsnorm/Standardizer/count'].value))

不幸的是,打印的结果是:

/obsnorm/Standardizer/count <HDF5 dataset "count": shape (), type "<f8">
count value: 512364.0
/obsnorm/Standardizer/count <HDF5 dataset "count": shape (), type "<f8">
count value: 512364.0

换句话说,在赋值之前,计数的类型是f8或float64。强制转换后,类型是浮动64。

1 回复 | 直到 9 年前

ComputerScientist 9 年前

正如hpaulj在评论中所建议的那样,我决定简单地重新创建一个重复的HDF5文件,除了生成类型为的数据集之外 f4 (与float32相同),我能够实现我的编码目标。

伪代码如下:

import h5py
import numpy as np

# Open the original file jointly with new file, with `float32` at the end.
with h5py.File(oldfile, 'r') as f, h5py.File(newfile[:-3]+'_float32.h5', 'w') as newf:
    # `head` is some directory structure
    # Create groups to follow the same directory structure
    newf.create_group(head)

    # When it comes time to create a dataset, make the cast here.
    newdata = (f[head+'/name_here'].value).astype('float32')
    newf.create_dataset(head+'/name_here', data=newdata, dtype='f4')

    # Proceed for all other datasets.

推荐文章

Megrez7 · C#ToArray转换合并为一行,导致数组元素更改

1 年前

bairog · 从按属性筛选的对象数组字典中创建值数组

1 年前

Anka HanÄ±m · 关于结构和动态数组地址的问题

1 年前

Swapnil Supekar · 将二维阵列转换为一维阵列,并将一维阵列粘贴到另一张图纸上

1 年前

Lorenzo Bertolaccini · 在Angular项目中通过对话框后,在控制台中显示但在HTML中不显示的数据数组

1 年前

MaSc. H. · 大小与阵列大小

1 年前

Geremia · 2D NumPy数组+1D数组?

1 年前

MARTIN · 交换第一个和最后一个单词,反转所有中间的字符

1 年前

Paul Williams · 迭代数组时输出有问题

1 年前

Oliver Morgan · PostgreSQL错误:无法累积不同维度的数组

1 年前