1 /*=========================================================================
4 Module: $RCSfile: test4DSplitter.cxx,v $
6 Date: $Date: 2011/03/30 14:50:44 $
7 Version: $Revision: 1.1 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
19 /* Raisons ne pas utiliser itkImageSeriesReader:
21 On Wed, Feb 16, 2011 at 11:51 AM, Roger Bramon Feixas <rogerbramon@gmail.com>
23 I'm developing with ITK 3.20 + GDCM 2.0.17 + VTK 5.6 and I've noticed
24 itkImageSeriesReader is ~2x slower than vtkGDCMImageReader (from GDCM2).
25 I compared both codes and I think the difference is the extra copy which
26 itkImageSeriesReader makes from ImageFileReader's output to its own output
27 (ImageSeriesReader::GenerateData() line 393). I've commented it just to test
28 the time needed without making the copy and now both readers take more or less
32 //#include "gdcmDocEntry.h"
33 //#include "gdcmDicomDir.h"
34 //#include "gdcmDicomDirPatient.h"
36 #include "gdcmFileHelper.h"
37 //#include "gdcmDirList.h"
38 //#include "gdcmDebug.h"
39 //#include "gdcmArgMgr.h"
41 #include "gdcmSerieHelper.h"
42 #include <vtkImageData.h>
44 #include "vtkGdcmReader.h"
47 #include "vtkGdcm4DSplitter.h"
53 typedef std::map<std::string, GDCM_NAME_SPACE::File*> SortedFiles;
55 // Due to a deep C++ trouble in 'SortOnTag', we use 'SortOnUserFunction'
56 bool myCompareFunction(GDCM_NAME_SPACE::File *file1, GDCM_NAME_SPACE::File *file2)
58 return atof(file1->GetEntryString(0x0018,0x1060).c_str()) < atof(file2->GetEntryString(0x0018,0x1060).c_str());
61 int main(int argc, char *argv[])
65 std::cout << "... inside " << argv[0] << std::endl;
69 std::string strDirName("/home/jpr/Desktop/Patients_Emilie/Patient.3T/AUB Jos/AUBERTIN JOSEPH/PROSTATE - 305629373/DYN7INJDYN6_901");
72 // ----- Begin Processing -----
74 unsigned short int grSplit;
75 unsigned short int elSplit;
77 unsigned short int grSort;
78 unsigned short int elSort;
80 // Pour un directory '4D'
81 // en sortie, chaque vtkImageData contiendra une coupe au cours du temps
84 std::vector<vtkImageData*> *output;
87 vtkGdcm4DSplitter *spl = new vtkGdcm4DSplitter();
88 spl->setDirName(strDirName);
89 spl->setRecursive(false);
90 spl->setSplitOnPosition();
92 // Time triger : 0018|1060
96 // ==> Big troubles with SortOnTag
97 //spl->setSortOnTag(grSort, elSort);
98 //spl->setSortConvertToFloat(true);
100 // ==> use SortOnUserFunction !
101 spl->setSortOnUserFunction(myCompareFunction);
103 std::cout << "Everything set" << std::endl;
106 std::cout << "GO() done, status " << res << std::endl;
109 std::cout << "plantage!" << std::endl;
112 output = spl->GetImageDataVector();
115 std::cout << "--------------------------------" << std::endl;
116 std::cout << "Vector size " << output->size() << std::endl;
117 std::cout << "--------------------------------" << std::endl;
119 std::vector<vtkImageData*>::iterator it;
120 // Print the first one (why not?)
122 (*it)->PrintSelf(std::cout, vtkIndent(2));
125 // Pour un directory '4D'
126 // en sortie, chaque vtkImageData contiendra un volume au cours du temps.