vtkPlaneWidget 实现实时切割

#include "vtkDICOMImageReader.h"

#include "vtkFixedPointVolumeRayCastMapper.h"

#include "vtkColorTransferFunction.h"

#include "vtkPiecewiseFunction.h"

#include "vtkVolumeProperty.h"

#include "vtkVolume.h"
#include "vtkRenderer.h"

#include "vtkRenderWindow.h"

#include "vtkRenderWindowInteractor.h"
#include "vtkImageMapToColors.h"

#include "vtkPolyDataMapper.h"

#include "vtkActor.h"

#include "vtkCutter.h"

#include "vtkPlane.h"

#include "vtkImageActor.h"

#include "vtkImageData.h"

#include "vtkImageCast.h"

#include "vtkPointData.h"

#include "vtkImageShiftScale.h"

#include "vtkPolyDataToImageStencil.h"

#include "vtkImageStencil.h"

#include "vtkLookupTable.h"

#include "vtkImagePlaneWidget.h"

#include "vtkCellArray.h"

#include "vtkType.h"

#include "vtkPlaneWidget.h"

#include "vtkProperty.h"

#include "vtkCommand.h"

#include "vtkProbeFilter.h"

#include "vtkImplicitPlaneWidget.h"

#include "vtkClipPolyData.h"

class vtkWidgetCall : public vtkCommand

{

public:

static vtkWidgetCall *New()

{

return new vtkWidgetCall;

}

public:

virtual void Execute(vtkObject *caller, unsigned long eventId, void *callData)

{

vtkPlaneWidget *pWidget = vtkPlaneWidget::SafeDownCast(caller);
pWidget->GetPlane(pPlane);
}

public:
vtkPlane *pPlane ;

};

void main()

{
vtkRenderer *pRender1 = vtkRenderer::New();

vtkRenderer *pRender2 = vtkRenderer::New();

vtkRenderWindow *pWin = vtkRenderWindow::New();

vtkRenderWindowInteractor *pInt = vtkRenderWindowInteractor::New();

pWin->SetSize(600, 600);

pWin->AddRenderer(pRender1);

pWin->AddRenderer(pRender2);

pRender1->SetViewport(0, 0, .5, 1);

pRender2->SetViewport(0.5, 0, 1, 1);

pRender2->SetBackground(.3, .3, .3);

pInt->SetRenderWindow(pWin);

vtkDICOMImageReader* pReader = vtkDICOMImageReader::New();

pReader->SetDirectoryName("E:/ExportedImages/2");

pReader->SetDataByteOrderToLittleEndian();

pReader->Update(); //读取
double range[2];

pReader->GetOutput()->GetScalarRange(range);

vtkImageShiftScale *m_pShift = vtkImageShiftScale::New();

m_pShift->SetShift(-1.0*range[0]);

m_pShift->SetScale(255.0/(range[1]-range[0]));

m_pShift->SetOutputScalarTypeToUnsignedChar();

m_pShift->SetInputConnection(pReader->GetOutputPort());

m_pShift->ReleaseDataFlagOff();

m_pShift->Update();
vtkFixedPointVolumeRayCastMapper *pMapper = vtkFixedPointVolumeRayCastMapper::New();

pMapper->SetInputConnection(pReader->GetOutputPort()); //映射:第一步
vtkColorTransferFunction* pColor = vtkColorTransferFunction::New(); //第二步:color

pColor->AddRGBSegment(0, 1, 1, 1, 255, 1, 1, 1);
vtkPiecewiseFunction* pPiceFun = vtkPiecewiseFunction::New();

pPiceFun->AddSegment(0, 0, 3000, 1);//第三步: vtkPiecwise分段
pMapper->SetBlendModeToMaximumIntensity();
vtkVolumeProperty* pProperty = vtkVolumeProperty::New();//设置属性

pProperty->SetColor(pColor);
pProperty->SetScalarOpacity(pPiceFun);

pProperty->SetInterpolationTypeToLinear();

pProperty->ShadeOff();
vtkVolume* pVolume = vtkVolume::New();

pVolume->SetProperty(pProperty);

pVolume->SetMapper(pMapper); //生成体对象
pRender1->AddVolume(pVolume); //增加体数据到渲染器中渲染

vtkPlaneWidget *pWidget = vtkPlaneWidget::New();

pWidget->SetInteractor(pInt);

pWidget->SetInput(m_pShift->GetOutput());

pWidget->SetResolution(10);

pWidget->GetPlaneProperty()->SetColor(.9, .4, .4);

pWidget->GetHandleProperty()->SetColor(0, .4, .7);

pWidget->GetHandleProperty()->SetLineWidth(1.5);

pWidget->NormalToYAxisOn();

pWidget->SetRepresentationToSurface();

pWidget->SetCenter(m_pShift->GetOutput()->GetCenter());

pWidget->PlaceWidget();

pWidget->On();

vtkPlane *pPlane = vtkPlane::New();

pPlane->SetOrigin(pReader->GetOutput()->GetCenter());

pPlane->SetNormal(1, 0, 0);

vtkCutter *pCut = vtkCutter::New();

pCut->SetCutFunction(pPlane);

pCut->SetInputConnection(m_pShift->GetOutputPort());

pCut->Update();

vtkPolyDataMapper *selectMapper = vtkPolyDataMapper::New();

selectMapper->SetInputConnection(pCut->GetOutputPort());
vtkActor *pActor = vtkActor::New();

pActor->SetMapper(selectMapper);

pRender2->AddActor(pActor);

vtkWidgetCall *pCall = vtkWidgetCall::New();

pCall->pPlane = pPlane;

pWidget->AddObserver(vtkCommand::EndInteractionEvent, pCall);
pInt->Initialize();

pInt->Start();

}

效果



本文来源：http://hi.baidu.com/vived/item/476b8d80365761efe596e060