bmpvolumerendering.cpp

vtk中一些类的具体定义方法

#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkPiecewiseFunction.h"
#include "vtkColorTransferFunction.h"
#include "vtkVolumeProperty.h"
#include "vtkVolumeRayCastCompositeFunction.h"
#include "vtkVolumeRayCastMapper.h"
#include "vtkVolume.h"
#include "vtkVolume16Reader.h"
#include "vtkCamera.h"
#include "vtkDICOMImageReader.h"
#include "vtkImageCast.h"
#include "vtkTIFFReader.h"


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

		vtkTIFFReader *tifReader = vtkTIFFReader::New();
		tifReader->SetDataByteOrderToLittleEndian();
		tifReader->SetFilePrefix("E://Data/MRI/HEAD");
		tifReader->SetDataExtent(0, 128, 1,128, 0, 53);
		tifReader->SetDataOrigin(0, 0, 0);
		tifReader->SetDataSpacing(3.2,3.2,4.3);
	

		vtkPiecewiseFunction *opacityTransferFunction = vtkPiecewiseFunction::New();
		opacityTransferFunction->AddPoint(0, 0.0);
		opacityTransferFunction->AddPoint(600, 0.3);
		opacityTransferFunction->AddPoint(700, 0.78);
		opacityTransferFunction->AddPoint(800, 0.89);
		opacityTransferFunction->AddPoint(900, 0.95);
		opacityTransferFunction->AddPoint(1000, 5.0);
		opacityTransferFunction->ClampingOn();
		

	

		vtkColorTransferFunction *colorTransferFunction = vtkColorTransferFunction::New();
			colorTransferFunction->AddRGBPoint(0.0, 0.0, 0.0, 0.0);
			colorTransferFunction->AddRGBPoint(500.0, 0.1, 0.1, 0.1);
			colorTransferFunction->AddRGBPoint(600.0, 0.3, 0.3, 0.3);
			colorTransferFunction->AddRGBPoint(700.0, 0.68, 0.68, 0.68);
			colorTransferFunction->AddRGBPoint(800.0, 0.8, 0.8, 0.8);
			colorTransferFunction->AddRGBPoint(900.0, 0.91, 0.91, 0.91);
			colorTransferFunction->AddRGBPoint(1000.0, 5.0, 5.0, 5.0);


		
		

		vtkVolumeProperty *volumeProperty = vtkVolumeProperty::New();
   	 	volumeProperty->SetColor(colorTransferFunction);
		volumeProperty->SetScalarOpacity(opacityTransferFunction);
		volumeProperty->ShadeOn();
		volumeProperty->SetInterpolationTypeToLinear();

		

		vtkVolumeRayCastCompositeFunction *compositeFunction = vtkVolumeRayCastCompositeFunction::New();
		vtkVolumeRayCastMapper *volumeMapper = vtkVolumeRayCastMapper::New();
		volumeMapper->SetVolumeRayCastFunction(compositeFunction);
		volumeMapper->SetInput(tifReader->GetOutput());
		volumeMapper->SetSampleDistance(0.5);


		vtkVolume *volume = vtkVolume::New();
		volume->SetMapper(volumeMapper);
		volume->SetProperty(volumeProperty);
     
		vtkCamera *camera = vtkCamera::New();
		camera->SetViewUp(0,0,1);
        camera->SetPosition(0,2,0);
		camera->SetFocalPoint(0,0,0);
		camera->Dolly(0.8);

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

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

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


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

		return 0;
		
}