src/creaImageIOSimpleView.cpp

   1 /*
   2 # ---------------------------------------------------------------------
   3 #
   4 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
   5 #                        pour la Santé)
   6 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
   7 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
   8 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
   9 #
  10 #  This software is governed by the CeCILL-B license under French law and
  11 #  abiding by the rules of distribution of free software. You can  use,
  12 #  modify and/ or redistribute the software under the terms of the CeCILL-B
  13 #  license as circulated by CEA, CNRS and INRIA at the following URL
  14 #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  15 #  or in the file LICENSE.txt.
  16 #
  17 #  As a counterpart to the access to the source code and  rights to copy,
  18 #  modify and redistribute granted by the license, users are provided only
  19 #  with a limited warranty  and the software's author,  the holder of the
  20 #  economic rights,  and the successive licensors  have only  limited
  21 #  liability.
  22 #
  23 #  The fact that you are presently reading this means that you have had
  24 #  knowledge of the CeCILL-B license and that you accept its terms.
  25 # ------------------------------------------------------------------------
  26 */
  27
  28 #include "boost/filesystem/operations.hpp"
  29 #include "boost/filesystem/fstream.hpp"
  30
  31 #include "creaImageIOSimpleView.h"
  32 #include "creaImageIOUnMosaicVtkImageData.h"
  33
  34 namespace creaImageIO
  35 {
  36         bool SimpleView::readFile(      std::vector<std::string> i_filenames,                                                                   // in information names
  37                                                                                 std::vector<std::string> i_attr,                                                                                        // in information attributes names
  38                                                                                 std::vector<std::map <std::string , std::string> > &i_imgAttr,  // out information attributes names-values
  39                                                                                 std::vector<vtkImageData *> &i_img)                                                                     // out information vtkImageData
  40         {
  41 printf("EED SimpleView::readFile 0\n");
  42
  43                   //Mosaic Attribute for Simens
  44                         i_attr.push_back("D0019_100a");  // Number of images in mosaic (Simens)
  45                         i_attr.push_back("D0018_0088");  // Spacing Between Slices
  46
  47
  48
  49                         creaImageIOUnMosaicVtkImageData Mosaic;
  50                         vtkImageData* tmpImg1;
  51                         vtkImageData* tmpImg2;
  52                         int numberOfImagesInMosaic;
  53                         double spc[3];
  54                         double spcZ;
  55
  56                         bool bresult, bfinal = true;
  57                         ImageReader *mReader = new ImageReader();
  58                         std::vector<std::string>::iterator it = i_filenames.begin();
  59                         for (; it != i_filenames.end(); it++)
  60                         {
  61                                 bresult = mReader->CanRead((*it).c_str());
  62                                 if(bresult)
  63                                 {
  64                                         std::map <std::string , std::string> mapAttr;
  65                                         mReader->getAttributes( (*it).c_str() , mapAttr, i_attr ); //  filename, outMapAttr-NameValue, inVectAttributeName
  66                                         i_imgAttr.push_back( mapAttr );
  67
  68 //EED borrame
  69 //                                      i_img.push_back( mReader->ReadImage( (*it).c_str() ) );
  70
  71 printf("EED SimpleView::readFile 1\n");
  72
  73                                         // EED 18 octt 2014
  74                                         std::map <std::string , std::string>::iterator itMapAttr;
  75                                         std::string numberOfImagesInMosaicStr="";
  76                                         std::string spaceBetweenSliceStr="";
  77
  78                                         itMapAttr = mapAttr.find("D0019_100a");
  79                                         if ( itMapAttr!=mapAttr.end() ) numberOfImagesInMosaicStr = itMapAttr->second;
  80
  81 //                                      std::string numberOfImagesInMosaicStr = mapAttr.find("D0019_100a")->second;
  82 printf("EED SimpleView::readFile 1.1\n");
  83                                         itMapAttr = mapAttr.find("D0018_0088");
  84                                         if ( itMapAttr!=mapAttr.end() ) spaceBetweenSliceStr = itMapAttr->second;
  85
  86 //                                      std::string spaceBetweenSliceStr = mapAttr.find("D0018_0088")->second;
  87 printf("EED SimpleView::readFile 1.2\n");
  88                                         if (numberOfImagesInMosaicStr=="")
  89                                         {
  90 printf("EED SimpleView::readFile 2\n");
  91                                                 i_img.push_back( mReader->ReadImage( (*it).c_str() ) );
  92                                         } else {
  93 printf("EED SimpleView::readFile 3\n");
  94                                                 numberOfImagesInMosaic = atoi(numberOfImagesInMosaicStr.c_str());
  95                                                 spcZ = atoi(spaceBetweenSliceStr.c_str());
  96                                                 tmpImg1 = mReader->ReadImage( (*it).c_str() ) ;
  97                                                 tmpImg1->GetSpacing( spc );
  98                                                 tmpImg2 = Mosaic.unMosaic( tmpImg1 , numberOfImagesInMosaic );
  99                                                 tmpImg2->SetSpacing( spc[0], spc[1], spcZ );
 100                                                 i_img.push_back( tmpImg2 );
 101                                         }// elemet
 102
 103                                 } else {
 104                                         printf("ERROR. Impossible to read file %s\n", (*it).c_str() );
 105                                         bfinal = false;
 106                                 } // if
 107                         } // for
 108                         delete mReader;
 109 printf("EED SimpleView::readFile 4\n");
 110
 111                         return bfinal;
 112         }
 113
 114
 115         bool SimpleView::readDirectory( const std::string i_pathname,                                                                                   // in information names
 116                                                                                                 std::vector<std::string> i_attr,                                                                                        // in information attributes names
 117                                                                                                 std::vector<std::map <std::string , std::string> > &i_imgAttr,  // out information attributes names-values
 118                                                                                                 std::vector<vtkImageData *> &i_imgs)                                                                    // out information vtkImageData
 119         {
 120                         bool bresult = true;
 121                         ImageReader *mReader = new ImageReader();
 122                         std::vector<std::string> names;
 123                         bresult = boost::filesystem::exists( i_pathname );
 124                         if (bresult)
 125                         {
 126                                 boost::filesystem::directory_iterator itr(i_pathname);
 127                                 boost::filesystem::directory_iterator end_itr;
 128                                 for(;itr != end_itr; ++itr)
 129                                 {
 130                                         if (!boost::filesystem::is_directory(itr->status()))
 131                                         {
 132                                                 if( mReader->CanRead(itr->path().string()) )
 133                                                 {
 134                                                         names.push_back(itr->path().string());
 135                                                 }
 136                                         }
 137                                 }
 138                                 std::sort (names.begin(), names.end()); // make sure names are in lexicographical order
 139                                 int lgr = (int)names.size();
 140
 141                                 for(int i=0; i<lgr; i++)
 142                                 {
 143                std::cout << names[i] << std::endl;
 144
 145                                         std::map <std::string , std::string> mapAttr;
 146                                         mReader->getAttributes( names[i] , mapAttr, i_attr ); //  filename, outMapAttr-NameValue, inVectAttributeName
 147                                         i_imgAttr.push_back( mapAttr );
 148
 149                                         //UnMosaic step..
 150
 151                                    i_imgs.push_back( mReader->ReadImage(names[i]) );
 152                                 }
 153                         }
 154                         return bresult;
 155         }
 156 }