使用vtk，求取两个牙齿咬合的polydata

要使用 VTK 计算两个牙齿的咬合，可以按照以下步骤进行：

1. 加载两个牙齿的 STL 模型文件，使用 vtkSTLReader 类读取模型数据，并分别存储为 vtkPolyData 对象。

vtkSmartPointer<vtkSTLReader> reader1 = vtkSmartPointer<vtkSTLReader>::New();
reader1->SetFileName("tooth1.stl");
reader1->Update();
vtkSmartPointer<vtkPolyData> polyData1 = reader1->GetOutput();

vtkSmartPointer<vtkSTLReader> reader2 = vtkSmartPointer<vtkSTLReader>::New();
reader2->SetFileName("tooth2.stl");
reader2->Update();
vtkSmartPointer<vtkPolyData> polyData2 = reader2->GetOutput();

2. 对每个牙齿的 vtkPolyData 对象进行预处理，包括计算表面法向量、扫描线三角形插值等操作。

vtkSmartPointer<vtkPolyDataNormals> normals1 = vtkSmartPointer<vtkPolyDataNormals>::New();
normals1->SetInputData(polyData1);
normals1->ComputePointNormalsOn();
normals1->ComputeCellNormalsOff();
normals1->Update();

vtkSmartPointer<vtkPolyDataNormals> normals2 = vtkSmartPointer<vtkPolyDataNormals>::New();
normals2->SetInputData(polyData2);
normals2->ComputePointNormalsOn();
normals2->ComputeCellNormalsOff();
normals2->Update();

vtkSmartPointer<vtkCleanPolyData> cleaner1 = vtkSmartPointer<vtkCleanPolyData>::New();
cleaner1->SetInputData(normals1->GetOutput());
cleaner1->Update();

vtkSmartPointer<vtkCleanPolyData> cleaner2 = vtkSmartPointer<vtkCleanPolyData>::New();
cleaner2->SetInputData(normals2->GetOutput());
cleaner2->Update();

3. 计算两个牙齿的交点，使用 vtkIntersectionPolyDataFilter 类进行计算。

vtkSmartPointer<vtkIntersectionPolyDataFilter> intersectionFilter = vtkSmartPointer<vtkIntersectionPolyDataFilter>::New();
intersectionFilter->SetInputData(0, cleaner1->GetOutput());
intersectionFilter->SetInputData(1, cleaner2->GetOutput());
intersectionFilter->Update();
vtkSmartPointer<vtkPolyData> intersectionPolyData = intersectionFilter->GetOutput();

4. 最后，将交点的 vtkPolyData 对象保存为 STL 文件或其他格式的文件。

vtkSmartPointer<vtkSTLWriter> stlWriter = vtkSmartPointer<vtkSTLWriter>::New();
stlWriter->SetFileName("intersection.stl");
stlWriter->SetInputData(intersectionPolyData);
stlWriter->Update();
stlWriter->Write();

以上是使用 VTK 计算两个牙齿咬合的大致步骤，具体实现可能需要根据实际情况进行调整。