creaImageIOGimmickView.cpp

Go to the documentation of this file.

#include <creaImageIOGimmickView.h>
#include <creaImageIOSystem.h>
#include "boost/filesystem.hpp"

#if defined(USE_GDCM)
#include <gdcmGlobal.h>
#include <gdcmSerieHelper.h>
#include <vtkGdcmReader.h>
#endif

#if defined(USE_GDCM2)
#include <vtkGDCMImageReader.h>
#include "gdcmSystem.h"
#include "gdcmCryptographicMessageSyntax.h"
#include "gdcmUIDGenerator.h"
#include "gdcmAnonymizer.h"
#include "gdcmGlobal.h"
#endif


namespace fs = boost::filesystem;
namespace creaImageIO
{

  class ImageExtent
  {
  public:
    ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
        {
                sscanf(x.c_str(),"%d",&mExtent[0]);
                sscanf(y.c_str(),"%d",&mExtent[1]);
                sscanf(z.c_str(),"%d",&mExtent[2]);
                sscanf(t.c_str(),"%d",&mExtent[3]);
                if(x==""){mExtent[0]=1;}
                if(y==""){mExtent[1]=1;}
                if(z==""){mExtent[2]=1;}
                if(t==""){mExtent[3]=1;}

                if (mExtent[3]>1) mDim=4;
                else if (mExtent[2]>1) mDim=3;
                else if (mExtent[1]>1) mDim=2;
                else if (mExtent[0]>1) mDim=1;
                else mDim=0;
        }
        

    void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }

    bool IsCompatible( const ImageExtent& );

    void Add ( const ImageExtent& );
                
    int Get(int i) { return mExtent[i]; }
    
    void SetDimension(int dim) { mDim=dim; }

    int GetDimension() { return mDim; }

