X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=f75d3212267441b0ee75e2a0ffe221d00ad9066a;hb=d55f025b18f68066a52b8f33c2dc6481e82c2580;hp=1d718bce26d6fb5d24498d02bbe01e5f3bddf5dc;hpb=039781919d2f8aded60740ec0e1f283bc5d410bb;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index 1d718bc..f75d321 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,110 +14,117 @@ - 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" -#include "vtkImageData.h" + +// clitk #include "clitkCommon.h" + +// vtk +#include +#include +#include +#include +#include +#include + +// std #include //-------------------------------------------------------------------- -vvImage::vvImage() { - mVtkImages.resize(0); - mTimeSpacing = 1; - mTimeOrigin = 0; +vvImage::vvImage() +{ + Init(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -vvImage::~vvImage() { - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - if (mVtkImages[i] != NULL) - mVtkImages[i]->Delete(); - } +void vvImage::Init() +{ + mTimeSpacing = 1; + mTimeOrigin = 0; + mImageDimension = 0; + mrange[0] = std::numeric_limits::max();;//min + mrange[1] = std::numeric_limits::min();;//max } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -void vvImage::SetImage(std::vector images) { - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - if (mVtkImages[i] != NULL) - mVtkImages[i]->Delete(); - } - mVtkImages.resize(0); - for (unsigned int i = 0; i < images.size(); i++) { - mVtkImages.push_back(images[i]); - } +vvImage::~vvImage() +{ + Reset(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -void vvImage::AddImage(vtkImageData* image) { - mVtkImages.push_back(image); +void vvImage::Reset() +{ + mVtkImages.clear(); + Init(); } //-------------------------------------------------------------------- //-------------------------------------------------------------------- -void vvImage::Init() { - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - if (mVtkImages[i] != NULL) - mVtkImages[i]->Delete(); - } - mVtkImages.resize(0); +void vvImage::AddVtkImage(vtkImageData* input) +{ + // RP: 20/12/2011 + // Note that we're simply adding a new image to the vector. + // mItkToVtkConverters is therefore not being updated, but + // up to here it's not being used anyway... + mImageDimension = 0; +#if VTK_MAJOR_VERSION <= 5 + int* extent = input->GetWholeExtent(); +#else + int* extent = input->GetInformation()->Get(vtkDataObject::DATA_EXTENT()); +#endif + if (extent[4] != extent[5]) + mImageDimension = 3; + else if (extent[3] != extent[4]) + mImageDimension = 2; + else if (extent[0] != extent[1]) + mImageDimension = 1; + + mVtkImages.push_back(input); } -//-------------------------------------------------------------------- +//-------------------------------------------------------------------- //-------------------------------------------------------------------- -int vvImage::GetNumberOfSpatialDimensions() { - int dim=GetNumberOfDimensions(); - if (IsTimeSequence()) - return dim-1; - else - return dim; +int vvImage::GetNumberOfSpatialDimensions() +{ + 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; - } +int vvImage::GetNumberOfDimensions() const +{ + if (mVtkImages.size()) { + if (IsTimeSequence()) + return mImageDimension+1; + else + return mImageDimension; + } 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]; - } + range[0]=mrange[0]; + range[1]=mrange[1]; } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -std::string vvImage::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"; @@ -127,74 +134,69 @@ std::string vvImage::GetScalarTypeAsString() { } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -int vvImage::GetNumberOfScalarComponents() { +int vvImage::GetNumberOfScalarComponents() +{ return mVtkImages[0]->GetNumberOfScalarComponents(); } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -int vvImage::GetScalarSize() { +int vvImage::GetScalarSize() +{ return mVtkImages[0]->GetScalarSize(); } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -std::vector vvImage::GetSpacing() { +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]); - } + 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 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]); - } + 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 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]); - } + 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 vvImage::GetActualMemorySize() +{ unsigned long size = 0; for (unsigned int i = 0; i < mVtkImages.size(); i++) { - size += mVtkImages[i]->GetActualMemorySize(); + size += mVtkImages[i]->GetActualMemorySize()*1024; } return size; } @@ -202,37 +204,39 @@ unsigned long vvImage::GetActualMemorySize() { //-------------------------------------------------------------------- -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); } //-------------------------------------------------------------------- - //-------------------------------------------------------------------- -bool vvImage::IsScalarTypeInteger(int t) { +bool vvImage::IsScalarTypeInteger(int t) +{ if ((t == VTK_BIT) || (t == VTK_CHAR) || (t == VTK_UNSIGNED_CHAR) || @@ -243,12 +247,77 @@ bool vvImage::IsScalarTypeInteger(int t) { (t == VTK_LONG) || (t == VTK_UNSIGNED_LONG)) { return true; - } - else { + } else { return false; } } //-------------------------------------------------------------------- +//-------------------------------------------------------------------- +const std::vector< vtkSmartPointer >& vvImage::GetTransform() +{ + return this->mTransform; +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +void vvImage::InitializeTransform() +{ + for (int i = 0; i < mTransform.size(); i++) + mTransform[i]->Identity(); +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +std::vector< std::vector > vvImage::GetDirection() +{ + int dim = this->GetNumberOfDimensions(); + vtkSmartPointer transform = vtkSmartPointer::New(); + transform->Identity(); + for (int i = 0; i < mTransform.size(); i++) + transform->Concatenate(this->mTransform[i]); + + vtkSmartPointer matrix = transform->GetMatrix(); + matrix->Invert(); + std::vector > direction0; + for (int i = 0; i < dim; i++) { + if (i != 3) { + std::vector direction1; + for (int j = 0; j < dim; j++) { + if (j != 3) +#if VTK_MAJOR_VERSION <= 6 + direction1.push_back((*matrix)[i][j]); +#else + direction1.push_back((*matrix).GetElement(i,j)); +#endif + } + direction0.push_back(direction1); + } + } + return direction0; +} +//-------------------------------------------------------------------- + + + +//-------------------------------------------------------------------- +bool vvImage::HaveSameSizeAndSpacingThan(vvImage * other) +{ + bool same = true; + for(int i=0; iGetSize()[i]) same = false; + if ( GetSpacing()[i] != other->GetSpacing()[i]) same = false; + } + return same; +} +//-------------------------------------------------------------------- +itk::MetaDataDictionary* vvImage::GetFirstMetaDataDictionary() +{ + return mDictionary[0]; +} +//-------------------------------------------------------------------- + #endif // VVIMAGE_CXX