Program: gdcm
Module: $RCSfile: gdcmFile.cxx,v $
Language: C++
- Date: $Date: 2004/10/12 04:35:46 $
- Version: $Revision: 1.139 $
+ Date: $Date: 2004/10/18 12:49:22 $
+ Version: $Revision: 1.145 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
#include "gdcmFile.h"
#include "gdcmDebug.h"
-#include "gdcmPixelConvert.h"
namespace gdcm
{
{
ImageDataSize = ImageDataSizeRaw;
}
+
+ PixelConverter.GrabInformationsFromHeader( HeaderInternal );
}
SaveInitialValues();
}
// 0028|1202 [US] [Green Palette Color Lookup Table Data]
// 0028|1203 [US] [Blue Palette Color Lookup Table Data]
- // Number of "Bits Allocated"
int numberBitsAllocated = HeaderInternal->GetBitsAllocated();
+ // Number of "Bits Allocated" is fixed to 16 when:
+ // - it is not defined (i.e. it's value is 0)
+ // - it's 12, since we will expand the image to 16 bits (see
+ // PixelConvert::ConvertDecompress12BitsTo16Bits() )
if ( ( numberBitsAllocated == 0 ) || ( numberBitsAllocated == 12 ) )
{
numberBitsAllocated = 16;
return ImageDataSize;
}
- // from Lut R + Lut G + Lut B
- uint8_t *newDest = new uint8_t[ImageDataSize];
- uint8_t *a = (uint8_t *)destination;
- uint8_t *lutRGBA = HeaderInternal->GetLUTRGBA();
-
- if ( lutRGBA )
+ FILE* fp = HeaderInternal->OpenFile();
+ if ( PixelConverter.BuildRGBImage( fp ) )
{
- int j;
- // move Gray pixels to temp area
- memmove(newDest, destination, ImageDataSizeRaw);
- for (size_t i=0; i<ImageDataSizeRaw; ++i)
- {
- // Build RGB Pixels
- j = newDest[i]*4;
- *a++ = lutRGBA[j];
- *a++ = lutRGBA[j+1];
- *a++ = lutRGBA[j+2];
- }
- delete[] newDest;
+ memmove( destination,
+ (void*)PixelConverter.GetRGB(),
+ PixelConverter.GetRGBSize() );
- // now, it's an RGB image
- // Lets's write it in the Header
+ // now, it's an RGB image
+ // Lets's write it in the Header
- // FIXME : Better use CreateOrReplaceIfExist ?
+ // FIXME : Better use CreateOrReplaceIfExist ?
- std::string spp = "3"; // Samples Per Pixel
- HeaderInternal->SetEntryByNumber(spp,0x0028,0x0002);
- std::string rgb = "RGB "; // Photometric Interpretation
- HeaderInternal->SetEntryByNumber(rgb,0x0028,0x0004);
- std::string planConfig = "0"; // Planar Configuration
- HeaderInternal->SetEntryByNumber(planConfig,0x0028,0x0006);
+ std::string spp = "3"; // Samples Per Pixel
+ HeaderInternal->SetEntryByNumber(spp,0x0028,0x0002);
+ std::string rgb = "RGB "; // Photometric Interpretation
+ HeaderInternal->SetEntryByNumber(rgb,0x0028,0x0004);
+ std::string planConfig = "0"; // Planar Configuration
+ HeaderInternal->SetEntryByNumber(planConfig,0x0028,0x0006);
}
- else // GetLUTRGBA() failed
+ else
{
+ // PixelConverter.BuildRGBImage( fp ) failed probably because
+ // PixelConverter.GetLUTRGBA() failed:
// (gdcm-US-ALOKA-16.dcm), contains Segmented xxx Palette Color
// that are *more* than 65535 long ?!?
// No idea how to manage such an image !
std::string photomInterp = "MONOCHROME1 "; // Photometric Interpretation
HeaderInternal->SetEntryByNumber(photomInterp,0x0028,0x0004);
}
+ HeaderInternal->CloseFile();
/// \todo Drop Palette Color out of the Header?
return ImageDataSize;
* @return On success, the number of bytes actually copied. Zero on
* failure e.g. MaxSize is lower than necessary.
*/
-size_t File::GetImageDataIntoVectorRaw (void* destination, size_t maxSize)
+void File::GetImageDataIntoVectorRaw (void* destination, size_t maxSize)
{
// we save the initial values of the following
// in order to be able to restore the header in a disk-consistent state
{
dbg.Verbose(0, "File::GetImageDataIntoVector: pixel data bigger"
"than caller's expected MaxSize");
- return (size_t)0;
+ return;
}
- ReadPixelData( destination );
-
- // Number of Bits Allocated for storing a Pixel
- int numberBitsAllocated = HeaderInternal->GetBitsAllocated();
- if ( numberBitsAllocated == 0 )
- {
- numberBitsAllocated = 16;
- }
-
- // Number of Bits actually used
- int numberBitsStored = HeaderInternal->GetBitsStored();
- if ( numberBitsStored == 0 )
- {
- numberBitsStored = numberBitsAllocated;
- }
-
- // High Bit Position
- int highBitPosition = HeaderInternal->GetHighBitPosition();
- if ( highBitPosition == 0 )
- {
- highBitPosition = numberBitsAllocated - 1;
- }
-
- bool signedPixel = HeaderInternal->IsSignedPixelData();
-
- PixelConvert::ConvertReorderEndianity(
- (uint8_t*) destination,
- ImageDataSize,
- numberBitsStored,
- numberBitsAllocated,
- HeaderInternal->GetSwapCode(),
- signedPixel );
-
- PixelConvert::ConvertReArrangeBits(
- (uint8_t*) destination,
- ImageDataSize,
- numberBitsStored,
- numberBitsAllocated,
- highBitPosition );
-
-#ifdef GDCM_DEBUG
- FILE* DebugFile;
- DebugFile = fopen( "SpuriousFile.RAW", "wb" );
- fwrite( PixelConvertor.GetUncompressed(),
- PixelConvertor.GetUncompressedsSize(),
- 1, DebugFile );
- fclose( DebugFile );
-#endif //GDCM_DEBUG
-
-// SPLIT ME
-//////////////////////////////////
-// Deal with the color
-
- // Monochrome pictures don't require color intervention
- if ( HeaderInternal->IsMonochrome() )
- {
- return ImageDataSize;
- }
-
- // Planar configuration = 0 : Pixels are already RGB
- // Planar configuration = 1 : 3 planes : R, G, B
- // Planar configuration = 2 : 1 gray Plane + 3 LUT
-
- // Well ... supposed to be !
- // See US-PAL-8-10x-echo.dcm: PlanarConfiguration=0,
- // PhotometricInterpretation=PALETTE COLOR
- // and heuristic has to be found :-(
-
- int planConf = HeaderInternal->GetPlanarConfiguration();
-
- // Planar configuration = 2 ==> 1 gray Plane + 3 LUT
- // ...and...
- // whatever the Planar Configuration might be, "PALETTE COLOR "
- // implies that we deal with the palette.
- if ( ( planConf == 2 ) || HeaderInternal->IsPaletteColor() )
- {
- return ImageDataSize;
- }
-
- // When planConf is 0, pixels are allready in RGB
+ FILE* fp = HeaderInternal->OpenFile();
+ PixelConverter.ReadAndDecompressPixelData( fp );
+ HeaderInternal->CloseFile();
+ memmove( destination,
+ (void*)PixelConverter.GetDecompressed(),
+ PixelConverter.GetDecompressedSize() );
- if ( planConf == 1 )
+ if ( ! PixelConverter.IsDecompressedRGB() )
{
- uint8_t* newDest = new uint8_t[ImageDataSize];
- // Warning : YBR_FULL_422 acts as RGB
- if ( HeaderInternal->IsYBRFull() )
- {
- ConvertYcBcRPlanesToRGBPixels((uint8_t*)destination, newDest);
- }
- else
- {
- ConvertRGBPlanesToRGBPixels((uint8_t*)destination, newDest);
- }
- memmove(destination, newDest, ImageDataSize);
- delete[] newDest;
+ return;
}
///////////////////////////////////////////////////
HeaderInternal->SetEntryByNumber(photInt,0x0028,0x0004);
HeaderInternal->SetEntryByNumber(planConfig,0x0028,0x0006);
- return ImageDataSize;
-}
-
-/**
- * \brief Convert (Y plane, cB plane, cR plane) to RGB pixels
- * \warning Works on all the frames at a time
- */
-void File::ConvertYcBcRPlanesToRGBPixels(uint8_t* source,
- uint8_t* destination)
-{
- // to see the tricks about YBR_FULL, YBR_FULL_422,
- // YBR_PARTIAL_422, YBR_ICT, YBR_RCT have a look at :
- // ftp://medical.nema.org/medical/dicom/final/sup61_ft.pdf
- // and be *very* affraid
- //
- int l = HeaderInternal->GetXSize() * HeaderInternal->GetYSize();
- int nbFrames = HeaderInternal->GetZSize();
-
- uint8_t* a = source;
- uint8_t* b = source + l;
- uint8_t* c = source + l + l;
- double R, G, B;
-
- /// \todo : Replace by the 'well known' integer computation
- /// counterpart. Refer to
- /// http://lestourtereaux.free.fr/papers/data/yuvrgb.pdf
- /// for code optimisation.
-
- for (int i = 0; i < nbFrames; i++)
- {
- for (int j = 0; j < l; j++)
- {
- R = 1.164 *(*a-16) + 1.596 *(*c -128) + 0.5;
- G = 1.164 *(*a-16) - 0.813 *(*c -128) - 0.392 *(*b -128) + 0.5;
- B = 1.164 *(*a-16) + 2.017 *(*b -128) + 0.5;
-
- if (R < 0.0) R = 0.0;
- if (G < 0.0) G = 0.0;
- if (B < 0.0) B = 0.0;
- if (R > 255.0) R = 255.0;
- if (G > 255.0) G = 255.0;
- if (B > 255.0) B = 255.0;
-
- *(destination++) = (uint8_t)R;
- *(destination++) = (uint8_t)G;
- *(destination++) = (uint8_t)B;
- a++;
- b++;
- c++;
- }
- }
-}
-
-/**
- * \brief Convert (Red plane, Green plane, Blue plane) to RGB pixels
- * \warning Works on all the frames at a time
- */
-void File::ConvertRGBPlanesToRGBPixels(uint8_t* source,
- uint8_t* destination)
-{
- int l = HeaderInternal->GetXSize() * HeaderInternal->GetYSize() * HeaderInternal->GetZSize();
-
- uint8_t* a = source;
- uint8_t* b = source + l;
- uint8_t* c = source + l + l;
-
- for (int j = 0; j < l; j++)
- {
- *(destination++) = *(a++);
- *(destination++) = *(b++);
- *(destination++) = *(c++);
- }
+ return;
}
/**
return true;
}
-//-----------------------------------------------------------------------------
-// Private
/**
- * \brief Read pixel data from disk (optionaly decompressing) into the
- * caller specified memory location.
- * @param destination where the pixel data should be stored.
- *
+ * \brief Access to the underlying \ref PixelConvertver RGBA LUT
*/
-bool File::ReadPixelData(void* destination)
+uint8_t* File::GetLutRGBA()
{
- FILE* fp = HeaderInternal->OpenFile();
-
- if ( !fp )
- {
- return false;
- }
- if ( fseek(fp, HeaderInternal->GetPixelOffset(), SEEK_SET) == -1 )
- {
- HeaderInternal->CloseFile();
- return false;
- }
-
- if ( HeaderInternal->GetBitsAllocated() == 12 )
- {
- PixelConvert::ConvertDecompress12BitsTo16Bits(
- (uint8_t*)destination,
- HeaderInternal->GetXSize(),
- HeaderInternal->GetYSize(),
- fp);
- HeaderInternal->CloseFile();
- return true;
- }
-
- // ---------------------- Uncompressed File
- if ( !HeaderInternal->IsDicomV3() ||
- HeaderInternal->IsImplicitVRLittleEndianTransferSyntax() ||
- HeaderInternal->IsExplicitVRLittleEndianTransferSyntax() ||
- HeaderInternal->IsExplicitVRBigEndianTransferSyntax() ||
- HeaderInternal->IsDeflatedExplicitVRLittleEndianTransferSyntax() )
- {
- size_t ItemRead = fread(destination, HeaderInternal->GetPixelAreaLength(), 1, fp);
- HeaderInternal->CloseFile();
- if ( ItemRead != 1 )
- {
- return false;
- }
- else
- {
- return true;
- }
- }
-
- // ---------------------- Run Length Encoding
- if ( HeaderInternal->IsRLELossLessTransferSyntax() )
- {
- bool res = PixelConvert::ReadAndDecompressRLEFile(
- destination,
- HeaderInternal->GetXSize(),
- HeaderInternal->GetYSize(),
- HeaderInternal->GetZSize(),
- HeaderInternal->GetBitsAllocated(),
- &(HeaderInternal->RLEInfo),
- fp );
- HeaderInternal->CloseFile();
- return res;
- }
-
- // --------------- SingleFrame/Multiframe JPEG Lossless/Lossy/2000
- int numberBitsAllocated = HeaderInternal->GetBitsAllocated();
- if ( ( numberBitsAllocated == 0 ) || ( numberBitsAllocated == 12 ) )
- {
- numberBitsAllocated = 16;
- }
-
- bool res = PixelConvert::ReadAndDecompressJPEGFile(
- (uint8_t*)destination,
- HeaderInternal->GetXSize(),
- HeaderInternal->GetYSize(),
- HeaderInternal->GetBitsAllocated(),
- HeaderInternal->GetBitsStored(),
- HeaderInternal->GetSamplesPerPixel(),
- HeaderInternal->GetPixelSize(),
- HeaderInternal->IsJPEG2000(),
- HeaderInternal->IsJPEGLossless(),
- &(HeaderInternal->JPEGInfo),
- fp );
- HeaderInternal->CloseFile();
- return res;
+ return PixelConverter.GetLutRGBA();
}
} // end namespace gdcm