medical3.cxx

VTK入门样例(种类齐全)

#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkVolume16Reader.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkOutlineFilter.h"
#include "vtkCamera.h"
#include "vtkStripper.h"
#include "vtkLookupTable.h"
#include "vtkImageDataGeometryFilter.h"
#include "vtkMarchingCubes.h"
#include "vtkImageShrink3D.h"
#include "vtkProperty.h"
#include "vtkDICOMImageReader.h"

void main( int argc, char *argv[] )
{
	// create the renderer stuff
	vtkRenderer *aRenderer = vtkRenderer::New();
	vtkRenderWindow *renWin = vtkRenderWindow::New();
	renWin->AddRenderer(aRenderer);
	vtkRenderWindowInteractor *iren =
	vtkRenderWindowInteractor::New();
	iren->SetRenderWindow(renWin);

	// read the volume
  vtkDICOMImageReader   *reader =  vtkDICOMImageReader::New();
	reader->SetDataByteOrderToLittleEndian();
	reader->SetDirectoryName("e://ct/head40/");
	reader->SetDataOrigin(0.0, 0.0, 0.0);
	reader->Update();
	int extent[6];
	reader->GetDataExtent(extent);
    double spacing [3];  
	reader->GetDataSpacing(spacing); 
	int xx,yy,zz;
	xx=(int)((extent[1]+1)*spacing[0]);
	yy=(int)((extent[3]+1)*spacing[1]);
	zz=(int)((extent[5]+1)*spacing[2]);
cout<<xx<<"  "<<yy<<"  "<<zz<<"  ";
  vtkImageShrink3D *shrink=vtkImageShrink3D::New();
		shrink->SetInput((vtkDataObject *)reader->GetOutput());
		shrink->SetShrinkFactors(4,4,1);

	// extract the skin
	vtkMarchingCubes *skinExtractor = vtkMarchingCubes::New();
	skinExtractor->SetInputConnection(shrink->GetOutputPort());
	skinExtractor->SetValue(0, -500);
	vtkStripper *skinStripper = vtkStripper::New();
	skinStripper->SetInput(skinExtractor->GetOutput());
	vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New();
	skinMapper->SetInput(skinStripper->GetOutput());
	skinMapper->ScalarVisibilityOff();
	vtkActor *skin = vtkActor::New();
	skin->SetMapper(skinMapper);
	skin->GetProperty()->SetDiffuseColor(1, .49, .25);
	skin->GetProperty()->SetSpecular(.3);
	skin->GetProperty()->SetSpecularPower(20);

	// extract the bone
	vtkMarchingCubes *boneExtractor = vtkMarchingCubes::New();
	boneExtractor->SetInputConnection(shrink->GetOutputPort());
	boneExtractor->SetValue(0, 200);
	vtkStripper *boneStripper = vtkStripper::New();
	boneStripper->SetInput(boneExtractor->GetOutput());
	vtkPolyDataMapper *boneMapper = vtkPolyDataMapper::New();
	boneMapper->SetInput(boneStripper->GetOutput());
	boneMapper->ScalarVisibilityOff();
	vtkActor *bone = vtkActor::New();
	bone->SetMapper(boneMapper);
	bone->GetProperty()->SetDiffuseColor(1, 1, .9412);

	// get an outline
	vtkOutlineFilter *outlineData = vtkOutlineFilter::New();
	outlineData->SetInputConnection(shrink->GetOutputPort());
	vtkPolyDataMapper *mapOutline = vtkPolyDataMapper::New();
	mapOutline->SetInput(outlineData->GetOutput());
	vtkActor *outline = vtkActor::New();
	outline->SetMapper(mapOutline);
	outline->GetProperty()->SetColor(0,0,0);

	// create a b/w lookup table
	vtkLookupTable *bwLut = vtkLookupTable::New();
	bwLut->SetTableRange (0, 2000);
	bwLut->SetSaturationRange (0, 0);
	bwLut->SetHueRange (0, 0);
	bwLut->SetValueRange (0, 1);
	// create a hue lookup table
	vtkLookupTable *hueLut = vtkLookupTable::New();
	hueLut->SetTableRange (0, 2000);
	hueLut->SetHueRange (0, 1);
	hueLut->SetSaturationRange (1, 1);
	hueLut->SetValueRange (1, 1);
	// create a saturation lookup table
	vtkLookupTable *satLut = vtkLookupTable::New();
	satLut->SetTableRange (0, 2000);
	satLut->SetHueRange (.6, .6);
	satLut->SetSaturationRange (0, 1);
	satLut->SetValueRange (1, 1);

	// sagittal
	vtkImageDataGeometryFilter *saggitalSection =
	vtkImageDataGeometryFilter::New();
	saggitalSection->SetExtent (65,65, 0,yy, 0, zz);
	saggitalSection->SetInputConnection(shrink->GetOutputPort());
	vtkPolyDataMapper *saggitalMapper = vtkPolyDataMapper::New();
	saggitalMapper->SetInput(saggitalSection->GetOutput());
	saggitalMapper->ScalarVisibilityOn();
	saggitalMapper->SetScalarRange (-1024, 2000);
	saggitalMapper->SetLookupTable (bwLut);
	vtkActor *sagittal = vtkActor::New();
	sagittal->SetMapper(saggitalMapper);

	// axial
	vtkImageDataGeometryFilter *axialSection =
	vtkImageDataGeometryFilter::New();
	axialSection->SetExtent (0,xx, 0,yy, 22, 22);
	axialSection->SetInputConnection(shrink->GetOutputPort());
	vtkPolyDataMapper *axialMapper = vtkPolyDataMapper::New();
	axialMapper->SetInput(axialSection->GetOutput());
	axialMapper->ScalarVisibilityOn();
	axialMapper->SetScalarRange (-1024, 2000);
	axialMapper->SetLookupTable (hueLut);
	vtkActor *axial = vtkActor::New();
	axial->SetMapper(axialMapper);

	// coronal
	vtkImageDataGeometryFilter *coronalSection =
	vtkImageDataGeometryFilter::New();
	coronalSection->SetExtent (0,xx, 65, 65, 0, zz);
	coronalSection->SetInputConnection(shrink->GetOutputPort());
	vtkPolyDataMapper *coronalMapper =
	vtkPolyDataMapper::New();
	coronalMapper->SetInput(coronalSection->GetOutput());
	coronalMapper->ScalarVisibilityOn();
	coronalMapper->SetScalarRange (-1024, 2000);
	coronalMapper->SetLookupTable (satLut);
	vtkActor *coronal = vtkActor::New();
	coronal->SetMapper(coronalMapper);

	// create a camera with the correct view up
	vtkCamera *aCamera = vtkCamera::New();
	aCamera->SetViewUp (0, 0, -1);
	aCamera->SetPosition (0, 1, 0);
	aCamera->SetFocalPoint (0, 0, 0);
	// now, tell the renderer our actors
	aRenderer->AddActor(outline);
	aRenderer->AddActor(sagittal);
	aRenderer->AddActor(axial);
	aRenderer->AddActor(coronal);
	aRenderer->AddActor(skin);
	aRenderer->AddActor(bone);

	// turn off bone for this examples
//	bone->VisibilityOff();
	// set skin to semi-transparent
	skin->GetProperty()->SetOpacity(0.5);
	bone->GetProperty()->SetOpacity(0.3);
	aRenderer->SetActiveCamera(aCamera);
	aRenderer->ResetCamera ();
	aCamera->Dolly(1.5);
	aRenderer->SetBackground(1,1,1);
	aRenderer->ResetCameraClippingRange();
	// interact with data
	renWin->SetSize(600, 600);
	renWin->Render();
//	SAVEIMAGE( renWin );
	iren->Start();
}