- free(newDest);
-
- // now, it's an RGB image
- // Lets's write it in the Header
-
- // CreateOrReplaceIfExist ?
-
- std::string spp = "3"; // Samples Per Pixel
- Header->SetPubElValByNumber(spp,0x0028,0x0002);
- std::string rgb= "RGB "; // Photometric Interpretation
- Header->SetPubElValByNumber(rgb,0x0028,0x0004);
- std::string planConfig = "0"; // Planar Configuration
- Header->SetPubElValByNumber(planConfig,0x0028,0x0006);
-
- } else {
- // need to make RGB Pixels (?)
- // from grey Pixels (?!)
- // and Gray Lut (!?!)
- // or Segmented xxx Palette Color Lookup Table Data and so on
-
- // Well . I'll wait till I find such an image
-
- // Oops! I get one (gdcm-US-ALOKA-16.dcm)
- // No idea how to manage it
- // It seems that *no Dicom Viewer* has any idea :-(
- // Segmented xxx Palette Color are *more* than 65535 long ?!?
-
- std::string rgb= "MONOCHROME1 "; // Photometric Interpretation
- Header->SetPubElValByNumber(rgb,0x0028,0x0004);
-
- }
-
- // TODO : Drop Palette Color out of the Header?
- return lgrTotale;
-}
-
-/**
- * \ingroup gdcmFile
- * \brief Allocates necessary memory, copies the pixel data
- * (image[s]/volume[s]) to newly allocated zone.
- * Transforms YBR pixels into RGB pixels if any
- Transforms 3 planes R, G, B into a single RGB Plane
- DOES NOT transform Grey plane + 3 Palettes into a RGB Plane
- * @return Pointer to newly allocated pixel data.
- * \ NULL if alloc fails
- */
-void * gdcmFile::GetImageDataRaw (void) {
- if (Header->HasLUT())
- lgrTotale /= 3; // TODO Let gdcmHeadar user a chance
- // to get the right value
- // Create a member lgrTotaleRaw ???
- PixelData = (void *) malloc(lgrTotale);
- if (PixelData)
- GetImageDataIntoVectorRaw(PixelData, lgrTotale);
- return(PixelData);
-}
-
-/**
- * \ingroup gdcmFile
- * \brief Copies at most MaxSize bytes of pixel data to caller's
- * memory space.
- * \warning This function was designed to avoid people that want to build
- * a volume from an image stack to need first to get the image pixels
- * and then move them to the volume area.
- * It's absolutely useless for any VTK user since vtk chooses
- * to invert the lines of an image, that is the last line comes first
- * (for some axis related reasons?). Hence he will have
- * to load the image line by line, starting from the end.
- * VTK users hace to call GetImageData
- * \warning DOES NOT transform the Grey Plane + Palette Color (if any)
- * into a single RGB Pixels Plane
- * the (VTK) user will manage the palettes
- *
- * @param destination Address (in caller's memory space) at which the
- * pixel data should be copied
- * @param MaxSize Maximum number of bytes to be copied. When MaxSize
- * is not sufficient to hold the pixel data the copy is not
- * executed (i.e. no partial copy).
- * @return On success, the number of bytes actually copied. Zero on
- * failure e.g. MaxSize is lower than necessary.
- */
-size_t gdcmFile::GetImageDataIntoVectorRaw (void* destination, size_t MaxSize) {
-
- int nb, nbu, highBit, signe;
- std::string str_nbFrames, str_nb, str_nbu, str_highBit, str_signe;
-
- if ( lgrTotale > MaxSize ) {
- dbg.Verbose(0, "gdcmFile::GetImageDataIntoVector: pixel data bigger"
- "than caller's expected MaxSize");
- return (size_t)0;
- }
-
- (void)ReadPixelData(destination);
-
- // Nombre de Bits Alloues pour le stockage d'un Pixel
- str_nb = Header->GetPubElValByNumber(0x0028,0x0100);
- if (str_nb == GDCM_UNFOUND ) {
- nb = 16;
- } else {
- nb = atoi(str_nb.c_str() );
- }
-
- // Nombre de Bits Utilises
- str_nbu=Header->GetPubElValByNumber(0x0028,0x0101);
- if (str_nbu == GDCM_UNFOUND ) {
- nbu = nb;
- } else {
- nbu = atoi(str_nbu.c_str() );
- }
-
- // Position du Bit de Poids Fort
- str_highBit=Header->GetPubElValByNumber(0x0028,0x0102);
- if (str_highBit == GDCM_UNFOUND ) {
- highBit = nb - 1;
- } else {
- highBit = atoi(str_highBit.c_str() );
- }
- // Pixel sign
- // 0 = Unsigned
- // 1 = Signed
- str_signe=Header->GetPubElValByNumber(0x0028,0x0103);
- if (str_signe == GDCM_UNFOUND ) {
- signe = 0; // default is unsigned
- } else {
- signe = atoi(str_signe.c_str() );
- }
-
- // re arange bytes inside the integer
- if (nb != 8)
- SwapZone(destination, Header->GetSwapCode(), lgrTotale, nb);
-
- // to avoid pb with some xmedcon breakers images
- if (nb==16 && nbu<nb && signe==0) {
- int l = (int)lgrTotale / (nb/8);
- guint16 *deb = (guint16 *)destination;
- for(int i = 0; i<l; i++) {
- if(*deb == 0xffff)
- *deb=0;
- deb++;
- }
- }