VTK读取序列的Dicom医学图片，用Marchingcube进行重建，并保存为obj文件
本文作者：StubbornHuang
版权声明：本文为站长原创文章，如果转载请注明原文链接！
原文标题：VTK读取序列的Dicom医学图片，用Marchingcube进行重建，并保存为obj文件
原文链接：https://www.stubbornhuang.com/200/
发布于：2019年11月4日 23:14
修改于：2020年1月3日 09:31
1 代码版本
代码针对于VTK版本5.10，新版本仅供参考。
2 VTK5.10版本示例代码
#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();
}
3 VTK6.0以上版本代码
第三节中代码是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();
}
4 结果示例
如果您觉得对您有帮助，可以请站长喝一杯咖啡哦！
记得在赞赏备注里写上您的昵称
您可在本站资助名单中查看你的打赏记录哦！
 支付宝扫一扫 |
 微信扫一扫 |
金额随意，礼清义重
当前分类随机文章推荐
- VTK读取一个TXT文档中的三维点坐标显示三维点云
- VTK读取序列的Dicom医学图片，用Marchingcube进行重建，并保存为obj文件
- VTK能干什么？VTK大部分功能的细节简介，VTK能打开的文件格式
- VTK以批量三维点坐标为中心（点云）绘制球体，可用于标识特征点或者是化学分子
- VTK - 冠脉重建点匹配坐标数据下载
全站随机文章推荐
- C++11 - std::bind简要介绍以及可绑定函数的几种形式总结
- Python - 运算符/ or // or %的含义和区别
- Duilib - 字体
- 资源分享 - Computer Graphics Principles and Practice (3rd edition) 英文高清PDF下载
- 资源分享 - The Magic of Computer Graphics - Landmarks in Rendering 英文高清PDF下载
- WordPress - 增加百度统计代码
- 资源分享 - GPU Pro 360 - Guide to Rendering 英文高清PDF下载
- FFmpeg - RGB图像编码为h264出现垂直旋转的问题
- 资源分享 - Python深度学习：基于PyTorch (吴茂贵著) 高清PDF下载
- Python3爬虫 - 下载反盗链图片的方式
纵使晴明无雨色，入云深处亦沾衣。
牛逼！！
谢谢你