/*========================================================================= Program: bbtk Module: $RCSfile: bbvtkConcatImages.cxx,v $ Language: C++ Date: $Date: 2010/04/08 14:39:45 $ Version: $Revision: 1.1 $ =========================================================================*/ /* --------------------------------------------------------------------- * Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale) * Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux * * This software is governed by the CeCILL-B license under French law and * abiding by the rules of distribution of free software. You can use, * modify and/ or redistribute the software under the terms of the CeCILL-B * license as circulated by CEA, CNRS and INRIA at the following URL * http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html * or in the file LICENSE.txt. * * As a counterpart to the access to the source code and rights to copy, * modify and redistribute granted by the license, users are provided only * with a limited warranty and the software's author, the holder of the * economic rights, and the successive licensors have only limited * liability. * * The fact that you are presently reading this means that you have had * knowledge of the CeCILL-B license and that you accept its terms. * ------------------------------------------------------------------------ */ /** * \file * \brief */ #ifdef _USE_VTK_ #include "bbvtkConcatImages.h" #include "bbvtkPackage.h" namespace bbvtk { BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,ConcatImages) BBTK_BLACK_BOX_IMPLEMENTATION(ConcatImages,bbtk::AtomicBlackBox); //--------------------------------------------------------------------- void ConcatImages::bbUserSetDefaultValues() { //std::cout << "-------- entree ds ConcatImages::bbUserSetDefaultValues()\n" << std::endl; // bbSetInputIn(NULL); /// \TODO fix it! IN is a std::vector // JPR mConcat = NULL; bbSetOutputOut(NULL); } //--------------------------------------------------------------------- void ConcatImages::bbUserInitializeProcessing() { //std::cout << "-------- entree ds ConcatImages::bbUserInitalizeProcessing()\n" << std::endl; //bbUserFinalizeProcessing(); mConcat = vtkImageData::New(); // Alloc depends on bbGetInputIn().size() } //--------------------------------------------------------------------- void ConcatImages::bbUserFinalizeProcessing() { //std::cout << "-------- entree ds ConcatImages::bbUserFinalizeProcessing()\n" << std::endl; // WTF? we never enter here // JPR bbUserFinalizeProcessing() JPR if (mConcat!=NULL) { // mConcat->Delete(); // mConcat=NULL; } bbSetOutputOut(mConcat); } //--------------------------------------------------------------------- /// \TODO :clean this dirty stuff : /// - receives a std::vector /// - exports a vtkImageData* /// - supposes *all* the images are consistent (same type, same number of components, same pixel type) /// - if images in the vector are already 3D, exports only the first one. ?!? // JPR /// void ConcatImages::Process() { int dim[3]; int nb = (unsigned int)bbGetInputIn().size(); if (nb == 0) {} // ??? JPR bbGetInputIn()[0]->GetDimensions(dim); if (dim[2] > 1 || nb == 1) // dim[2] > 1 : Hopeless for vtk : the first file contains already a 3D object // nb == 1 : only one image, nothing to do. { mConcat->Delete(); mConcat = bbGetInputIn()[0]; //mConcat->PrintSelf(std::cout, vtkIndent(2)); } else { // We suppose *all* the images within the directory are consistent (same sizeS, same pixel Type?...) vtkImageData * firstImage = bbGetInputIn()[0]; //std::cout << "--------PrintSelf firstImage " << std::endl; // firstImage->PrintSelf(std::cout, vtkIndent(2)); int nbComponents = firstImage->GetNumberOfScalarComponents(); mConcat->SetNumberOfScalarComponents(nbComponents); mConcat->SetScalarType(firstImage->GetScalarType( )); firstImage->GetDimensions(dim); dim[2]=nb; mConcat->SetDimensions(dim); mConcat->SetExtent(0, dim[0]-1, 0, dim[1]-1, 0, dim[2]-1); mConcat->SetSpacing(firstImage->GetSpacing()); mConcat->AllocateScalars(); mConcat->Update(); int index_image, index_lignes, index_colonnes, index_components; for(index_image=0; index_imageSetScalarComponentFromDouble(index_lignes, index_colonnes, index_image, index_components,(bbGetInputIn()[index_image])->GetScalarComponentAsDouble(index_colonnes,index_lignes,0,index_components)); } } } } } // Devrait etre dans bbUserFinalizeProcessing() ? On n'y entre jamais // JPR bbSetOutputOut(mConcat); } }//namespace bbtk #endif // _USE_VTK_