Testing/TestChangeHeader.cxx

   1 #include <iostream>
   2 //#include "gdcm.h"
   3 #include "gdcmHeader.h"
   4 #include "gdcmFile.h"
   5
   6
   7         // ecriture d'un fichier DICOM à partir d'un dcmHeader correct.
   8         // et des pixels d'une autre image
   9
  10
  11 int main(int argc, char* argv[])
  12 {
  13         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                 printf ("usage :\n testChangeEntete nomFichierPourEntete nomFichierPourDonnées \n");
  23                 exit(0);
  24         }
  25
  26         premier = argv[1];
  27         f1 = new gdcmFile(premier);
  28         f1->LoadElements();
  29
  30         deuxieme = argv[2];
  31         f2 = new gdcmFile(deuxieme);
  32         f2->LoadElements();
  33         printf("apres LoadElements\n");
  34
  35         //f1->PrintPubElVal();
  36
  37         // On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
  38
  39         f1->ReplaceOrCreateByNumber(0x0028, 0x0008, f2->GetPubElValSet().GetElValueByNumber(0x0028, 0x0008));// nb Frames
  40         f1->ReplaceOrCreateByNumber(0x0028, 0x0010, f2->GetPubElValSet().GetElValueByNumber(0x0028, 0x0010));// nbLig
  41         f1->ReplaceOrCreateByNumber(0x0028, 0x0011, f2->GetPubElValSet().GetElValueByNumber(0x0028, 0x0011));// nbCol
  42
  43         // sans doute d'autres à mettre à jour...
  44
  45         // TODO : une routine qui recoit une liste de couples (gr,el), et qui fasse le boulot.
  46
  47
  48         dataSize = f2->GetImageDataSize();
  49         printf ("dataSize %d\n",dataSize);
  50         imageData= f2->GetImageData();
  51
  52         f1->SetImageData(imageData,dataSize);
  53
  54         // ou, plus joli:
  55         //f1->SetImageData(f2->GetImageData(),f2->GetImageDataSize());
  56
  57         sprintf(resultat, "%s.vol", deuxieme.c_str());
  58         printf ("WriteDCM\n");
  59         f1->WriteDcm(resultat);
  60 }
  61
  62
  63