  private:
    int mExtent[4];
    int mDim;
  };

  //======================================================================

  //======================================================================
  // CTor
  GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
    : mGimmick(gimmick),
          mReader(threads)
  {
    GimmickDebugMessage(1,"GimmickView::GimmickView"
                        <<std::endl);
        // Anciently started the threads ...
    // Threads now automatically start at first image request
    //mReader.Start();

  }
  //======================================================================

  //======================================================================
  GimmickView::~GimmickView()
  {
    GimmickDebugMessage(1,"GimmickView::~GimmickView"
                        <<std::endl);
  }
   //======================================================================
 
  //======================================================================  
  void GimmickView::Initialize()
  {
        mReaderStarted=false;
  }
  //======================================================================

  //======================================================================
  void GimmickView::Finalize()
  {
  }

  //======================================================================

  //======================================================================
  void GimmickView::CreateTreeViews()
  {
    GimmickMessage(2,"Creating the tree views"<<std::endl);
    Gimmick::TreeHandlerMapType::const_iterator i;
    for (i = mGimmick->GetTreeHandlerMap().begin();
         i!= mGimmick->GetTreeHandlerMap().end();
         ++i)
      {
        this->CreateTreeView(i->second);
      }
  }

  void GimmickView::CreateSingleTreeView(std::string &i_name)
  {
        this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
     
  }


  //======================================================================

  //======================================================================
  void GimmickView::UpdateTreeViewLevel(const std::string& t, int l)
  {
    TreeViewMapType::iterator i;
    i = GetTreeViewMap().find(t);
    if ( i == GetTreeViewMap().end() )
      {
        GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
                     <<t<<"' is not in TreeViewMap");
      }
    i->second->UpdateLevel(l);
  }

  //======================================================================
  void GimmickView::ResetExtent()
  {
          if(mImageExtent!=0)
          {
                  mImageExtent.reset();
          }
          valid=true;
  }

  //======================================================================

  //======================================================================
  bool ImageExtent::IsCompatible(const ImageExtent& ie)
  {
          bool compatible=true;
          ImageExtent * extent= (ImageExtent*)&ie;
          if((*extent).Get(0)!=Get(0)
                 || (*extent).Get(1)!=Get(1))
          {
                  compatible=false;
          }
          return compatible;
  }

  //======================================================================

  //======================================================================
  void ImageExtent::Add(const ImageExtent& ie)
  {
          ImageExtent * extent= (ImageExtent*)&ie;
          mExtent[2]+=(*extent).Get(2);
          if(mExtent[2]>1)
          {
          SetDimension(3);
          }
  }

  //======================================================================
  bool GimmickView::NoValidateSelected ()
  {
        GimmickDebugMessage(2,"Validating selected"<<std::endl);
        std::string mMessage;
        mMessage="Cannot have 0 images selected!";
        valid=false;
        modifyValidationSignal(valid);
        SetMessage(mMessage);
        return valid;
  }

  //======================================================================
  bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
  {
        GimmickDebugMessage(2,"Validating selected"<<std::endl);
        std::string mMessage;
        
        if(sel==0)
        {
                mMessage="Cannot have 0 images selected!";
                valid=false;
        }
        else
        {
                        boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
                                                                        (*sel).GetAttribute("D0028_0011"),
                                                                        (*sel).GetAttribute("D0028_0012"), 
                                                                        ""));

        if(mImageExtent==0)
        {
                mImageExtent=ie;
                if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
                {
                        valid=false;
                }
                else
                {
                        std::stringstream out;
                        out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
                        mMessage = out.str();
                        mImageExtent->SetDimension(2);
                        valid=true;
                }
        }
        else
        {
                if(mImageExtent->IsCompatible(*ie))
                {
                        if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
                        {
                                std::stringstream out;
                                out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
                                mMessage=out.str();
                                valid=false;
                        }
                        else if(max_dim<3)
                        {
                                std::stringstream out;
                                out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
                                mMessage=out.str();
                                valid=false;
                        }
                        else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
                        {
                                std::stringstream out;
                                out << "Cannot build the selection as it would result in a ";
                                out << mImageExtent->GetDimension();
                                out << "D image, and the minimum is ";
                                out << min_dim;
                                out << "D!";
                                mMessage=out.str();
                                valid=false;
                        }
                        else
                        {
                                mImageExtent->Add(*ie);
                                std::stringstream out;
                                out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
                                mMessage = out.str();
                        }
                }
                else
                {
                        mMessage="The selected images are not compatible.";
                        valid=false;
                }
          }
        }

        modifyValidationSignal(valid);
        SetMessage(mMessage);
        return valid;
  }

  //======================================================================
  void GimmickView::modifyValidationSignal(bool ivalid)
  {

          mValidationSignal(ivalid);
  }

  void GimmickView::stopReader()
  {
          mReader.Stop();
  }

   //======================================================================
