X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=9a6fc8ff489f39c1d661b1a55940a46253d3a2c3;hb=6ab0361c289035b4ba8487fe074ea9b268341da9;hp=2a47a99d6463989ff9157009aa6057f87f509b0a;hpb=48d8569538b93bf923de2f8d5c4b9db51640978f;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index 2a47a99..9a6fc8f 100644 --- a/common/vvImage.cxx +++ b/common/vvImage.cxx @@ -1,9 +1,9 @@ /*========================================================================= Program: vv http://www.creatis.insa-lyon.fr/rio/vv - Authors belong to: + Authors belong to: - University of LYON http://www.universite-lyon.fr/ - - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr + - Léon Bérard cancer center http://www.centreleonberard.fr - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr This software is distributed WITHOUT ANY WARRANTY; without even @@ -14,228 +14,196 @@ - BSD See included LICENSE.txt file - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html -======================================================================-====*/ +===========================================================================**/ #ifndef VVIMAGE_CXX #define VVIMAGE_CXX + +// vv #include "vvImage.h" + +// clitk #include "clitkCommon.h" -#include + +// vtk #include -#include +#include + +// std #include //-------------------------------------------------------------------- -vvImage::vvImage() { +vvImage::vvImage():mTransform(vtkSmartPointer::New()) +{ Init(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -void vvImage::Init() { - mTimeSpacing = 1; - mTimeOrigin = 0; - - if (CLITK_EXPERIMENTAL) - { - mVtkImageReslice = vtkSmartPointer::New(); - transform = vtkSmartPointer::New(); - mVtkImageReslice->SetInterpolationModeToLinear(); - mVtkImageReslice->AutoCropOutputOn(); - mVtkImageReslice->SetBackgroundColor(-1000,-1000,-1000,1); - mVtkImageReslice->SetResliceTransform(transform); - } -} -//-------------------------------------------------------------------- - -//-------------------------------------------------------------------- -vvImage::~vvImage() { - Reset(); +void vvImage::Init() +{ + mTimeSpacing = 1; + mTimeOrigin = 0; + mImageDimension = 0; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -void vvImage::Reset() { - if (CLITK_EXPERIMENTAL) - { - for (unsigned int i = 0; i < mVtkImages.size(); i++) - mVtkImages[i] = (vtkImageData*)mVtkImageReslice->GetInput(i); - mVtkImageReslice->RemoveAllInputs(); - } - for (unsigned int i = 0; i < mVtkImages.size(); i++) - mVtkImages[i]->Delete(); - mVtkImages.resize(0); - Init(); +vvImage::~vvImage() +{ + Reset(); } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -void vvImage::SetImage(std::vector< vtkImageData* > images) { - Reset(); - for (unsigned int i = 0; i < images.size(); i++) - AddImage(images[i]); +void vvImage::Reset() +{ + mVtkImages.clear(); + Init(); } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -void vvImage::AddImage(vtkImageData* image) { - if (CLITK_EXPERIMENTAL) - { - mVtkImageReslice->SetInput(mVtkImages.size(), image); - mVtkImageReslice->Update(); - mVtkImages.push_back( mVtkImageReslice->GetOutput( mVtkImages.size() ) ); - } - else - mVtkImages.push_back(image); +int vvImage::GetNumberOfSpatialDimensions() +{ + return mImageDimension; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -int vvImage::GetNumberOfSpatialDimensions() { - int dim=GetNumberOfDimensions(); +int vvImage::GetNumberOfDimensions() const +{ + if (mVtkImages.size()) { if (IsTimeSequence()) - return dim-1; + return mImageDimension+1; else - return dim; + return mImageDimension; + } + return 0; } //-------------------------------------------------------------------- -// -//-------------------------------------------------------------------- -int vvImage::GetNumberOfDimensions() const { - if (mVtkImages.size()) - { - int dimension = 2; - int extent[6]; - mVtkImages[0]->GetWholeExtent(extent); - if (extent[5] - extent[4] >= 1) - dimension++; - if (mVtkImages.size() > 1) - dimension++; - return dimension; - } - return 0; -} + //-------------------------------------------------------------------- void vvImage::GetScalarRange(double* range) { - assert(mVtkImages.size()); - double * temp = mVtkImages[0]->GetScalarRange(); - range[0]=temp[0];range[1]=temp[1]; - for (unsigned int i=1;iGetScalarRange(); - if (temp[0] < range[0]) range[0]=temp[0]; - if (temp[1] > range[1]) range[1]=temp[1]; - } + assert(mVtkImages.size()); + double * temp = mVtkImages[0]->GetScalarRange(); + range[0]=temp[0]; + range[1]=temp[1]; + for (unsigned int i=1; iGetScalarRange(); + if (temp[0] < range[0]) range[0]=temp[0]; + if (temp[1] > range[1]) range[1]=temp[1]; + } } +//-------------------------------------------------------------------- //-------------------------------------------------------------------- -std::string vvImage::GetScalarTypeAsString() { - return mVtkImages[0]->GetScalarTypeAsString(); +std::string vvImage::GetScalarTypeAsITKString() +{ + // WARNING VTK pixel type different from ITK Pixel type + std::string vtktype = mVtkImages[0]->GetScalarTypeAsString(); + if (vtktype == "unsigned char") return "unsigned_char"; + if (vtktype == "unsigned short") return "unsigned_short"; + if (vtktype == "unsigned int") return "unsigned_int"; + return vtktype; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -int vvImage::GetNumberOfScalarComponents() { - return mVtkImages[0]->GetNumberOfScalarComponents(); +int vvImage::GetNumberOfScalarComponents() +{ + return mVtkImages[0]->GetNumberOfScalarComponents(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -int vvImage::GetScalarSize() { - return mVtkImages[0]->GetScalarSize(); +int vvImage::GetScalarSize() +{ + return mVtkImages[0]->GetScalarSize(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -std::vector vvImage::GetSpacing() { - std::vector spacing; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) - { - if (i == 3) - spacing.push_back(1); - else - spacing.push_back(mVtkImages[0]->GetSpacing()[i]); - if (i == 3) - spacing.push_back(mTimeSpacing); - else - spacing.push_back(mVtkImages[0]->GetSpacing()[i]); - } - return spacing; +std::vector vvImage::GetSpacing() +{ + std::vector spacing; + int dim = this->GetNumberOfDimensions(); + for (int i = 0; i < dim; i++) { + if (i == 3) + spacing.push_back(mTimeSpacing); + else + spacing.push_back(mVtkImages[0]->GetSpacing()[i]); + } + return spacing; } //-------------------------------------------------------------------- - -//-------------------------------------------------------------------- -std::vector vvImage::GetOrigin() const { - std::vector origin; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) - { - if (i == 3) - origin.push_back(0); - else - origin.push_back(mVtkImages[0]->GetOrigin()[i]); - if (i == 3) - origin.push_back(mTimeOrigin); - else - origin.push_back(mVtkImages[0]->GetOrigin()[i]); - } - return origin; +std::vector vvImage::GetOrigin() const +{ + std::vector origin; + int dim = this->GetNumberOfDimensions(); + for (int i = 0; i < dim; i++) { + if (i == 3) + origin.push_back(mTimeOrigin); + else + origin.push_back(mVtkImages[0]->GetOrigin()[i]); + } + return origin; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -std::vector vvImage::GetSize() { - std::vector size0; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) - { - if (i == 3) - size0.push_back(mVtkImages.size()); - else - size0.push_back(mVtkImages[0]->GetDimensions()[i]); - } - return size0; +std::vector vvImage::GetSize() +{ + std::vector size0; + int dim = this->GetNumberOfDimensions(); + for (int i = 0; i < dim; i++) { + if (i == 3) + size0.push_back(mVtkImages.size()); + else + size0.push_back(mVtkImages[0]->GetDimensions()[i]); + } + return size0; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -unsigned long vvImage::GetActualMemorySize() { - unsigned long size = 0; - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - size += mVtkImages[i]->GetActualMemorySize(); - } - return size; +unsigned long vvImage::GetActualMemorySize() +{ + unsigned long size = 0; + for (unsigned int i = 0; i < mVtkImages.size(); i++) { + size += mVtkImages[i]->GetActualMemorySize(); + } + return size; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -bool vvImage::IsTimeSequence() { +bool vvImage::IsTimeSequence() const +{ return mVtkImages.size()>1; } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -const std::vector& vvImage::GetVTKImages() { +const std::vector& vvImage::GetVTKImages() +{ return mVtkImages; } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -vtkImageData* vvImage::GetFirstVTKImageData() { - return mVtkImages[0]; +vtkImageData* vvImage::GetFirstVTKImageData() +{ + return mVtkImages[0]; } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -bool vvImage::IsScalarTypeInteger() { +bool vvImage::IsScalarTypeInteger() +{ assert(mVtkImages.size()> 0); int t = mVtkImages[0]->GetScalarType(); return IsScalarTypeInteger(t); @@ -243,7 +211,8 @@ bool vvImage::IsScalarTypeInteger() { //-------------------------------------------------------------------- //-------------------------------------------------------------------- -bool vvImage::IsScalarTypeInteger(int t) { +bool vvImage::IsScalarTypeInteger(int t) +{ if ((t == VTK_BIT) || (t == VTK_CHAR) || (t == VTK_UNSIGNED_CHAR) || @@ -254,103 +223,31 @@ bool vvImage::IsScalarTypeInteger(int t) { (t == VTK_LONG) || (t == VTK_UNSIGNED_LONG)) { return true; - } - else { + } else { return false; } } - //-------------------------------------------------------------------- //-------------------------------------------------------------------- - -void vvImage::SetVTKImageReslice(vtkImageReslice* reslice) { - mVtkImageReslice=reslice; - -} -//-------------------------------------------------------------------- - -//-------------------------------------------------------------------- -vtkAbstractTransform * vvImage::GetTransform() +vtkSmartPointer vvImage::GetTransform() { - return mVtkImageReslice->GetResliceTransform(); + return this->mTransform; } //-------------------------------------------------------------------- -//------------------------------------------------------------------- -void vvImage::SetRotateX(int xvalue) -{ - transform->PostMultiply(); - transform->Translate(-origin[0],-origin[1],-origin[2]); - transform->RotateX(xvalue); - transform->Translate(origin[0],origin[1],origin[2]); - mVtkImageReslice->Update(); -} - -//-------------------------------------------------------------------- -//------------------------------------------------------------------- -void vvImage::SetRotateY(int yvalue) -{ - transform->PostMultiply(); - transform->Translate(-origin[0],-origin[1],-origin[2]); - transform->RotateY(yvalue); - transform->Translate(origin[0],origin[1],origin[2]); - mVtkImageReslice->Update(); -} //-------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetRotateZ(int zvalue) +bool vvImage::HaveSameSizeAndSpacingThan(vvImage * other) { - transform->PostMultiply(); - transform->Translate(-origin[0],-origin[1],-origin[2]); - transform->RotateZ(zvalue); - transform->Translate(origin[0],origin[1],origin[2]); - mVtkImageReslice->Update(); -} -//-------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetTranslationX(int xvalue) -{ - transform->Translate(xvalue,0,0); - mVtkImageReslice->Update(); -} -//-------------------------------------------------------------------- - -//------------------------------------------------------------------- - -void vvImage::SetTranslationY(int yvalue) -{ - transform->Translate(0,yvalue,0); - mVtkImageReslice->Update(); + bool same = true; + for(int i=0; iGetSize()[i]) same = false; + if ( GetSpacing()[i] != other->GetSpacing()[i]) same = false; + } + return same; } //-------------------------------------------------------------------- -//------------------------------------------------------------------- -void vvImage::SetTranslationZ(int zvalue) -{ - transform->Translate(0,0,zvalue); - mVtkImageReslice->Update(); -} -//------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetOrigin(double value[3]) -{ - origin=new double[mVtkImageReslice->GetOutputDimensionality()]; - origin[0]=value[0]; - origin[1]=value[1]; - origin[2]=value[2]; -} - -//-------------------------------------------------------------------- -void vvImage::SetTransform(vtkAbstractTransform *transform) -{ - mVtkImageReslice->SetResliceTransform(transform); - mVtkImageReslice->Update(); -} -//-------------------------------------------------------------------- #endif // VVIMAGE_CXX