使用h5py访问HDF5文件中命名类型(Datatypes)内部信息的问题

核心思路 #

h5m是MOAB的网格格式，其中的命名类型（Named Datatypes）本质是HDF5的自定义类型（如数组、复合类型），不能直接迭代或下标访问，但可以通过h5py的类型API读取结构信息，无需重新打开已获取的类型对象。

具体操作步骤 #

1. 确认命名类型对象

首先确保你遍历到的是h5py.Datatype实例，可通过类型检查验证：

import h5py

def check_datatype(obj):
    if isinstance(obj, h5py.Datatype):
        print("找到有效命名类型")
        return True
    return False

2. 根据类型类别读取信息

HDF5命名类型分多种类别（数组、复合、枚举等），需针对性处理：

针对数组类型（如你遇到的('<f8', (6,))）

这类是固定长度数组，可通过低级API获取元素类型和维度：

# 假设my_datatype是已获取的h5py.Datatype对象
type_id = my_datatype.id
if type_id.get_class() == h5py.h5t.ARRAY:
    # 获取数组维度
    dims = h5py.h5t.get_array_dims(type_id)
    # 获取元素基础类型
    base_type_id = h5py.h5t.get_super(type_id)
    base_dtype = h5py.h5t.dtype_from_id(base_type_id)
    
    print(f"数组形状: {dims}")
    print(f"元素类型: {base_dtype}")

针对复合类型（多字段结构）

如果是包含多个字段的复合类型，可遍历所有字段获取名称和类型：

if type_id.get_class() == h5py.h5t.COMPOUND:
    field_count = type_id.get_nmembers()
    print(f"复合类型包含{field_count}个字段：")
    for idx in range(field_count):
        field_name = type_id.get_member_name(idx)
        field_type_id = type_id.get_member_type(idx)
        field_dtype = h5py.h5t.dtype_from_id(field_type_id)
        print(f"  {field_name}: {field_dtype}")

其他类型（枚举、字符串等）

可通过type_id.get_class()判断类型类别，再调用对应API：

• 枚举类型：type_id.get_nmembers()获取枚举值数量，type_id.get_member_value(idx)获取对应值
• 字符串类型：type_id.get_size()获取长度，type_id.get_cset()获取编码集

3. 修正错误的打开方式

你之前用h5py.h5t.open(my_group.id, my_datatype.id.encode())报错，是因为my_datatype.id已经是类型的ID，不需要重新打开。如果要通过名称获取类型，正确写法是：

# 从父组中通过类型名称获取
target_datatype = my_group.get_type("your_type_name")

遍历所有命名类型的工具函数 #

可以用以下代码快速定位文件中所有命名类型并打印信息：

def process_datatype(name, obj):
    if isinstance(obj, h5py.Datatype):
        print(f"\n===== 命名类型: {name} =====")
        print(f"基础dtype: {obj.dtype}")
        type_id = obj.id
        type_class = type_id.get_class()
        print(f"类型类别: {type_class}")
        
        # 根据类别处理
        if type_class == h5py.h5t.ARRAY:
            dims = h5py.h5t.get_array_dims(type_id)
            base_type = h5py.h5t.dtype_from_id(h5py.h5t.get_super(type_id))
            print(f"数组信息: 形状{dims}, 元素类型{base_type}")
        elif type_class == h5py.h5t.COMPOUND:
            field_count = type_id.get_nmembers()
            print(f"复合字段列表:")
            for idx in range(field_count):
                fname = type_id.get_member_name(idx)
                ftype = h5py.h5t.dtype_from_id(type_id.get_member_type(idx))
                print(f"  {fname}: {ftype}")

# 遍历文件
with h5py.File("your_file.h5m", "r") as f:
    f.visititems(process_datatype)

补充说明 #

如果h5m文件中的命名类型是MOAB自定义的网格元数据（如节点、单元标签），也可以考虑使用MOAB的Python绑定pymoab直接解析，能更贴合格式逻辑，但仅用h5py也能获取所有类型结构信息。

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