X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=9b9ca8289956b832b19e9dca0b6b47300208c6e1;hb=78ae141332b4f8d5146952b08dced58e4d248e06;hp=7d68a1438e70e4656a430d153ec0145b247225d2;hpb=758de1c14c7f16d924e1a2a575d52cd3e48ec503;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index 7d68a14..9b9ca82 100644 --- a/common/vvImage.cxx +++ b/common/vvImage.cxx @@ -1,189 +1,252 @@ -#ifndef VVIMAGE_CXX -#define VVIMAGE_CXX - /*========================================================================= + Program: vv http://www.creatis.insa-lyon.fr/rio/vv -Program: vv -Language: C++ -Author : Pierre Seroul (pierre.seroul@gmail.com) - -Copyright (C) 2008 -Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr -CREATIS-LRMN http://www.creatis.insa-lyon.fr - -This program is free software: you can redistribute it and/or modify -it under the terms of the GNU General Public License as published by -the Free Software Foundation, version 3 of the License. - -This program is distributed in the hope that it will be useful, -but WITHOUT ANY WARRANTY; without even the implied warranty of -MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the -GNU General Public License for more details. + Authors belong to: + - University of LYON http://www.universite-lyon.fr/ + - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr + - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr -You should have received a copy of the GNU General Public License -along with this program. If not, see . + This software is distributed WITHOUT ANY WARRANTY; without even + the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR + PURPOSE. See the copyright notices for more information. -=========================================================================*/ + It is distributed under dual licence + - 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 #include "vvImage.h" #include "vtkImageData.h" #include "clitkCommon.h" #include -//==================================================================== +//-------------------------------------------------------------------- vvImage::vvImage() { - mVtkImages.resize(0); + mVtkImages.resize(0); } -//==================================================================== +//-------------------------------------------------------------------- + -//==================================================================== +//-------------------------------------------------------------------- vvImage::~vvImage() { - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - if (mVtkImages[i] != NULL) - mVtkImages[i]->Delete(); - } + for (unsigned int i = 0; i < mVtkImages.size(); i++) { + if (mVtkImages[i] != NULL) + mVtkImages[i]->Delete(); + } } -//==================================================================== +//-------------------------------------------------------------------- + -//==================================================================== +//-------------------------------------------------------------------- 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]); - } + 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]); + } } -//==================================================================== +//-------------------------------------------------------------------- -//==================================================================== +//-------------------------------------------------------------------- void vvImage::AddImage(vtkImageData* image) { - mVtkImages.push_back(image); + mVtkImages.push_back(image); } -//==================================================================== +//-------------------------------------------------------------------- -//==================================================================== +//-------------------------------------------------------------------- void vvImage::Init() { - 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 < mVtkImages.size(); i++) { + if (mVtkImages[i] != NULL) + mVtkImages[i]->Delete(); + } + mVtkImages.resize(0); } -//==================================================================== +//-------------------------------------------------------------------- -//==================================================================== + +//-------------------------------------------------------------------- int vvImage::GetNumberOfSpatialDimensions() { - int dim=GetNumberOfDimensions(); - if (IsTimeSequence()) - return dim-1; - else - return dim; -} -//==================================================================== -// -//==================================================================== + int dim=GetNumberOfDimensions(); + if (IsTimeSequence()) + return dim-1; + else + return dim; +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- int vvImage::GetNumberOfDimensions() const { - if (mVtkImages.size()) + 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 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; + 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(); + 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]; + temp = mVtkImages[i]->GetScalarRange(); + 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(); + // 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(); + return mVtkImages[0]->GetNumberOfScalarComponents(); } -//==================================================================== +//-------------------------------------------------------------------- + -//==================================================================== +//-------------------------------------------------------------------- int vvImage::GetScalarSize() { - return mVtkImages[0]->GetScalarSize(); + return mVtkImages[0]->GetScalarSize(); } -//==================================================================== +//-------------------------------------------------------------------- + -//==================================================================== +//-------------------------------------------------------------------- std::vector vvImage::GetSpacing() { - std::vector spacing; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) + 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(1); + else + spacing.push_back(mVtkImages[0]->GetSpacing()[i]); } - return spacing; + return spacing; } -//==================================================================== +//-------------------------------------------------------------------- -//==================================================================== + +//-------------------------------------------------------------------- std::vector vvImage::GetOrigin() const { - std::vector origin; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) + 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(0); + else + origin.push_back(mVtkImages[0]->GetOrigin()[i]); } - return origin; + return origin; } -//==================================================================== +//-------------------------------------------------------------------- + -//==================================================================== +//-------------------------------------------------------------------- std::vector vvImage::GetSize() { - std::vector size0; - int dim = this->GetNumberOfDimensions(); - for (int i = 0; i < dim; i++) + 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]); + if (i == 3) + size0.push_back(mVtkImages.size()); + else + size0.push_back(mVtkImages[0]->GetDimensions()[i]); } - return size0; + 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 size = 0; + for (unsigned int i = 0; i < mVtkImages.size(); i++) { + size += mVtkImages[i]->GetActualMemorySize(); + } + return size; } -//==================================================================== +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +bool vvImage::IsTimeSequence() { + return mVtkImages.size()>1; +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +const std::vector& vvImage::GetVTKImages() { + return mVtkImages; +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +vtkImageData* vvImage::GetFirstVTKImageData() { + return mVtkImages[0]; +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +bool vvImage::IsScalarTypeInteger() { + assert(mVtkImages.size()> 0); + int t = mVtkImages[0]->GetScalarType(); + return IsScalarTypeInteger(t); +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +bool vvImage::IsScalarTypeInteger(int t) { + if ((t == VTK_BIT) || + (t == VTK_CHAR) || + (t == VTK_UNSIGNED_CHAR) || + (t == VTK_SHORT) || + (t == VTK_UNSIGNED_SHORT) || + (t == VTK_INT) || + (t == VTK_UNSIGNED_INT) || + (t == VTK_LONG) || + (t == VTK_UNSIGNED_LONG)) { + return true; + } + else { + return false; + } +} +//-------------------------------------------------------------------- + #endif // VVIMAGE_CXX