X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=9a6fc8ff489f39c1d661b1a55940a46253d3a2c3;hb=e1098dccba022441144bfc1c44a2b58db7dde34e;hp=a394c7608456b9d0de2f06f88ceeb3a7bfbd993c;hpb=a97c480d25285a6d143bd035d6b5fa0804a7d068;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index a394c76..9a6fc8f 100644 --- a/common/vvImage.cxx +++ b/common/vvImage.cxx @@ -3,7 +3,7 @@ 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,18 +14,25 @@ - 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(); } @@ -37,9 +44,7 @@ void vvImage::Init() { mTimeSpacing = 1; mTimeOrigin = 0; - - if (CLITK_EXPERIMENTAL) - mTransform = vtkSmartPointer::New(); + mImageDimension = 0; } //-------------------------------------------------------------------- @@ -53,76 +58,31 @@ vvImage::~vvImage() //-------------------------------------------------------------------- void vvImage::Reset() { - if (CLITK_EXPERIMENTAL) - for (unsigned int i = 0; i < mVtkImages.size(); i++) - mVtkImageReslice[i]->GetInput(0)->Delete(); - else - for (unsigned int i = 0; i < mVtkImages.size(); i++) - mVtkImages[i]->Delete(); - - mVtkImages.resize(0); - mVtkImageReslice.resize(0); - + mVtkImages.clear(); Init(); } //-------------------------------------------------------------------- - -//-------------------------------------------------------------------- -void vvImage::SetImage(std::vector< vtkImageData* > images) -{ - Reset(); - for (unsigned int i = 0; i < images.size(); i++) - AddImage(images[i]); -} -//-------------------------------------------------------------------- - - -//-------------------------------------------------------------------- -void vvImage::AddImage(vtkImageData* image) -{ - if (CLITK_EXPERIMENTAL) { - mVtkImageReslice.push_back(vtkSmartPointer::New()); - - mVtkImageReslice.back()->SetInterpolationModeToLinear(); - mVtkImageReslice.back()->AutoCropOutputOn(); - mVtkImageReslice.back()->SetBackgroundColor(-1000,-1000,-1000,1); - mVtkImageReslice.back()->SetResliceTransform(mTransform); - mVtkImageReslice.back()->SetInput(0, image); - mVtkImageReslice.back()->Update(); - - mVtkImages.push_back( mVtkImageReslice.back()->GetOutput(0) ); - } else - mVtkImages.push_back(image); -} -//-------------------------------------------------------------------- - //-------------------------------------------------------------------- int vvImage::GetNumberOfSpatialDimensions() { - int dim=GetNumberOfDimensions(); - if (IsTimeSequence()) - return dim-1; - else - return dim; + return mImageDimension; } //-------------------------------------------------------------------- -// + //-------------------------------------------------------------------- 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; + if (IsTimeSequence()) + return mImageDimension+1; + else + return mImageDimension; } return 0; } +//-------------------------------------------------------------------- + //-------------------------------------------------------------------- void vvImage::GetScalarRange(double* range) { @@ -136,9 +96,10 @@ void vvImage::GetScalarRange(double* range) if (temp[1] > range[1]) range[1]=temp[1]; } } +//-------------------------------------------------------------------- //-------------------------------------------------------------------- -std::string vvImage::GetScalarTypeAsString() +std::string vvImage::GetScalarTypeAsITKString() { // WARNING VTK pixel type different from ITK Pixel type std::string vtktype = mVtkImages[0]->GetScalarTypeAsString(); @@ -176,8 +137,6 @@ std::vector vvImage::GetSpacing() } return spacing; } -//-------------------------------------------------------------------- - //-------------------------------------------------------------------- std::vector vvImage::GetOrigin() const { @@ -221,7 +180,7 @@ unsigned long vvImage::GetActualMemorySize() //-------------------------------------------------------------------- -bool vvImage::IsTimeSequence() +bool vvImage::IsTimeSequence() const { return mVtkImages.size()>1; } @@ -235,7 +194,6 @@ const std::vector& vvImage::GetVTKImages() } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- vtkImageData* vvImage::GetFirstVTKImageData() { @@ -243,7 +201,6 @@ vtkImageData* vvImage::GetFirstVTKImageData() } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- bool vvImage::IsScalarTypeInteger() { @@ -270,82 +227,27 @@ bool vvImage::IsScalarTypeInteger(int t) return false; } } - //-------------------------------------------------------------------- //-------------------------------------------------------------------- -vtkTransform * vvImage::GetTransform() +vtkSmartPointer vvImage::GetTransform() { - return mTransform; + return this->mTransform; } //-------------------------------------------------------------------- -//------------------------------------------------------------------- -void vvImage::SetRotateX(int xvalue) -{ - mTransform->PostMultiply(); - mTransform->Translate(-origin[0],-origin[1],-origin[2]); - mTransform->RotateX(xvalue); - mTransform->Translate(origin[0],origin[1],origin[2]); -} //-------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetRotateY(int yvalue) +bool vvImage::HaveSameSizeAndSpacingThan(vvImage * other) { - mTransform->PostMultiply(); - mTransform->Translate(-origin[0],-origin[1],-origin[2]); - mTransform->RotateY(yvalue); - mTransform->Translate(origin[0],origin[1],origin[2]); -} -//-------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetRotateZ(int zvalue) -{ - mTransform->PostMultiply(); - mTransform->Translate(-origin[0],-origin[1],-origin[2]); - mTransform->RotateZ(zvalue); - mTransform->Translate(origin[0],origin[1],origin[2]); -} -//-------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetTranslationX(int xvalue) -{ - mTransform->Translate(xvalue,0,0); -} -//-------------------------------------------------------------------- - -//------------------------------------------------------------------- - -void vvImage::SetTranslationY(int yvalue) -{ - mTransform->Translate(0,yvalue,0); + 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) -{ - mTransform->Translate(0,0,zvalue); -} -//------------------------------------------------------------------- - -//------------------------------------------------------------------- -void vvImage::SetOrigin(double value[3]) -{ - origin[0]=value[0]; - origin[1]=value[1]; - origin[2]=value[2]; -} - -//-------------------------------------------------------------------- -void vvImage::SetTransform(vtkTransform *transform) -{ - mTransform = transform; -} -//-------------------------------------------------------------------- #endif // VVIMAGE_CXX