- * \ingroup gdcmHeader
- * \brief
- * @return Lookup Table RGB
- * \ when (0028,0004),Photometric Interpretation = [PALETTE COLOR ]
- * \ and (0028,1201),(0028,1202),(0028,1202) are found
- * \warning : hazardous ! Use better GetPubElValVoidAreaByNumber
- */
-void * gdcmHeader::GetLUTRGB(void) {
-// Not so easy : see
-// http://www.barre.nom.fr/medical/dicom2/limitations.html#Color%20Lookup%20Tables
-// and OT-PAL-8-face.dcm
-
- if (GetPubElValByNumber(0x0028,0x0004) == GDCM_UNFOUND) {
- dbg.Verbose(0, "gdcmHeader::GetLUTRGB: unfound Photometric Interpretation");
- return NULL;
- }
- void * LutR,*LutG,*LutB;
- int l;
-
- // Maybe, some day we get an image
- // that respects the definition ...
- // Let's consider no ones does.
-
- l= GetLUTLength();
- if(l==0)
- return (NULL);
- int nBits=GetLUTNbits();
-
- // a virer quand on aura trouve UNE image
- // qui correspond VRAIMENT à la definition !
- std::cout << "l " << l << " nBits " << nBits;
-
- l= l/(nBits/8);
-
- LutR =GetPubElValVoidAreaByNumber(0x0028,0x1201);
- LutG =GetPubElValVoidAreaByNumber(0x0028,0x1202);
- LutB =GetPubElValVoidAreaByNumber(0x0028,0x1203);
-
- // Warning : Any value for nBits as to be considered as 8
- // Any value for Length as to be considered as 256
- // That's DICOM ...
-
- // Just wait before removing the following code
- /*
- if (nBits == 16) {
- guint16 * LUTRGB, *rgb;
- LUTRGB = rgb = (guint16 *) malloc(3*l*sizeof( guint16));
- guint16 * r = (guint16 *)LutR;
- guint16 * g = (guint16 *)LutG;
- guint16 * b = (guint16 *)LutB;
- for(int i=0;i<l;i++) {
- *rgb++ = *r++;
- *rgb++ = *g++;
- *rgb++ = *b++;
- }
- return(LUTRGB);
- } else
-
- */ { // we assume it's always 8 Bits
- l=256; // we assume ...
- unsigned char * LUTRGB, *rgb;
- LUTRGB = rgb = (unsigned char *) malloc(3*l*sizeof( char));
- unsigned char * r = (unsigned char *)LutR;
- unsigned char * g = (unsigned char *)LutG;
- unsigned char * b = (unsigned char *)LutB;
- for(int i=0;i<l;i++) {
- //std::cout << "lut16 " << i << " : " << *r << " " << *g << " " << *b
- // << std::endl;
- printf("lut 8 %d : %d %d %d \n",i,*r,*g,*b);
- *rgb++ = *r++;
- *rgb++ = *g++;
- *rgb++ = *b++;
- }
- free(LutR); free(LutB); free(LutG);
- return(LUTRGB);
- }
+ * \brief Initialize a default DICOM header that should contain all the
+ * field require by other reader. DICOM standart does not
+ * explicitely defines thoses fields, heuristic has been choosen.
+ * This is not perfect as we are writting a CT image...
+ */
+void Header::InitializeDefaultHeader()
+{
+ typedef struct
+ {
+ const char *value;
+ uint16_t group;
+ uint16_t elem;
+ } DICOM_DEFAULT_VALUE;
+
+ std::string date = Util::GetCurrentDate();
+ std::string time = Util::GetCurrentTime();
+ std::string uid = Util::CreateUniqueUID();
+
+ static DICOM_DEFAULT_VALUE defaultvalue[] = {
+ { "76", 0x0002, 0x0000}, // MetaElementGroup Length // FIXME: how to recompute ?
+ { "1.2.840.10008.5.1.4.1.1.2", 0x0002, 0x0002}, // MediaStorageSOPInstanceUID (CT Image Storage)
+ { uid.c_str(), 0x0002, 0x0012}, // META Implementation Class UID
+ { "ISO_IR 100",0x0008, 0x0005}, // Specific Character Set
+ { "DERIVED\\SECONDARY\\OTHER\\GDCM", 0x0008, 0x0008}, // Image Type
+ { "1.2.840.10008.5.1.4.1.1.2", 0x0008, 0x0016}, // SOP Class UID
+ { "", 0x0008, 0x0050}, // Accession Number
+ { "GDCM", 0x0008, 0x0070}, // Manufacturer
+ { "GDCM", 0x0008, 0x0080}, // Institution Name
+ { "http://www-creatis.insa-lyon.fr/Public/Gdcm/", 0x0008, 0x0081}, // Institution Address
+ { "", 0x0008, 0x0090}, // Refering Physician Name
+ { "", 0x0008, 0x1030}, // Study Description
+ { "", 0x0008, 0x1050}, // Performing Physician's Name
+ { "", 0x0008, 0x1060}, // Name of Physician(s) Reading Study
+ { "", 0x0010, 0x0040}, // Patient's Sex
+ { uid.c_str(), 0x0020, 0x000d}, // StudyInstanceUID
+ { uid.c_str(), 0x0020, 0x000e}, // SeriesInstanceUID
+ { "", 0x0020, 0x0011}, // AcquisitionNumber
+ { "1\\0\\0\\0\\1\\0", 0x0020, 0x0037}, // Image Orientation Patient
+ { "1", 0x0028, 0x0002}, // Samples per pixel 1 or 3
+ { "MONOCHROME2",0x0028, 0x0004}, // photochromatic interpretation
+
+// Date and timeG
+
+ { date.c_str(), 0x0008, 0x0012 } , // Instance Creation Date
+ { time.c_str(), 0x0008, 0x0013 } , // Instance Creation Time
+ { date.c_str(), 0x0008, 0x0020 } , // Study Date
+ { date.c_str(), 0x0008, 0x0022 } , // Acquisition Date
+ { date.c_str(), 0x0008, 0x0023 } , // Content Date
+ { time.c_str(), 0x0008, 0x0030 } , // Study Time
+ { "CT", 0x0008, 0x0060 } , // Modality
+ { "GDCM", 0x0010, 0x0010 } , // Patient's Name
+ { "GDCMID", 0x0010, 0x0020 } , // Patient ID
+ { "1", 0x0020, 0x0010 } , // StudyID
+ { "1", 0x0020, 0x0011 } , // SeriesNumber
+ { "1.0", 0x0018, 0x0050 } , // slice Thickness
+ { "1.0", 0x0018, 0x0088 } , // space between slices
+ { "1.0\\1.0", 0x0028, 0x0030 } , // pixelSpacing
+ { "64", 0x0028, 0x0010 } , // nbRows
+ { "64", 0x0028, 0x0011 } , // nbCols
+ { "16", 0x0028, 0x0100 } , // BitsAllocated 8 or 16
+ { "12", 0x0028, 0x0101 } , // BitsStored 8 or 12
+ { "11", 0x0028, 0x0102 } , // HighBit 7 or 11
+ { "0", 0x0028, 0x0103 } , // Pixel Representation 0(unsigned) or 1(signed)
+ { "1000.0", 0x0028, 0x1051 } , // Window Width
+ { "500.0", 0x0028, 0x1050 } , // Window Center
+ { 0, 0, 0 }
+ };
+
+ // default value
+ // Special case this is the image (not a string)
+ GrPixel = 0x7fe0;
+ NumPixel = 0x0010;
+ ReplaceOrCreateByNumber(0, 0, GrPixel, NumPixel);
+
+ // All remaining strings:
+ unsigned int i = 0;
+ DICOM_DEFAULT_VALUE current = defaultvalue[i];
+ while( current.value )
+ {
+ std::string value = Util::DicomString( current.value ); //pad the string properly
+ ReplaceOrCreateByNumber(value, current.group, current.elem);
+ current = defaultvalue[++i];
+ }