2 # ---------------------------------------------------------------------
4 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
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
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.
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
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 # ------------------------------------------------------------------------
29 #include <creaImageIOGimmickView.h>
31 #include "boost/filesystem.hpp"
34 #include <gdcmGlobal.h>
35 #include <gdcmSerieHelper.h>
36 #include <gdcmFileHelper.h>
38 #include <vtkGdcmReader.h>
41 /*#if defined(USE_GDCM2)
42 #include <vtkGDCMImageReader.h>
43 #include "gdcmSystem.h"
44 #include "gdcmCryptographicMessageSyntax.h"
45 #include "gdcmUIDGenerator.h"
46 #include "gdcmAnonymizer.h"
47 #include "gdcmGlobal.h"
50 #pragma warning(disable: 4996)
53 namespace fs = boost::filesystem;
57 ///Class used to represent the actual state of the image selected and to perform comparisons on its values
61 ImageExtent(const std::string& x, const std::string& y, const std::string& z, const std::string& t)
63 sscanf(x.c_str(),"%d",&mExtent[0]);
64 sscanf(y.c_str(),"%d",&mExtent[1]);
65 sscanf(z.c_str(),"%d",&mExtent[2]);
66 sscanf(t.c_str(),"%d",&mExtent[3]);
67 if(x==""){mExtent[0]=1;}
68 if(y==""){mExtent[1]=1;}
69 if(z==""){mExtent[2]=1;}
70 if(t==""){mExtent[3]=1;}
72 if (mExtent[3]>1) mDim=4;
73 else if (mExtent[2]>1) mDim=3;
74 else if (mExtent[1]>1) mDim=2;
75 else if (mExtent[0]>1) mDim=1;
81 void Clear() { mExtent[0] = mExtent[1] = mExtent[2] = mExtent[3] = 1; }
83 ///Returns true if the two extents are compatible
84 bool IsCompatible( const ImageExtent& );
86 ///Adds the extent passed as a parameter to the current extent
87 void Add ( const ImageExtent& );
89 ///Returns the ieth position of the extent
90 int Get(int i) { return mExtent[i]; }
92 ///Returns the dimension of the current image
93 void SetDimension(int dim) { mDim=dim; }
95 ///Returns the dimension of the current image
96 int GetDimension() { return mDim; }
103 //======================================================================
105 //======================================================================
107 GimmickView::GimmickView(boost::shared_ptr<Gimmick> gimmick, int threads)
111 GimmickDebugMessage(1,"GimmickView::GimmickView"
113 // Anciently started the threads ...
114 // Threads now automatically start at first image request
118 //======================================================================
120 //======================================================================
122 GimmickView::~GimmickView()
124 GimmickDebugMessage(1,"GimmickView::~GimmickView"
127 //======================================================================
129 //======================================================================
130 /// Initializes the view :
131 /// Creates the TreeViews for all the TreeHandler of the Controller
133 void GimmickView::Initialize()
135 mReaderStarted=false;
137 //======================================================================
139 //======================================================================
141 void GimmickView::Finalize()
145 //======================================================================
147 //======================================================================
148 /// Create the tree views
149 void GimmickView::CreateTreeViews()
151 GimmickMessage(2,"Creating the tree views"<<std::endl);
152 Gimmick::TreeHandlerMapType::const_iterator i;
153 for (i = mGimmick->GetTreeHandlerMap().begin();
154 i!= mGimmick->GetTreeHandlerMap().end();
157 this->CreateTreeView(i->second);
161 /// Create a tree view with a given name
162 void GimmickView::CreateSingleTreeView(std::string &i_name)
164 this->CreateTreeView(mGimmick->GetTreeHandlerMap()[i_name]);
169 //======================================================================
171 //======================================================================
172 /// Updates the TreeView of given name from level l to bottom
173 /// (calls the virtual method TreeView::Update())
174 void GimmickView::UpdateTreeViewLevel(const std::string& t, int l)
176 TreeViewMapType::iterator i;
177 i = GetTreeViewMap().find(t);
178 if ( i == GetTreeViewMap().end() )
180 GimmickError("INTERNAL ERROR : GimmickView::UpdateTreeView : '"
181 <<t<<"' is not in TreeViewMap");
183 i->second->UpdateLevel(l);
186 //======================================================================
187 /// Clears the status and begins a new selection process
188 void GimmickView::ResetExtent()
192 mImageExtent.reset();
197 //======================================================================
199 //======================================================================
200 bool ImageExtent::IsCompatible(const ImageExtent& ie)
202 bool compatible=true;
203 ImageExtent * extent= (ImageExtent*)&ie;
204 if((*extent).Get(0)!=Get(0)
205 || (*extent).Get(1)!=Get(1))
212 //======================================================================
214 //======================================================================
215 void ImageExtent::Add(const ImageExtent& ie)
217 ImageExtent * extent= (ImageExtent*)&ie;
218 mExtent[2]+=(*extent).Get(2);
225 //======================================================================
226 /// No selected image
227 bool GimmickView::NoValidateSelected ()
229 GimmickDebugMessage(2,"Validating selected"<<std::endl);
230 std::string mMessage;
231 mMessage="Cannot have 0 images selected!";
233 modifyValidationSignal(valid);
234 SetMessage(mMessage);
238 //======================================================================
239 ///Validates the dimension compliance of the images with the maximum and
240 ///minimum given, and between their sizes
241 bool GimmickView::ValidateSelected (tree::Node* sel, int min_dim, int max_dim)
243 GimmickDebugMessage(2,"Validating selected"<<std::endl);
244 std::string mMessage;
248 mMessage="Cannot have 0 images selected!";
253 boost::shared_ptr<ImageExtent> ie=boost::shared_ptr<ImageExtent>(new ImageExtent((*sel).GetAttribute("D0028_0010"),
254 (*sel).GetAttribute("D0028_0011"),
255 (*sel).GetAttribute("D0028_0012"),
261 if((mImageExtent->Get(min_dim-1)<2)||(mImageExtent->Get(max_dim)>1))
267 std::stringstream out;
268 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
269 mMessage = out.str();
270 mImageExtent->SetDimension(2);
276 if(mImageExtent->IsCompatible(*ie))
278 if(mImageExtent->GetDimension()==max_dim && mImageExtent->Get(max_dim)>2)
280 std::stringstream out;
281 out<<"Cannot add this image to selection : would result in a "<<mImageExtent->GetDimension()+1<<"D image!";
287 std::stringstream out;
288 out<<"Selecting "<<mImageExtent->GetDimension()<<"D images is not allowed !";
292 else if(min_dim==3 && (ie->Get(2)+mImageExtent->Get(2))<2)
294 std::stringstream out;
295 out << "Cannot build the selection as it would result in a ";
296 out << mImageExtent->GetDimension();
297 out << "D image, and the minimum is ";
305 mImageExtent->Add(*ie);
306 std::stringstream out;
307 out << mImageExtent->GetDimension() << "D image " << mImageExtent->Get(0) << "x"<< mImageExtent->Get(1) << "x"<< mImageExtent->Get(2) <<" selected";
308 mMessage = out.str();
313 mMessage="The selected images are not compatible.";
319 modifyValidationSignal(valid);
320 SetMessage(mMessage);
324 //======================================================================
325 void GimmickView::modifyValidationSignal(bool ivalid)
328 mValidationSignal(ivalid);
331 void GimmickView::stopReader()
336 //======================================================================
337 ///Reads Images (Non Threaded)
338 void GimmickView::ReadImagesNotThreaded(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
343 #define GIMMICK_NO_IMAGE_SELECTION 0
344 #define GIMMICK_2D_IMAGE_SELECTION 2
345 #define GIMMICK_3D_IMAGE_SELECTION 3
346 #define GIMMICK_4D_IMAGE_SELECTION 4
353 // Create the output data
356 vtkImageData * out=vtkImageData::New();
357 out->ShallowCopy(mReader.GetImage(im.front()));
360 else if (im.size()>1) // Test inutile ? JPR
362 vtkImageData* first = mReader.GetImage( im.front());
366 std::vector<std::string>::iterator it;
367 for (it=im.begin(); it!=im.end(); ++it)
369 vtkImageData* out = vtkImageData::New();
370 out->ShallowCopy(mReader.GetImage(*it));
377 vtkImageData* out = vtkImageData::New();
378 // out->CopyStructure(first);
380 //first->GetExtent(ext); // JPR
381 //EED 2017-01-01 Migration VTK7
382 #if VTK_MAJOR_VERSION <= 5
383 first->GetWholeExtent(ext); // renvoie egalement 0,0 en Z // JPR
385 first->GetExtent(ext); // renvoie egalement 0,0 en Z // JPR
390 ext[5] = (int)im.size()-1;
394 ext[5] = ext[5] * (int)im.size()-1; // to deal with multiframes - JPR
398 // LG : TODO : Z Spacing ?
401 first->GetDimensions(dim);
402 out->SetDimensions(dim[0], dim[1], (int)im.size() );
404 //EED 2017-01-01 Migration VTK7
405 #if VTK_MAJOR_VERSION <= 5
406 out->SetScalarType(first->GetScalarType());
407 out->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
408 out->AllocateScalars();
411 out->AllocateScalars(first->GetScalarType(),first->GetNumberOfScalarComponents());
414 unsigned long imsize = dim[0] * dim[1];
415 imsize = imsize * dim[2] ; // deal with multiframes // JPR
419 // differents formats char , short, etc...
420 // differents components 1..3 ex. jpg ->RGB 3
421 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
424 // Order the file name vector
427 first->GetSpacing(spc);
429 // OrderTheFileNameVector is not here anymore.
430 // Try orderFilesWithZSpacing from OutputModel FCY
431 // spc[2]=OrderTheFileNameVector(im);
434 out->SetSpacing(spc);
437 std::vector<std::string>::iterator it;
439 for (it=im.begin(); it!=im.end(); ++it)
441 vtkImageData* cur = mReader.GetImage( (*it) );
442 memcpy(out->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
452 //======================================================================
456 //////////////////////////////////////////////////////////
457 // Test if the image is a multiple or single frame.
458 // For the moment only with the creation of vtkImageDta.
459 // TO DO with image size and dim!!!
460 //////////////////////////////////////////////////////////
461 bool GimmickView::isSingle(const std::string i_file)
464 vtkImageData* first = mReader.GetImage( i_file);
466 first->GetDimensions(dim);
477 //////////////////////////////////////////////////////////
478 // get Attributes values for a file
479 //////////////////////////////////////////////////////////
481 void GimmickView::getAttributes(const std::string i_file, std::map<std::string, std::string> &o_infos, OutputAttr i_attr)
483 if(i_attr.inside.size() >0)
485 mGimmick->GetAttributes(i_file,o_infos,i_attr);
487 if(i_attr.outside.size()>0)
489 mReader.getAttributes(i_file,o_infos, i_attr.outside);
493 //////////////////////////////////////////////////////////
494 // create an output structure with n entries = n output
495 //////////////////////////////////////////////////////////
496 void GimmickView::readImages1(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
500 i_attr.outside.push_back("D0019_100a"); // simens Number Of Images In Mosaic
503 int i,isize=i_attr.outside.size();
504 for (i=0;i<isize;i++)
506 printf("EED GimmickView::readImages1 A i_attr out %s \n", i_attr.outside[i].c_str());
511 std::vector<std::string>::iterator it;
512 for (it=im.begin(); it!=im.end(); ++it)
515 out.img = vtkImageData::New();
516 out.img->ShallowCopy(mReader.GetImage(*it));
517 printf("EED GimmickView::readImages1 A\n ");
520 getAttributes((*it),out.infos,i_attr);
521 printf("EED GimmickView::readImages1 B %s \n ", out.infos.find("D0019_100a")->second.c_str() );
523 o_output.push_back(out);
525 // If we want only one output information structure, we set it outside the loop
528 getAttributes(im.front(), o_output.front().infos, i_attr);
529 printf("EED GimmickView::readImages1 C %s \n ", o_output.front().infos.find("D0019_100a")->second.c_str() );
534 //////////////////////////////////////////////////////////
535 // create an output structure with n entries = 1 output
536 //////////////////////////////////////////////////////////
537 void GimmickView::readImages3(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
538 OutputAttr i_attr, double i_zspc)
541 vtkImageData* first = mReader.GetImage( im.front());
542 out.img = vtkImageData::New();
544 //EED 2017-01-01 Migration VTK7
545 #if VTK_MAJOR_VERSION <= 5
546 first->GetWholeExtent(ext); // send also 0,0 in Z
548 first->GetExtent(ext); // send also 0,0 in Z
552 ext[5] = (int)im.size()-1;
554 ext[5] = ext[5] * (int)im.size()-1; // to deal with multiframes
556 out.img->SetExtent(ext);
559 first->GetDimensions(dim);
560 first->GetSpacing(spac);
561 out.img->SetSpacing(spac);
562 out.img->SetDimensions(dim[0], dim[1], (int)im.size() );
566 //EED 2017-01-01 Migration VTK7
567 #if VTK_MAJOR_VERSION <= 5
568 out.img->SetScalarType(first->GetScalarType());
569 out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
570 out.img->AllocateScalars();
573 out.img->AllocateScalars(first->GetScalarType(), first->GetNumberOfScalarComponents());
577 unsigned long imsize = dim[0] * dim[1];
578 imsize = imsize * dim[2] ; // deal with multiframes here
579 // differents formats char , short, etc...
580 // differents components 1..3 ex. jpg ->RGB 3
581 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
582 // Order the file name vector already did with the OutputModel!!!
583 //!!!!out.img->SetSpacing(i_zspc);
585 std::vector<std::string>::iterator it;
586 for (it=im.begin(); it!=im.end(); ++it)
588 vtkImageData* cur = mReader.GetImage( (*it) );
589 memcpy(out.img->GetScalarPointer(0,0,slice), cur->GetScalarPointer(0,0,0), imsize);
592 getAttributes(im.front(),out.infos, i_attr);
593 o_output.push_back(out);
597 // TO DO NO VERY SURE : NEED TO BE TESTED
598 //////////////////////////////////////////////////////////
599 // create an output structure with n entries (T size) = T output (n size)
600 //////////////////////////////////////////////////////////
601 void GimmickView::readImages2(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
602 OutputAttr i_attr, double i_zspc)
604 vtkImageData* first = mReader.GetImage( im.front());
607 first->GetDimensions(dim);
608 first->GetSpacing(spac);
609 // differents formats char , short, etc...
610 // differents components 1..3 ex. jpg ->RGB 3
611 unsigned long imsize = dim[0] * dim[1];
612 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
614 // Order the file name vector already did with the OutputModel!!!
615 std::vector<std::string>::iterator it;
616 std::vector<OutStrGimmick>::iterator it_out = o_output.begin();
618 for (it=im.begin(); it!=im.end(); ++it)//, it_out ++)
620 vtkImageData* cur = mReader.GetImage( (*it) );
621 for (int slice= 0 ; slice <dim[2]; slice++)
624 out.img = vtkImageData::New();
625 out.img->SetSpacing(spac);
627 //EED 2017-01-01 Migration VTK7
628 #if VTK_MAJOR_VERSION <= 5
629 first->GetWholeExtent(ext); // send also 0,0 in Z
631 first->GetExtent(ext); // send also 0,0 in Z
634 out.img->SetExtent(ext);
635 out.img->SetDimensions(dim[0], dim[1], 1 );
637 //EED 2017-01-01 Migration VTK7
638 #if VTK_MAJOR_VERSION <= 5
639 out.img->SetScalarType(first->GetScalarType());
640 out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
641 out.img->AllocateScalars();
644 out.img->AllocateScalars(first->GetScalarType(),first->GetNumberOfScalarComponents());
648 memcpy(out.img->GetScalarPointer(0,0,0), cur->GetScalarPointer(0,0,slice), imsize);
649 o_output.push_back(out);
652 // getAttributes((*it),(*it_out).infos,i_attr);
656 getAttributes(im.front(), o_output.front().infos,i_attr);
661 //////////////////////////////////////////////////////////
662 // create an output structure with n entries (T size) = T + n output
663 //////////////////////////////////////////////////////////
664 void GimmickView::readImages4(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
667 std::vector<std::string>::iterator it;
668 std::vector<OutStrGimmick>::iterator it_out = o_output.begin();
669 vtkImageData* first = mReader.GetImage( im.front());
671 first->GetDimensions(dim);
673 for (int slice= 0 ; slice <dim[2]; slice++)
676 out.img = vtkImageData::New();
680 //EED 2017-01-01 Migration VTK7
681 #if VTK_MAJOR_VERSION <= 5
682 first->GetWholeExtent(ext); // send also 0,0 in Z
684 first->GetExtent(ext); // send also 0,0 in Z
689 out.img->SetExtent(ext);
690 first->GetSpacing(spac);
691 out.img->SetSpacing(spac);
692 out.img->SetDimensions(dim[0], dim[1], (int)im.size() );
694 //EED 2017-01-01 Migration VTK7
695 #if VTK_MAJOR_VERSION <= 5
696 out.img->SetScalarType(first->GetScalarType());
697 out.img->SetNumberOfScalarComponents(first->GetNumberOfScalarComponents());
698 out.img->AllocateScalars();
701 out.img->AllocateScalars(first->GetScalarType(), first->GetNumberOfScalarComponents());
704 unsigned long imsize = dim[0] * dim[1];
705 imsize = imsize * first->GetScalarSize() * first->GetNumberOfScalarComponents();
708 for (it=im.begin(); it!=im.end(); ++it, index ++)
710 vtkImageData* cur = mReader.GetImage( (*it) );
711 memcpy(out.img->GetScalarPointer(0,0,index), cur->GetScalarPointer(0,0,slice), imsize);
712 o_output.push_back(out);
715 if(!i_attr.mult) // No sense to take informations in all images
717 getAttributes(im.front(), o_output.front().infos,i_attr);
723 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
724 // Global function to read Images and create wished output (informations, multiple, type and size of output...)
725 // In function of type (file, vector) and size, the correct readImages function is selected
726 //////////////////////////////////////////////////////////////////////////////////////////////////////////////
727 void GimmickView::readImages(std::vector<OutStrGimmick>& o_output, std::vector<std::string> im,
728 OutputAttr i_attr, int i_dim, double i_zspc)
730 int size = (int)im.size();
734 } else if (size == 1) {
736 // Only one image : give it
737 // But take in count multiframe possibility
738 if ( isSingle(im.front()) || i_dim != 1)
740 readImages1(o_output,im, i_attr);
742 readImages2(o_output,im, i_attr,i_zspc);
745 // multiple or single frame
746 if ( isSingle(im.front()) )
748 //we deal with 2D images
752 readImages3(o_output,im, i_attr,i_zspc);
754 readImages1(o_output,im, i_attr);
757 // we deal with multiple frames n x (2D x T)
758 // Differents outputs are avaialable
761 // put all in one output
762 readImages3(o_output,im, i_attr,i_zspc);
763 } else if( i_dim == 2) {
764 // put in a vector of n x T (2D)
765 readImages2(o_output,im, i_attr,i_zspc);
766 } else if( i_dim == 3) {
767 // put in a vector of n (2D x T)
768 // No transformations.
769 readImages1(o_output,im, i_attr);
771 // put in a vector of T (2D x n)
772 readImages4(o_output,im, i_attr);
774 } //isSingle(im.front())
777 //EED UnMosaic step...
778 //How to verifie if is a mosaic file , with how many images inside??
784 void GimmickView::ReadImagesNotThreadedInVector(std::vector<vtkImageData*>& s, std::vector<std::string> im, int dimension)
786 // Create the output data
789 // Only one image : give it
790 vtkImageData* out = vtkImageData::New();
791 GimmickDebugMessage(3, "State Check: Full Filename: "
794 out->ShallowCopy(mReader.GetImage(im.front()));
797 else if (im.size()>1) // Test inutile ? JPR
799 /// \TODO fix unused variable 'first'
800 // vtkImageData* first = mReader.GetImage( im.front());
804 std::vector<std::string>::iterator it;
805 for (it=im.begin(); it!=im.end(); ++it)
807 vtkImageData* out = vtkImageData::New();
808 out->ShallowCopy(mReader.GetImage(*it));
813 // n *2D + T in a vector :
815 std::vector<std::string>::iterator it;
816 for (it=im.begin(); it!=im.end(); ++it)
818 vtkImageData* out = mReader.GetImage( (*it));
824 //======================================================================
826 //======================================================================
827 ///Requests the reading of an image
828 void GimmickView::RequestReading(tree::Node* n,
829 int prio, int selection_index, boost::shared_ptr<ImagePointerHolder> p)
836 ImageEventType t(n,selection_index);
838 mImageEventMap[n->GetAttribute("FullFileName")] = t;
839 mReader.Request(this,n->GetAttribute("FullFileName"),prio);
841 //======================================================================
843 //======================================================================
845 OnMultiThreadImageReaderEvent(const std::string& filename,
846 MultiThreadImageReaderUser::EventType e,
849 GimmickDebugMessage(7,
850 "MultiThreadImageReader event : "<<e<<std::endl);
853 if (filename.size()==0)
855 //What to do in this case?
857 GimmickDebugMessage(5,
858 "Pushing unknown image in queue"
860 mImageEventQueue.push_back(ImageEventType(image));*/
863 ImageEventTypeMap::iterator i;
864 //JCP 22-06-2009, test mImageEventMap.size() > 0
865 if(mImageEventMap.size()>0){
866 i = mImageEventMap.find(filename);
867 if (i!=mImageEventMap.end())
869 GimmickDebugMessage(5,
870 "Putting image of file '"<<filename<<"' on pointer"
872 ImageEventType ie(i->second);
874 ie.pointerHolder->Set(ie.image);
875 //mImageEventMap.erase(i);
881 std::string mess="ERROR: MultiThreadImageReader: Cannot read image in file ";
884 GimmickMessage(1,mess);
885 ImageEventTypeMap::iterator i;
886 i = mImageEventMap.find(filename);
887 if (i!=mImageEventMap.end())
889 ImageEventType ie(i->second);
891 ie.pointerHolder->Set(GetDefaultImage());
892 //mImageEventMap.erase(i);
896 else if (e==ImageUnloaded)
898 std::string mess="Unloaded image in file ";
901 GimmickMessage(1,mess);
902 ImageEventTypeMap::iterator i;
903 i = mImageEventMap.find(filename);
904 if (i!=mImageEventMap.end())
906 ImageEventType ie(i->second);
908 ie.pointerHolder->Set(GetDefaultImage());
909 //mImageEventMap.erase(i);
914 //====================================================================
916 //====================================================================
917 void GimmickView::ConnectValidationObserver(ValidationCallbackType callback)
919 mValidationSignal.connect(callback);
923 #if defined(USE_GDCM)
924 ////////////////////////////////////////////////////////////////////////
925 void GimmickView::Anonymize(std::vector<std::string> i_filenames, int type)
928 std::vector<GDCM_NAME_SPACE::FileHelper *> filesH;
929 std::vector<std::string> suid;
930 std::map<std::string, std::string> msuid;
931 std::string tempuid = GDCM_NAME_SPACE::Util::CreateUniqueUID();
933 std::vector<std::string>::iterator it = i_filenames.begin();
934 for(; it != i_filenames.end(); it++)
937 GDCM_NAME_SPACE::File *file;
938 file = GDCM_NAME_SPACE::File::New( );
939 file->SetLoadMode( GDCM_NAME_SPACE::LD_ALL );
940 file->SetFileName( (*it).c_str() );
944 std::cerr << "Sorry, " << (*it).c_str() <<" not a gdcm-readable "
945 << "DICOM / ACR File" <<std::endl;
949 std::cout << " ... is readable " << std::endl;
951 // We need a gdcm::FileHelper, since we want to load the pixels
952 GDCM_NAME_SPACE::FileHelper *fh = GDCM_NAME_SPACE::FileHelper::New(file);
955 //EED uint8_t *imageData = fh->GetImageData();
958 file->AddAnonymizeElement(0x0008, 0x0080, "*");
960 file->AddAnonymizeElement(0x0010, 0x0010, "*");
962 file->AddAnonymizeElement( 0x0010, 0x0020,"1515" );
963 // Study Instance UID
964 file->AddAnonymizeElement(0x0020, 0x000d, tempuid );
966 file->AddAnonymizeElement(0x0010, 0x2154, "3615" );
968 // Aware user will add here more fields to anonymize here
970 // The gdcm::File is modified in memory
972 file->AnonymizeFile();
976 fh->SetContentType(GDCM_NAME_SPACE::UNMODIFIED_PIXELS_IMAGE);
978 fh->WriteDcmExplVR(file->GetFileName() +".ano1" );
979 std::cout << i <<" End Anonymize" << std::endl;
980 file->ClearAnonymizeList();
986 #if defined(USE_GDCM2)
987 void GimmickView::Anonymize(std::vector<std::string> i_filenames, int type)
991 ////////////////////////////////////////////////////////////////////////
992 //void GimmickView::Anonymize(std::vector<std::string> i_filenames, int type)
995 // gdcm::FileMetaInformation::SetSourceApplicationEntityTitle( "gdcmanon" );
996 // gdcm::Global& g = gdcm::Global::GetInstance();
997 // //if( !resourcespath )
999 // // const char *xmlpathenv = getenv("GDCM_RESOURCES_PATH");
1000 // // if( xmlpathenv )
1002 // // // Make sure to look for XML dict in user explicitly specified dir first:
1003 // // xmlpath = xmlpathenv;
1004 // // resourcespath = 1;
1007 // // if( resourcespath )
1009 // // // xmlpath is set either by the cmd line option or the env var
1010 // // if( !g.Prepend( xmlpath.c_str() ) )
1012 // // std::cerr << "Specified Resources Path is not valid: " << xmlpath << std::endl;
1016 // // All set, then load the XML files:
1017 // if( !g.LoadResourcesFiles() )
1021 // const gdcm::Defs &defs = g.GetDefs(); (void)defs;
1024 // // only read the env var if no explicit cmd line option
1025 // // maybe there is an env var defined... let's check
1026 // const char *rootuid_env = getenv("GDCM_ROOT_UID");
1027 // if( rootuid_env )
1030 // root = rootuid_env;
1035 // // root is set either by the cmd line option or the env var
1036 // if( !gdcm::UIDGenerator::IsValid( root.c_str() ) )
1038 // std::cerr << "specified Root UID is not valid: " << root << std::endl;
1041 // gdcm::UIDGenerator::SetRoot( root.c_str() );
1046 // // Get private key/certificate
1047 // gdcm::CryptographicMessageSyntax cms;
1050 // if( !GetRSAKeys(cms, rsa_path.c_str(), cert_path.c_str() ) )
1054 // cms.SetCipherType( ciphertype );
1057 // // Setup gdcm::Anonymizer
1058 // gdcm::Anonymizer anon;
1060 // anon.SetCryptographicMessageSyntax( &cms );
1064 // for(unsigned int i = 0; i < nfiles; ++i)
1066 // const char *in = filenames[i].c_str();
1067 // const char *out = outfilenames[i].c_str();
1068 // if( !AnonymizeOneFileDumb(anon, in, out, empty_tags, remove_tags, replace_tags_value) )
1070 // //std::cerr << "Could not anonymize: " << in << std::endl;
1077 // // for(unsigned int i = 0; i < nfiles; ++i)
1079 // // const char *in = filenames[i].c_str();
1080 // // const char *out = outfilenames[i].c_str();
1081 // // if( !AnonymizeOneFile(anon, in, out) )
1083 // // //std::cerr << "Could not anonymize: " << in << std::endl;
1091 } // EO namespace creaImageIO