VTK 表面重建技术（三维点云曲面重建）

最新推荐文章于 2024-04-28 10:30:53 发布

恋恋西风

最新推荐文章于 2024-04-28 10:30:53 发布

阅读量2.1k

点赞数 2

分类专栏： VTK 文章标签： VTK Reconstruction

本文链接：https://blog.csdn.net/q610098308/article/details/128107500

版权

VTK 专栏收录该内容

123 篇文章 76 订阅

订阅专栏

环境：windows 10 vs2017 VTK 8.2

1.点云重建

虽然Delaunay三角剖分算法可以实现网格曲面重建，但是其应用主要在二维剖分，在三维空间网格生成中遇到了问题。因为在三维点云曲面重建中，Delaunay条件不在满足，不仅基于最大最小角判断的对角线交换准则不在成立，而且基于外接圆判据的Delaunay三角化也不能保证网格质量。

VTKSurfaceReconstructionFilter则实现了一种隐式曲面重建方法，即将曲面看做一个符号距离函数的等值面，曲面内外的距离值得符号相反，而零等值面即为所求的曲面。该方法需要对点云数据进行网格划分，然后估算每个点的切平面和方向，并以每个点与最近的切平面距离来近似表面距离。这样即可得到一个符号距离的体数据。这样，我们就可以利用VTKContourFilter来提取零等值面即可得到相应的网格。

利用人脸点云数据进行人脸网格曲面重建实例如下：

2. 重建注意参数调整，可能会崩溃不同数据

#include "vtkAutoInit.h"
VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkInteractionStyle);
VTK_MODULE_INIT(vtkRenderingFreeType)

#include <vtkSmartPointer.h>
#include <vtkPolyDataReader.h>
#include <vtkPolyData.h>
#include <vtkSurfaceReconstructionFilter.h>
#include <vtkContourFilter.h>
#include <vtkVertexGlyphFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkCamera.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkProperty.h>
#include <vtkPLYReader.h>

int main()
{
     std::string filename = "E:/test.ply";    
    vtkSmartPointer<vtkPLYReader> reader =
        vtkSmartPointer<vtkPLYReader>::New();
    reader->SetFileName(filename.c_str());
    reader->Update();
 

    vtkSmartPointer<vtkPolyData> points =
        vtkSmartPointer<vtkPolyData>::New();
    points->SetPoints(reader->GetOutput()->GetPoints()); //获得网格模型中的几何数据：点集

    vtkSmartPointer<vtkSurfaceReconstructionFilter> surf =
        vtkSmartPointer<vtkSurfaceReconstructionFilter>::New();
    surf->SetInputData(points);
    surf->SetNeighborhoodSize(5);
    surf->SetSampleSpacing(0.5);
    surf->Update();

    vtkSmartPointer<vtkContourFilter> contour =
        vtkSmartPointer<vtkContourFilter>::New();
    contour->SetInputConnection(surf->GetOutputPort());
    contour->SetValue(0, 0.0);
    contour->Update();
    //
    vtkSmartPointer <vtkVertexGlyphFilter> vertexGlyphFilter =
        vtkSmartPointer<vtkVertexGlyphFilter>::New();
    vertexGlyphFilter->AddInputData(points);
    vertexGlyphFilter->Update();
    vtkSmartPointer<vtkPolyDataMapper> pointMapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
    pointMapper->SetInputData(vertexGlyphFilter->GetOutput());
    pointMapper->ScalarVisibilityOff();

    vtkSmartPointer<vtkActor> pointActor =
        vtkSmartPointer<vtkActor>::New();
    pointActor->SetMapper(pointMapper);
    pointActor->GetProperty()->SetColor(1, 0, 0);
    pointActor->GetProperty()->SetPointSize(4);

    vtkSmartPointer<vtkPolyDataMapper> contourMapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
    contourMapper->SetInputData(contour->GetOutput());
    vtkSmartPointer<vtkActor> contourActor =
        vtkSmartPointer<vtkActor>::New();
    contourActor->SetMapper(contourMapper);
    ///
    double pointView[4] = { 0, 0, 0.5, 1 };
    double contourView[4] = { 0.5, 0, 1, 1 };

    vtkSmartPointer<vtkRenderer> pointRender =
        vtkSmartPointer<vtkRenderer>::New();
    pointRender->AddActor(pointActor);
    pointRender->SetViewport(pointView);
    pointRender->SetBackground(1, 1, 1);

    vtkSmartPointer<vtkRenderer> contourRender =
        vtkSmartPointer<vtkRenderer>::New();
    contourRender->AddActor(contourActor);
    contourRender->SetViewport(contourView);
    contourRender->SetBackground(0, 1, 0);

    pointRender->GetActiveCamera()->SetPosition(0, -1, 0);
    pointRender->GetActiveCamera()->SetFocalPoint(0, 0, 0);
    pointRender->GetActiveCamera()->SetViewUp(0,0,1);
    pointRender->GetActiveCamera()->Azimuth(30);
    pointRender->GetActiveCamera()->Elevation(30);
    pointRender->ResetCamera();
    contourRender->SetActiveCamera(pointRender->GetActiveCamera());

    vtkSmartPointer<vtkRenderWindow> rw =
        vtkSmartPointer<vtkRenderWindow>::New();
    rw->AddRenderer(pointRender);
    rw->AddRenderer(contourRender);
    rw->SetSize(640, 320);
    rw->SetWindowName("3D Surface Reconstruction ");
    rw->Render();

    vtkSmartPointer<vtkRenderWindowInteractor> rwi =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    rwi->SetRenderWindow(rw);
    rwi->Initialize();
    rwi->Start();

    return 0;
}