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 main(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                 printf ("usage :\n testChangeEntete nomFichierPourEntete nomFichierPourDonnées \n");
  23                 exit(0);
  24         }
  25
  26         premier = argv[1];
  27         f1 = new gdcmFile(premier);
  28
  29         deuxieme = argv[2];
  30         f2 = new gdcmFile(deuxieme);
  31
  32         //f1->PrintPubElVal();
  33
  34         // On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
  35
  36         std::string nbFrames = f2->GetPubElValByNumber(0x0028, 0x0008);
  37         if(nbFrames != "gdcm::Unfound") {
  38            f1->ReplaceOrCreateByNumber( nbFrames, 0x0028, 0x0008);
  39         }
  40
  41         f1->ReplaceOrCreateByNumber( f2->GetPubElValByNumber(0x0028, 0x0010),
  42           0x0028, 0x0010);// nbLig
  43         f1->ReplaceOrCreateByNumber( f2->GetPubElValByNumber(0x0028, 0x0011),
  44           0x0028, 0x0011);// nbCol
  45
  46
  47         // sans doute d'autres à mettre à jour...
  48
  49         // TODO : rajouter une valeur par defaut.
  50         // TODO : une routine qui recoit une liste de couples (gr,el), et qui fasse le boulot.
  51
  52
  53         dataSize = f2->GetImageDataSize();
  54         printf ("dataSize %d\n",dataSize);
  55         imageData= f2->GetImageData();
  56
  57         // TODO : ne devrait-on pas fusionner ces 2 fonctions ?
  58         f1->SetImageData(imageData,dataSize);
  59         f1->SetImageDataSize(dataSize);
  60
  61         f1->PrintPubElVal();
  62
  63         // ou, plus joli:
  64         //f1->SetImageData(f2->GetImageData(),f2->GetImageDataSize());
  65
  66         std::string s0 =f2->GetPubElValByNumber(0x7fe0, 0x0000);
  67         std::string s10=f2->GetPubElValByNumber(0x7fe0, 0x0010);
  68         printf("lgr 7fe0, 0000 %s\n",s0.c_str());
  69         printf("lgr 7fe0, 0010 %s\n",s10.c_str());
  70
  71         sprintf(resultat, "%s.vol", deuxieme.c_str());
  72         printf ("WriteDCM\n");
  73         //f1->WriteDcmImplVR(resultat);
  74         f1->WriteAcr(resultat);
  75 }
  76
  77
  78