//#include "gdcm.h"
#include "gdcmHeader.h"
#include "gdcmFile.h"
+#include <string>
// ecriture d'un fichier DICOM à partir d'un dcmHeader correct.
// et des pixels d'une autre image
-int main(int argc, char* argv[])
-{
- string premier, deuxieme;
- char resultat[200];
+int testChangeEntete(int argc, char* argv[]) {
+ std::string premier, deuxieme;
+ char resultat[200];
- gdcmFile *f1, *f2;
+ gdcmFile *f1, *f2;
//gdcmHeader *e1, *e2;
- void* imageData;
- int dataSize;
+ void* imageData;
+ int dataSize;
- if (argc < 3) {
- printf ("usage :\n testChangeEntete nomFichierPourEntete nomFichierPourDonnées \n");
- exit(0);
- }
+ if (argc < 3) {
+ std::cerr << "usage :" << std::endl <<
+ argv[0] << " nomFichierPourEntete nomFichierPourDonnées" << std::endl;
+ return 0;
+ }
- premier = argv[1];
- f1 = new gdcmFile(premier);
+ premier = argv[1];
+ f1 = new gdcmFile(premier);
- deuxieme = argv[2];
- f2 = new gdcmFile(deuxieme);
+ deuxieme = argv[2];
+ f2 = new gdcmFile(deuxieme);
//f1->PrintPubElVal();
// On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
- f1->ReplaceOrCreateByNumber(0x0028, 0x0008,
- f2->GetPubElValByNumber(0x0028, 0x0008));// nb Frames
- f1->ReplaceOrCreateByNumber(0x0028, 0x0010,
- f2->GetPubElValByNumber(0x0028, 0x0010));// nbLig
- f1->ReplaceOrCreateByNumber(0x0028, 0x0011,
- f2->GetPubElValByNumber(0x0028, 0x0011));// nbCol
+ std::string nbFrames = f2->GetHeader()->GetEntryByNumber(0x0028, 0x0008);
+ if(nbFrames != "gdcm::Unfound") {
+ f1->GetHeader()->ReplaceOrCreateByNumber( nbFrames, 0x0028, 0x0008);
+ }
+
+ f1->GetHeader()->ReplaceOrCreateByNumber(f2->GetHeader()->GetEntryByNumber(0x0028, 0x0010),
+ 0x0028, 0x0010);// nbLig
+ f1->GetHeader()->ReplaceOrCreateByNumber( f2->GetHeader()->GetEntryByNumber(0x0028, 0x0011),
+ 0x0028, 0x0011);// nbCol
+
// sans doute d'autres à mettre à jour...
- // TODO : une routine qui recoit une liste de couples (gr,el), et qui fasse le boulot.
+ // TODO : rajouter une valeur par defaut.
+ // TODO : une routine qui recoit une liste de couples (gr,el),
+ // et qui fasse le boulot.
- dataSize = f2->GetImageDataSize();
- printf ("dataSize %d\n",dataSize);
- imageData= f2->GetImageData();
+ dataSize = f2->GetImageDataSize();
+ printf ("dataSize %d\n",dataSize);
+ imageData= f2->GetImageData();
+
+ // TODO : ne devrait-on pas fusionner ces 2 fonctions ?
+ f1->SetImageData(imageData,dataSize);
+ f1->GetHeader()->SetImageDataSize(dataSize);
- f1->SetImageData(imageData,dataSize);
+ f1->GetHeader()->PrintEntry();
- // ou, plus joli:
- //f1->SetImageData(f2->GetImageData(),f2->GetImageDataSize());
+ std::string s0 =f2->GetHeader()->GetEntryByNumber(0x7fe0, 0x0000);
+ std::string s10=f2->GetHeader()->GetEntryByNumber(0x7fe0, 0x0010);
+ printf("lgr 7fe0, 0000 %s\n",s0.c_str());
+ printf("lgr 7fe0, 0010 %s\n",s10.c_str());
+
+ sprintf(resultat, "%s.vol", deuxieme.c_str());
+ printf ("WriteDCM\n");
+ //f1->WriteDcmExplVR("itk-gdcm-ex.dcm");
+ //f1->WriteDcmImplVR(resultat);
+ f1->WriteAcr(resultat);
- sprintf(resultat, "%s.vol", deuxieme.c_str());
- printf ("WriteDCM\n");
- f1->WriteDcmImplVR(resultat);
+ return 0;
}