VTK读取序列的Dicom医学图片，用Marchingcube进行重建，并保存为obj文件

最新推荐文章于 2022-10-07 13:25:37 发布
最新推荐文章于 2022-10-07 13:25:37 发布
阅读量9.3k 收藏 28
点赞数 3
分类专栏： VTK
本文为博主原创文章，未经博主允许不得转载，如需转载请先得到博主的同意，如需疑问，请联系[email protected]，也可以加入计算机图形图像群526867211，以及访问我的个人站点：www.stubbornhuang.com，谢谢。
本文链接：https://blog.csdn.net/hw140701/article/details/51682645
版权 
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkDICOMImageReader.h"
#include "vtkMarchingCubes.h"
#include "vtkStripper.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkSmartPointer.h"
#include "vtkProperty.h"
#include "vtkCamera.h"
#include "vtkOutlineFilter.h"
#include "vtkOBJExporter.h"

void main()
{
	vtkSmartPointer<vtkRenderer>ren=vtkSmartPointer<vtkRenderer>::New();             //设置绘制者(绘制对象指针)
	vtkSmartPointer<vtkRenderWindow>renWin=vtkSmartPointer<vtkRenderWindow>::New();         //设置绘制窗口
	renWin->AddRenderer(ren);                                 //将绘制者加入绘制窗口

	vtkSmartPointer<vtkRenderWindowInteractor>iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();//设置绘制交互操作窗口的
	iren->SetRenderWindow(renWin);                             //将绘制窗口添加到交互窗口

	vtkSmartPointer<vtkDICOMImageReader>Reader=vtkSmartPointer<vtkDICOMImageReader>::New(); //创建读取dicom图片指针对象
	Reader->SetDirectoryName("E:\\dcm\\dcm"); //设置医学图像文件夹路径
	Reader->SetDataByteOrderToLittleEndian();
	Reader->Update();
	cout<<"读取数据完成"<<endl;


	//
	vtkSmartPointer<vtkMarchingCubes>marchingcube=vtkSmartPointer<vtkMarchingCubes>::New();       //建立一个Marching Cubes 算法的对象	
	//marchingcube->SetInput((vtkDataSet *)Reader->GetOutput());      //获得所读取的数据
	marchingcube->SetInputConnection(Reader->GetOutputPort());      //第二种读取数据的方法
	marchingcube->SetValue(0,140);                                  //提取出灰度值为45的东西

	vtkSmartPointer<vtkStripper>Stripper=vtkSmartPointer<vtkStripper>::New();                   //建立三角带对象
	Stripper->SetInput( marchingcube->GetOutput() );              //将生成的三角片连接成三角带

	vtkSmartPointer<vtkPolyDataMapper>Mapper=vtkSmartPointer<vtkPolyDataMapper>::New();   //建立一个数据映射对象
	Mapper->SetInput(Stripper->GetOutput() );            //将三角带映射为几何数据
	Mapper->ScalarVisibilityOff();

	vtkSmartPointer<vtkActor>actor=vtkSmartPointer<vtkActor>::New();                            //建立一个代表皮肤的演员
	actor->SetMapper(Mapper);                                  //获得皮肤几何数据的属性
	actor->GetProperty()->SetDiffuseColor(1, .49, .25);            //设置皮肤颜色的属性//(1, .49, .25)
	//actor->GetProperty()->SetDiffuseColor(1, 1, .25);
	actor->GetProperty()->SetSpecular(0.3);                         //设置反射率
	actor->GetProperty()->SetSpecularPower(20);                    //设置反射光强度 
	actor->GetProperty()->SetOpacity(1.0);                 
	actor->GetProperty()->SetColor(1,0,0);                 //设置角色的颜色
	actor->GetProperty()->SetRepresentationToWireframe();


	vtkSmartPointer<vtkCamera>aCamera=vtkSmartPointer<vtkCamera>::New();     // 创建摄像机
	aCamera->SetViewUp ( 0, 0, -1 );        //设置相机的“上”方向
	aCamera->SetPosition ( 0, 1, 0 );       //位置：世界坐标系，设置相机位置
	aCamera->SetFocalPoint( 0, 0, 0 );     //焦点：世界坐标系，控制相机方向
	aCamera->ComputeViewPlaneNormal();     //重置视平面法向，基于当前的位置和焦点。否则会出现斜推剪切效果

	///
	vtkSmartPointer<vtkOutlineFilter>outlinefilter=vtkSmartPointer<vtkOutlineFilter>::New();
	outlinefilter->SetInputConnection(Reader->GetOutputPort());

	vtkSmartPointer<vtkPolyDataMapper>outlineMapper=vtkSmartPointer<vtkPolyDataMapper>::New();
	outlineMapper->SetInputConnection(outlinefilter->GetOutputPort());

	vtkSmartPointer<vtkActor>OutlineActor=vtkSmartPointer<vtkActor>::New();
	OutlineActor->SetMapper(outlineMapper);
	OutlineActor->GetProperty()->SetColor(0,0,0);
	// 告诉绘制者将要在绘制窗口中进行显示的演员
	ren->AddActor(actor);
	ren->AddActor(OutlineActor);
	
	ren->SetActiveCamera(aCamera);               设置渲染器的相机
	ren->ResetCamera();
	aCamera->Dolly(1.5);                //使用Dolly()方法沿着视平面法向移动相机，实现放大或缩小可见角色物体
	ren->SetBackground(1,1,1);               //设置背景颜色
	//ren->ResetCameraClippingRange();
	renWin->SetSize(1000, 600);
	renWin->Render();
	iren->Initialize();
	iren->Start(); 
	vtkSmartPointer<vtkOBJExporter> porter=vtkSmartPointer<vtkOBJExporter>::New();
	porter->SetFilePrefix("E:\\PolyDataWriter.obj");
	porter->SetInput(renWin);
	porter->Write();

}

 

