First, I implemented this in Java, but here's what I did. I created a data structure that is an array of sets. The length of the array is the same as the number of planes in the direction I was drawing on (for me that was the XY plane). The elements in the array were sets because I can have multiple contours on each plane. For each contour I used a vtkImageActorPointPlacer. This causes the contours to always be drawn above the image. Then, every time the XY plane that was being viewed changed, I set every contour in my dataset to Off(), and turned On only those contours that were on the plane being viewed. This gives the effect of contours being associated with specific planes.<br>
<br><div class="gmail_quote">On Sat, Oct 16, 2010 at 11:58 PM, <span dir="ltr"><<a href="mailto:shengwen.guo@gmail.com">shengwen.guo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear Mr.Jonathan Morra<br>
I have a question about the contours for 3D data as same as it at <a href="http://vtk.1045678.n5.nabble.com/Question-about-contour-widget-td3208404.html#a3211419" target="_blank">http://vtk.1045678.n5.nabble.com/Question-about-contour-widget-td3208404.html#a3211419</a><br>
Would you please help me to solve this problem or send me a sample code?<br>
Thanks,<br>
Shengwen<br>
<br>
My code is as following:<br>
<br>
<br>
<br>
///////////////////////////////////////////<br>
// TestContours.cxx<br>
// lars<br>
////////////////////////////////////////////<br>
<br>
#include "vtkImageData.h"<br>
#include "vtkVolume16Reader.h"<br>
#include "vtkImageShiftScale.h"<br>
#include "vtkImageViewer2.h"<br>
#include "vtkPoints.h"<br>
#include "vtkPolyData.h"<br>
#include "vtkRenderer.h"<br>
#include "vtkRenderWindow.h"<br>
#include "vtkRenderWindowInteractor.h"<br>
#include "vtkSliderRepresentation.h"<br>
#include "vtkSliderRepresentation2D.h"<br>
#include "vtkSliderWidget.h"<br>
#include "vtkCommand.h"<br>
#include "vtkContourWidget.h"<br>
// Add the following header files by shengwen@GZ 10/17/2010<br>
<br>
#include "vtkOrientedGlyphContourRepresentation.h"<br>
#include "vtkContourWidget.h"<br>
#include "vtkImageActorPointPlacer.h"<br>
#include "vtkTestUtilities.h"<br>
#include "vtkBoundedPlanePointPlacer.h"<br>
<br>
class vtkSliderCallback : public vtkCommand<br>
{<br>
public:<br>
static vtkSliderCallback *New()<br>
{ return new vtkSliderCallback; }<br>
void Delete()<br>
{ delete this; }<br>
void SetImageViewer(vtkImageViewer2 *viewer)<br>
{ m_Viewer = viewer; }<br>
virtual void Execute(vtkObject *caller, unsigned long ul, void* vd)<br>
{<br>
vtkSliderWidget *slider = (vtkSliderWidget *)caller;<br>
vtkSliderRepresentation *sliderRepres = (vtkSliderRepresentation *)slider->GetRepresentation();<br>
int pos = (int)sliderRepres->GetValue();<br>
<br>
m_Viewer->SetSlice(pos);<br>
}<br>
protected:<br>
vtkImageViewer2 *m_Viewer;<br>
};<br>
<br>
int main(int argc, char *argv[])<br>
{<br>
char* fname = vtkTestUtilities::ExpandDataFileName(argc, argv, "Data/headsq/quarter");<br>
//char* fname = "/home/philipp/libs/vtk/VTKData/Data/headsq/quarter";<br>
<br>
vtkVolume16Reader* v16 = vtkVolume16Reader::New();<br>
v16->SetDataDimensions(64, 64);<br>
v16->SetDataByteOrderToLittleEndian();<br>
v16->SetImageRange(1, 93);<br>
v16->SetDataSpacing(3.2, 3.2, 1.5);<br>
v16->SetFilePrefix(fname);<br>
v16->ReleaseDataFlagOn();<br>
v16->SetDataMask(0x7fff);<br>
v16->Update();<br>
//delete[] fname;<br>
<br>
double range[2];<br>
v16->GetOutput()->GetScalarRange(range);<br>
<br>
vtkImageShiftScale* shifter = vtkImageShiftScale::New();<br>
shifter->SetShift(-1.0*range[0]);<br>
shifter->SetScale(255.0/(range[1]-range[0]));<br>
shifter->SetOutputScalarTypeToUnsignedChar();<br>
shifter->SetInputConnection(v16->GetOutputPort());<br>
shifter->ReleaseDataFlagOff();<br>
shifter->Update();<br>
<br>
<br>
vtkImageViewer2 *ImageViewer = vtkImageViewer2::New();<br>
ImageViewer->SetInput(shifter->GetOutput());<br>
ImageViewer->SetColorLevel(127);<br>
ImageViewer->SetColorWindow(255);<br>
<br>
vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();<br>
ImageViewer->SetupInteractor(iren);<br>
<br>
ImageViewer->Render();<br>
ImageViewer->GetRenderer()->ResetCamera();<br>
<br>
ImageViewer->Render();<br>
<br>
int *dims = v16->GetOutput()->GetDimensions();<br>
std::cout << dims[0] << "," << dims[1] << "," << dims[2] << "\n";<br>
<br>
vtkSliderRepresentation2D *SliderRepres = vtkSliderRepresentation2D::New();<br>
int min = ImageViewer->GetSliceMin();<br>
int max = ImageViewer->GetSliceMax();<br>
SliderRepres->SetMinimumValue(min);<br>
SliderRepres->SetMaximumValue(max);<br>
SliderRepres->SetValue((int)((min + max) / 2));<br>
SliderRepres->SetTitleText("Slice");<br>
SliderRepres->GetPoint1Coordinate()->SetCoordinateSystemToNormalizedDisplay();<br>
SliderRepres->GetPoint1Coordinate()->SetValue(0.3, 0.05);<br>
SliderRepres->GetPoint2Coordinate()->SetCoordinateSystemToNormalizedDisplay();<br>
SliderRepres->GetPoint2Coordinate()->SetValue(0.7, 0.05);<br>
SliderRepres->SetSliderLength(0.02);<br>
SliderRepres->SetSliderWidth(0.03);<br>
SliderRepres->SetEndCapLength(0.01);<br>
SliderRepres->SetEndCapWidth(0.03);<br>
SliderRepres->SetTubeWidth(0.005);<br>
SliderRepres->SetLabelFormat("%3.0lf");<br>
SliderRepres->SetTitleHeight(0.02);<br>
SliderRepres->SetLabelHeight(0.02);<br>
<br>
vtkSliderWidget *SliderWidget = vtkSliderWidget::New();<br>
SliderWidget->SetInteractor(iren);<br>
SliderWidget->SetRepresentation(SliderRepres);<br>
SliderWidget->KeyPressActivationOff();<br>
SliderWidget->SetAnimationModeToAnimate();<br>
SliderWidget->SetEnabled(true);<br>
<br>
vtkSliderCallback *SliderCb = vtkSliderCallback::New();<br>
SliderCb->SetImageViewer(ImageViewer);<br>
SliderWidget->AddObserver(vtkCommand::InteractionEvent, SliderCb);<br>
<br>
ImageViewer->SetSlice((int)SliderRepres->GetValue());<br>
<br>
// vtkContourWidget *ContourWidget = vtkContourWidget::New();<br>
// Change new codes<br>
<br>
vtkOrientedGlyphContourRepresentation *rep = vtkOrientedGlyphContourRepresentation::New();<br>
vtkContourWidget *ContourWidget = vtkContourWidget::New();<br>
ContourWidget->SetRepresentation(rep);<br>
// Add<br>
vtkBoundedPlanePointPlacer *placer = vtkBoundedPlanePointPlacer::New();<br>
placer->SetProjectionNormalToZAxis();<br>
//placer.SetProjectionPosition(panel.getImageViewer().GetImageActor().GetCenter()[2]);<br>
rep->SetPointPlacer(placer);<br>
//<br>
<br>
vtkImageActorPointPlacer * imageActorPointPlacer = vtkImageActorPointPlacer::New();<br>
imageActorPointPlacer->SetImageActor(ImageViewer->GetImageActor());<br>
rep->SetPointPlacer(imageActorPointPlacer);<br>
<br>
imageActorPointPlacer->Delete();<br>
rep->Delete();<br>
<br>
// end of changing new codes<br>
<br>
<br>
ContourWidget->SetInteractor(iren);<br>
ContourWidget->SetEnabled(true);<br>
ContourWidget->ProcessEventsOn();<br>
<br>
iren->Start();<br>
<br>
return -1;<br>
}<br>
<br>
<br>
</blockquote></div><br>