如何使用PyVista实现两个平面的相交并可视化？

错误原因 #

你使用的pv.Plane是有限范围的多边形平面（默认是边长为1的正方形），而非数学意义上的无限平面。尽管两个平面在数学上相交，但默认的小尺寸+低分辨率网格导致VTK的相交算法无法检测到有效交点，从而报错。

解决方案 #

方案1：调整平面参数，让有限平面范围重叠 #

通过增大平面尺寸、提高网格分辨率，确保两个有限平面的区域有重叠，让相交算法能识别交点：

import pyvista as pv

# 创建尺寸更大、分辨率更高的平面
p1 = pv.Plane(center=(0, 0, 0), direction=(0, 0, 1), size=5, resolution=10)
p2 = pv.Plane(center=(0, 0, 0), direction=(1, 0, 0), size=5, resolution=10)

# 计算相交结果
intersection = p1.intersection(p2)

# 可视化
plotter = pv.Plotter()
plotter.add_mesh(p1, color='blue', opacity=0.5)
plotter.add_mesh(p2, color='red', opacity=0.5)
if intersection:
    plotter.add_mesh(intersection, color='green', line_width=5)
plotter.show()

方案2：直接计算无限平面的交线（更高效准确） #

如果需要的是数学上无限平面的交线，可通过平面方程推导后直接创建交线：

• 平面1（方向(0,0,1)）方程：z=0
• 平面2（方向(1,0,0)）方程：x=0
• 交线为x=0且z=0，即Y轴方向的直线

代码示例：

import pyvista as pv

# 创建交线（可按需调整线段长度）
intersection_line = pv.Line(pointa=(0, -5, 0), pointb=(0, 5, 0))

# 可视化
plotter = pv.Plotter()
plotter.add_mesh(pv.Plane(center=(0,0,0), direction=(0,0,1)), color='blue', opacity=0.5)
plotter.add_mesh(pv.Plane(center=(0,0,0), direction=(1,0,0)), color='red', opacity=0.5)
plotter.add_mesh(intersection_line, color='green', line_width=5)
plotter.show()

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