ENH: Adding a CT-PET matcher

diff --git a/Example/CMakeLists.txt b/Example/CMakeLists.txt
index 3b5d5a1f9797e77cf6584d483d35dcac4d4f4038..8d906038f9e5ed9af568066c31fb7be8e119c7ee 100644 (file)
--- a/Example/CMakeLists.txt
+++ b/Example/CMakeLists.txt
@@ -42,6 +42,7 @@ SET(EXAMPLE_SOURCES
   Volume2Dicom
   WriteDicomSimple
   WriteRead
+  exCTPET
 )
 
 FOREACH(name ${EXAMPLE_SOURCES})

diff --git a/Example/exCTPET.cxx b/Example/exCTPET.cxx
new file mode 100644 (file)
index 0000000..e52598b
--- /dev/null
+++ b/Example/exCTPET.cxx
@@ -0,0 +1,131 @@
+/*=========================================================================
+                                                                                
+  Program:   gdcm
+  Module:    $RCSfile: exCTPET.cxx,v $
+  Language:  C++
+  Date:      $Date: 2006/01/02 21:54:23 $
+  Version:   $Revision: 1.1 $
+                                                                                
+  Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
+  l'Image). All rights reserved. See Doc/License.txt or
+  http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
+                                                                                
+     This software is distributed WITHOUT ANY WARRANTY; without even
+     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+     PURPOSE.  See the above copyright notices for more information.
+                                                                                
+=========================================================================*/
+#include "gdcmFile.h"
+#include "gdcmSerieHelper.h"
+#include "gdcmDirList.h"
+#include "gdcmUtil.h"
+#include "gdcmDataEntry.h"
+
+int main(int argc, char *argv[])
+{
+  if(argc < 2 )
+    {
+    std::cerr << "RTFM" << std::endl;
+    return 1;
+    }
+
+  // Get the directory name
+  const char *directory = argv[1];
+#if 0
+  gdcm::SerieHelper *s;
+   std::cout << "Dir Name :[" << directory << "]" << std::endl;
+
+   s = gdcm::SerieHelper::New();
+   s->SetLoadMode(gdcm::LD_NOSHADOW | gdcm::LD_NOSEQ);
+   // 0008|0060 [CS] [Modality] [CT]   
+   gdcm::TagKey modality(0x0008,0x0060);
+   s->AddRestriction(modality, "CT", gdcm::GDCM_EQUAL); // Keep only files where
+                                              // restriction is true
+   s->SetDirectory(directory, true); // true : recursive exploration
+
+//   std::cout << " ---------------------------------------- Finish parsing :["
+//             << directory << "]" << std::endl;
+//
+//   s->Print();
+//   std::cout << " ---------------------------------------- Finish printing (1)"
+//             << std::endl;
+
+   int nbFiles;
+   // For all the Single SerieUID Files Sets of the gdcm::Serie
+   gdcm::FileList *l = s->GetFirstSingleSerieUIDFileSet();
+   while (l)
+   { 
+      nbFiles = l->size() ;
+      //if ( l->size() > 3 ) // Why not ? Just an example, for testing
+      {
+         //std::cout << "Sort list : " << nbFiles << " long" << std::endl;
+         s->OrderFileList(l);  // sort the list
+      }
+      l = s->GetNextSingleSerieUIDFileSet();
+   } 
+  s->Delete();
+#endif
+
+  // Open a file A
+  // Open another file B
+  // Same Serie/Study ?
+  // No -> return
+  // Yes -> continue
+  // Same Frame of Reference
+  // No -> Return
+  // Yes -> continue
+  // A is CT and B is PT (PET)
+  // No -> return
+  // Yes -> continue
+  // Same Image Position (no string comparison !!)
+  // Even if floating point comparison are dangerous, string comp will not work in general case
+  // No -> Return
+  // Yes: We found a match
+
+  const char filename1[] = "/tmp/PETCT/case1/WHOLEB001_CT001.dcm";
+  const char filename2[] = "/tmp/PETCT/case1/CT_PET001_CT100.dcm";
+  gdcm::File *fileRef = gdcm::File::New();
+  fileRef->SetFileName( filename2 );
+  fileRef->Load();
+  // 0008 0060 CS 1 Modality
+  std::string modalityRef = fileRef->GetEntryString(0x0008,0x0060);
+  if( modalityRef == gdcm::GDCM_UNFOUND ) return 1;
+  if ( !gdcm::Util::DicomStringEqual(modalityRef, "CT") ) return 1;
+    // 0020 000d UI 1 Study Instance UID
+    // 0020 000e UI REL Series Instance UID
+  std::string series_uid_ref = fileRef->GetEntryString(0x0020, 0x000e);
+  // 0020 0052 UI 1 Frame of Reference UID
+  std::string frame_uid_ref = fileRef->GetEntryString(0x0020, 0x0052);
+  // 0020 0032 DS 3 Image Position (Patient)
+  gdcm::DataEntry *imagePosRef = fileRef->GetDataEntry(0x0020,0x0032);
+  assert( imagePosRef->GetValueCount() == 3 );
+
+  gdcm::DirList dirList( directory, true );
+  const gdcm::DirListType filenames = dirList.GetFilenames();
+  gdcm::DirListType::const_iterator it = filenames.begin();
+  gdcm::File *file = gdcm::File::New();
+  for( ; it != filenames.end(); ++it)
+    {
+    file->SetFileName( *it );
+    file->Load();
+    std::string modality   = file->GetEntryString(0x0008,0x0060);
+    // This is a dual modality: modality should be *different*
+    if( modality == modalityRef ) continue;
+    if ( !gdcm::Util::DicomStringEqual(modality, "PT") ) continue;
+    std::string series_uid = file->GetEntryString(0x0020, 0x000e);
+    // Not same series !
+    if( series_uid == series_uid_ref ) continue;
+    std::string frame_uid = file->GetEntryString(0x0020, 0x0052);
+    if( frame_uid_ref != frame_uid ) continue;
+    gdcm::DataEntry *imagePos = file->GetDataEntry(0x0020,0x0032);
+    assert( imagePos->GetValueCount() == 3 );
+    if( imagePos->GetValue(0) == imagePosRef->GetValue(0)
+      && imagePos->GetValue(1) == imagePosRef->GetValue(1)
+      && imagePos->GetValue(2) == imagePosRef->GetValue(2) )
+      {
+      std::cerr << "We found a match: " << *it << std::endl;
+      }
+    }
+  return 0;
+}
+