vtk-ExtractSurface（官例）

最新推荐文章于 2024-05-22 00:15:00 发布
最新推荐文章于 2024-05-22 00:15:00 发布
阅读量1.1k 收藏 1
点赞数
分类专栏： VTK
原文链接：https://vtk.org/Wiki/VTK/Examples/Cxx/Points/ExtractSurface
版权

转载：https://vtk.org/Wiki/VTK/Examples/Cxx/Points/ExtractSurface
链接：https://blog.csdn.net/hopefullight/article/details/79157144

//拟合类的算法以Hugues Hoppe提出的Signed Distance Function和Poisson方法为代表。所得曲面并不完全与数据点重合，而是进行了一定的合并和简化，并对噪声有一定鲁棒性。
//经典的Signed Distance Function重建算法主要流程如下：
//1、对每个数据点，搜索其邻域数据点，使用特征向量方法计算法线和法平面
//2、由于法线的方向可有正负两个不确定，故对全局使用最小生成树近似计算法线朝向
//3、以立方体素(voxel)为单位，沿着曲面滑动延展，计算每个格顶点到法平面的距离
//4、使用Marching Cube算法从立方体素中(插值地)提取三角面
//vtkExtractSurface使用了类似的方法，并可以利用Range的可见性信息以确定法线方向。
#include <vtkSmartPointer.h>
#include <vtkPLYReader.h>
#include <vtkXMLPolyDataReader.h>
#include <vtkOBJReader.h>
#include <vtkSTLReader.h>
#include <vtkPointSource.h>

#include <vtkPCANormalEstimation.h>
#include <vtkSignedDistance.h>
#include <vtkExtractSurface.h>
#include <vtkPointData.h>

#include <vtkPolyDataMapper.h>
#include <vtkProperty.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkCamera.h>

#include <vtksys/SystemTools.hxx>

static vtkSmartPointer<vtkPolyData> ReadPolyData(const char *fileName);