以上代码是VTK5.10版本所适用的，有一些小伙伴发现使用高版本如VTK6.0，VTK7.0的时候，数据传输会发现错误，是因为高版本的改变了数据的传输方式，需要改变数据连接方式，如将

 

Mapper->SetInput(Stripper->GetOutput() );            //将三角带映射为几何数据

改成

 

 

Mapper->SetInputConnection(Stripper->GetOutputPort() );            //将三角带映射为几何数据


就可以了，以下附上高版本的修改代码

 

 

#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkDICOMImageReader.h"
#include "vtkMarchingCubes.h"
#include "vtkStripper.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkSmartPointer.h"
#include "vtkProperty.h"
#include "vtkCamera.h"
#include "vtkOutlineFilter.h"
#include "vtkOBJExporter.h"

void main()
{
	vtkSmartPointer<vtkRenderer>ren=vtkSmartPointer<vtkRenderer>::New();             //设置绘制者(绘制对象指针)
	vtkSmartPointer<vtkRenderWindow>renWin=vtkSmartPointer<vtkRenderWindow>::New();         //设置绘制窗口
	renWin->AddRenderer(ren);                                 //将绘制者加入绘制窗口

	vtkSmartPointer<vtkRenderWindowInteractor>iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();//设置绘制交互操作窗口的
	iren->SetRenderWindow(renWin);                             //将绘制窗口添加到交互窗口

	vtkSmartPointer<vtkDICOMImageReader>Reader=vtkSmartPointer<vtkDICOMImageReader>::New(); //创建读取dicom图片指针对象
	Reader->SetDirectoryName("C:\\Users\\Administrator\\Desktop\\医学CT图片\\COU IV"); //设置医学图像文件夹路径
	Reader->SetDataByteOrderToLittleEndian();
	Reader->Update();
	cout<<"读取数据完成"<<endl;


	//
	vtkSmartPointer<vtkMarchingCubes>marchingcube=vtkSmartPointer<vtkMarchingCubes>::New();       //建立一个Marching Cubes 算法的对象	
	//marchingcube->SetInput((vtkDataSet *)Reader->GetOutput());      //获得所读取的数据
	marchingcube->SetInputConnection(Reader->GetOutputPort());      //第二种读取数据的方法
	marchingcube->SetValue(0,140);                                  //提取出灰度值为45的东西

	vtkSmartPointer<vtkStripper>Stripper=vtkSmartPointer<vtkStripper>::New();                   //建立三角带对象
	Stripper->SetInputConnection( marchingcube->GetOutputPort());              //将生成的三角片连接成三角带

	vtkSmartPointer<vtkPolyDataMapper>Mapper=vtkSmartPointer<vtkPolyDataMapper>::New();   //建立一个数据映射对象
	Mapper->SetInputConnection(Stripper->GetOutputPort() );            //将三角带映射为几何数据
	Mapper->ScalarVisibilityOff();

	vtkSmartPointer<vtkActor>actor=vtkSmartPointer<vtkActor>::New();                            //建立一个代表皮肤的演员
	actor->SetMapper(Mapper);                                  //获得皮肤几何数据的属性
	actor->GetProperty()->SetDiffuseColor(1, .49, .25);            //设置皮肤颜色的属性//(1, .49, .25)
	//actor->GetProperty()->SetDiffuseColor(1, 1, .25);
	actor->GetProperty()->SetSpecular(0.3);                         //设置反射率
	actor->GetProperty()->SetSpecularPower(20);                    //设置反射光强度 
	actor->GetProperty()->SetOpacity(1.0);                 
	actor->GetProperty()->SetColor(1,0,0);                 //设置角色的颜色
	actor->GetProperty()->SetRepresentationToWireframe();


	vtkSmartPointer<vtkCamera>aCamera=vtkSmartPointer<vtkCamera>::New();     // 创建摄像机
	aCamera->SetViewUp ( 0, 0, -1 );        //设置相机的“上”方向
	aCamera->SetPosition ( 0, 1, 0 );       //位置：世界坐标系，设置相机位置
	aCamera->SetFocalPoint( 0, 0, 0 );     //焦点：世界坐标系，控制相机方向
	aCamera->ComputeViewPlaneNormal();     //重置视平面法向，基于当前的位置和焦点。否则会出现斜推剪切效果

	///
	vtkSmartPointer<vtkOutlineFilter>outlinefilter=vtkSmartPointer<vtkOutlineFilter>::New();
	outlinefilter->SetInputConnection(Reader->GetOutputPort());

	vtkSmartPointer<vtkPolyDataMapper>outlineMapper=vtkSmartPointer<vtkPolyDataMapper>::New();
	outlineMapper->SetInputConnection(outlinefilter->GetOutputPort());

	vtkSmartPointer<vtkActor>OutlineActor=vtkSmartPointer<vtkActor>::New();
	OutlineActor->SetMapper(outlineMapper);
	OutlineActor->GetProperty()->SetColor(0,0,0);
	// 告诉绘制者将要在绘制窗口中进行显示的演员
	ren->AddActor(actor);
	ren->AddActor(OutlineActor);
	
	ren->SetActiveCamera(aCamera);               设置渲染器的相机
	ren->ResetCamera();
	aCamera->Dolly(1.5);                //使用Dolly()方法沿着视平面法向移动相机，实现放大或缩小可见角色物体
	ren->SetBackground(1,1,1);               //设置背景颜色
	//ren->ResetCameraClippingRange();
	renWin->SetSize(1000, 600);
	renWin->Render();
	iren->Initialize();
	iren->Start(); 
	vtkSmartPointer<vtkOBJExporter> porter=vtkSmartPointer<vtkOBJExporter>::New();
	porter->SetFilePrefix("E:\\PolyDataWriter.obj");
	porter->SetInput(renWin);
	porter->Write();

}


 

