使用pyvista.PolyDataFilters.boolean_union()时内核崩溃问题求助

内核崩溃的可能原因 #

布尔运算（boolean_union）对输入网格的拓扑和几何质量要求极高，你的情况大概率是以下问题导致：

• 网格拓扑缺陷：通过delaunay_3d()生成的体网格，提取表面后可能存在非流形边、重复面/点，或者内部残留的无效几何结构，这些都会让布尔运算的底层算法出现计算异常，直接引发内核崩溃。
• 数值精度问题：两个建筑的接触面虽然看起来共面，但浮点坐标的微小误差可能导致算法无法正确识别重合区域，进而陷入计算死循环或内存溢出。

解决步骤 #

1. 替换网格生成方式：手动构建表面网格 #

放弃delaunay_3d()的方式，它是为散乱点生成体网格设计的，对规则立方体来说完全没必要，反而会引入冗余结构。直接用PyVista的Cube类构建建筑的表面网格，保证拓扑结构干净。

示例代码：

import numpy as np
import pyvista as pv

# 定义建筑1的边界（min_x, max_x, min_y, max_y, min_z, max_z）
building1_bounds = [52.56, 102.47, 12.13, 42.43, 0, 6]
building1 = pv.Cube(bounds=building1_bounds).triangulate()

# 定义建筑2的边界
building2_bounds = [52.56, 94.67, 0, 12.13, 0, 6]
building2 = pv.Cube(bounds=building2_bounds).triangulate()

2. 预处理原有网格：清理冗余点和面 #

如果坚持用原有网格生成方式，先对网格做清理操作，消除拓扑缺陷：

# 合并重合点（设置容差处理浮点误差）
building1_clean = building1_triangulated.clean(tolerance=1e-6)
building2_clean = building2_triangulated.clean(tolerance=1e-6)

# 删除重复面
building1_clean.remove_duplicate_cells()
building2_clean.remove_duplicate_cells()

3. 执行布尔运算 #

使用清理后的网格执行布尔并集，同时指定容差参数帮助算法处理数值误差：

# 执行布尔并集，设置容差
union_mesh = building2_clean.boolean_union(building1_clean, tolerance=1e-6)

# 可视化结果
union_mesh.plot(show_edges=True, opacity=0.5, color="gray")

4. 备选方案：手动合并共面区域 #

如果布尔运算还是有问题，可以手动处理两个建筑的共面部分：

# 找到建筑1中y=12.13的接触面
building1_contact_faces = building1.extract_cells(building1.points[:,1] == 12.13)
# 找到建筑2中y=12.13的接触面
building2_contact_faces = building2.extract_cells(building2.points[:,1] == 12.13)

# 删除接触面的面
building1_no_contact = building1.remove_cells(building1_contact_faces.cell_ids)
building2_no_contact = building2.remove_cells(building2_contact_faces.cell_ids)

# 合并剩余网格
final_mesh = building1_no_contact.merge(building2_no_contact, tolerance=1e-6)
final_mesh.plot(show_edges=True, opacity=0.5, color="gray")

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