int main(int argc, char *argv[]) {
    vtkSmartPointer<vtkSTLReader> reader =
        vtkSmartPointer<vtkSTLReader>::New();
    reader->SetFileName("../data/test.stl");
    reader->Update();

    vtkSmartPointer<vtkPolyData> polyData =
        vtkSmartPointer<vtkPolyData>::New();
    polyData = reader->GetOutput();

    double bounds[6];
    polyData->GetBounds(bounds);
    double range[3];
    for(int i = 0; i < 3; ++i) {
        range[i] = bounds[2 * i + 1] - bounds[2 * i];
    }

    int sampleSize = polyData->GetNumberOfPoints() * .00005;
    if(sampleSize < 10) {
        sampleSize = 10;
    }
    std::cout << "Sample size is: " << sampleSize << std::endl;
    // Do we need to estimate normals?
    vtkSmartPointer<vtkSignedDistance> distance =
        vtkSmartPointer<vtkSignedDistance>::New();
    if(polyData->GetPointData()->GetNormals()) {
        std::cout << "Using normals from input file" << std::endl;
        distance->SetInputData(polyData);
    } else {
        std::cout << "Estimating normals using PCANormalEstimation" << std::endl;
        vtkSmartPointer<vtkPCANormalEstimation> normals =
            vtkSmartPointer<vtkPCANormalEstimation>::New();
        normals->SetInputData(polyData);
        normals->SetSampleSize(sampleSize);
        normals->SetNormalOrientationToGraphTraversal();
        normals->FlipNormalsOn();
        distance->SetInputConnection(normals->GetOutputPort());
    }
    std::cout << "Range: "
              << range[0] << ", "
              << range[1] << ", "
              << range[2] << std::endl;
    int dimension = 256;
    double radius = range[0] * .02;
    radius = range[0] / static_cast<double>(dimension) * 3; // ~3 voxels
    std::cout << "Radius: " << radius << std::endl;

    distance->SetRadius(radius);
    distance->SetDimensions(dimension, dimension, dimension);
    distance->SetBounds(
        bounds[0] - range[0] * .1,
        bounds[1] + range[0] * .1,
        bounds[2] - range[1] * .1,
        bounds[3] + range[1] * .1,
        bounds[4] - range[2] * .1,
        bounds[5] + range[2] * .1);

    vtkSmartPointer<vtkExtractSurface> surface =
        vtkSmartPointer<vtkExtractSurface>::New();
    surface->SetInputConnection(distance->GetOutputPort());
    surface->SetRadius(radius * .99);
    surface->Update();

    vtkSmartPointer<vtkPolyDataMapper> surfaceMapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
    surfaceMapper->SetInputConnection(surface->GetOutputPort());

    vtkSmartPointer<vtkActor> surfaceActor =
        vtkSmartPointer<vtkActor>::New();
    surfaceActor->SetMapper(surfaceMapper);
    surfaceActor->GetProperty()->SetColor(1.0000, 0.3882, 0.2784);

    // Create graphics stuff
    //
    vtkSmartPointer<vtkRenderer> ren1 =
        vtkSmartPointer<vtkRenderer>::New();
    ren1->SetBackground(.3, .4, .6);

    vtkSmartPointer<vtkRenderWindow> renWin =
        vtkSmartPointer<vtkRenderWindow>::New();
    renWin->AddRenderer(ren1);
    renWin->SetSize(512, 512);

    vtkSmartPointer<vtkRenderWindowInteractor> iren =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    iren->SetRenderWindow(renWin);

    // Add the actors to the renderer, set the background and size
    //
    ren1->AddActor(surfaceActor);

    // Generate an interesting view
    //
    ren1->ResetCamera();
    ren1->GetActiveCamera()->Azimuth(120);
    ren1->GetActiveCamera()->Elevation(30);
    ren1->GetActiveCamera()->Dolly(1.0);
    ren1->ResetCameraClippingRange();

    iren->Initialize();
    iren->Start();

    return EXIT_SUCCESS;
}
//下面这个是为读取不同形式的数据准备的
static vtkSmartPointer<vtkPolyData> ReadPolyData(const char *fileName) {
    vtkSmartPointer<vtkPolyData> polyData;
    std::string extension = vtksys::SystemTools::GetFilenameExtension(std::string(fileName));
    if(extension == ".ply") {
        vtkSmartPointer<vtkPLYReader> reader =
            vtkSmartPointer<vtkPLYReader>::New();
        reader->SetFileName(fileName);
        reader->Update();
        polyData = reader->GetOutput();
    } else if(extension == ".vtp") {
        vtkSmartPointer<vtkXMLPolyDataReader> reader =
            vtkSmartPointer<vtkXMLPolyDataReader>::New();
        reader->SetFileName(fileName);
        reader->Update();
        polyData = reader->GetOutput();
    } else if(extension == ".obj") {
        vtkSmartPointer<vtkOBJReader> reader =
            vtkSmartPointer<vtkOBJReader>::New();
        reader->SetFileName(fileName);
        reader->Update();
        polyData = reader->GetOutput();

    } else if(extension == ".stl") {
        vtkSmartPointer<vtkSTLReader> reader =
            vtkSmartPointer<vtkSTLReader>::New();
        reader->SetFileName(fileName);
        reader->Update();
        polyData = reader->GetOutput();
    } else {
        vtkSmartPointer<vtkPointSource> points =
            vtkSmartPointer<vtkPointSource>::New();
        points->SetNumberOfPoints(100000);
        points->SetRadius(10.0);
        points->SetCenter(vtkMath::Random(-100, 100),
                          vtkMath::Random(-100, 100),
                          vtkMath::Random(-100, 100));
        points->SetDistributionToShell();
        points->Update();
        polyData = points->GetOutput();
    }
    return polyData;
}

在这里插入图片描述