如果您觉得这篇博文有用，请访问我的个人站：http://www.stubbornhuang.com，更多博文干货等着您。

 


 

HW140701
关注
  • 3
    点赞
  • 踩
  • 28
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 64
    评论
专栏目录
一个命令行工具，转换DICOM卷到3d表面网格(obj, stl或ply)。几个网格处理程序可以启用，如网格减少，平滑或清洁。适用于Linux, OSX和Windows。- AOT-AG / DicomToMesh
01-27
一个命令行工具，转换DICOM卷到3d表面网格(obj, stl或ply)。几个网格处理程序可以启用，如网格减少，平滑或清洁。适用于Linux, OSX和Windows。- AOT-AG / DicomToMesh-源码
Marching Cubes算法在医学三维重建中的应用
qq_33598781的博客
06-01 2109
概述 　　之前的博文已经完整的介绍了三维图像数据和三角形网格数据。在实际应用中，利用遥感硬件或者各种探测仪器，可以获得表征现实世界中物体的三维图像。比如利用CT机扫描人体得到人体断层扫描图像，就是一个表征人体内部组织器官形状的一个三维图像。其中的感兴趣的组织器官通过体素的颜色和背景加以区别。如下图的人体足骨扫描图像。医生通过观察这样的图像可以分析病人足骨的特征，从而对症下药。 　　这类应用在计算机领域叫做科学可视化。由于本文主要不是讨论可视化这个大的命题，所以只是简要的讲述一下三维可视化的两大类实现
VS2008、QT及VTK实现DICOM图像三维重建之二：VS2008下VTK的安装使用
flywithmj的专栏
01-05 3139
网上有很多VS2008下VTK安装配置的文章，这里也就不再赘述，主要提醒大家注意一下几点： 1，选用VTK版本和CMake版本时，一定要注意CMake是否支持当前版本VTK及两者是否均支持VS2008。在我的工程中，我使用VTK-5.6.1，CMake版本时当时最新的2.8.6版本。Visual studio版本无疑就是VS2008了。   2，在用Cmake生成VS2008工程时，首先选择
DICOM三维重建：基于VTK实现三种面绘制算法
一直在笑的博客
10-21 4797
版本：VS2019、VTK9.0.1 #define vtkRenderingCore_AUTOINIT 3(vtkInteractionStyle,vtkRenderingFreeType,vtkRenderingOpenGL2) #define vtkRenderingVolume_AUTOINIT 1(vtkRenderingVolumeOpenGL2) #include "vtkDICOMImageReader.h" #include "vtkRenderWindowInteractor.h" #i
VTK序列图像的读取
东灵工作室
05-07 1万+
医学图像处理的应用程序中，经常会碰到读取一个序列图像的操作。比如CT、MR等所成的图像都是一个切面一个切面地存储的，医学图像处理程序要处理这些数据，第一步当然是把这些数据从磁盘等外部存储介质中导入内存。         利用VTK可以读取多种格式的图像文件，支持读取单个的二维图像(比如*.BMP、*.JPEG、*.PNG等)或者三维图像文件(*.VTK、*.mhd、*.mha等)，也支持序列图像
itk 读取dicom，用vtk重建显示
simple的博客
12-30 696
using PixelType = float; constexpr unsigned int Dimension = 3; using ImageType = itk::Image< PixelType, Dimension >; using ReaderType = itk::ImageSeriesReader< ImageType >; using ImageIOType = itk::GDCMImageIO; using NamesGeneratorType = .
VTK实现三维模型的导出保存，STL、OBJ和PLY等格式
qq_41060452的博客
05-14 8734
基于之前的博客中实现的三维模型，将其导出保存为3D格式的文件，生成的文件可以使用通用的3D浏览器进行查看。附录的代码实现中各变量的命名使用可以查照之前的博客。 1.STL 保存 STL 文件格式，最早由美国的一家公司提出，是计算机实现三维数据打印的 一种通用格式文件。由于该格式文件存储方式简单，具有规范的接口格式，是现 阶段快速成型系统中被应用最多的一种标准文件格式，可以说市面上对其的应用 是十分广泛的，尤其是在医学影像和文物保护领域。STL 文件是通过很多小三角 形面片进行堆积来实现三维模型表面的构造的一
QT C++ vtk 加载二维dicom 三维渲染 3DMPR
10-07 1688
vtk 容积协议 效果 ；QT VTK dicomtk 3DMPR 等
ITK、VTK 读取 DiCOM 文件信息，并实现交互预览
小张Python
07-20 4883
DICOM 文件介绍 DICOM 全名叫 Digital Imaging and Communications in Medicine(医学数字图像与通信)，是医学图像和相关信息的过节标准，也定义了临床所需数据医学图像格式； DICOM 数据格式包含大量属性，例如病人名字、ID、年龄等基本信息，这些信息从储存在 “头文件” 中，这种储存方式避免了病人与图像误配； 一个病人样本 包含多个 DICOM 子文件，单个 DICOM 对象只包含具有一帧像素数据，帧与帧之间以序列形式存在，可以来回切换形成类似于电影放映
【医学影像数据处理】nii、npz、npy、dcm、mhd 的数据格式互转，及多目标分割处理汇总
热门推荐
钱多多先森
04-22 2万+
医学领域进行数据处理的DICOM、.nii、.npy、.mha、.npz，jpeg、png等格式。
VTK.js dicom图像三维重建
04-19
VTK.js 实现网页版dicom图像三维重建
读取并可视化DICOM图像（VTK & OpenCV） - louishao的博客 - CSDN博客1
08-04
本人处理的CT图的格式是DICOM格式，这个用OpenCV不能直接读取，我采用了VTK进行读取，然后用OpenCV处理的思路（当然网上也有人说可以直接用ITK读
在vtk中保存为obj文件格式
09-17
在vtk中将三维模型的文件格式保存为obj文件格式
Python vtk读取并显示dicom文件示例
01-20
因为做项目的原因，所以接触到了医学图像dicom文件。vtk刚开始看，这里仅仅只是其最简单的读取显示功能。此处用到了vtk库，可自行百度安装方法。 下面附上代码： from vtk import * # reader the dicom file reader...
VTK读取30张Dicom图片三维重建
04-20
利用VTK三维重建软件与VS2010结合，读取格式为Dicom的系列图像，并能显示三维重建结果，语言为C++语言(Using VTK 3D reconstruction software and VS2010 combination, read the format of Dicom series of images, ...
VTK读取30张Dicom图片三维重建.rar_C ，DICOM_dicom重建软件_三维重建 图片_图像 三维重建_图片 三
07-15
利用VTK三维重建软件与VS2010结合，读取格式为Dicom的系列图像，并能显示三维重建结果，语言为C++语言
C++VTK+MFC读取医学dicom图像及GPU三维重建(VTK MFC GPU)
04-19
VTK+MFC读取医学dicom图像及GPU三维重建(VTK MFC GPU)
三维模型重建（1）：关于三维模型重建的一些简介
HW140701的博客
07-07 1万+
三维模型重建顾名思义就是从一些数据（点云，图片，三维轮廓线等）重建出物体的三维逼真的三维模型，在其重建的过程中针对不同的数据的的三维重建会有不同的处理算法，如针对点云数据的三维重建有很多种重建方法,如基于Delaunay三角化，Voronoi图，隐式曲面等方法，另外在三维模型的重建过程（MarchingCube，RayCast,网格构建等）以及三维模型生成后的后处理算法（三维网格简化，三维网格加密...
VTK读取一个TXT文档中的三维点坐标显示三维点云
HW140701的博客
02-20 1万+
VTK读取一个TXT文档中的三维点坐标就可以显示三维点云，txt文档中的格式为 X坐标 Y坐标 Z坐标 如下所示： 附上代码如下： #include <iostream> #include <vector> #include "vtkActor.h" #include "vtkRenderer.h" #include "vtkRenderWindo...
c语言 二维数组和指针数组实现十二个月份字典序排序
最新发布
03-28
#include <stdio.h> #include <string.h> int main() { char months[12][10] = {"January", "February", "March", "April", "May", "June", "July", "August", "September", "October", "November", "December"}; char *ptr[12]; int i, j; // 将指针数组指向二维数组中的每个字符串 for (i = 0; i < 12; i++) { ptr[i] = months[i]; } // 冒泡排序 for (i = 0; i < 11; i++) { for (j = 0; j < 11 - i; j++) { if (strcmp(ptr[j], ptr[j+1]) > 0) { char *temp = ptr[j]; ptr[j] = ptr[j+1]; ptr[j+1] = temp; } } } // 输出排序后的月份 for (i = 0; i < 12; i++) { printf("%s\n", ptr[i]); } return 0; }

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 64
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

HW140701

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值