VTK 超声图像体绘\BMP图像体绘

#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 “vtkBMPReader.h”
#include “vtkFixedPointVolumeRayCastMapper.h”
#include “vtkStructuredPointsWriter.h”

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

vtkBMPReader *reader = vtkBMPReader::New();	
	reader->SetDataExtent(0,127,0,127,1,51);
	reader->SetFilePrefix("F://超声图像");
	reader->SetFilePattern("%s%d.bmp");
	reader->SetDataSpacing (1, 1, 1.5);//像素间的间隔
	reader->SetNumberOfScalarComponents (1);
	reader->SetAllow8BitBMP(8);
	reader->SwapBytesOn(); 
	reader->SetReleaseDataFlag(1);
	reader->Update();


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

/*
vtkStructuredPointsWriter *wSP=vtkStructuredPointsWriter::New();
wSP->SetInput((vtkDataObject *)readerImageCast->GetOutput());
wSP->SetFileName(“E://CT/mr.VTK”);
wSP->Write();
wSP->Delete();
return 1;
/**/
vtkPiecewiseFunction *opacityTransferFunction = vtkPiecewiseFunction::New();
opacityTransferFunction->AddPoint(0, 0);
opacityTransferFunction->AddPoint(40, 0);
opacityTransferFunction->AddPoint(50, 0.0);
opacityTransferFunction->AddPoint(70, 0.0);
opacityTransferFunction->AddPoint(80, 0.1);
opacityTransferFunction->AddPoint(100, 0.15);
opacityTransferFunction->AddPoint(150, 0.1);
opacityTransferFunction->AddPoint(180, 0.05);
opacityTransferFunction->AddPoint(200, 0);

vtkColorTransferFunction *colorTransferFunction = vtkColorTransferFunction::New();
	colorTransferFunction->AddRGBPoint(0.0, 0.0, 0.0, 0.0);
	colorTransferFunction->AddRGBPoint(50,0.9,0.85,0.8);
	colorTransferFunction->AddRGBPoint(60,1,1,1);
	colorTransferFunction->AddRGBPoint(80,1,0,0);
	colorTransferFunction->AddRGBPoint(100,1,1,0);
	colorTransferFunction->AddRGBPoint(130,0,1,1);
	colorTransferFunction->AddRGBPoint(150,1,0,0.85);
	colorTransferFunction->AddRGBPoint(200,0,0,0);

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();
vtkVolumeRayCastMapper *volumeMapper = vtkVolumeRayCastMapper::New();
volumeMapper->SetVolumeRayCastFunction(compositeFunction);
volumeMapper->SetInputConnection(readerImageCast->GetOutputPort());
*/
vtkFixedPointVolumeRayCastMapper *volumeMapper = vtkFixedPointVolumeRayCastMapper::New();
volumeMapper->SetInput(readerImageCast->GetOutput());
vtkVolume *volume = vtkVolume::New();
volume->SetMapper(volumeMapper);
volume->SetProperty(volumeProperty);
volume->RotateX(130);
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;

}