绘制四面体(切割后封闭)

最新推荐文章于 2025-03-11 15:18:36 发布

努力减肥的小胖子5

最新推荐文章于 2025-03-11 15:18:36 发布

阅读量686

点赞数

CC 4.0 BY-SA版权

分类专栏： VTK

本文链接：https://blog.csdn.net/yuxing55555/article/details/103476129

VTK 专栏收录该内容

58 篇文章

订阅专栏

本文介绍使用VTK库创建四面体并进行平面剪切的过程。通过具体的C++代码示例，展示了如何定义四面体的顶点和单元，设置平面的原点和法线，应用剪切操作并可视化结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

参考url：https://lorensen.github.io/VTKExamples/site/Cxx/GeometricObjects/Cell3DDemonstration

效果图：
在这里插入图片描述

实现：


vtkSmartPointer<vtkUnstructuredGrid> MakeTetrahedron(){
    // Make a tetrahedron.
    int numberOfVertices = 4;

    vtkSmartPointer<vtkCellArray> cellArray =
      vtkSmartPointer<vtkCellArray>::New();

    vtkSmartPointer< vtkPoints > points =
      vtkSmartPointer< vtkPoints > :: New();
    points->InsertNextPoint(0, 0, 0);
    points->InsertNextPoint(1, 0, 0);
    points->InsertNextPoint(0, 1, 0);
    points->InsertNextPoint(0, 0, 1);

    points->InsertNextPoint(3, 0, 0);
    points->InsertNextPoint(4, 0, 0);
    points->InsertNextPoint(3, 1, 0);
    points->InsertNextPoint(3, 0, 1);
    {
        vtkSmartPointer<vtkTetra> tetra =
          vtkSmartPointer<vtkTetra>::New();
        for (int i = 0; i < numberOfVertices; ++i)
        {
            tetra->GetPointIds()->SetId(i, i);
        }
        cellArray->InsertNextCell(tetra);
    }
    {
        vtkSmartPointer<vtkTetra> tetra =
          vtkSmartPointer<vtkTetra>::New();
        for (int i = 0; i < numberOfVertices; ++i)
        {
            tetra->GetPointIds()->SetId(i, i+4);
        }
        cellArray->InsertNextCell(tetra);
    }

    vtkSmartPointer<vtkUnstructuredGrid> unstructuredGrid =
      vtkSmartPointer<vtkUnstructuredGrid>::New();
    unstructuredGrid->SetPoints(points);
    unstructuredGrid->SetCells(VTK_TETRA, cellArray);
    //unstructuredGrid->InsertNextCell(VTK_TETRA, 4, )

    return unstructuredGrid;
}


int main(int argc, char *argv[])
{
    vtkSmartPointer<vtkUnstructuredGrid> test = MakeTetrahedron();
    vtkSmartPointer<vtkPlane> plane =
      vtkSmartPointer<vtkPlane>::New();
    plane->SetOrigin(0,0,0.5);
    plane->SetNormal(0.0, 0.0, -1.0);

    vtkSmartPointer<vtkClipDataSet> clipper =
      vtkSmartPointer<vtkClipDataSet>::New();
    clipper->SetInputData(test);
    clipper->SetClipFunction(plane);
    //clipper->SetValue(0);
    clipper->Update();

    vtkSmartPointer<vtkDataSetMapper> mapper =
        vtkSmartPointer<vtkDataSetMapper>::New();
    mapper->SetInputData(clipper->GetOutput());

    vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
    actor->SetMapper(mapper);

    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(actor);

    vtkSmartPointer<vtkRenderWindow> renderWindow =
        vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);

    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    renderWindowInteractor->SetRenderWindow(renderWindow);

    vtkSmartPointer<vtkInteractorStyleTrackballCamera> style =
        vtkSmartPointer<vtkInteractorStyleTrackballCamera>::New();
    renderWindowInteractor->SetInteractorStyle(style);
    renderWindowInteractor->Start();

    return 0;
}