int main() {
    boost::progress_timer pt;
    std::string filename = "../points.txt";
    std::ifstream filestream(filename.c_str()); //文件流
    std::string line;
    vtkNew<vtkPoints> m_Points;
    vtkNew<vtkCellArray> vertices ;
    double x, y, z;
    int i = 0;
    while (std::getline(filestream, line)) { //整行读取文件
        ///double x, y, z;
        std::stringstream linestream;
        linestream << line;
        //cout << line << endl;
        linestream >> x >> y >> z;
        m_Points->InsertNextPoint(x, y, z);  //新读取的数据赋予点的几何结构
        vertices->InsertNextCell(1);     //_加入细胞顶点信息----用于渲染点集
        vertices->InsertCellPoint(i++);
    }
    filestream.close();  //关闭文件流操作
    //好，文件的读取到这就完成了

    vtkNew<vtkPolyData> polyData ;
    polyData->SetPoints(m_Points);获得网格模型中的几何数据：点集
    polyData->SetPolys(vertices);
    //下面就是点云曲面的重建
    double bounds[6];
    polyData->GetBounds(bounds);

    double range[3];
    for (int i = 0; i < 3; ++i) {
        range[i] = bounds[2 * i + 1] - bounds[2 * i];
    }
    int sampleSize = polyData->GetNumberOfPoints() * .00005;//0.00005
    if (sampleSize < 10) {
        sampleSize = 10;//sampleSize =   10
    }
    std::cout << "Sample size is: " << sampleSize << std::endl;
    // 法向量
    vtkSmartPointer<vtkSignedDistance> distance =
        vtkSmartPointer<vtkSignedDistance>::New();
    if (polyData->GetPointData()->GetNormals()) {
        std::cout << "Using normals from input file" << std::endl;
        distance->SetInputData(polyData);
    } else {
        std::cout << "Estimating normals using PCANormalEstimation" << std::endl;
        vtkSmartPointer<vtkPCANormalEstimation> normals =
            vtkSmartPointer<vtkPCANormalEstimation>::New();
        normals->SetInputData(polyData);
        normals->SetSampleSize(sampleSize);
        normals->SetNormalOrientationToGraphTraversal();
        normals->FlipNormalsOn();
        distance->SetInputConnection(normals->GetOutputPort());
    }
    std::cout << "Range: "
              << range[0] << ", "
              << range[1] << ", "
              << range[2] << std::endl;
    int dimension = 300;//看情况作调整
    double radius = range[0] * .02;
    radius = range[0] / static_cast<double>(dimension) * 3; // ~3 voxels
    std::cout << "Radius: " << radius << std::endl;

    distance->SetRadius(radius);
    distance->SetDimensions(dimension, dimension, dimension);
    distance->SetBounds(
        bounds[0] - range[0] * .1,
        bounds[1] + range[0] * .1,
        bounds[2] - range[1] * .1,
        bounds[3] + range[1] * .1,
        bounds[4] - range[2] * .1,
        bounds[5] + range[2] * .1);

    vtkSmartPointer<vtkExtractSurface> surface =
        vtkSmartPointer<vtkExtractSurface>::New();
    surface->SetInputConnection(distance->GetOutputPort());
    surface->SetRadius(radius * 2.0); 参数自己看着调整
    surface->Update();

    vtkSmartPointer<vtkPolyDataMapper> surfaceMapper =
        vtkSmartPointer<vtkPolyDataMapper>::New();
    surfaceMapper->SetInputConnection(surface->GetOutputPort());

    vtkSmartPointer<vtkActor> surfaceActor =
        vtkSmartPointer<vtkActor>::New();
    surfaceActor->SetMapper(surfaceMapper);
    surfaceActor->GetProperty()->SetColor(1.0000, 0.3882, 0.2784);
    //--------------------------------------渲染管线可视化-----------------------------------------------
    vtkSmartPointer<vtkRenderer> ren1 =
        vtkSmartPointer<vtkRenderer>::New();
    ren1->SetBackground(.3, .4, .6);

    vtkSmartPointer<vtkRenderWindow> renWin =
        vtkSmartPointer<vtkRenderWindow>::New();
    renWin->AddRenderer(ren1);
    renWin->SetSize(512, 512);

    vtkSmartPointer<vtkRenderWindowInteractor> iren =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
    iren->SetRenderWindow(renWin);
    renWin->Render();
    renWin->SetWindowName("vtkExtractSurface");
    ren1->AddActor(surfaceActor);
    ren1->ResetCamera();
    ren1->GetActiveCamera()->Azimuth(120);
    ren1->GetActiveCamera()->Elevation(30);
    ren1->GetActiveCamera()->Dolly(1.0);
    ren1->ResetCameraClippingRange();

    iren->Initialize();
    iren->Start();
    cout << "模型渲染完毕" << endl;
    return EXIT_SUCCESS;
}
榴莲小怪兽
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【VTKExamples::Points】第六期 ExtractSurface
qq_40041064的博客
03-23 159
本文分享VTK样例ExtractSurface，并解析接口vtkSignedDistance & vtkExtractSurface，希望对各位小伙伴有所帮助！半径较小时半径较大时该样例演示了如何通过点云数据提取曲面。i < 3;++i)else//第一种参数设置*///第二种参数设置//////
VTK: vtkClipClosedSurface/vtkExtractSurface
qq_40041064的博客
10-28 842
vtkClipClosedSurface
VTK 建模方法：Extrusion
最新发布
ProgramNovice的博客
05-22 896
VTK 建模方法：Extrusion
VTK 获取物体外轮廓 ExtractOutsideSurface
q610098308的专栏
02-19 797
有时我们需要从多个表面体中，提取 最后层哪个，简单的方法 可以使用这种：
VTK 曲面构建+颜色映射
weixin_34216196的博客
04-01 1059
Vtk，（visualization toolkit）是一个开源的免费软件系统，主要用于三维计算机图形学、图像处理和可视化。 VTK编程中主要用到的几个对象 vtkRenderer ，vtkRenderWindow，vtkActor，vtkMapper，其渲染场景如下图所示（图片来自东灵工作室博客） 在VTK的封装类中有一个专门为用...
VTK曲面重建技术(CPR)
Ericohe的博客
04-27 2716
#include <vtkAutoInit.h> VTK_MODULE_INIT(vtkInteractionStyle) VTK_MODULE_INIT(vtkRenderingFreeType) VTK_MODULE_INIT(vtkRenderingOpenGL2) VTK_MODULE_INIT(vtkRenderingVolumeOpenGL2) #include <vtkActor.h> #include <vtkCamera.h> #includ...
VTK-9.1.0-cp38-cp38-win_amd64.zip
11-27
在本例中，`VTK-9.1.0-cp38-cp38-win_amd64.whl` 是VTK库的预编译版本，专为Python 3.8和Windows x64环境设计。 安装这个VTK库的步骤非常简单。首先，确保你的系统已经安装了Python 3.8及pip（Python的包管理器）。...
vtk-v8.0.0.zip
09-08
VTK-8.0.0是该库的一个重要版本，它包含了丰富的数据处理和图形渲染功能。在给定的“vtk-v8.0.0.zip”压缩包中，用户可以找到编译好的VTK库文件，用于与其他项目集成，特别是与PCL（Point Cloud Library）1.8.1版本...
VTK-7.1.1.zip VTK-7.1.1.源码
06-05
VTK-7.1.1.zip 是该库的一个特定版本，包含了VTK 7.1.1的完整源码，使得开发者能够深入理解其内部工作原理，进行定制化开发或优化。 在VTK 7.1.1这个版本中，我们可以期待一系列关键特性和改进。首先，VTK是一个C++...
VTK-8.2.0-cp36-cp36m-win_amd64.zip
11-27
VTK-8.2.0-cp36-cp36m-win_amd64.zip 是一个专为基于Windows x64平台设计的Python 3.6版本的VTK库压缩包。它包含了Python接口的VTK库，允许开发者在Python环境中进行高效的数据可视化。 此压缩包中的“VTK-8.2.0-cp...
VTK-8.1.0.zip
02-02
VTK-8.1.0.zip 是一个包含了VTK版本8.1.0的压缩包，提供了完整的源代码、文档和其他相关资源，允许用户在自己的项目中集成和自定义VTK的功能。 在VTK 8.1.0中，有几个关键的知识点值得深入探讨： 1. **模块化设计*...
基于VTK的面绘制
10-09
很好的基于VTK开发库的医学图像三维重建，这里用了VTK里的面绘制方法就行了三维重建，对初学者很有益，望交流学习！
VTK体中的曲面提取(其中包含曲线绘制算法)
yuxing55555的博客
05-18 3422
在体中提取曲面 1.绘制曲线 demo：https://lorensen.github.io/VTKExamples/site/Cxx/PolyData/ParametricSpline/ // Create a vtkPoints object and store the points in it vtkSmartPointer&amp;amp;lt;vtkPoints&amp;amp;gt; points ...
VTK 实现 曲面重建（CPR）
Beyond欣
05-25 8118
VTK 实现 曲面重建（CPR） 原理 CPR±+Curved+Planar+Reformation 下载: https://www.jianguoyun.com/p/DXY5HOoQgYiuCRia-PED (访问密码 : r0ga0d) 曲面重建CPR步骤 输入点拟合样条曲线（利用vtkSplineFilter） 计算样条曲线各点的法向量(弗莱纳公式) 根据点和法向量截取图片(vtkProbeFilter) 每张图片拼接起来(vtkImageAppend) VTK实现 常用功能vtk肯定提供现成办法
VTK初学一，vtkDelaunay2D创建球冠曲面
weixin_30306905的博客
12-03 177
#ifndef INITIAL_OPENGL #define INITIAL_OPENGL #include <vtkAutoInit.h> VTK_MODULE_INIT(vtkRenderingOpenGL) VTK_MODULE_INIT(vtkInteractionStyle) #endif #include <iostream> using na...
利用VTK手动切割体数据，曲面切割
我才睡醒的博客
01-16 8032
vtk手动切割CT成像的体数据 支持Windows系统和Linux系统，前提是有VTK环境
vtk三维点云隐式曲面重建：vtkSurfaceConstructionFilter类
HW140701的博客
10-12 8937
vtk三维点云重建vtkSurfaceConstructionFilter vtkSurfaceReconstructionFilter获取假设为位于实体3D对象表面上的点列表。 计算与表面距离的带符号度量，并在常规网格上采样。 然后可以将网格轮廓化为零以提取表面，使用了三维点云隐式曲面重构技术。 #include "vtkRenderer.h" #include "vtkRender...
VTK 三角剖分 Delaunay2D（五）—— 折线裁剪封闭曲面
GIS探索之路
03-18 2697
引言 vtkClipClosedSurface使用一组平面集合剪切闭合曲面。 输入vtkPolyData闭合曲面和vtkPlaneCollection一组剪裁平面。 输出vtkPolyData新多边形闭合曲面。 描述 通过画折线的方式来裁剪封闭的多边形曲面。 效果如下： 实现 模拟交互式生成折线 //裁剪多段线 double pt0[3] = { -2.0, 2.0, 2.0}...
VTK-3 三维点云曲面重建
C++ 爱好者 make_it_simple888
12-02 4873
表面重建技术（三维点云曲面重建） https://blog.csdn.net/shenziheng1/article/details/54880915 多分辨率策略（模型细化的三种方法） https://blog.csdn.net/shenziheng1/article/details/54861856 VTK 点云重建 https://blog.csdn.net/weixin_3850011...
VTK-m用户指南：构建高性能可视化
"VTKmUsersGuide.pdf 是一个关于VTK-m的用户指南，版本1.2.0，由Kenneth Moreland等人编写，并由Kitware Inc.出版。该指南介绍了VTK-m，一个用于科学可视化的现代并行库，特别强调了在高性能计算中的应用。VTK-m旨在...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值