3 #include "gdcmHeader.h"
8 // ecriture d'un fichier DICOM à partir d'un dcmHeader correct.
9 // et des pixels d'une autre image
12 int TestChangeHeader(int argc, char* argv[]) {
13 std::string premier, deuxieme;
17 //gdcmHeader *e1, *e2;
22 std::cerr << "usage :" << std::endl <<
23 argv[0] << " nomFichierPourEntete nomFichierPourDonnées" << std::endl;
28 f1 = new gdcmFile(premier);
31 f2 = new gdcmFile(deuxieme);
33 //f1->PrintPubElVal();
35 // On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
37 std::string nbFrames = f2->GetHeader()->GetEntryByNumber(0x0028, 0x0008);
38 if(nbFrames != "gdcm::Unfound") {
39 f1->GetHeader()->ReplaceOrCreateByNumber( nbFrames, 0x0028, 0x0008);
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
48 // sans doute d'autres à mettre à jour...
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.
55 dataSize = f2->GetImageDataSize();
56 printf ("dataSize %d\n",dataSize);
57 imageData= f2->GetImageData();
59 // TODO : ne devrait-on pas fusionner ces 2 fonctions ?
60 f1->SetImageData(imageData,dataSize);
61 f1->GetHeader()->SetImageDataSize(dataSize);
63 f1->GetHeader()->Print();
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());
70 sprintf(resultat, "%s.vol", deuxieme.c_str());
71 printf ("WriteDCM\n");
72 //f1->WriteDcmExplVR("itk-gdcm-ex.dcm");
73 //f1->WriteDcmImplVR(resultat);
74 f1->WriteAcr(resultat);