-/*=========================================================================
+/*
+ # ---------------------------------------------------------------------
+ #
+ # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
+ # pour la SantÈ)
+ # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
+ # Previous Authors : Laurent Guigues, Jean-Pierre Roux
+ # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
+ #
+ # This software is governed by the CeCILL-B license under French law and
+ # abiding by the rules of distribution of free software. You can use,
+ # modify and/ or redistribute the software under the terms of the CeCILL-B
+ # license as circulated by CEA, CNRS and INRIA at the following URL
+ # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+ # or in the file LICENSE.txt.
+ #
+ # As a counterpart to the access to the source code and rights to copy,
+ # modify and redistribute granted by the license, users are provided only
+ # with a limited warranty and the software's author, the holder of the
+ # economic rights, and the successive licensors have only limited
+ # liability.
+ #
+ # The fact that you are presently reading this means that you have had
+ # knowledge of the CeCILL-B license and that you accept its terms.
+ # ------------------------------------------------------------------------ */
+
+
+/*=========================================================================
Program: bbtk
Module: $RCSfile: bbvtkImagePlanes.cxx,v $
Language: C++
- Date: $Date: 2010/05/06 09:05:32 $
- Version: $Revision: 1.33 $
+ Date: $Date: 2012/11/16 08:51:58 $
+ Version: $Revision: 1.38 $
=========================================================================*/
-/* ---------------------------------------------------------------------
-
-* Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale)
-* Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux
-*
-* This software is governed by the CeCILL-B license under French law and
-* abiding by the rules of distribution of free software. You can use,
-* modify and/ or redistribute the software under the terms of the CeCILL-B
-* license as circulated by CEA, CNRS and INRIA at the following URL
-* http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
-* or in the file LICENSE.txt.
-*
-* As a counterpart to the access to the source code and rights to copy,
-* modify and redistribute granted by the license, users are provided only
-* with a limited warranty and the software's author, the holder of the
-* economic rights, and the successive licensors have only limited
-* liability.
-*
-* The fact that you are presently reading this means that you have had
-* knowledge of the CeCILL-B license and that you accept its terms.
-* ------------------------------------------------------------------------ */
-
/**
* \file
* \brief
#include "bbstdRelay.h"
#include "vtkObjectFactory.h"
+#include "vtkImageFlip.h"
+#include "vtkImageReslice.h"
+#include "vtkImageChangeInformation.h"
namespace bbstd
{
};
//================================================================
- //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.33 $");
+ //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.38 $");
//================================================================
image = NULL;
_transform = NULL;
_matrix = NULL;
+
+ bbSetInputInterpolation(1);
}
vtkImagePlaneWidget* planeWidget3Pts = GetPlaneWidget('3', 0, 1, 1, picker);
planeWidget3Pts->SetLookupTable(planeWidgetX->GetLookupTable());
+ // EED MPR view orientation correction..
+ vtkImageFlip *flipYFilter = vtkImageFlip::New();
+ flipYFilter->SetFilteredAxis(1); // flip y axis
+ flipYFilter->SetInput( planeWidgetX->GetResliceOutput() );
+ flipYFilter->Update();
+
+ vtkImageChangeInformation *image = vtkImageChangeInformation::New();
+ image->SetInput( planeWidgetY->GetResliceOutput() );
+ image->SetOutputSpacing( 1,1,1 );
+ image->CenterImageOn();
+ image->Update();
+ _imageTransform = vtkTransform::New();
+ vtkImageReslice *slicer =vtkImageReslice::New();
+ slicer->SetInput( image->GetOutput() );
+ slicer->SetInformationInput( image->GetOutput() );
+ slicer->SetResliceTransform( _imageTransform );
+ slicer->SetOutputOrigin(0 , 0 , 0 );
+ slicer->SetInterpolationModeToNearestNeighbor();
+ slicer->Update();
+ vtkImageChangeInformation *imageResult = vtkImageChangeInformation::New();
+ imageResult->SetInput( slicer->GetOutput() );
+ double spc[3];
+ planeWidgetY->GetResliceOutput()->GetSpacing(spc);
+ imageResult->SetOutputSpacing( spc[1], spc[0], spc[2] );
+ imageResult->SetOutputOrigin( 0,0,0 );
+
+
bbSetOutputPlaneX(planeWidgetX);
bbSetOutputPlaneY(planeWidgetY);
bbSetOutputPlaneZ(planeWidgetZ);
bbSetOutputPlane3Pts(planeWidget3Pts);
- bbSetOutputImageX(planeWidgetX->GetResliceOutput());
- bbSetOutputImageY(planeWidgetY->GetResliceOutput());
+ bbSetOutputImageX( flipYFilter->GetOutput() ); // EED MPR view orientation correction..
+ bbSetOutputImageY( imageResult->GetOutput() ); // EED MPR view orientation correction..
bbSetOutputImageZ(planeWidgetZ->GetResliceOutput());
bbSetInputInteractor(0);
//bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput());
- picker->UnRegister(NULL);
+ if(picker != 0)
+ picker->UnRegister(NULL);
mVtkCallback = VtkCallbackType::New();
mVtkCallback->SetBlackBox(this);
planeWidgetX->AddObserver(vtkCommand::InteractionEvent,mVtkCallback);
planeWidgetY->AddObserver(vtkCommand::InteractionEvent,mVtkCallback);
- planeWidgetZ->AddObserver(vtkCommand::InteractionEvent,mVtkCallback);
+ planeWidgetZ->AddObserver(vtkCommand::InteractionEvent,mVtkCallback);
+
}
//---------------------------------------------------------------------
void ImagePlanes::bbUserFinalizeProcessing()
{
-
+
if (bbGetOutputPlaneX())
{
//---------------------------------------------------------------------
void ImagePlanes::Process()
{
+
if (bbGetInputIn()!=0)
{
- int xMin, xMax, yMin, yMax, zMin, zMax;
- bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
-
+ int dim[3];
+ int ext[6];
+ bbGetOutputPlaneX()->GetResliceOutput()->GetWholeExtent(ext);
+ dim[0] = ext[1]-ext[0]+1;
+ dim[1] = ext[3]-ext[2]+1;
+ dim[2] = ext[5]-ext[4]+1;
+ _imageTransform->Identity();
+ _imageTransform->PostMultiply();
+ _imageTransform->Translate( (int)(-(dim[0]/2)) , (int)(-(dim[1]/2)) ,0);
+ _imageTransform->RotateZ(90);
+
+
if ( image != bbGetInputIn()){//bbGetInputStatus("In") != bbtk::UPTODATE ){
// Input image has changed : reinitialize planes
image = bbGetInputIn();
// Initial values : center of the volume (in real world, not in pixels!)
+ int xMin, xMax, yMin, yMax, zMin, zMax;
+ bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
double xSpacing, ySpacing, zSpacing;
bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing);
{
bbGetOutputPlaneZ()->SetWindowLevel(bbGetInputWindowLevel()[0],
bbGetInputWindowLevel()[1]);
- }
- else
- {
+ } else {
double *range = image->GetScalarRange();
bbGetOutputPlaneZ()->SetWindowLevel(range[1]-range[0],
0.5*(range[1]+range[0]));
- }
+ } // windowlevel
updateInteractor();
- }
+
+ dim[0] = xMax-xMin+1;
+ dim[1] = yMax-yMin+1;
+ dim[2] = zMax-zMin+1;
+ _imageTransform->Identity();
+ _imageTransform->PostMultiply();
+ _imageTransform->Translate( (int)(-(dim[0]/2)*(1/xSpacing)) , (int)(-(dim[2]/2)*(1/zSpacing)) ,0);
+ _imageTransform->RotateZ(90);
+
+ } // image
+
// UPDATE DES SORTIES
+ bbGetOutputPlaneX()->SetResliceInterpolate( bbGetInputInterpolation() );
+ bbGetOutputPlaneY()->SetResliceInterpolate( bbGetInputInterpolation() );
+ bbGetOutputPlaneZ()->SetResliceInterpolate( bbGetInputInterpolation() );
+
+
bbGetOutputPlaneX()->GetResliceOutput()->Update();
bbGetOutputPlaneY()->GetResliceOutput()->Update();
bbGetOutputPlaneZ()->GetResliceOutput()->Update();
-
+
std::vector<int> pointsx = bbGetInputPointsX();
std::vector<int> pointsy = bbGetInputPointsY();
std::vector<int> pointsz = bbGetInputPointsZ();
if (_imageReslicer==NULL){
_imageReslicer = vtkImageReslice::New();
_imageReslicer->SetOutputDimensionality(2);
- _imageReslicer->SetInterpolationModeToLinear();
_transform = vtkTransform::New();
_matrix = vtkMatrix4x4::New();
}
+ _imageReslicer->SetInterpolationMode( bbGetInputInterpolation() );
_imageReslicer->SetInput( bbGetInputIn() );
_imageReslicer->SetInformationInput(bbGetInputIn());
//fill out the information with the created vectors and using the spacing of the image
//set the transformation out to be used by other bbBoxes
bbSetOutputTransform3Pts((vtkLinearTransform*)_transform);
- }
- }
+ } // pointsx pointsy pointsz
+ } // bbGetInputIn
}
void ImagePlanes::updateInteractor(){