vtk获取三角面片的坐标

在 VTK 中，获取三角面片的坐标需要先获取三角面片的顶点坐标，然后再根据顶点坐标计算得到三角面片的坐标。可以使用 vtkTriangleFilter 过滤器获取三角面片，然后使用 vtkCellArray 类的 GetCell 方法获取三角面片的顶点索引，再根据索引获取顶点坐标。以下是示例代码：

#include <vtkSmartPointer.h>
#include <vtkTriangleFilter.h>
#include <vtkPolyData.h>
#include <vtkCellArray.h>
#include <vtkPoints.h>

int main()
{
    // 创建一个立方体
    vtkSmartPointer<vtkCubeSource> cubeSource = vtkSmartPointer<vtkCubeSource>::New();

    // 使用三角面片过滤器获取立方体的三角面片
    vtkSmartPointer<vtkTriangleFilter> triangleFilter = vtkSmartPointer<vtkTriangleFilter>::New();
    triangleFilter->SetInputConnection(cubeSource->GetOutputPort());
    triangleFilter->Update();

    // 获取三角面片的顶点坐标
    vtkPolyData* polyData = triangleFilter->GetOutput();
    vtkCellArray* cells = polyData->GetPolys();
    vtkPoints* points = polyData->GetPoints();

    // 遍历三角面片，获取顶点坐标
    for (vtkIdType i = 0; i < cells->GetNumberOfCells(); i++)
    {
        vtkCell* cell = cells->GetCell(i);
        vtkIdType* pointIds = cell->GetPointIds()->GetPointer(0);
        double* p0 = points->GetPoint(pointIds[0]);
        double* p1 = points->GetPoint(pointIds[1]);
        double* p2 = points->GetPoint(pointIds[2]);

        // 输出顶点坐标
        cout << "Triangle " << i << ":\n";
        cout << "(" << p0[0] << ", " << p0[1] << ", " << p0[2] << ")\n";
        cout << "(" << p1[0] << ", " << p1[1] << ", " << p1[2] << ")\n";
        cout << "(" << p2[0] << ", " << p2[1] << ", " << p2[2] << ")\n";
    }

    return 0;
}

这段代码会输出立方体的所有三角面片的顶点坐标。你可以根据自己的需求修改代码，只输出你需要的三角面片的坐标。