X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=inline;f=vv%2FvvImageContour.cxx;h=173bde0e900e9f80b28cdf5c85bf87d906b8afb7;hb=5a7da4aedae5c204bc55c187717193e5950f9a44;hp=27703963c345be05e09715957a59423295319d38;hpb=dd8a36ea0968ce26b9f25166c1c62354cbb4717b;p=clitk.git diff --git a/vv/vvImageContour.cxx b/vv/vvImageContour.cxx index 2770396..173bde0 100644 --- a/vv/vvImageContour.cxx +++ b/vv/vvImageContour.cxx @@ -36,6 +36,7 @@ vvImageContour::vvImageContour() mDisplayModeIsPreserveMemory = true; SetPreserveMemoryModeEnabled(true); mPreviousOrientation = -1; + mDepth = 1.0; } //------------------------------------------------------------------------------ @@ -267,6 +268,7 @@ void vvImageContour::CreateNewActor(int numImage) { clipper->SetInput(mHiddenImage->GetVTKImages()[0]); else clipper->SetInput(mSlicer->GetImage()->GetVTKImages()[numImage]); + squares->SetInput(clipper->GetOutput()); squaresMapper->SetInput(squares->GetOutput()); squaresMapper->ScalarVisibilityOff(); @@ -290,8 +292,7 @@ void vvImageContour::UpdateActor(vtkActor * actor, vtkMarchingSquares * squares, vtkImageClip * clipper, double threshold, int orientation, int slice) { - - // Set parameter for the MarchigSquare + // Set parameter for the MarchigSquare squares->SetValue(0, threshold); // Get image extent @@ -311,15 +312,16 @@ void vvImageContour::UpdateActor(vtkActor * actor, s = s-mHiddenImage->GetFirstVTKImageData()->GetOrigin()[orientation]; // from corner second image s = s/mHiddenImage->GetFirstVTKImageData()->GetSpacing()[orientation]; // in voxel - if (s == floor(s)) { - extent2[orientation*2] = extent2[orientation*2+1] = (int)floor(s); - } else { - extent2[orientation*2] = (int)floor(s); - extent2[orientation*2+1] = extent2[orientation*2]; - } + // Rint to the closest slice + extent2[orientation*2+1] = extent2[orientation*2] = (int)lrint(s); // Do not display a contour if there is no contour on this slice - if (extent2[orientation*2+1] > extent3[orientation*2+1]) { + // DD(extent2[orientation*2+1]); + // DD(extent3[orientation*2+1]); + // DD(extent2[orientation*2]); + // DD(extent3[orientation*2]); + if ((extent2[orientation*2+1] > extent3[orientation*2+1]) || + (extent2[orientation*2] < extent3[orientation*2])) { actor->VisibilityOff(); return; } @@ -333,56 +335,14 @@ void vvImageContour::UpdateActor(vtkActor * actor, clipper->SetOutputWholeExtent(extent2[0],extent2[1],extent2[2], extent2[3],extent2[4],extent2[5]); - //std::cout << mTSlice << " " << mSlice << " " << extent2[0] << " " << extent2[1] << " " << extent2[2] << " " << extent2[3] << " " << extent2[4] << " " << extent2[5] << std::endl; - if (mHiddenImageIsUsed) delete extent2; // Move the actor to be visible double position[3] = {0, 0, 0}; - position[orientation] = -1; + position[orientation] = -mDepth; actor->SetPosition(position); -// switch (orientation) { -// case 0: -// actor->SetPosition(-1,0,0); -// /* -// // DD(mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[0]); -// if (mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[0] > slice) { -// actor->SetPosition(1,0,0); -// } else { -// actor->SetPosition(-1,0,0); -// }*/ -// break; -// case 1: -// actor->SetPosition(0,-1,0); -// /* -// // DD(mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[1]); -// if (mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[1] > slice) { -// actor->SetPosition(0,1,0); -// } else { -// actor->SetPosition(0,-1,0); -// } -// */ -// break; -// case 2: -// actor->SetPosition(0,0,-1); -// /* -// DD(mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[2]); -// if (mSlicer->GetRenderer()->GetActiveCamera()->GetPosition()[2] > slice) { -// DD("1"); -// actor->SetPosition(0,0,1); -// } else { -// DD("-1"); -// actor->SetPosition(0,0,-1); -// } -// */ -// break; -// } - mapper->Update(); - - - } //------------------------------------------------------------------------------