轉(zhuǎn)載于http://apps.hi.baidu.com/share/detail/7831982?感謝原作者！ 1.用vtkDICOMImageReader #include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h" #include "vtkImageActor.h"
#include "vtkDICOMImageReader.h"
#include "vtkImageShiftScale.h"
int main ()
{
???vtkDICOMImageReader *dcmReader = vtkDICOMImageReader::New();
???dcmReader->SetDirectoryName("E:\medicalimages\vhm");
???dcmReader->SetDataByteOrderToLittleEndian();
?double x,y,z;
?dcmReader->GetDataSpacing(x,y,z);
?dcmReader->SetDataSpacing(x,y,z); vtkImageShiftScale *shifter =vtkImageShiftScale::New();
?shifter->SetInputConnection(dcmReader->GetOutputPort());
?shifter->SetShift(70);
?shifter->SetScale(0.5);
?shifter->SetOutputScalarTypeToUnsignedChar(); vtkImageActor *actor =vtkImageActor::New();
?actor->SetInput(shifter->GetOutput()); vtkRenderer *aRender = vtkRenderer::New();//設(shè)置繪制類
?aRender->AddActor(actor); vtkRenderWindow*renWin=vtkRenderWindow::New();//設(shè)置繪制窗口
?renWin->SetSize(500, 500);//設(shè)置背景顏色和繪制窗口大小
?renWin->AddRenderer(aRender);//裝載繪制類 vtkRenderWindowInteractor*iRen=vtkRenderWindowInteractor::New();//設(shè)置繪制窗口的交互
?iRen->SetRenderWindow(renWin);//裝載繪制窗口
??renWin->Render();//窗口進(jìn)行繪制
?iRen->Initialize();?
?iRen->Start(); //初始化并進(jìn)行交互繪制 dcmReader->Delete();
?shifter->Delete();
?actor->Delete();
?aRender->Delete();
???renWin->Delete();
?iRen->Delete();?
????return 0;
}?
2.用vtkVolume16Reader,不限定DCM文件,不過文件后綴要以.1 .2.3...來命名,工作有得做,而且不能讀太多圖片,vtkArray承載不了太多數(shù)據(jù).而且圖片會錯位
????????????vtkVolume16Readerv16 = new vtkVolume16Reader();
????????????v16.SetDataDimensions(512,512);
????????????v16.SetDataByteOrderToLittleEndian();
????????????v16.SetFilePrefix("F:\LiverGUI\SE4\IM");
????????????v16.SetImageRange(1,100);
????????????v16.SetDataSpacing(1.0,1.0, 1.0);?

3.讀raw data...用vtkImageData,讀取用vtkImageReader更快!具體查vtk手冊吧!
??????private void volumeRendering(String filePath, int[] dims, int[]shrinkFactor,vtkRenderWindow renWin)
????????{
????????????//接口說明
????????????//1.filePath:dat文件名
????????????//2.dims:體數(shù)據(jù)的三個維度所構(gòu)成的數(shù)組
????????????//3.shrinkFactor:三個元素的數(shù)組,代表長寬 高 的縮放比例
????????????//4.vtkRenderWindow,窗口中添加的vtkFormsWindowControl通過GetRenderWindow方法獲得的變量
????????????//Readdat data
????????????FileStreamfs = null;
????????????byte[]point = new byte[2];
????????????if(!File.Exists(filePath))
????????????{
????????????????MessageBox.Show("RawData doesn't exist!");
????????????????return;
????????????}
????????????fs= new FileStream(filePath, FileMode.Open, FileAccess.Read,FileShare.Read);
????????????inti, j, k;
????????????short[,,]rawData= new short[dims[0], dims[1], dims[2]];
????????????for(i = 0; i < dims[0]; i++)
????????????{
????????????????for(j = 0; j < dims[1]; j++)
????????????????{
????????????????????for(k = 0; k < dims[2]; k++)
????????????????????{
????????????????????????point[0]= (byte)fs.ReadByte();
????????????????????????point[1]= (byte)fs.ReadByte();
????????????????????????rawData= BitConverter.ToInt16(point, 0);
????????????????????}
????????????????}
????????????}
????????????//VolumeRendeing
????????????vtkImageDataid = new vtkImageData();
????????????id.SetDimensions(dims[0],dims[1], dims[2]);
????????????id.SetScalarTypeToShort();
????????????id.SetNumberOfScalarComponents(1);
????????????id.AllocateScalars();
????????????vtkImageShrink3Dmask = new vtkImageShrink3D();
????????????mask.SetInput(id);
????????????mask.SetShrinkFactors(shrinkFactor[0],shrinkFactor[1], shrinkFactor[2]);
????????????unsafe
????????????{
????????????????short*a = (short*)(id.GetScalarPointer().ToPointer());
????????????????for(i = 0; i < dims[0]; i++)
????????????????{
????????????????????for(j = 0; j < dims[1]; j++)
????????????????????{
????????????????????????for(k = 0; k < dims[2]; k++)
????????????????????????{
????????????????????????????*a++= rawData;
????????????????????????}
????????????????????}
????????????????}
????????????}
????????????vtkContourFiltercf = new vtkContourFilter();
????????????cf.SetInputConnection(mask.GetOutputPort());
????????????vtkPolyDataMappermapper = new vtkPolyDataMapper();
????????????mapper.SetInputConnection(cf.GetOutputPort());
????????????mapper.SetScalarRange(-2048,2047);
????????????vtkActoractor = new vtkActor();
????????????actor.SetMapper(mapper);
????????????vtkCameraaCamera = new vtkCamera();
????????????aCamera.SetViewUp(0,5, 0);
????????????aCamera.SetPosition(5,0, 0);
????????????aCamera.SetFocalPoint(0,0, 0);
????????????aCamera.ComputeViewPlaneNormal();
????????????vtkRendererren = new vtkRenderer();
????????????ren.AddActor(actor);
????????????ren.ResetCamera();
????????????aCamera.Dolly(1.5);
????????????ren.ResetCameraClippingRange();
????????????renWin.AddRenderer(ren);
????????????ren.SetBackground(0.0f,0.0f, 0.0f);
????????}
????}

總結(jié)

以上是生活随笔為你收集整理的vtk读取文件并显示的几种方法的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。