VTK 16位图像体绘

#include “vtkRenderer.h”
#include “vtkRenderWindow.h”
#include “vtkRenderWindowInteractor.h”
#include “vtkPiecewiseFunction.h”
#include “vtkColorTransferFunction.h”
#include “vtkVolumeProperty.h”
#include “vtkVolumeRayCastIsosurfaceFunction.h”
#include “vtkVolumeRayCastCompositeFunction.h”
#include “vtkVolumeRayCastMapper.h”
#include “vtkVolume.h”
#include “vtkImageCast.h”
#include “vtkVolume16Reader.h”

int main()
{
vtkRenderer *aRender = vtkRenderer::New();
vtkRenderWindow *renWin = vtkRenderWindow::New();
renWin->AddRenderer(aRender);
vtkRenderWindowInteractor *iRen = vtkRenderWindowInteractor::New();
iRen->SetRenderWindow(renWin);

vtkVolume16Reader *reader = vtkVolume16Reader::New();
//初始化16位图像类型,前缀等
reader->SetDataDimensions (256,256);//设数据范围大小
reader->SetImageRange (1,30);//读取图像的文件名 .1–.123
reader->SetDataByteOrderToLittleEndian();//字节格式(字节顺序),如果是littleEndian要倒过来
reader->SetFilePrefix (“E://CT”);//设置文件目录
reader->SetDataSpacing (1, 1, 6);//像素间的间隔

vtkImageCast *readerImageCast = vtkImageCast::New();
	readerImageCast->SetInput((vtkDataObject *)reader->GetOutput());
	readerImageCast->SetOutputScalarTypeToUnsignedShort();

// readerImageCast->ClampOverflowOn();

vtkPiecewiseFunction *opacityTransferFunction = vtkPiecewiseFunction::New();
	opacityTransferFunction->AddPoint(200, 0.0);
    opacityTransferFunction->AddPoint(590, 0.0);
	opacityTransferFunction->AddPoint(1200, 0.0);
    opacityTransferFunction->AddPoint(1300, 0.0);
    opacityTransferFunction->AddPoint(2000, 0.4);
    opacityTransferFunction->AddPoint(4096, 1);
	opacityTransferFunction->ClampingOn();	
	
vtkColorTransferFunction *colorTransferFunction = vtkColorTransferFunction::New();
  colorTransferFunction->AddRGBPoint(0.0, 0.5, 0.0, 0.0);
  colorTransferFunction->AddRGBPoint(500.0, 1.0, 0.5, 0.23);//皮肤
  colorTransferFunction->AddRGBPoint(1280, 1, 1, 0.96);//骨骼
  colorTransferFunction->AddRGBPoint(2550.0, 1, 1, 1);

vtkVolumeProperty *volumeProperty = vtkVolumeProperty::New();
  volumeProperty->SetColor(colorTransferFunction);
  volumeProperty->SetScalarOpacity(opacityTransferFunction);
  volumeProperty->ShadeOn();
  volumeProperty->SetInterpolationTypeToLinear();
  volumeProperty->SetAmbient(0.2);
  volumeProperty->SetDiffuse(0.9);
  volumeProperty->SetSpecular(0.2);
  volumeProperty->SetSpecularPower(10);	

vtkVolumeRayCastCompositeFunction *compositeFunction = vtkVolumeRayCastCompositeFunction::New();
vtkVolumeRayCastIsosurfaceFunction *RayIso=vtkVolumeRayCastIsosurfaceFunction::New();
  RayIso->SetIsoValue  (12000) ;
  
vtkVolumeRayCastMapper *volumeMapper = vtkVolumeRayCastMapper::New();
  volumeMapper->SetVolumeRayCastFunction(RayIso);
  volumeMapper->SetInputConnection(reader->GetOutputPort());


vtkVolume *volume = vtkVolume::New();
	  volume->SetMapper(volumeMapper);
	  volume->SetProperty(volumeProperty);
	  volume->RotateX(130);

// volume->RotateY(35);
// volume->RotateZ(-35);

aRender->AddVolume(volume);
aRender->SetBackground(1, 1, 1);

renWin->SetSize(600, 600);
renWin->Render();

iRen->Initialize();
renWin->Render();
iRen->Start();

aRender->Delete();
renWin->Delete();
iRen->Delete();
reader->Delete();
opacityTransferFunction->Delete();
colorTransferFunction->Delete();
volumeMapper->Delete();
volumeProperty->Delete();
compositeFunction->Delete();
volume->Delete();

return 0;

}