/*========================================================================= Program: gdcm Module: $RCSfile: gdcmRLEFramesInfo.cxx,v $ Language: C++ Date: $Date: 2005/01/31 06:17:22 $ Version: $Revision: 1.10 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details. This software is distributed WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the above copyright notices for more information. =========================================================================*/ #include "gdcmRLEFramesInfo.h" #include "gdcmDebug.h" namespace gdcm { RLEFramesInfo::~RLEFramesInfo() { for(RLEFrameList::iterator it = Frames.begin(); it != Frames.end(); ++it) { delete (*it); } Frames.clear(); } /** * \brief Print self. * @param indent Indentation string to be prepended during printing. * @param os Stream to print to. */ void RLEFramesInfo::Print( std::ostream &os, std::string indent ) { os << std::endl; os << indent << "----------------- RLE frames --------------------------------" << std::endl; os << indent << "Total number of Frames : " << Frames.size() << std::endl; int frameNumber = 0; for(RLEFrameList::iterator it = Frames.begin(); it != Frames.end(); ++it) { os << indent << " frame number :" << frameNumber++ << std::endl; (*it)->Print( os, indent + " " ); } } void RLEFramesInfo::AddFrame(RLEFrame *frame) { Frames.push_back(frame); } RLEFrame *RLEFramesInfo::GetFirstFrame() { ItFrames = Frames.begin(); if (ItFrames != Frames.end()) return *ItFrames; return NULL; } RLEFrame *RLEFramesInfo::GetNextFrame() { gdcmAssertMacro (ItFrames != Frames.end()); ++ItFrames; if (ItFrames != Frames.end()) return *ItFrames; return NULL; } /** * \brief Reads from disk the Pixel Data of 'Run Length Encoded' * Dicom encapsulated file and decompress it. * @param fp already open File Pointer * at which the pixel data should be copied * @return Boolean */ bool RLEFramesInfo::DecompressRLEFile( std::ifstream *fp , uint8_t *raw, int xSize, int ySize, int zSize, int bitsAllocated ) { uint8_t *subRaw = raw; long rawSegmentSize = xSize * ySize; // Loop on the frame[s] for(RLEFrameList::iterator it = Frames.begin(); it != Frames.end(); ++it) { subRaw = (*it)->ReadAndDecompressRLEFrame( subRaw, rawSegmentSize, fp); } if ( bitsAllocated == 16 ) { // Try to deal with RLE 16 Bits ConvertRLE16BitsFromRLE8Bits( raw, xSize, ySize, zSize ); } return true; } /** * \brief Try to deal with RLE 16 Bits. * We assume the RLE has already been parsed and loaded in * Raw (through \ref ReadAndDecompressJPEGFile ). * We here need to make 16 Bits Pixels from Low Byte and * High Byte 'Planes'...(for what it may mean) * @return Boolean */ bool RLEFramesInfo::ConvertRLE16BitsFromRLE8Bits( uint8_t* raw, int xSize, int ySize,int numberOfFrames ) { size_t pixelNumber = xSize * ySize; size_t rawSize = xSize * ySize * numberOfFrames; // We assumed Raw contains the decoded RLE pixels but as // 8 bits per pixel. In order to convert those pixels to 16 bits // per pixel we cannot work in place within Raw and hence // we copy it in a safe place, say copyRaw. uint8_t* copyRaw = new uint8_t[rawSize * 2]; memmove( copyRaw, raw, rawSize * 2 ); uint8_t* x = raw; uint8_t* a = copyRaw; uint8_t* b = a + pixelNumber; for ( int i = 0; i < numberOfFrames; i++ ) { for ( unsigned int j = 0; j < pixelNumber; j++ ) { *(x++) = *(b++); *(x++) = *(a++); } } delete[] copyRaw; /// \todo check that operator new []didn't fail, and sometimes return false return true; } } // end namespace gdcm