在python中使用VTK求三维空间中两个柱面的交点

import vtk

cylinder1 = vtk.vtkCylinderSource()
cylinder1.SetRadius(1.0)  # 设置柱面1的半径
cylinder1.SetHeight(3.0)  # 设置柱面1的高度
cylinder1.SetResolution(100)  # 设置柱面1的分辨率

cylinder2 = vtk.vtkCylinderSource()
cylinder2.SetRadius(0.8)  # 设置柱面2的半径
cylinder2.SetHeight(4.0)  # 设置柱面2的高度
cylinder2.SetResolution(100)  # 设置柱面2的分辨率

mapper1 = vtk.vtkPolyDataMapper()
mapper1.SetInputConnection(cylinder1.GetOutputPort())

actor1 = vtk.vtkActor()
actor1.SetMapper(mapper1)

mapper2 = vtk.vtkPolyDataMapper()
mapper2.SetInputConnection(cylinder2.GetOutputPort())

actor2 = vtk.vtkActor()
actor2.SetMapper(mapper2)

renderer = vtk.vtkRenderer()
renderer.AddActor(actor1)
renderer.AddActor(actor2)

window = vtk.vtkRenderWindow()
window.AddRenderer(renderer)

interactor = vtk.vtkRenderWindowInteractor()
interactor.SetRenderWindow(window)

window.Render()
interactor.Start()

cylinder1_polydata = cylinder1.GetOutput()
cylinder2_polydata = cylinder2.GetOutput()

intersection = vtk.vtkIntersectionPolyDataFilter()
intersection.SetInputData(0, cylinder1_polydata)
intersection.SetInputData(1, cylinder2_polydata)
intersection.Update()

points = intersection.GetOutput().GetPoints()
num_points = points.GetNumberOfPoints()

if num_points > 0:
    for i in range(num_points):
        point = points.GetPoint(i)
        print("交点坐标{}：{}".format(i+1, point))
else:
    print("两个柱面无交点")