#include <iostream>
#include "gdcm.h"


	// 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];
	
	gdcmFile  *f1, *f2;
	//gdcmHeader *e1, *e2;
	void* imageData;
	int dataSize;
		
	if (argc < 3) {
		printf ("usage :\n testChangeEntete nomFichierPourEntete nomFichierPourDonnées \n");
		exit(0);	
	}

	premier = argv[1];
	f1 = new gdcmFile(premier);
	f1->LoadElements();
	
	deuxieme = argv[2];
	f2 = new gdcmFile(deuxieme);
	f2->LoadElements();
	printf("apres LoadElements\n");
	
	//f1->PrintPubElVal();
	
	// On suppose que les champs DICOM du 2ieme fichier existent *effectivement*
	
	f1->ReplaceOrCreateByNumber(0x0028, 0x0008, f2->GetPubElValSet().GetElValueByNumber(0x0028, 0x0008));// nb Frames
	f1->ReplaceOrCreateByNumber(0x0028, 0x0010, f2->GetPubElValSet().GetElValueByNumber(0x0028, 0x0010));// nbLig
	f1->ReplaceOrCreateByNumber(0x0028, 0x0011, f2->GetPubElValSet().GetElValueByNumber(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.

				
	dataSize = f2->GetImageDataSize();
	printf ("dataSize %d\n",dataSize);
	imageData= f2->GetImageData();
	
	f1->SetImageData(imageData,dataSize);
	
	// ou, plus joli:
	//f1->SetImageData(f2->GetImageData(),f2->GetImageDataSize());	

	sprintf(resultat, "%s.vol", deuxieme.c_str());
	printf ("WriteDCM\n");
	f1->WriteDcm(resultat);		
}