]> Creatis software - gdcm.git/blob - Testing/TestDicomDir.cxx
Let's take into account (again) LoadMode !
[gdcm.git] / Testing / TestDicomDir.cxx
1 /*=========================================================================
2                                                                                 
3   Program:   gdcm
4   Module:    $RCSfile: TestDicomDir.cxx,v $
5   Language:  C++
6   Date:      $Date: 2005/10/25 14:52:30 $
7   Version:   $Revision: 1.42 $
8                                                                                 
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.
12                                                                                 
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.
16                                                                                 
17 =========================================================================*/
18 #include "gdcmDocEntry.h"
19 #include "gdcmDataEntry.h"
20 #include "gdcmDicomDir.h"
21 #include "gdcmDicomDirPatient.h"
22 #include "gdcmDicomDirStudy.h"
23 #include "gdcmDicomDirSerie.h"
24 #include "gdcmDicomDirImage.h"
25 #include "gdcmTS.h"
26
27 #include <iostream>
28 #include <fstream>
29
30 // check *all* the dicom elements (gdcm::DocEntry)
31 // of this gdcm::DicomDirObject
32 int CompareSQItem(gdcm::SQItem *pa1, gdcm::SQItem *pa2 )
33 {
34    gdcm::DocEntry *e1;
35    gdcm::DocEntry *e2;
36
37    e2 = pa2->GetFirstEntry();
38    while (!e2)
39    {
40       // locate the corresponding element in 'source' file 
41       e1 = pa1->GetDocEntry( e2->GetGroup(),e2->GetElement() );
42
43       // an element doesn't exist in origin file
44       if (!e1)
45       {
46        std::cout << "DicomDir element " << std::hex 
47                  << e2->GetGroup() << "," <<e2->GetElement() << std::endl;
48        return 1; 
49       }
50       // skip SeqEntries (I don't want to deal with 'recursion pbs' here)
51       if ( !dynamic_cast<gdcm::DataEntry *>(e1) ||
52            !dynamic_cast<gdcm::DataEntry *>(e2) )
53          continue;
54
55       // a value is read as GDCM_UNFOUND 
56       if ( ((gdcm::DataEntry *)e1)->GetString() == gdcm::GDCM_UNFOUND )
57       {
58          std::cout << "for gdcm source DicomDir : element (" << std::hex 
59                    << e1->GetGroup() << "," <<e1->GetElement() 
60                    << ") has values [" << gdcm::GDCM_UNFOUND << "]"
61                    << std::endl;
62          return 1;
63       }
64
65       // values differ in source file and destination file
66       if ( ((gdcm::DataEntry *)e1)->GetString() != 
67            ((gdcm::DataEntry *)e2)->GetString() )
68       {
69  
70          // serious trouble : values differ in source and destination file
71          std::cout << "for gdcm DicomDir element (" << std::hex 
72                    << e2->GetGroup() << "," <<e2->GetElement() 
73                    << ") values differ [" 
74                    << ((gdcm::DataEntry *)e1)->GetString() << "] != [" 
75                    << ((gdcm::DataEntry *)e2)->GetString() << "]"
76                    << std::endl;
77           return 1;
78       }
79    }
80    return 0;
81
82  
83 int TestDicomDir(int argc, char *argv[])
84 {  
85    gdcm::DicomDir *dicomdir;
86    
87    gdcm::DicomDirPatient *pa1;
88    gdcm::DicomDirStudy   *st1;
89    gdcm::DicomDirSerie   *se1;
90    gdcm::DicomDirImage   *im1;
91
92    gdcm::TSKey v;
93     
94    std::string file; 
95    if (argc > 1) 
96       file = argv[1];    
97    else 
98    {
99       file += GDCM_DATA_ROOT;
100       file += "/DICOMDIR";
101    }
102
103    dicomdir = gdcm::DicomDir::New();
104    dicomdir->SetFileName(file);
105    dicomdir->Load();
106    if (argc > 2) 
107    {
108       int level = atoi(argv[2]);   
109       dicomdir->SetPrintLevel(level);
110    }
111
112    // Test if the DicomDir is readable
113    if( !dicomdir->IsReadable() )
114    {
115       std::cout<<"          DicomDir '"<<file
116                <<"' is not readable"<<std::endl
117                <<"          ...Failed"<<std::endl;
118
119       dicomdir->Delete();
120       return 1;
121    }
122    else
123    {
124       std::cout<<"          DicomDir '"<<file
125                <<"' is readable"<<std::endl;
126    }
127
128    // Test if the DicomDir contains any patient
129    if( !dicomdir->GetFirstPatient() )
130    {
131       std::cout<<"          DicomDir '"<<file
132                <<" has no patient"<<std::endl
133                <<"          ...Failed"<<std::endl;
134
135       dicomdir->Delete();
136       return 1;
137    }
138
139    // step by step structure full exploitation
140    std::cout << std::endl << std::endl  
141              << " = PATIENT/STUDY/SERIE/IMAGE List ============================" 
142              << std::endl<< std::endl;
143
144    pa1 = dicomdir->GetFirstPatient(); 
145    while ( pa1 ) 
146    {  // we process all the PATIENT of this DICOMDIR 
147       std::cout << pa1->GetEntryString(0x0010, 0x0010) << std::endl; // Patient's Name
148
149       st1 = pa1->GetFirstStudy();
150       while ( st1 ) 
151       { // we process all the STUDY of this patient
152          std::cout << "--- "<< st1->GetEntryString(0x0008, 0x1030) // Study Description
153          << std::endl;  
154          std::cout << " Stud.ID:[" << st1->GetEntryString(0x0020, 0x0010) // Study ID
155          << "]"; 
156
157          se1 = st1->GetFirstSerie();
158          while ( se1 ) 
159          { // we process all the SERIES of this study
160             std::cout << "--- --- "<< se1->GetEntryString(0x0008, 0x103e) << std::endl;      // Serie Description
161             std::cout << " Ser.nb:["         <<  se1->GetEntryString(0x0020, 0x0011);        // Series number
162             std::cout << "] Mod.:["          <<  se1->GetEntryString(0x0008, 0x0060) << "]"; // Modality
163
164             im1 = se1->GetFirstImage();
165             while ( im1 ) { // we process all the IMAGE of this serie
166                std::cout << "--- --- --- "<< im1->GetEntryString(0x0004, 0x1500) << std::endl; // File name
167                im1 = se1->GetNextImage();   
168             }
169             se1 = st1->GetNextSerie();   
170          }
171          st1 = pa1->GetNextStudy();
172       } 
173       pa1 = dicomdir->GetNextPatient();
174    }  
175
176    std::cout << std::endl << std::endl  
177              << " = DICOMDIR full content ====================================" 
178              << std::endl<< std::endl;
179   // dicomdir->Print();
180
181    // ------------------------- second stage ---------------------------
182     
183    // Write on disc what we read
184    dicomdir->Write("NewDICOMDIR");
185
186    std::cout << std::endl << std::endl  
187              << "NewDICOMDIR written on disc =================================" 
188              << std::endl<< std::endl;
189   // Read what we wrote  
190    gdcm::DicomDir *d2 = gdcm::DicomDir::New();
191    d2->SetFileName("NewDICOMDIR");
192    d2->Load();
193    if (!d2->IsReadable())
194    {
195       std::cout << std::endl << std::endl  
196                 << "Read NewDicomDir from disc failed ========================" 
197                 << std::endl<< std::endl;
198       return 1;
199    }
200    std::cout << std::endl << std::endl  
201              << "NewDICOMDIR successfully read from disc =================================" 
202              << std::endl<< std::endl;
203   
204    gdcm::DicomDirPatient *pa2;
205    gdcm::DicomDirStudy   *st2;
206    gdcm::DicomDirSerie   *se2;
207    gdcm::DicomDirImage   *im2;
208
209    pa1 = dicomdir->GetFirstPatient(); 
210    pa2 = d2->GetFirstPatient(); 
211
212    if (!d2)
213    {
214       std::cout << "NewDICOMDIR contains no Patient ?!?" << std::endl;
215       dicomdir->Delete();
216       d2->Delete();
217       return 1;
218    }
219    
220    while ( pa1 && pa2 ) 
221    {  // we process all the PATIENT of this DICOMDIR
222
223       if ( CompareSQItem(pa2,pa1) == 1 )
224       {
225          dicomdir->Delete();
226          d2->Delete();
227          return 1;
228       }
229   
230       // just to allow human reader to be sure ...
231       std::cout << pa2->GetEntryString(0x0010, 0x0010) 
232                 << std::endl; // Patient's Name
233  
234       st1 = pa1->GetFirstStudy();
235       st2 = pa2->GetFirstStudy();
236
237       while ( st1 && st2 )   
238       {
239          if ( CompareSQItem(st2,st1) == 1 )
240          {
241             dicomdir->Delete();
242             d2->Delete();
243             return 1;
244          }
245
246          // just to allow human reader to be sure ...
247          std::cout << "--- "<< st2->GetEntryString(0x0008, 0x1030);
248          // << std::endl;    // Study Description
249          std::cout << " Stud.ID:["          
250                    << st2->GetEntryString(0x0020, 0x0010)
251                    << "]" << std::endl; // Study ID
252   
253          se1 = st1->GetFirstSerie();
254          se2 = st2->GetFirstSerie();
255
256          while ( se1 && se2 ) 
257          { // we process all the SERIE of this study
258             if ( CompareSQItem(se2,se1) == 1 )
259               return 1;
260
261             std::cout << "--- --- " << se2->GetEntryString(0x0008, 0x103e);      // Serie Description
262             std::cout << " Ser.nb:["<< se2->GetEntryString(0x0020, 0x0011);        // Series number
263             std::cout << "] Mod.:[" << se2->GetEntryString(0x0008, 0x0060) << "]" << std::endl; // Modality
264             im1 = se1->GetFirstImage();
265             im2 = se2->GetFirstImage();
266
267             while ( im1 && im2 ) // we process all the IMAGE of this serie
268             {
269                if ( CompareSQItem(im2,im1) == 1 )
270                {
271                   dicomdir->Delete();
272                   d2->Delete();
273                   return 1; 
274                }
275
276                im1 = se1->GetNextImage();   
277                im2 = se2->GetNextImage();   
278             }
279             se1 = st1->GetNextSerie();   
280             se2 = st2->GetNextSerie();   
281          }
282          st1 = pa1->GetNextStudy();
283          st2 = pa2->GetNextStudy();
284       }
285       pa1 = dicomdir->GetNextPatient();
286       pa2 = dicomdir->GetNextPatient();
287    }
288    
289    std::cout << std::flush;
290    dicomdir->Delete();
291    d2->Delete();
292
293    return 0;
294 }