X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=821de16ebb35c34067b2488ba93be5161ba9599f;hb=4c78eb9adeaad50d0885769f21fecdfd4804f0c0;hp=c0c3aa61c683a734854a21b90c17024bfa73b191;hpb=931a42358442f4ee4f314613c991c838d4b4e3b7;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index c0c3aa6..821de16 100644 --- a/common/vvImage.cxx +++ b/common/vvImage.cxx @@ -1,192 +1,309 @@ +/*========================================================================= + Program: vv http://www.creatis.insa-lyon.fr/rio/vv + + 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 + + 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 "clitkCommon.h" +#include +#include +#include -/*========================================================================= +#include +#include +#include -Program: vv -Module: $RCSfile: vvImage.cxx,v $ -Language: C++ -Date: $Date: 2010/01/06 13:32:01 $ -Version: $Revision: 1.1 $ -Author : Pierre Seroul (pierre.seroul@gmail.com) +//-------------------------------------------------------------------- +vvImage::vvImage():mGeneralTransform(vtkSmartPointer::New()), + mTransform (vtkSmartPointer::New()), + mGridTransform (vtkSmartPointer::New()) +{ + mGeneralTransform->Concatenate(mTransform); + mGeneralTransform->Concatenate(mGridTransform); + Init(); +} +//-------------------------------------------------------------------- -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. +//-------------------------------------------------------------------- +void vvImage::Init() +{ + mTimeSpacing = 1; + mTimeOrigin = 0; +} +//-------------------------------------------------------------------- -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. +//-------------------------------------------------------------------- +vvImage::~vvImage() +{ + DD("vvImage got deleted!"); + Reset(); +} +//-------------------------------------------------------------------- -You should have received a copy of the GNU General Public License -along with this program. If not, see . +//-------------------------------------------------------------------- +void vvImage::Reset() +{ + for (unsigned int i = 0; i < mVtkImages.size(); i++) + mVtkImageReslice[i]->GetInput()->Delete(); -=========================================================================*/ + mVtkImages.resize(0); + mVtkImageReslice.resize(0); -#include "vvImage.h" -#include "vtkImageData.h" -#include "clitkCommon.h" -#include + Init(); +} +//-------------------------------------------------------------------- -//==================================================================== -vvImage::vvImage() { - mVtkImages.resize(0); -} -//==================================================================== - -//==================================================================== -vvImage::~vvImage() { - 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]); - } -} -//==================================================================== - -//==================================================================== -void vvImage::AddImage(vtkImageData* 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); -} -//==================================================================== - -//==================================================================== -int vvImage::GetNumberOfSpatialDimensions() { - int dim=GetNumberOfDimensions(); - if (IsTimeSequence()) - return dim-1; - else - return dim; -} -//==================================================================== -// -//==================================================================== -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::SetImage(std::vector< vtkImageData* > images) +{ + Reset(); + for (unsigned int i = 0; i < images.size(); i++) + AddImage(images[i]); +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +void vvImage::AddImage(vtkImageData* image) +{ + mVtkImageReslice.push_back(vtkSmartPointer::New()); + + mVtkImageReslice.back()->SetInterpolationModeToLinear(); + mVtkImageReslice.back()->AutoCropOutputOn(); + mVtkImageReslice.back()->SetBackgroundColor(-1000,-1000,-1000,1); + mVtkImageReslice.back()->SetResliceTransform(mGeneralTransform); + mVtkImageReslice.back()->SetInput(0, image); + mVtkImageReslice.back()->Update(); + + mVtkImages.push_back( mVtkImageReslice.back()->GetOutput(0) ); +} +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +int vvImage::GetNumberOfSpatialDimensions() +{ + int dim=GetNumberOfDimensions(); + if (IsTimeSequence()) + return dim-1; + else + return dim; +} +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +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]; - } -} - -//==================================================================== -std::string vvImage::GetScalarTypeAsString() { - return mVtkImages[0]->GetScalarTypeAsString(); -} -//==================================================================== - -//==================================================================== -int vvImage::GetNumberOfScalarComponents() { - return mVtkImages[0]->GetNumberOfScalarComponents(); -} -//==================================================================== - -//==================================================================== -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]); - } - 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]); - } - 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; -} -//==================================================================== - -//==================================================================== -unsigned long vvImage::GetActualMemorySize() { - unsigned long size = 0; - for (unsigned int i = 0; i < mVtkImages.size(); i++) { - size += mVtkImages[i]->GetActualMemorySize(); - } - return size; -} -//==================================================================== + 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::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::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(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(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; +} +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +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() +{ + 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; + } +} + +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +vtkSmartPointer vvImage::GetTransform() +{ + return mTransform; +} +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +vtkSmartPointer vvImage::GetGridTransform() +{ + return mGridTransform; +} +//-------------------------------------------------------------------- + +//-------------------------------------------------------------------- +void vvImage::UpdateReslice() +{ + for (unsigned int i=0; iUpdate(); + } +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +vtkImageData * CopyAndCastToFloatFrom(vtkImageData * input) +{ + vtkSmartPointer p = vtkSmartPointer::New(); + p->SetExtent(input->GetExtent ()); // Only first ! could not be 4D + p->SetScalarTypeToFloat(); + p->AllocateScalars(); + p->CopyAndCastFrom(input, input->GetExtent()); + return p; +} +//-------------------------------------------------------------------- + #endif // VVIMAGE_CXX