Gmsh二进制网格格式解析及Python转ASCII方案咨询

背景 #

用Gmsh Python API生成大型网格，ASCII格式文件可达千兆字节，因此改用二进制格式保存（体积更小、读写更快）。现在需要在Python中读取二进制网格文件，转换为可解析的ASCII格式，兼容现有依赖ASCII的代码库，要求内存中转换且速度优于直接读取ASCII文件。

具体问题 #

1. Gmsh保存的二进制.msh文件采用何种编码/格式？
2. 如何在Python中读取二进制网格文件并转换为与ASCII导出完全一致的格式？

测试代码 #

以下是可切换ASCII/二进制保存的简化Gmsh教程代码：

import gmsh

binary = True
argv = ["", "-bin"] if binary else []
gmsh.initialize(argv=argv)

gmsh.model.add("t1")

lc = 1e-2
gmsh.model.geo.addPoint(0, 0, 0, lc, 1)
gmsh.model.geo.addPoint(.1, 0, 0, lc, 2)
gmsh.model.geo.addPoint(.1, .3, 0, lc, 3)
p4 = gmsh.model.geo.addPoint(0, .3, 0, lc)

gmsh.model.geo.addLine(1, 2, 1)
gmsh.model.geo.addLine(3, 2, 2)
gmsh.model.geo.addLine(3, p4, 3)
gmsh.model.geo.addLine(4, 1, p4)

gmsh.model.geo.addCurveLoop([4, 1, -2, 3], 1)
gmsh.model.geo.addPlaneSurface([1], 1)

gmsh.model.geo.synchronize()

gmsh.model.addPhysicalGroup(1, [1, 2, 4], 5)
gmsh.model.addPhysicalGroup(2, [1], name="My surface")

gmsh.model.mesh.generate(2)
gmsh.write("t1_binary.msh" if binary else "t1_ascii.msh")
gmsh.finalize()

尝试过的方案及报错 #

直接以二进制模式读取文件后调用.decode('utf-8')，在$Entities行之后解码失败，报错：

UnicodeDecodeError: 'utf-8' codec can't decode byte 0x9a in position 72: invalid start byte

Linux下执行file t1_binary.msh仅返回t1_binary.msh: data，无法获取格式细节。

补充说明 #

• 已查看gmsh.write()源码，未找到明确的二进制格式定义
• 要求转换后的ASCII内容与binary=False时保存的t1_ascii.msh完全一致

解决方案 #

问题1：Gmsh二进制文件格式说明 #

Gmsh的二进制.msh文件是自定义混合格式：

• 文件开头的$MeshFormat、$Entities等标记行是ASCII文本
• 标记后的实际数据段（节点、单元等）采用二进制存储，使用小端字节序，数据类型对应Gmsh的版本定义（比如3.0版本默认用8字节双精度浮点数、8字节整数）

问题2：Python中二进制转ASCII的高效方法 #

最可靠且高效的方式是直接利用Gmsh Python API完成内存中的格式转换，无需手动解析二进制：

import gmsh

# 初始化Gmsh并读取二进制网格
gmsh.initialize()
gmsh.open("t1_binary.msh")

# 直接在内存中导出为ASCII格式的字符串
ascii_mesh_content = gmsh.write_string()

# 此时ascii_mesh_content就是与t1_ascii.msh完全一致的ASCII内容
# 可直接传递给依赖ASCII格式的代码，或写入文件验证
with open("converted_ascii.msh", "w") as f:
    f.write(ascii_mesh_content)

gmsh.finalize()

为什么不推荐手动解析？

手动解析需要严格遵循Gmsh的.msh格式规范，处理不同版本的格式差异（比如2.2和4.0版本的二进制存储逻辑不同），实现复杂且容易出错，远不如直接调用API高效可靠。

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