vtkimageviewer 和vtkcontourwidget交互提取感兴趣区域！

 #include "vtkImageData.h" 

#include "vtkVolume16Reader.h" 

#include "vtkImageShiftScale.h" 

#include "vtkImageViewer2.h" 

#include "vtkPoints.h" 

#include "vtkPolyData.h" 

#include "vtkRenderer.h" 

#include "vtkRenderWindow.h" 

#include "vtkRenderWindowInteractor.h" 

#include "vtkSliderRepresentation.h" 

#include "vtkSliderRepresentation2D.h" 

#include "vtkSliderWidget.h" 

#include "vtkCommand.h" 

#include "vtkContourWidget.h" 

#include "vtkOrientedGlyphContourRepresentation.h"

#include "vtkImageActorPointPlacer.h"

 

class vtkSliderCallback : public vtkCommand 

{ 

public: 

 static vtkSliderCallback *New() 

 { return new vtkSliderCallback; } 

 void Delete() 

 { delete this; } 

 void SetImageViewer(vtkImageViewer2 *viewer) 

 { m_Viewer =  viewer; } 

 virtual void Execute(vtkObject *caller, unsigned long ul, void* vd) 

 { 

 vtkSliderWidget *slider = (vtkSliderWidget *)caller; 

 vtkSliderRepresentation *sliderRepres = (vtkSliderRepresentation *)slider->GetRepresentation();       

 int pos = (int)sliderRepres->GetValue(); 

 

 m_Viewer->SetSlice(pos); 

 } 

protected: 

 vtkImageViewer2 *m_Viewer; 

}; 

 

int main(int argc, char *argv[]) 

{ 

 char* fname = "D:/VTK/vtkdata/Data/headsq/quarter";   

 

 vtkVolume16Reader* v16 = vtkVolume16Reader::New(); 

 v16->SetDataDimensions(64, 64); 

 v16->SetDataByteOrderToLittleEndian(); 

 v16->SetImageRange(1, 93); 

 v16->SetDataSpacing(3.2, 3.2, 1.5); 

 v16->SetFilePrefix(fname); 

 v16->ReleaseDataFlagOn(); 

 v16->SetDataMask(0x7fff); 

 v16->Update(); 

 //delete[] fname; 

 

 double range[2]; 

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

 

 vtkImageShiftScale* shifter = vtkImageShiftScale::New(); 

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

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

 shifter->SetOutputScalarTypeToUnsignedChar(); 

 shifter->SetInputConnection(v16->GetOutputPort()); 

 shifter->ReleaseDataFlagOff(); 

 shifter->Update(); 

 

 

 vtkImageViewer2 *ImageViewer = vtkImageViewer2::New(); 

 ImageViewer->SetInput(shifter->GetOutput()); 

 ImageViewer->SetColorLevel(127); 

 ImageViewer->SetColorWindow(255); 

 

 vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New(); 

 ImageViewer->SetupInteractor(iren); 

 

 ImageViewer->Render(); 

 ImageViewer->GetRenderer()->ResetCamera(); 

 

 ImageViewer->Render();     

 

 int *dims = v16->GetOutput()->GetDimensions(); 

 std::cout << dims[0] << "," << dims[1] << "," << dims[2] << "\n"; 

 

 vtkSliderRepresentation2D *SliderRepres = vtkSliderRepresentation2D::New(); 

 int min = ImageViewer->GetSliceMin(); 

 int max = ImageViewer->GetSliceMax(); 

 SliderRepres->SetMinimumValue(min); 

 SliderRepres->SetMaximumValue(max); 

 SliderRepres->SetValue((int)((min + max) / 2)); 

 SliderRepres->SetTitleText("Slice"); 

 SliderRepres->GetPoint1Coordinate()->SetCoordinateSystemToNormalizedDisplay(); 

 SliderRepres->GetPoint1Coordinate()->SetValue(0.3, 0.05); 

 SliderRepres->GetPoint2Coordinate()->SetCoordinateSystemToNormalizedDisplay(); 

 SliderRepres->GetPoint2Coordinate()->SetValue(0.7, 0.05); 

 SliderRepres->SetSliderLength(0.02); 

 SliderRepres->SetSliderWidth(0.03); 

 SliderRepres->SetEndCapLength(0.01); 

 SliderRepres->SetEndCapWidth(0.03); 

 SliderRepres->SetTubeWidth(0.005); 

 SliderRepres->SetLabelFormat("%3.0lf"); 

 SliderRepres->SetTitleHeight(0.02); 

 SliderRepres->SetLabelHeight(0.02); 

 

 vtkSliderWidget *SliderWidget = vtkSliderWidget::New(); 

 SliderWidget->SetInteractor(iren); 

 SliderWidget->SetRepresentation(SliderRepres); 

 SliderWidget->KeyPressActivationOff(); 

 SliderWidget->SetAnimationModeToAnimate(); 

 SliderWidget->SetEnabled(true); 

 

 vtkSliderCallback *SliderCb = vtkSliderCallback::New(); 

 SliderCb->SetImageViewer(ImageViewer); 

 SliderWidget->AddObserver(vtkCommand::InteractionEvent, SliderCb);   

 

 ImageViewer->SetSlice((int)SliderRepres->GetValue()); 

 

 vtkContourWidget *ContourWidget = vtkContourWidget::New(); 

 

 vtkOrientedGlyphContourRepresentation *rep = 

 vtkOrientedGlyphContourRepresentation::New(); 

 ContourWidget->SetRepresentation(rep); 

 

 vtkImageActorPointPlacer * imageActorPointPlacer = vtkImageActorPointPlacer::New(); 

 imageActorPointPlacer->SetImageActor(ImageViewer->GetImageActor()); 

 rep->SetPointPlacer(imageActorPointPlacer); 

 

 imageActorPointPlacer->Delete(); 

 rep->Delete(); 

 

 ContourWidget->SetInteractor(iren); 

 ContourWidget->SetEnabled(true); 

 ContourWidget->ProcessEventsOn(); 

 

 iren->Start(); 

 

 return -1; 

}

 

 

本文来源：http://www.vislab.cn/bbs/viewthread.php?tid=1467&extra=page%3D7