X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmRLE.cxx;h=4387ebcfd78bb1460a1b434c9c886a81399184ab;hb=49484255e3674c038cbf2259a8c3325c0039cd41;hp=bb993d4828e74053fb20b2daf926987fa4ab7aca;hpb=33413fc6fb41f4fd7efbce32083f943f252e8b6f;p=gdcm.git diff --git a/src/gdcmRLE.cxx b/src/gdcmRLE.cxx index bb993d48..4387ebcf 100644 --- a/src/gdcmRLE.cxx +++ b/src/gdcmRLE.cxx @@ -1,35 +1,22 @@ - +// gdcmRLE.cxx +//----------------------------------------------------------------------------- #include #include "gdcmFile.h" #include /* to declare isprint() */ #define str2num(str, typeNum) *((typeNum *)(str)) -static int _gdcm_read_RLE_fragment (char ** image_buffer, - long lengthToDecode, - long uncompressedSegmentSize, - FILE* fp); -// static because nothing but gdcm_read_RLE_file may call it - -// ---------------------------------------------------------------------------- +//----------------------------------------------------------------------------- /** * \ingroup gdcmFile * \brief Reads a 'Run Length Encoded' Dicom encapsulated file + * @param fp already open File Pointer * @param image_buffer destination Address (in caller's memory space) * at which the pixel data should be copied * - * @return int acts as a Boolean + * @return Boolean */ - -// This is a debug version. -// Forget the printf as they will be removed -// as soon as the last Heuristics are checked - -// pb with RLE 16 Bits : - - -int -gdcmFile::gdcm_read_RLE_file (void * image_buffer) { +bool gdcmFile::gdcm_read_RLE_file (FILE *fp,void * image_buffer) { long fragmentBegining; // for ftell, fseek char * im = (char *)image_buffer; @@ -38,21 +25,21 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { guint32 nbRleSegments; guint32 RleSegmentOffsetTable[15]; guint16 ItemTagGr,ItemTagEl; - uncompressedSegmentSize=GetXSize()*GetYSize(); + uncompressedSegmentSize=Header->GetXSize()*Header->GetYSize(); ftellRes=ftell(fp); // Basic Offset Table with Item Value // Item Tag fread(&ItemTagGr,2,1,fp); // Reading (fffe):Basic Offset Table Item Tag Gr fread(&ItemTagEl,2,1,fp); // Reading (e000):Basic Offset Table Item Tag El - if(GetSwapCode()) { - ItemTagGr=SwapShort(ItemTagGr); - ItemTagEl=SwapShort(ItemTagEl); + if(Header->GetSwapCode()) { + ItemTagGr=Header->SwapShort(ItemTagGr); + ItemTagEl=Header->SwapShort(ItemTagEl); } // Item Length ftellRes=ftell(fp); fread(&ln,4,1,fp); - if(GetSwapCode()) - ln=SwapLong(ln); // Basic Offset Table Item Lentgh + if(Header->GetSwapCode()) + ln=Header->SwapLong(ln); // Basic Offset Table Item Lentgh if (ln != 0) { // What is it used for ?? char * BasicOffsetTableItemValue= (char *)malloc(ln+1); @@ -66,9 +53,9 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { ftellRes=ftell(fp); fread(&ItemTagGr,2,1,fp); // Reading (fffe) : Item Tag Gr fread(&ItemTagEl,2,1,fp); // Reading (e000) : Item Tag El - if(GetSwapCode()) { - ItemTagGr=SwapShort(ItemTagGr); - ItemTagEl=SwapShort(ItemTagEl); + if(Header->GetSwapCode()) { + ItemTagGr=Header->SwapShort(ItemTagGr); + ItemTagEl=Header->SwapShort(ItemTagEl); } // while 'Sequence Delimiter Item' (fffe,e0dd) not found @@ -76,20 +63,20 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { // Parse fragments of the current Fragment (Frame) ftellRes=ftell(fp); fread(&fragmentLength,4,1,fp); - if(GetSwapCode()) - fragmentLength=SwapLong(fragmentLength); // length + if(Header->GetSwapCode()) + fragmentLength=Header->SwapLong(fragmentLength); // length //------------------ scanning (not reading) fragment pixels fread(&nbRleSegments,4,1,fp); // Reading : Number of RLE Segments - if(GetSwapCode()) - nbRleSegments=SwapLong(nbRleSegments); + if(Header->GetSwapCode()) + nbRleSegments=Header->SwapLong(nbRleSegments); for(int k=1; k<=15; k++) { // Reading RLE Segments Offset Table ftellRes=ftell(fp); fread(&RleSegmentOffsetTable[k],4,1,fp); - if(GetSwapCode()) - RleSegmentOffsetTable[k]=SwapLong(RleSegmentOffsetTable[k]); + if(Header->GetSwapCode()) + RleSegmentOffsetTable[k]=Header->SwapLong(RleSegmentOffsetTable[k]); } if (nbRleSegments>1) { @@ -97,7 +84,7 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { RleSegmentLength[k]=RleSegmentOffsetTable[k+1]-RleSegmentOffsetTable[k]; ftellRes=ftell(fp); fragmentBegining=ftell(fp); - _gdcm_read_RLE_fragment (&im, RleSegmentLength[k],uncompressedSegmentSize,fp); + gdcm_read_RLE_fragment (&im, RleSegmentLength[k],uncompressedSegmentSize,fp); fseek(fp,fragmentBegining,SEEK_SET); fseek(fp,RleSegmentLength[k],SEEK_CUR); } @@ -105,7 +92,7 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { RleSegmentLength[nbRleSegments] = fragmentLength - RleSegmentOffsetTable[nbRleSegments]; ftellRes=ftell(fp); fragmentBegining=ftell(fp); - _gdcm_read_RLE_fragment (&im, RleSegmentLength[nbRleSegments],uncompressedSegmentSize, fp); + gdcm_read_RLE_fragment (&im, RleSegmentLength[nbRleSegments],uncompressedSegmentSize, fp); fseek(fp,fragmentBegining,SEEK_SET); fseek(fp,RleSegmentLength[nbRleSegments],SEEK_CUR); @@ -114,19 +101,19 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { ftellRes=ftell(fp); fread(&ItemTagGr,2,1,fp); // Reading (fffe) : Item Tag Gr fread(&ItemTagEl,2,1,fp); // Reading (e000) : Item Tag El - if(GetSwapCode()) { - ItemTagGr=SwapShort(ItemTagGr); - ItemTagEl=SwapShort(ItemTagEl); + if(Header->GetSwapCode()) { + ItemTagGr=Header->SwapShort(ItemTagGr); + ItemTagEl=Header->SwapShort(ItemTagEl); } } - if (GetBitsAllocated()==16) { // try to deal with RLE 16 Bits + if (Header->GetBitsAllocated()==16) { // try to deal with RLE 16 Bits im = (char *)image_buffer; // need to make 16 Bits Pixels from Low Byte and Hight Byte 'Planes' - int l = GetXSize()*GetYSize(); - int nbFrames = GetZSize(); + int l = Header->GetXSize()*Header->GetYSize(); + int nbFrames = Header->GetZSize(); char * newDest = (char*) malloc(l*nbFrames*2); char *x = newDest; @@ -143,24 +130,14 @@ gdcmFile::gdcm_read_RLE_file (void * image_buffer) { free(newDest); } - return (1); + return(true); } - /* -------------------------------------------------------------------- */ - // - // RLE LossLess Fragment - // - /* -------------------------------------------------------------------- */ - - // static because nothing but gdcm_read_RLE_file can call it - // DO NOT doxygen ! - -static int -_gdcm_read_RLE_fragment (char ** areaToRead, - long lengthToDecode, - long uncompressedSegmentSize, - FILE* fp) { +// ---------------------------------------------------------------------------- +// RLE LossLess Fragment +int gdcmFile::gdcm_read_RLE_fragment(char **areaToRead, long lengthToDecode, + long uncompressedSegmentSize, FILE *fp) { long ftellRes; int count; long numberOfOutputBytes=0; @@ -168,7 +145,6 @@ _gdcm_read_RLE_fragment (char ** areaToRead, ftellRes =ftell(fp); while(numberOfOutputBytes