Program: gdcm
Module: $RCSfile: ToInTag.cxx,v $
Language: C++
- Date: $Date: 2006/05/31 16:23:25 $
- Version: $Revision: 1.1 $
+ Date: $Date: 2007/07/04 15:40:51 $
+ Version: $Revision: 1.16 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
* according to their Patient/Study/Serie/Image characteristics
* - fills a single level Directory with *all* the files,
* converted into a Brucker-like Dicom, Intags compliant
- *
*/
-typedef std::map<std::string, gdcm::File*> SortedFiles;
+typedef std::map<std::string, GDCM_NAME_SPACE::File*> SortedFiles;
int main(int argc, char *argv[])
{
START_USAGE(usage)
- " \n ToInTag :\n ",
+ " \n ToInTag :\n ",
" - explores recursively the given directory, ",
" - keeps the requested series/ drops the unrequested series ",
" - orders the gdcm-readable found Files according to their ",
" - fills a single level (*) Directory with *all* the files, ",
" converted into a Brucker-like Dicom, InTags compliant ",
" (*) actually : creates as many directories as Patients ",
- " -that shouldn't appear, but being carefull is better ! -",
+ " -that shouldn't appear, but being carefull is better!- ",
" or ",
" - fills a tree-like structure of directories as : ",
" - Patient ",
" -- Serie ",
" --- Position ",
- " Images are (sorted by Trigger Time / ",
- " Encoding Direction (Row, Column) ",
- " use : ",
+ " ---- Images (sorted by Trigger Time / ",
+ " Encoding Direction (Row, Column) ",
+ " ",
+ " Note : when (0008|1090) [Model Name ] equals 'TrioTim ' : ",
+ " - (0008|103e)[Series Description ] is checked for ",
+ " '90' (-> COL) or '0' (-> ROW) ",
+ " - (0x0020, 0x000e) [Series Instance UID] is NOT dealt with, ",
+ " since row an col tagging are in 2 different Series ",
+ " DO NOT supply a directory holding different exams ",
+ " for the same Slice level! ",
+ " uses : ",
" 0x0021, 0x1020 : 'SLICE INDEX' ",
" 0x0021, 0x1040 : 'FRAME INDEX' ",
" 0x0020, 0x0012 : 'SESSION INDEX' (Acquisition Number) ",
+ " ",
" usage: ",
- " PhilipsToBrucker dirin=rootDirectoryName ",
+ " ----- ",
+ " ToInTag dirin=rootDirectoryName ",
" dirout=outputDirectoryName ",
" { [keep= list of seriesNumber to process] ",
" | [drop= list of seriesNumber to ignore] } ",
" [taggrid] [skel] ",
" [input = {ACR|DCM}] ",
" [extent=image suffix (.IMA, .NEMA, .DCM, ...)] ",
- " [listonly] [split] ",
+ " [listonly] [split] [rubout] ",
" [noshadowseq][noshadow][noseq] [verbose] [debug] ",
" ",
" dirout : will be created if doesn't exist ",
" he gives the list of 'SeriesNumber' (tag 0020|0011) ",
" SeriesNumber are short enough to be human readable ",
" e.g : 1030,1035,1043 ",
- " taggrid : user knows all the images are 'grid' -ie : not 'col', not 'raw'",
+ " taggrid : user knows all the images are 'grid' -ie: not 'col', not 'raw'-",
" extent : DO NOT forget the leading '.' ! ",
" skel: name skeleton eg : patName_1.nema -> skel=patName_ ",
" split: creates a tree-like structure of directories as : ",
" - Patient ",
" -- Serie ",
" --- Position ",
- " Images are (sorted by Trigger Time / ",
- " Encoding Direction (Row, Column) ",
+ " ---- Images (sorted by Trigger Time / ",
+ " Encoding Direction (Row, Column) ",
+ " rubout : user asks to rubout burnt-in image number ",
" noshadowseq: user doesn't want to load Private Sequences ",
" noshadow : user doesn't want to load Private groups (odd number) ",
" noseq : user doesn't want to load Sequences ",
" debug : *developer* wants to run the program in 'debug mode' ",
FINISH_USAGE
+ std::cout << "inside ToInTag" << std::endl;
+
// ----- Initialize Arguments Manager ------
- gdcm::ArgMgr *am = new gdcm::ArgMgr(argc, argv);
+ GDCM_NAME_SPACE::ArgMgr *am = new GDCM_NAME_SPACE::ArgMgr(argc, argv);
if (argc == 1 || am->ArgMgrDefined("usage"))
{
const char *dirNameout;
dirNameout = am->ArgMgrGetString("dirout",".");
- int loadMode = gdcm::LD_ALL;
+ int loadMode = GDCM_NAME_SPACE::LD_ALL;
if ( am->ArgMgrDefined("noshadowseq") )
- loadMode |= gdcm::LD_NOSHADOWSEQ;
+ loadMode |= GDCM_NAME_SPACE::LD_NOSHADOWSEQ;
else
{
if ( am->ArgMgrDefined("noshadow") )
- loadMode |= gdcm::LD_NOSHADOW;
+ loadMode |= GDCM_NAME_SPACE::LD_NOSHADOW;
if ( am->ArgMgrDefined("noseq") )
- loadMode |= gdcm::LD_NOSEQ;
+ loadMode |= GDCM_NAME_SPACE::LD_NOSEQ;
}
if (am->ArgMgrDefined("debug"))
- gdcm::Debug::DebugOn();
+ GDCM_NAME_SPACE::Debug::DebugOn();
- bool verbose = am->ArgMgrDefined("verbose");
- bool split = am->ArgMgrDefined("split");
- bool listonly = am->ArgMgrDefined("listonly");
-
+ int verbose = am->ArgMgrDefined("verbose");
+ int split = am->ArgMgrDefined("split");
+ int listonly = am->ArgMgrDefined("listonly");
+
+ bool rubout = ( 0 != am->ArgMgrDefined("rubout") );
+
int nbSeriesToKeep;
int *seriesToKeep = am->ArgMgrGetListOfInt("keep", &nbSeriesToKeep);
int nbSeriesToDrop;
// ----- Begin Processing -----
- if ( ! gdcm::DirList::IsDirectory(dirNamein) )
+ if ( ! GDCM_NAME_SPACE::DirList::IsDirectory(dirNamein) )
{
std::cout << "KO : [" << dirNamein << "] is not a Directory." << std::endl;
- exit(0);
+ return 0;
+
}
else
{
std::cout << "Check for output directory :[" << dirNameout << "]."
<<std::endl;
- if ( ! gdcm::DirList::IsDirectory(dirNameout) ) // dirout not found
+ if ( ! GDCM_NAME_SPACE::DirList::IsDirectory(dirNameout) ) // dirout not found
{
std::string strDirNameout(dirNameout); // to please gcc 4
systemCommand = "mkdir " +strDirNameout; // create it!
if (verbose)
std::cout << systemCommand << std::endl;
system (systemCommand.c_str());
- if ( ! gdcm::DirList::IsDirectory(dirNameout) ) // be sure it worked
+ if ( ! GDCM_NAME_SPACE::DirList::IsDirectory(dirNameout) ) // be sure it worked
{
std::cout << "KO : not a dir : [" << dirNameout << "] (creation failure ?)" << std::endl;
- exit(0);
+ return 0;
+
}
else
{
}
std::string strDirNamein(dirNamein);
- gdcm::DirList dirList(strDirNamein, true); // get recursively the list of files
+ GDCM_NAME_SPACE::DirList dirList(strDirNamein, true); // get recursively the list of files
if (listonly)
{
std::cout << "------------List of found files ------------" << std::endl;
dirList.Print();
+ std::cout << std::endl;
}
- gdcm::DirListType fileNames;
+ GDCM_NAME_SPACE::DirListType fileNames;
fileNames = dirList.GetFilenames();
- gdcm::SerieHelper *s; // Needed only to may use SerieHelper::AddSeriesDetail()
- s = gdcm::SerieHelper::New();
+ GDCM_NAME_SPACE::SerieHelper *s; // Needed to use SerieHelper::AddSeriesDetail()
+ s = GDCM_NAME_SPACE::SerieHelper::New();
std::string token = "%%%"; // Hope it's enough!
/*
s->Print();
*/
- gdcm::File *f;
- gdcm::FileHelper *fh;
+ GDCM_NAME_SPACE::File *f;
+ GDCM_NAME_SPACE::FileHelper *fh;
std::vector<std::string> tokens;
std::vector<std::string> tokensForFileName;
+
+ // For Siemens pb, we need Manufacturer's Model Name
+
+ GDCM_NAME_SPACE::DirListType::iterator it1 = fileNames.begin();
+ f = GDCM_NAME_SPACE::File::New();
+ f->SetLoadMode(GDCM_NAME_SPACE::LD_ALL);
+ f->SetFileName( *it1 );
+ f->Load();
+ std::string modelName = f->GetEntryString(0x0008,0x1090);
+ f->Delete();
+
/*
std::cout << "---------------Print Unique Series identifiers---------"
<< std::endl;
std::string uniqueSeriesIdentifier;
- for (gdcm::DirListType::iterator it) = fileNames.begin();
+ for (GDCM_NAME_SPACE::DirListType::iterator it) = fileNames.begin();
it != fileNames.end();
++it)
{
std::cout << "File Name : " << *it << std::endl;
- f = gdcm::File::New();
- f->SetLoadMode(gdcm::LD_ALL);
+ f = GDCM_NAME_SPACE::File::New();
+ f->SetLoadMode(GDCM_NAME_SPACE::LD_ALL);
f->SetFileName( *it );
f->Load();
SortedFiles sf;
s->AddSeriesDetail(0x0010, 0x0010, false); // Patient's Name
- s->AddSeriesDetail(0x0020, 0x000e, false); // Series Instance UID
+ // for Siemens TrioTim, don't deal with 'Series Instance UID'
+ if ( !GDCM_NAME_SPACE::Util::DicomStringEqual(modelName,"TrioTim") )
+ s->AddSeriesDetail(0x0020, 0x000e, false); // Series Instance UID
+ else
+ s->AddSeriesDetail(0x9999, 0x9999, false); // dirty trick to ignore 'Series Instance UID'
s->AddSeriesDetail(0x0020, 0x0032, false); // Image Position (Patient)
s->AddSeriesDetail(0x0018, 0x1060, true); // Trigger Time (true: convert to keep numerical order)
s->AddSeriesDetail(0x0018, 0x1312, false); // In-plane Phase Encoding Direction
-
- for (gdcm::DirListType::iterator it = fileNames.begin();
+ s->AddSeriesDetail(0x0008, 0x103e, false); // Series Description (special Siemens ...)
+
+ //uint8_t *imageData; // Useless : pixels will not be loaded
+ // (images are overwritten)
+
+ for (GDCM_NAME_SPACE::DirListType::iterator it = fileNames.begin();
it != fileNames.end();
++it)
{
- f = gdcm::File::New();
+ f = GDCM_NAME_SPACE::File::New();
f->SetLoadMode(loadMode);
f->SetFileName( *it );
f->Load();
-
std::string strSeriesNumber;
int seriesNumber;
f->Delete();
continue;
}
- }
+ }
userFileIdentifier=s->CreateUserDefinedFileIdentifier(f);
tokens.clear();
- gdcm::Util::Tokenize (userFileIdentifier, tokens, token);
+ GDCM_NAME_SPACE::Util::Tokenize (userFileIdentifier, tokens, token);
int imageNum; // Within FileName
char newName[1024];
- if ( tokens[3] == "gdcmUnfound") // sometimes Trigger Time is not found. CreateUserDefinedFileIdentifier is not aware of the pb.
+ if ( tokens[3] == GDCM_NAME_SPACE::GDCM_UNFOUND) // sometimes Trigger Time is not found. CreateUserDefinedFileIdentifier is not aware of the pb.
{
///this is a trick to build up a lexicographical compliant name :
/// eg : fich001.ima vs fich100.ima as opposed to fich1.ima vs fich100.ima
- std::string name = gdcm::Util::GetName( *it );
+ std::string name = GDCM_NAME_SPACE::Util::GetName( *it );
if (hasSkel)
{
- gdcm::Util::Tokenize (name, tokensForFileName, skel);
+ GDCM_NAME_SPACE::Util::Tokenize (name, tokensForFileName, skel);
imageNum = atoi ( tokensForFileName[0].c_str() );
// probabely we could write something much more complicated using C++ !
sprintf (newName, "%s%06d%s", skel, imageNum, extent);
userFileIdentifier = tokens[0] + token + tokens[1] + token + tokens[2] + token
- + tokens[3] + token + tokens[4] + token;
+ + tokens[3] + token + tokens[4] + token + tokens[5] + token;
}
if (verbose)
- std::cout << "[" << userFileIdentifier << "]" << std::endl;
+ std::cout << "[" << userFileIdentifier << "] : " << *it << std::endl;
// storing in a map ensures automatic sorting !
sf[userFileIdentifier] = f;
}
+
+ if (verbose)
+ std::cout << " " << std::endl;
std::string fullFilename, lastFilename;
std::string previousPatientName, currentPatientName;
std::string previousImagePosition, currentImagePosition;
std::string previousPhaseEncodingDirection, currentPhaseEncodingDirection;
std::string previousTriggerTime, currentTriggerTime;
+
+ std::string currentStudyUID;
+ std::string seriesDescription;
std::string writeDir, currentWriteDir;
std::string currentPatientWriteDir, currentSerieWriteDir,
std::string fullWriteFilename;
std::string strExtent(extent);
- writeDir = gdcm::Util::NormalizePath(dirNameout);
+ writeDir = GDCM_NAME_SPACE::Util::NormalizePath(dirNameout);
SortedFiles::iterator it2;
previousPatientName = "";
previousTriggerTime = "";
int sliceIndex = 0; // Is incremented *at the beginning* of processing
- int frameIndex = 1;
+ int frameIndex;
+ if (taggrid)
+ frameIndex = 0;
+ else
+ frameIndex = 1;
+
int flag = 0;
- gdcm::File *currentFile;
+ GDCM_NAME_SPACE::File *currentFile;
+
+ std::string defaultStudyUID = GDCM_NAME_SPACE::Util::CreateUniqueUID();
+ std::string defaultSerieUID;
for (it2 = sf.begin() ; it2 != sf.end(); ++it2)
{
currentFile = it2->second;
fullFilename = currentFile->GetFileName();
- lastFilename = gdcm::Util::GetName( fullFilename );
- std::cout << "Rewrite [" <<lastFilename << "]" << std::endl;
+ lastFilename = GDCM_NAME_SPACE::Util::GetName( fullFilename );
+ if (verbose)
+ std::cout <<" ------------------------------------------------------------------------------"
+ << std::endl << " Deal with [" << it2->first << "] : ["<<fullFilename << "]"
+ << std::endl;
tokens.clear();
- gdcm::Util::Tokenize (it2->first, tokens, token);
+ GDCM_NAME_SPACE::Util::Tokenize (it2->first, tokens, token);
currentPatientName = tokens[0];
currentSerieInstanceUID = tokens[1];
currentImagePosition = tokens[2];
currentTriggerTime = tokens[3];
- currentPhaseEncodingDirection = tokens[4];
+ currentPhaseEncodingDirection = tokens[4];
+ seriesDescription = tokens[5]; // For Siemens pb
if ( currentImagePosition[0] == '-')
currentImagePosition[0] = 'M';
if ( currentImagePosition[0] == '+')
currentImagePosition[0] = 'P';
+
+ // Add a default ImagePositionPatient to avoid confusion at post processing time
+ if ( currentFile->GetEntryString(0x0020,0x0032) == GDCM_NAME_SPACE::GDCM_UNFOUND &&
+ currentFile->GetEntryString(0x0020,0x0030) == GDCM_NAME_SPACE::GDCM_UNFOUND )
+ {
+ currentFile->InsertEntryString("0.\\0.\\0.",0x0020, 0x0032, "DS" );
+ }
+
+ // Add a default ImagePositionPatient to avoid confusion at post processing time
+ if ( currentFile->GetEntryString(0x0020,0x0037) == GDCM_NAME_SPACE::GDCM_UNFOUND &&
+ currentFile->GetEntryString(0x0020,0x0035) == GDCM_NAME_SPACE::GDCM_UNFOUND )
+ {
+ currentFile->InsertEntryString("1.\\0.\\0.\\0.\\1.\\0.",0x0020, 0x0037, "DS" );
+ }
if (previousPatientName != currentPatientName)
- {
+ {
+ if ( currentFile->GetEntryString(0x0020,0x000d) == GDCM_NAME_SPACE::GDCM_UNFOUND)
+ {
+ if (verbose)
+ std::cout << "--- new Study UID created" << std::endl;
+ defaultStudyUID = GDCM_NAME_SPACE::Util::CreateUniqueUID();
+ currentFile->InsertEntryString(defaultStudyUID, 0x0020, 0x000d, "UI" );
+ }
+
previousPatientName = currentPatientName;
if (verbose)
std::cout << "==== new Patient [" << currentPatientName << "]" << std::endl;
previousPhaseEncodingDirection = ""; //currentPhaseEncodingDirection;
currentPatientWriteDir = writeDir + currentPatientName;
- //if ( ! gdcm::DirList::IsDirectory(currentPatientWriteDir) )
+ //if ( ! GDCM_NAME_SPACE::DirList::IsDirectory(currentPatientWriteDir) )
{
systemCommand = "mkdir " + currentPatientWriteDir;
if (verbose)
}
}
+ if ( GDCM_NAME_SPACE::Util::DicomStringEqual(modelName,"TrioTim") ) // for Siemens TrioTim , don't deal with 'Series Instance UID'
+
if (previousSerieInstanceUID != currentSerieInstanceUID)
{
if (verbose)
std::cout << "==== === new Serie [" << currentSerieInstanceUID << "]"
<< std::endl;
+
+ if ( currentFile->GetEntryString(0x0020,0x000e) == GDCM_NAME_SPACE::GDCM_UNFOUND)
+ {
+ if (verbose)
+ std::cout << "--- --- new Serie UID created" << std::endl;
+ defaultSerieUID = GDCM_NAME_SPACE::Util::CreateUniqueUID();
+ currentFile->InsertEntryString(defaultSerieUID, 0x0020, 0x000e, "UI" );
+ }
+
if (split)
{
- currentSerieWriteDir = currentPatientWriteDir + gdcm::GDCM_FILESEPARATOR
+ currentSerieWriteDir = currentPatientWriteDir + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
+ currentSerieInstanceUID;
systemCommand = "mkdir " + currentSerieWriteDir;
system (systemCommand.c_str());
previousImagePosition = ""; //currentImagePosition;
previousPhaseEncodingDirection = ""; //currentPhaseEncodingDirection;
}
-
+ // end of modelName != "TrioTim "
+
if (previousImagePosition != currentImagePosition)
{
frameIndex = 1;
<< std::endl;
if (split)
{
- currentPositionWriteDir = currentSerieWriteDir + gdcm::GDCM_FILESEPARATOR
+ currentPositionWriteDir = currentSerieWriteDir + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
+ currentImagePosition;
systemCommand = "mkdir " + currentPositionWriteDir;
system (systemCommand.c_str());
sliceIndex = 1; // only *one* slice in a given directory
else
sliceIndex += 1;
- }
+ }
// We don't split on Row/Column!
/*
if (split)
{
currentPhaseEncodingDirectionWriteDir = currentPositionWriteDir
- + gdcm::GDCM_FILESEPARATOR
+ + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
+ currentPhaseEncodingDirection;
systemCommand = "mkdir " + currentPhaseEncodingDirectionWriteDir;
system (systemCommand.c_str());
previousPhaseEncodingDirection = currentPhaseEncodingDirection;
}
*/
-
- if (verbose)
- std::cout << "--- --- --- --- --- " << (it2->second)->GetFileName()
- << std::endl;
- if ( gdcm::Debug::GetDebugFlag())
+ if ( GDCM_NAME_SPACE::Debug::GetDebugFlag())
std::cout << "--- --- --- --- --- " << it2->first << " "
<< (it2->second)->GetFileName() << " "
- << gdcm::Util::GetName( fullFilename ) << std::endl;
+ << GDCM_NAME_SPACE::Util::GetName( fullFilename ) << std::endl;
// Transform the image to be 'Brucker-Like'
// ----------------------------------------
// Deal with 0x0020, 0x0012 : 'SESSION INDEX' (Acquisition Number)
std::string chSessionIndex;
// CLEANME
- if (taggrid)
+
+ if (taggrid)
+ {
chSessionIndex = "1";
+ }
else
{
- if (currentPhaseEncodingDirection == "COL" || currentPhaseEncodingDirection == "COL " || currentPhaseEncodingDirection == " COL")
- chSessionIndex = "1";
- else if (currentPhaseEncodingDirection == "ROW" || currentPhaseEncodingDirection == "ROW "|| currentPhaseEncodingDirection == " ROW")
- chSessionIndex = "2";
- else
+ /* for SIEMENS MRI :
+ D 0008|1090 [LO] [Manufacturer's Model Name ] [Triotim ]
+ we have to deal with :
+
+ D 0008|103e [LO] [Series Description ] [fl2d9_line PA 15 90deg] or anything that contains '90' !
+ D 0008|103e [LO] [Series Description ] [fl2d9_line PA 15 0deg ]
+ (everything is flagged as 'ROW')
+ */
+
+ if ( GDCM_NAME_SPACE::Util::DicomStringEqual(modelName,"TrioTim") )
{
- std::cout << "====================== PhaseEncodingDirection "
- << " neither COL nor ROW (?!?) : [ "
- << currentPhaseEncodingDirection << "]" << std::endl;
- chSessionIndex = "1";
+ if (seriesDescription.find("90", 0) != std::string::npos)
+ chSessionIndex = "1"; // 90 deg -> COL
+ else if (seriesDescription.find("0", 0)!= std::string::npos)
+ chSessionIndex = "2"; // 0 deg -> ROW
+ else
+ {
+ std::cout << "====================== seriesDescription doesn't contain"
+ << " neither '90' nor '0' (?!?) : ["
+ << seriesDescription << "]" << std::endl;
+ chSessionIndex = "1";
+ }
+ }
+ else // for all other 'normal' cases
+ {
+ if (currentPhaseEncodingDirection == "COL" || currentPhaseEncodingDirection == "COL " || currentPhaseEncodingDirection == " COL")
+ chSessionIndex = "1";
+ else if (currentPhaseEncodingDirection == "ROW" || currentPhaseEncodingDirection == "ROW "|| currentPhaseEncodingDirection == " ROW")
+ chSessionIndex = "2";
+ else
+ {
+ std::cout << "====================== PhaseEncodingDirection "
+ << " neither COL nor ROW (?!?) : [ "
+ << currentPhaseEncodingDirection << "]" << std::endl;
+ chSessionIndex = "1";
+ }
}
}
- if (currentFile->IsVRCoherent(0x0020) == 1 )
+
+ if (currentFile->IsVRCoherent(0x0020) == 1 )
currentFile->InsertEntryString(chSessionIndex, 0x0020, 0x0012, " ");
- else
+ else
currentFile->InsertEntryString(chSessionIndex, 0x0020, 0x0012, "IS");
// Deal with 0x0021, 0x1020 : 'SLICE INDEX'
stringVR = "IS";
currentFile->InsertEntryString(strChSliceIndex, 0x0021, 0x1020, stringVR);
- currentFile->InsertEntryString(chFrameIndex, 0x0021, 0x1040, stringVR);
+ currentFile->InsertEntryString(chFrameIndex, 0x0021, 0x1040, stringVR);
+
+ if (verbose) {
+ std::cout << "0x0021, 0x1020 : strChSliceIndex " << strChSliceIndex << std::endl;
+ std::cout << "0x0021, 0x1040 : chFrameIndex " << chFrameIndex << std::endl;
+ std::cout << "0x0020, 0x0012 : chSessionIndex " << chSessionIndex << std::endl;
+ }
+
+ std::string strImagePositionPatient = currentFile->GetEntryString(0x0020, 0x0032 );
+ if (strImagePositionPatient == GDCM_NAME_SPACE::GDCM_UNFOUND)
+ {
+ if (verbose)
+ std::cout << "Duplicate ImagePosition into ImagePositionPatient" << std::endl;
+ currentFile->InsertEntryString(currentFile->GetEntryString(0x0020, 0x0030), 0x0020, 0x0032, "DS" );
+ }
+ std::string strImageOrientationPatient = f->GetEntryString(0x0020, 0x0037 );
+ if (strImageOrientationPatient == GDCM_NAME_SPACE::GDCM_UNFOUND)
+ {
+ if (verbose)
+ std::cout << "Duplicate ImageOrientation into ImageOrientationPatient" << std::endl;
+ currentFile->InsertEntryString(currentFile->GetEntryString(0x0020, 0x0035), 0x0020, 0x0037, "DS" );
+ }
+
if (taggrid)
frameIndex++;
else
if (split)
- //fullWriteFilename = currentPhaseEncodingDirectionWriteDir + gdcm::GDCM_FILESEPARATOR
+ //fullWriteFilename = currentPhaseEncodingDirectionWriteDir + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
// + lastFilename + strExtent;
- fullWriteFilename = currentPositionWriteDir + gdcm::GDCM_FILESEPARATOR
+ fullWriteFilename = currentPositionWriteDir + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
+ lastFilename + strExtent;
else
- fullWriteFilename = currentPatientWriteDir + gdcm::GDCM_FILESEPARATOR
+ fullWriteFilename = currentPatientWriteDir + GDCM_NAME_SPACE::GDCM_FILESEPARATOR
+ lastFilename + strExtent;
-
- /*
- systemCommand = "cp " + fullFilename + " " + fullWriteFilename;
- std::cout << systemCommand << std::endl;
- system ( systemCommand.c_str() );
- */
// Load the pixels in RAM.
- fh = gdcm::FileHelper::New(currentFile);
- fh->GetImageDataRaw(); // Don't convert (Gray Pixels + LUT) into (RGB pixels) ?!?
- fh->SetWriteTypeToDcmExplVR();
- // We didn't modify pixels -> keep unchanged the following :
+ fh = GDCM_NAME_SPACE::FileHelper::New(currentFile);
+ uint8_t *imageData = fh->GetImageDataRaw(); // Don't convert (Gray Pixels + LUT) into (RGB pixels) ?!?
+ fh->SetWriteTypeToDcmExplVR();
+
+ if (rubout) {
+ // Put to Black the burnt-in number.
+ nX = currentFile->GetXSize();
+ nY = currentFile->GetYSize();
+ for(int y=nY-15; y<nY; y++)
+ for(int x=nX/3; x<nX/2+50; x++)
+ imageData[ y*nX*2 + x ] = 0;
+ }
+
+ // We didn't make any computation on the pixels -> keep unchanged the following :
// 'Media Storage SOP Class UID' (0x0002,0x0002)
// 'SOP Class UID' (0x0008,0x0016)
// 'Image Type' (0x0008,0x0008)
- // 'Conversion Type' (0x0008,0x0064)
- fh->SetContentType(gdcm::UNMODIFIED_PIXELS_IMAGE);
+ // 'Conversion Type' (0x0008,0x0064)
+ fh->SetContentType(GDCM_NAME_SPACE::UNMODIFIED_PIXELS_IMAGE);
+
if (!fh->Write(fullWriteFilename))
{
std::cout << "Fail to write :[" << fullWriteFilename << "]"
}
fh->Delete();
}
+ return 0;
}
+