1 /*=========================================================================
4 Module: $RCSfile: ReWrite.cxx,v $
6 Date: $Date: 2005/08/30 15:13:05 $
7 Version: $Revision: 1.11 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
19 #include "gdcmFileHelper.h"
20 #include "gdcmDebug.h"
22 #include "gdcmArgMgr.h"
26 int main(int argc, char *argv[])
30 " Re write a full gdcm-readable Dicom image ",
31 " (usefull when the file header is not very straight). ",
33 " usage: ReWrite filein=inputFileName fileout=anonymizedFileName ",
34 " [mode=write mode] [noshadow] [noseq][debug] ",
36 " mode = a (ACR), x (Explicit VR Dicom), r (RAW : only pixels) ",
37 " noshadowseq: user doesn't want to load Private Sequences ",
38 " noshadow : user doesn't want to load Private groups (odd number)",
39 " noseq : user doesn't want to load Sequences ",
40 " rgb : user wants to tranform LUT (if any) to RGB pixels ",
41 " debug : user wants to run the program in 'debug mode' ",
44 // ----- Initialize Arguments Manager ------
45 gdcm::ArgMgr *am = new gdcm::ArgMgr(argc, argv);
47 if (argc == 1 || am->ArgMgrDefined("usage"))
49 am->ArgMgrUsage(usage); // Display 'usage'
53 char *fileName = am->ArgMgrWantString("filein",usage);
54 if ( fileName == NULL )
56 std::cout << "'filein= ...' is mandatory" << std::endl;
61 char *outputFileName = am->ArgMgrWantString("fileout",usage);
62 if ( outputFileName == NULL )
64 std::cout << "'fileout= ...' is mandatory" << std::endl;
69 char *mode = am->ArgMgrGetString("mode",(char *)"X");
71 int loadMode = gdcm::LD_ALL;
72 if ( am->ArgMgrDefined("noshadowseq") )
73 loadMode |= gdcm::LD_NOSHADOWSEQ;
76 if ( am->ArgMgrDefined("noshadow") )
77 loadMode |= gdcm::LD_NOSHADOW;
78 if ( am->ArgMgrDefined("noseq") )
79 loadMode |= gdcm::LD_NOSEQ;
82 bool rgb = ( 0 != am->ArgMgrDefined("RGB") );
84 if (am->ArgMgrDefined("debug"))
85 gdcm::Debug::DebugOn();
87 // if unused Params we give up
88 if ( am->ArgMgrPrintUnusedLabels() )
90 am->ArgMgrUsage(usage);
95 delete am; // we don't need Argument Manager any longer
97 // ----------- End Arguments Manager ---------
99 gdcm::File *f = new gdcm::File();
100 f->SetLoadMode( loadMode );
101 f->SetFileName( fileName );
102 bool res = f->Load();
109 if (!f->IsReadable())
111 std::cerr << "Sorry, not a Readable DICOM / ACR File" <<std::endl;
116 gdcm::FileHelper *fh = new gdcm::FileHelper(f);
122 dataSize = fh->GetImageDataSize();
123 imageData = fh->GetImageData(); // somewhat important... can't remember
124 fh->SetWriteModeToRGB();
128 dataSize = fh->GetImageDataRawSize();
129 imageData = fh->GetImageDataRaw();// somewhat important... can't remember
130 fh->SetWriteModeToRaw();
133 if ( imageData == 0 ) // to avoid warning
135 std::cout << "Was unable to read pixels " << std::endl;
137 std::cout <<std::endl <<" dataSize " << dataSize << std::endl;
138 int nX,nY,nZ,sPP,planarConfig;
139 std::string pixelType, transferSyntaxName;
143 std::cout << " DIMX=" << nX << " DIMY=" << nY << " DIMZ=" << nZ << std::endl;
145 pixelType = f->GetPixelType();
146 sPP = f->GetSamplesPerPixel();
147 planarConfig = f->GetPlanarConfiguration();
149 std::cout << " pixelType=" << pixelType
150 << " SampleserPixel=" << sPP
151 << " PlanarConfiguration=" << planarConfig
152 << " PhotometricInterpretation="
153 << f->GetEntryValue(0x0028,0x0004)
156 int numberOfScalarComponents=f->GetNumberOfScalarComponents();
157 std::cout << "NumberOfScalarComponents " << numberOfScalarComponents
159 transferSyntaxName = f->GetTransferSyntaxName();
160 std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]"
167 // Writting an ACR file
168 // from a full gdcm readable File
170 std::cout << "WriteACR" << std::endl;
171 fh->WriteAcr(outputFileName);
174 case 'D' : // Not documented in the 'usage', because the method
175 case 'd' : // is known to be bugged.
176 // Writting a DICOM Implicit VR file
177 // from a full gdcm readable File
179 std::cout << "WriteDCM Implicit VR" << std::endl;
180 fh->WriteDcmImplVR(outputFileName);
185 // writting a DICOM Explicit VR
186 // from a full gdcm readable File
188 std::cout << "WriteDCM Explicit VR" << std::endl;
189 fh->WriteDcmExplVR(outputFileName);
194 // Writting a Raw File,
196 std::cout << "WriteRaw" << std::endl;
197 fh->WriteRawData(outputFileName);