void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
{
        stopReader();
/* remember!

#define GIMMICK_NO_IMAGE_SELECTION 0
#define GIMMICK_2D_IMAGE_SELECTION 2
#define GIMMICK_3D_IMAGE_SELECTION 3
#define GIMMICK_4D_IMAGE_SELECTION 4

#define NATIVE 0
#define _2D    2
#define _3D    3

*/
        // Create the output data
        if (im.size()==1) 
        {
                vtkImageData * out=vtkImageData::New();
                out->ShallowCopy(mReader.GetImage(im.front()));
                s.push_back(out);
        }
    else if (im.size()>1) // Test inutile ? JPR
        {
                vtkImageData* first = mReader.GetImage( im.front());
                if (dimension == 2)
                {
                 // n3D
                    std::vector<std::string>::iterator it;
                        for (it=im.begin(); it!=im.end(); ++it) 
                        {
                                vtkImageData* out = vtkImageData::New();
                                out->ShallowCopy(mReader.GetImage(*it));
                                s.push_back(out);
                        }
                }         
                else 
                {
                        // n*2D to 3D
                        vtkImageData* out = vtkImageData::New();
//                      out->CopyStructure(first);      
                        out->SetScalarType(first->GetScalarType());
                        out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
                        int ext[6];
                        //first->GetExtent(ext);  // JPR
                        first->GetWholeExtent(ext);  // renvoie egalement 0,0 en Z // JPR

                        if(ext[5] == 0)
                        {
                                ext[5] = im.size()-1;
                        }
                        else
                        {
                                ext[5] = ext[5] * im.size()-1; // to deal with multiframes - JPR
                        }
                        out->SetExtent(ext);

                        // LG : TODO : Z Spacing  ?

                        int dim[3];
                        first->GetDimensions(dim);

                        out->SetDimensions(dim[0], dim[1], im.size() );
                        out->AllocateScalars();
                        out->Update();

                        unsigned long imsize = dim[0] * dim[1];
                        imsize = imsize * dim[2] ;  // deal with multiframes // JPR


//EED 03-11-2009
                        // differents formats char , short, etc...
                        // differents components 1..3  ex. jpg ->RGB 3
                        imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();


                        // Order the file name vector

                        double spc[3];
                        first->GetSpacing(spc);

                        // OrderTheFileNameVector is not here anymore.
                        // Try orderFilesWithZSpacing from OutputModel FCY
                        // spc[2]=OrderTheFileNameVector(im);   
                        spc[2] =1;

                        out->SetSpacing(spc);

                        int slice = 0;
                        std::vector<std::string>::iterator it;

                        for (it=im.begin(); it!=im.end(); ++it) 
                        {
                                vtkImageData* cur = mReader.GetImage( (*it) );
                                memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
                                slice++;
                        }       
                        s.push_back(out);

                }  // dimension == 3

        }  // size >1

}
  //======================================================================



        // Test if the image is a multiple or single frame. 
        // For the moment only with the creation of vtkImageDta.
        // TO DO with image size and dim!!!
        bool GimmickView::isSingle(const std::string i_file)
        {
                bool bres = true;
                vtkImageData* first = mReader.GetImage( i_file);
                int dim[3];
                first->GetDimensions(dim);
                if (dim[2] > 1)
                {
                        bres = false;
                }
                else
                {
                }
                return bres;
        }

        // get Attributes values for a file
        
        void GimmickView::getAttributes(const std::string i_file, std::map<std::string, std::string> &o_infos, OutputAttr i_attr)
        {
                if(i_attr.inside.size() >0)
                {
                        mGimmick->GetAttributes(i_file,o_infos,i_attr);
                }
                if(i_attr.outside.size()>0)
                {
                        mReader.getAttributes(i_file,o_infos, i_attr.outside);
                }
        }

        // create an output structure with n entries = n output
        void GimmickView::readImages1(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
                         OutputAttr i_attr)
        {
                std::vector<std::string>::iterator it;
                for (it=im.begin(); it!=im.end(); ++it)
                {
                        OutStrGimmick out;
                        out.img = vtkImageData::New();
                        out.img->ShallowCopy(mReader.GetImage(*it));
                        if(i_attr.mult)
                                getAttributes((*it),out.infos,i_attr);
                        o_output.push_back(out);
                }
                // If we want only one output information structure, we set it outside the loop
                if(!i_attr.mult)
                {
                        getAttributes(im.front(), o_output.front().infos, i_attr);
                }
                
        }

        // create an output structure with n entries = 1 output
        void GimmickView::readImages3(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
                         OutputAttr i_attr, double i_zspc)
        {
                OutStrGimmick out;
                vtkImageData* first = mReader.GetImage( im.front());
                out.img  = vtkImageData::New();
                out.img->SetScalarType(first->GetScalarType());
                out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
                int ext[6];
                first->GetWholeExtent(ext);  // send also 0,0 in Z 
                if(ext[5] == 0)
                {
                        ext[5] = im.size()-1;
                }
                else
                {
                        ext[5] = ext[5] * im.size()-1; // to deal with multiframes 
                }
                out.img->SetExtent(ext);
                int dim[3];
                first->GetDimensions(dim);
                out.img->SetDimensions(dim[0], dim[1], im.size() );
                out.img->AllocateScalars();
                out.img->Update();
                unsigned long imsize = dim[0] * dim[1];
                imsize = imsize * dim[2] ;  // deal with multiframes here
                // differents formats char , short, etc...
                // differents components 1..3  ex. jpg ->RGB 3
                imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
                // Order the file name vector already did with the OutputModel!!!
                int slice = 0;
                std::vector<std::string>::iterator it;
                for (it=im.begin(); it!=im.end(); ++it) 
                {
                        vtkImageData* cur = mReader.GetImage( (*it) );
                        memcpy(out.img->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
                        slice++;
                }       
                getAttributes(im.front(),out.infos, i_attr);
                o_output.push_back(out);
        }


        // TO DO  NO VERY SURE : NEED TO BE TESTED
        // create an output structure with n entries (T size) = T output (n size)
        void GimmickView::readImages2(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
                         OutputAttr i_attr, double i_zspc)
        {
                vtkImageData* first = mReader.GetImage( im.front());
                int dim[3];
                first->GetDimensions(dim);
                // differents formats char , short, etc...
                // differents components 1..3  ex. jpg ->RGB 3
                unsigned long imsize = dim[0] * dim[1];
                imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();

                // Order the file name vector already did with the OutputModel!!!
                std::vector<std::string>::iterator it;
                std::vector<OutStrGimmick>::iterator it_out = o_output.begin();

                for (it=im.begin(); it!=im.end(); ++it, it_out += dim[2])
                {
                        vtkImageData* cur = mReader.GetImage( (*it) );
                        for (int slice= 0 ; slice <dim[2]; slice++)
                        {
                                OutStrGimmick out;
                                out.img = vtkImageData::New();
                                out.img->SetScalarType(first->GetScalarType());

                                out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
                                int ext[6];
                                first->GetWholeExtent(ext);  // send also 0,0 in Z 
                                ext[5] = 0;
                                out.img->SetExtent(ext);

                                out.img->SetDimensions(dim[0], dim[1], 1 );
                                out.img->AllocateScalars();
                                out.img->Update();
                                memcpy(out.img->GetScalarPointer(0,0,0), cur->GetScalarPointer(0,0,slice), imsize);
                                o_output.push_back(out);
                        }
                        if(i_attr.mult)
                                getAttributes((*it),(*it_out).infos,i_attr);
                }
                if(!i_attr.mult)
                {
                        getAttributes(im.front(), o_output.front().infos,i_attr);
                }

        }

        // create an output structure with n entries (T size) = T + n output
        void GimmickView::readImages4(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
                         OutputAttr i_attr)
        {
                std::vector<std::string>::iterator it;
                std::vector<OutStrGimmick>::iterator it_out = o_output.begin();
                vtkImageData* first = mReader.GetImage( im.front());
                int dim[3];
                first->GetDimensions(dim);
                        
                for (int slice= 0 ; slice <dim[2]; slice++)
                {
                        OutStrGimmick out;
                        out.img = vtkImageData::New();
                        out.img->SetScalarType(first->GetScalarType());
                        out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
                        
                        int ext[6];
                        first->GetWholeExtent(ext);  // send also 0,0 in Z 
                        ext[5] = 0;
                        out.img->SetExtent(ext);
                        
                        out.img->SetDimensions(dim[0], dim[1], im.size() );
                        out.img->AllocateScalars();
                        out.img->Update();
                        unsigned long imsize = dim[0] * dim[1];
                        imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
                        int index = 0;
        
                        for (it=im.begin(); it!=im.end(); ++it, index ++)
                        {
                                vtkImageData* cur = mReader.GetImage( (*it) );
                                memcpy(out.img->GetScalarPointer(0,0,index), cur->GetScalarPointer(0,0,slice), imsize);
                                o_output.push_back(out);
                        }
                }
                if(!i_attr.mult) // No sense to take informations in all images
                {
                        getAttributes(im.front(), o_output.front().infos,i_attr);
                }

        }


// Global function to read Images and create wished output (informations, multiple, type and size of output...)
// In function of type (file, vector) and size, the correct readImages function is selected
        void GimmickView::readImages(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
                         OutputAttr i_attr, int i_dim, double i_zspc)
        {
                int size = im.size();
                if ( size == 0)
                {
                        return;
                }
                else if (size == 1)
                {
                        // Simplest case
                        // Only one image : give it
                        // But take in count multiframe possibility
                        if ( isSingle(im.front()) || i_dim != 1)
                        {
                                 readImages1(o_output,im, i_attr);
                        }
                        else
                        {
                                readImages2(o_output,im, i_attr,i_zspc);
                        }
                        
                }
                else
                {
                        // multiple or single frame
                        if ( isSingle(im.front()) )
                        {
                                //we deal with 2D images
                                if(i_dim == 1)
                                {
                                        // 2D to 3D
                                        readImages3(o_output,im, i_attr,i_zspc);
                                }
                                else
                                {
                                        readImages1(o_output,im, i_attr);
                                }
                        }
                        else
                        {
                                // we deal with multiple frames n x (2D x T)
                                // Differents outputs are avaialable
                                if(i_dim == 1)
                                {
                                        // put all in one output
                                        readImages3(o_output,im, i_attr,i_zspc);

                                }
                                else if( i_dim == 2)
                                {
                                        // put in a vector of n x T (2D)
                                        readImages2(o_output,im, i_attr,i_zspc);
                                }
                                else if( i_dim == 3)
                                {
                                        // put in a vector of n (2D x T)
                                        // No transformations.
                                        readImages1(o_output,im, i_attr);
                                }
                                else
                                {
                                        // put in a vector of T (2D x n)
                                        readImages4(o_output,im, i_attr);
                                }
                        }
                }
        }



void GimmickView::ReadImagesNotThreadedInVector(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
{
        // Create the output data
        if (im.size()==1)
        {
                // Only one image : give it
                vtkImageData* out = vtkImageData::New();
                GimmickDebugMessage(3, "State Check: Full Filename: "
                                                <<im.front()
                                                <<std::endl);
                out->ShallowCopy(mReader.GetImage(im.front()));
                s.push_back( out );
        }
        else if (im.size()>1) // Test inutile ? JPR
        {
                vtkImageData* first = mReader.GetImage( im.front());
                if (dimension == 2)
                {
                 // n3D
                  std::vector<std::string>::iterator it;
                        for (it=im.begin(); it!=im.end(); ++it)
                        {
                                vtkImageData* out = vtkImageData::New();
                                out->ShallowCopy(mReader.GetImage(*it));
                                s.push_back(out);
                        }
                }
                else
                {
                        // n2D to 3D // NO!
                        // n *2D + T in a vector :
                        
                        std::vector<std::string>::iterator it;
                        for (it=im.begin(); it!=im.end(); ++it) 
                        {
                                vtkImageData* out = mReader.GetImage( (*it));
                                s.push_back(out);
                        }
                }
        }
}
  //======================================================================

  //======================================================================
  void GimmickView::RequestReading(tree::Node* n, 
          int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
  {
          if(!mReaderStarted)
          {
                mReader.Start();
                mReaderStarted=true;
          }
    ImageEventType t(n,selection_index);
        t.pointerHolder = p;
    mImageEventMap[n->GetAttribute("FullFileName")] = t;    
    mReader.Request(this,n->GetAttribute("FullFileName"),prio);
  }
  //======================================================================

  //======================================================================
  void GimmickView::
  OnMultiThreadImageReaderEvent(const std::string& filename,
                                MultiThreadImageReaderUser::EventType e,
                                vtkImageData* image)
  {
    GimmickDebugMessage(7,
                        "MultiThreadImageReader event : "<<e<<std::endl);
        if (e==ImageLoaded)
        {
                if (filename.size()==0)
                {
                  //What to do in this case?
                  /*
                        GimmickDebugMessage(5,
                                        "Pushing unknown image in queue"
                                        <<std::endl);
                        mImageEventQueue.push_back(ImageEventType(image));*/
                        return;
                }
                ImageEventTypeMap::iterator i;
//JCP 22-06-2009, test mImageEventMap.size() > 0
                if(mImageEventMap.size()>0){
                        i = mImageEventMap.find(filename);
                        if (i!=mImageEventMap.end())
                        {
                                GimmickDebugMessage(5,
                                                "Putting image of file '"<<filename<<"' on pointer"
                                                <<std::endl);
                                ImageEventType ie(i->second);
                                ie.image = image;
                                ie.pointerHolder->Set(ie.image);
                                //mImageEventMap.erase(i);
                        }
                }    
        }
        else if (e==Error)
        {
                std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
                mess+=filename;
                mess+="\n";
                GimmickMessage(1,mess);
                ImageEventTypeMap::iterator i;
                i = mImageEventMap.find(filename);
                if (i!=mImageEventMap.end())
                {
                ImageEventType ie(i->second);
                ie.image = image;
                ie.pointerHolder->Set(GetDefaultImage());
                //mImageEventMap.erase(i);
                }
        }

        else if (e==ImageUnloaded)
        {
                std::string mess="Unloaded image in file ";
                mess+=filename;
                mess+="\n";
                GimmickMessage(1,mess);
                                ImageEventTypeMap::iterator i;
                i = mImageEventMap.find(filename);
                if (i!=mImageEventMap.end())
                {
                ImageEventType ie(i->second);
                ie.image = image;
                ie.pointerHolder->Set(GetDefaultImage());
                //mImageEventMap.erase(i);
                }
        }
  }

  //====================================================================

  //====================================================================
  void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
  {
    mValidationSignal.connect(callback);
  }       
  


//void GimmickView::Anonymize(std::vector<std::string> i_filenames, int type)
//{
//
//      gdcm::FileMetaInformation::SetSourceApplicationEntityTitle( "gdcmanon" );
//    gdcm::Global& g = gdcm::Global::GetInstance();
//      //if( !resourcespath )
// //   {
//      //      const char *xmlpathenv = getenv("GDCM_RESOURCES_PATH");
//      //      if( xmlpathenv )
//      //        {
//      //        // Make sure to look for XML dict in user explicitly specified dir first:
//      //        xmlpath = xmlpathenv;
//      //        resourcespath = 1;
//      //        }
// //   }
// // if( resourcespath )
// //   {
// //   // xmlpath is set either by the cmd line option or the env var
// //   if( !g.Prepend( xmlpath.c_str() ) )
// //     {
// //     std::cerr << "Specified Resources Path is not valid: " << xmlpath << std::endl;
// //     return 1;
// //     }
// //   }
//  // All set, then load the XML files:
//  if( !g.LoadResourcesFiles() )
//  {
//    return ;
//  }
//  const gdcm::Defs &defs = g.GetDefs(); (void)defs;
//  if( !rootuid )
//    {
//    // only read the env var if no explicit cmd line option
//    // maybe there is an env var defined... let's check
//    const char *rootuid_env = getenv("GDCM_ROOT_UID");
//    if( rootuid_env )
//      {
//      rootuid = 1;
//      root = rootuid_env;
//      }
//    }
//  if( rootuid )
//    {
//    // root is set either by the cmd line option or the env var
//    if( !gdcm::UIDGenerator::IsValid( root.c_str() ) )
//      {
//      std::cerr << "specified Root UID is not valid: " << root << std::endl;
//      return 1;
//      }
//    gdcm::UIDGenerator::SetRoot( root.c_str() );
//    }
//
//      if(type == 0)
//      {
//                 // Get private key/certificate
//           gdcm::CryptographicMessageSyntax cms;
//              if( !dumb_mode )
//              {
//                      if( !GetRSAKeys(cms, rsa_path.c_str(), cert_path.c_str() ) )
//                      {
//                              return 1;
//                      }
//                      cms.SetCipherType( ciphertype );
//              }
//
//      // Setup gdcm::Anonymizer
//      gdcm::Anonymizer anon;
//              if( !dumb_mode )
//              anon.SetCryptographicMessageSyntax( &cms );
//
//              if( dumb_mode )
//              {
//                      for(unsigned int i = 0; i < nfiles; ++i)
//              {
//              const char *in  = filenames[i].c_str();
//              const char *out = outfilenames[i].c_str();
//              if( !AnonymizeOneFileDumb(anon, in, out, empty_tags, remove_tags, replace_tags_value) )
//        {
//                      //std::cerr << "Could not anonymize: " << in << std::endl;
//                      return 1;
//              }
//      }
//    }
//  //else
//  //  {
//  //  for(unsigned int i = 0; i < nfiles; ++i)
//  //    {
//  //    const char *in  = filenames[i].c_str();
//  //    const char *out = outfilenames[i].c_str();
//  //    if( !AnonymizeOneFile(anon, in, out) )
//  //      {
//  //      //std::cerr << "Could not anonymize: " << in << std::endl;
//  //      return 1;
//  //      }
//  //    }
//  //  }
//      }
//}
//
} // EO namespace creaImageIO