Program: gdcm
Module: $RCSfile: TestCopyRescaleDicom.cxx,v $
Language: C++
- Date: $Date: 2004/12/07 16:55:56 $
- Version: $Revision: 1.1 $
+ Date: $Date: 2005/01/21 11:40:54 $
+ Version: $Revision: 1.8 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
-#include "gdcmHeader.h"
#include "gdcmFile.h"
+#include "gdcmFileHelper.h"
#include "gdcmValEntry.h"
#include "gdcmBinEntry.h"
//////////////// Step 1:
std::cout << " 1...";
- gdcm::Header *originalH = new gdcm::Header( filename );
- gdcm::Header *copyH = new gdcm::Header( );
+ gdcm::File *originalH = new gdcm::File( filename );
+ gdcm::File *copyH = new gdcm::File( );
//First of all copy the header field by field
// Warning :Accessor gdcmElementSet::GetEntry() should not exist
// It was commented out by Mathieu, that was a *good* idea
- // (the user does NOT have to know the way we implemented the Header !)
+ // (the user does NOT have to know the way we implemented the File !)
// Waiting for a 'clean' solution, I keep the method ...JPRx
//////////////// Step 2:
std::cout << "2...";
- originalH->Initialize();
- gdcm::DocEntry* d = originalH->GetNextEntry();
-
+ // Copy of the header content
+ gdcm::DocEntry* d = originalH->GetFirstEntry();
while(d)
{
- if ( gdcm::ValEntry* v = dynamic_cast<gdcm::ValEntry*>(d) )
+ if ( gdcm::BinEntry* b = dynamic_cast<gdcm::BinEntry*>(d) )
+ {
+ copyH->ReplaceOrCreate(
+ b->GetBinArea(),
+ b->GetLength(),
+ b->GetGroup(),
+ b->GetElement(),
+ b->GetVR() );
+ }
+ else if ( gdcm::ValEntry* v = dynamic_cast<gdcm::ValEntry*>(d) )
{
- copyH->ReplaceOrCreateByNumber(
+ copyH->ReplaceOrCreate(
v->GetValue(),
v->GetGroup(),
v->GetElement(),
v->GetVR() );
}
+ else
+ {
+ // We skip pb of SQ recursive exploration
+ }
d=originalH->GetNextEntry();
}
- gdcm::File *original = new gdcm::File( originalH );
- gdcm::File *copy = new gdcm::File( copyH );
+ gdcm::FileHelper *original = new gdcm::FileHelper( originalH );
+ gdcm::FileHelper *copy = new gdcm::FileHelper( copyH );
size_t dataSize = original->GetImageDataSize();
- uint8_t* imageData = original->GetImageData();
- size_t rescaleSize = dataSize / 2;
-#ifdef VERSION_1
- uint8_t *rescaleImage = new uint8_t[dataSize];
-#endif
-
- //2 case now, if 16 bits input then try to downscale to 8bits just for
- //fun:
-
- const std::string & bitsStored = originalH->GetEntryByNumber(0x0028,0x0101);
-// std::cout << "BitsStored: " << bitsStored << std::endl;
+ size_t rescaleSize;
+ uint8_t *rescaleImage;
+ const std::string & bitsStored = originalH->GetEntry(0x0028,0x0101);
if( bitsStored == "16" )
{
- uint16_t* imageData16 = (uint16_t*)original->GetImageData();
- copyH->ReplaceOrCreateByNumber( "8", 0x0028, 0x0100); // BitsAllocated
- copyH->ReplaceOrCreateByNumber( "8", 0x0028, 0x0101); // BitsStored
- copyH->ReplaceOrCreateByNumber( "7", 0x0028, 0x0102); // HighBit
- copyH->ReplaceOrCreateByNumber( "0", 0x0028, 0x0103); //Pixel Representation
+ std::cout << "Rescale...";
+ copyH->ReplaceOrCreate( "8", 0x0028, 0x0100); // BitsAllocated
+ copyH->ReplaceOrCreate( "8", 0x0028, 0x0101); // BitsStored
+ copyH->ReplaceOrCreate( "7", 0x0028, 0x0102); // HighBit
+ copyH->ReplaceOrCreate( "0", 0x0028, 0x0103); //Pixel Representation
// We assume the value were from 0 to uint16_t max
+ rescaleSize = dataSize / 2;
+ rescaleImage = new uint8_t[dataSize];
+
+ uint16_t* imageData16 = (uint16_t*)original->GetImageData();
for(unsigned int i=0; i<rescaleSize; i++)
{
-#ifdef VERSION_1
rescaleImage[i] = imageData16[i]/256;
-#else
- imageData[i] = imageData16[i]/256;
-#endif
}
-#ifdef VERSION_1
- copy->SetImageData(rescaleImage, rescaleSize);
-#endif
}
else
{
- copy->SetImageData(imageData, dataSize);
+ std::cout << "Same...";
+ rescaleSize = dataSize;
+ rescaleImage = new uint8_t[dataSize];
+ memcpy(rescaleImage,original->GetImageData(),dataSize);
}
+ copy->SetImageData(rescaleImage, rescaleSize);
+
//////////////// Step 3:
std::cout << "3...";
copy->SetWriteModeToRGB();
delete copy;
delete originalH;
delete copyH;
+ delete[] rescaleImage;
return 1;
}
//////////////// Step 4:
std::cout << "4...";
- copy = new gdcm::File( output );
+ copy = new gdcm::FileHelper( output );
//Is the file written still gdcm parsable ?
- if ( !copy->GetHeader()->IsReadable() )
+ if ( !copy->GetFile()->IsReadable() )
{
std::cout << " Failed" << std::endl
<< " " << output << " not readable" << std::endl;
delete original;
delete originalH;
+ delete[] rescaleImage;
return 1;
}
size_t dataSizeWritten = copy->GetImageDataSize();
uint8_t* imageDataWritten = copy->GetImageData();
- if (dataSize != dataSizeWritten)
+ if (originalH->GetXSize() != copy->GetFile()->GetXSize() ||
+ originalH->GetYSize() != copy->GetFile()->GetYSize() ||
+ originalH->GetZSize() != copy->GetFile()->GetZSize())
+ {
+ std::cout << "Failed" << std::endl
+ << " X Size differs: "
+ << "X: " << originalH->GetXSize() << " # "
+ << copy->GetFile()->GetXSize() << " | "
+ << "Y: " << originalH->GetYSize() << " # "
+ << copy->GetFile()->GetYSize() << " | "
+ << "Z: " << originalH->GetZSize() << " # "
+ << copy->GetFile()->GetZSize() << std::endl;
+ delete original;
+ delete copy;
+ delete originalH;
+ delete[] rescaleImage;
+
+ return 1;
+ }
+
+ if (rescaleSize != dataSizeWritten)
{
std::cout << " Failed" << std::endl
<< " Pixel areas lengths differ: "
delete original;
delete copy;
delete originalH;
+ delete[] rescaleImage;
return 1;
}
- if (int res = memcmp(imageData, imageDataWritten, dataSize) !=0)
+ if (int res = memcmp(rescaleImage, imageDataWritten, rescaleSize) !=0)
{
(void)res;
std::cout << " Failed" << std::endl
delete original;
delete copy;
delete originalH;
+ delete[] rescaleImage;
return 1;
}
delete original;
delete copy;
delete originalH;
-#ifdef VERSION_1
delete[] rescaleImage;
-#endif
return 0;
}
<< std::endl;
std::cout << " apply the following to each filename.xxx: "
<< std::endl;
- std::cout << " step 1: parse the image (as gdcmHeader) and call"
+ std::cout << " step 1: parse the image (as gdcmFile) and call"
<< " IsReadable(). After that, call GetImageData() and "
<< "GetImageDataSize() "
<< std::endl;