Testing/TestChangeHeader.cxx

   1 #include <iostream>
   2 //#include "gdcm.h"
   3 #include "gdcmHeader.h"
   4 #include "gdcmFile.h"
   5 #include <string>
   6
   7
   8         // ecriture d'un fichier DICOM à partir d'un dcmHeader correct.
   9         // et des pixels d'une autre image
  10
  11
  12 int testChangeEntete(int argc, char* argv[]) {
  13         std::string premier, deuxieme;
  14         char resultat[200];
  15
  16         gdcmFile  *f1, *f2;
  17         //gdcmHeader *e1, *e2;
  18         void* imageData;
  19         int dataSize;
  20
  21         if (argc < 3) {
  22                 std::cerr << "usage :" << std::endl <<
  23       argv[0] << " nomFichierPourEntete nomFichierPourDonnées" << std::endl;
  24                 return 0;
  25         }
  26
  27         premier = argv[1];
  28         f1 = new gdcmFile(premier);
  29
  30         deuxieme = argv[2];
  31         f2 = new gdcmFile(deuxieme);
  32
  33         //f1->PrintPubElVal();
  34
  35         // On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
  36
  37         std::string nbFrames = f2->GetHeader()->GetEntryByNumber(0x0028, 0x0008);
  38         if(nbFrames != "gdcm::Unfound") {
  39            f1->GetHeader()->ReplaceOrCreateByNumber( nbFrames, 0x0028, 0x0008);
  40         }
  41
  42         f1->GetHeader()->ReplaceOrCreateByNumber(f2->GetHeader()->GetEntryByNumber(0x0028, 0x0010),
  43           0x0028, 0x0010);// nbLig
  44         f1->GetHeader()->ReplaceOrCreateByNumber( f2->GetHeader()->GetEntryByNumber(0x0028, 0x0011),
  45           0x0028, 0x0011);// nbCol
  46
  47
  48         // sans doute d'autres à mettre à jour...
  49
  50         // TODO : rajouter une valeur par defaut.
  51         // TODO : une routine qui recoit une liste de couples (gr,el),
  52         //                        et qui fasse le boulot.
  53
  54
  55         dataSize = f2->GetImageDataSize();
  56         printf ("dataSize %d\n",dataSize);
  57         imageData= f2->GetImageData();
  58
  59         // TODO : ne devrait-on pas fusionner ces 2 fonctions ?
  60         f1->SetImageData(imageData,dataSize);
  61         f1->GetHeader()->SetImageDataSize(dataSize);
  62
  63         f1->GetHeader()->PrintEntry();
  64
  65         std::string s0 =f2->GetHeader()->GetEntryByNumber(0x7fe0, 0x0000);
  66         std::string s10=f2->GetHeader()->GetEntryByNumber(0x7fe0, 0x0010);
  67         printf("lgr 7fe0, 0000 %s\n",s0.c_str());
  68         printf("lgr 7fe0, 0010 %s\n",s10.c_str());
  69
  70         sprintf(resultat, "%s.vol", deuxieme.c_str());
  71         printf ("WriteDCM\n");
  72         //f1->WriteDcmImplVR(resultat);
  73         f1->WriteAcr(resultat);
  74
  75   return 0;
  76 }
  77
  78
  79