00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030 #ifndef DIYBRPXT_H
00031 #define DIYBRPXT_H
00032
00033 #include "dcmtk/config/osconfig.h"
00034
00035 #include "dcmtk/dcmimage/dicopxt.h"
00036 #include "dcmtk/dcmimgle/diinpx.h"
00037
00038
00039
00040
00041
00042
00045 template<class T1, class T2>
00046 class DiYBRPixelTemplate
00047 : public DiColorPixelTemplate<T2>
00048 {
00049
00050 public:
00051
00061 DiYBRPixelTemplate(const DiDocument *docu,
00062 const DiInputPixel *pixel,
00063 EI_Status &status,
00064 const unsigned long planeSize,
00065 const int bits,
00066 const OFBool rgb)
00067 : DiColorPixelTemplate<T2>(docu, pixel, 3, status)
00068 {
00069 if ((pixel != NULL) && (this->Count > 0) && (status == EIS_Normal))
00070 convert(OFstatic_cast(const T1 *, pixel->getData()) + pixel->getPixelStart(), planeSize, bits, rgb);
00071 }
00072
00075 virtual ~DiYBRPixelTemplate()
00076 {
00077 }
00078
00079
00080 private:
00081
00089 void convert(const T1 *pixel,
00090 const unsigned long planeSize,
00091 const int bits,
00092 const OFBool rgb)
00093 {
00094 if (Init(pixel))
00095 {
00096 const T1 offset = OFstatic_cast(T1, DicomImageClass::maxval(bits - 1));
00097
00098
00099 const unsigned long count = (this->InputCount < this->Count) ? this->InputCount : this->Count;
00100 if (rgb)
00101 {
00102 register T2 *r = this->Data[0];
00103 register T2 *g = this->Data[1];
00104 register T2 *b = this->Data[2];
00105 const T2 maxvalue = OFstatic_cast(T2, DicomImageClass::maxval(bits));
00106 DiPixelRepresentationTemplate<T1> rep;
00107 if (bits == 8 && !rep.isSigned())
00108 {
00109 Sint16 rcr_tab[256];
00110 Sint16 gcb_tab[256];
00111 Sint16 gcr_tab[256];
00112 Sint16 bcb_tab[256];
00113 const double r_const = 0.7010 * OFstatic_cast(double, maxvalue);
00114 const double g_const = 0.5291 * OFstatic_cast(double, maxvalue);
00115 const double b_const = 0.8859 * OFstatic_cast(double, maxvalue);
00116 register unsigned long l;
00117 for (l = 0; l < 256; ++l)
00118 {
00119 rcr_tab[l] = OFstatic_cast(Sint16, 1.4020 * OFstatic_cast(double, l) - r_const);
00120 gcb_tab[l] = OFstatic_cast(Sint16, 0.3441 * OFstatic_cast(double, l));
00121 gcr_tab[l] = OFstatic_cast(Sint16, 0.7141 * OFstatic_cast(double, l) - g_const);
00122 bcb_tab[l] = OFstatic_cast(Sint16, 1.7720 * OFstatic_cast(double, l) - b_const);
00123 }
00124 register Sint32 sr;
00125 register Sint32 sg;
00126 register Sint32 sb;
00127 if (this->PlanarConfiguration)
00128 {
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143 register const T1 *y = pixel;
00144 register const T1 *cb = y + planeSize;
00145 register const T1 *cr = cb + planeSize;
00146 register unsigned long i = count;
00147 while (i != 0)
00148 {
00149
00150 for (l = planeSize; (l != 0) && (i != 0); --l, --i, ++y, ++cb, ++cr)
00151 {
00152 sr = OFstatic_cast(Sint32, *y) + OFstatic_cast(Sint32, rcr_tab[OFstatic_cast(Uint32, *cr)]);
00153 sg = OFstatic_cast(Sint32, *y) - OFstatic_cast(Sint32, gcb_tab[OFstatic_cast(Uint32, *cb)]) - OFstatic_cast(Sint32, gcr_tab[OFstatic_cast(Uint32, *cr)]);
00154 sb = OFstatic_cast(Sint32, *y) + OFstatic_cast(Sint32, bcb_tab[OFstatic_cast(Uint32, *cb)]);
00155 *(r++) = (sr < 0) ? 0 : (sr > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sr);
00156 *(g++) = (sg < 0) ? 0 : (sg > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sg);
00157 *(b++) = (sb < 0) ? 0 : (sb > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sb);
00158 }
00159
00160 y += 2 * planeSize;
00161 cb += 2 * planeSize;
00162 cr += 2 * planeSize;
00163 }
00164 }
00165 else
00166 {
00167 register const T1 *p = pixel;
00168 register T1 y;
00169 register T1 cb;
00170 register T1 cr;
00171 register unsigned long i;
00172 for (i = count; i != 0; --i)
00173 {
00174 y = *(p++);
00175 cb = *(p++);
00176 cr = *(p++);
00177 sr = OFstatic_cast(Sint32, y) + OFstatic_cast(Sint32, rcr_tab[OFstatic_cast(Uint32, cr)]);
00178 sg = OFstatic_cast(Sint32, y) - OFstatic_cast(Sint32, gcb_tab[OFstatic_cast(Uint32, cb)]) - OFstatic_cast(Sint32, gcr_tab[OFstatic_cast(Uint32, cr)]);
00179 sb = OFstatic_cast(Sint32, y) + OFstatic_cast(Sint32, bcb_tab[OFstatic_cast(Uint32, cb)]);
00180 *(r++) = (sr < 0) ? 0 : (sr > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sr);
00181 *(g++) = (sg < 0) ? 0 : (sg > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sg);
00182 *(b++) = (sb < 0) ? 0 : (sb > OFstatic_cast(Sint32, maxvalue)) ? maxvalue : OFstatic_cast(T2, sb);
00183 }
00184 }
00185 }
00186 else
00187 {
00188 if (this->PlanarConfiguration)
00189 {
00190
00191
00192
00193
00194
00195
00196
00197
00198 register unsigned long l;
00199 register unsigned long i = count;
00200 register const T1 *y = pixel;
00201 register const T1 *cb = y + planeSize;
00202 register const T1 *cr = cb + planeSize;
00203 while (i != 0)
00204 {
00205
00206 for (l = planeSize; (l != 0) && (i != 0); --l, --i)
00207 {
00208 convertValue(*(r++), *(g++), *(b++), removeSign(*(y++), offset), removeSign(*(cb++), offset),
00209 removeSign(*(cr++), offset), maxvalue);
00210 }
00211
00212 y += 2 * planeSize;
00213 cb += 2 * planeSize;
00214 cr += 2 * planeSize;
00215 }
00216 }
00217 else
00218 {
00219 register const T1 *p = pixel;
00220 register T2 y;
00221 register T2 cb;
00222 register T2 cr;
00223 register unsigned long i;
00224 for (i = count; i != 0; --i)
00225 {
00226 y = removeSign(*(p++), offset);
00227 cb = removeSign(*(p++), offset);
00228 cr = removeSign(*(p++), offset);
00229 convertValue(*(r++), *(g++), *(b++), y, cb, cr, maxvalue);
00230 }
00231 }
00232 }
00233 } else {
00234 register const T1 *p = pixel;
00235 if (this->PlanarConfiguration)
00236 {
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248
00249
00250
00251 register unsigned long l;
00252 register unsigned long i = 0;
00253 while (i < count)
00254 {
00255
00256 const unsigned long iStart = i;
00257 for (int j = 0; j < 3; ++j)
00258 {
00259
00260 for (l = planeSize, i = iStart; (l != 0) && (i < count); --l, ++i)
00261 this->Data[j][i] = removeSign(*(p++), offset);
00262 }
00263 }
00264 }
00265 else
00266 {
00267 register int j;
00268 register unsigned long i;
00269 for (i = 0; i < count; ++i)
00270 for (j = 0; j < 3; ++j)
00271 this->Data[j][i] = removeSign(*(p++), offset);
00272 }
00273 }
00274 }
00275 }
00276
00279 inline void convertValue(T2 &red, T2 &green, T2 &blue, const T2 y, const T2 cb, const T2 cr, const T2 maxvalue)
00280 {
00281 double dr = OFstatic_cast(double, y) + 1.4020 * OFstatic_cast(double, cr) - 0.7010 * OFstatic_cast(double, maxvalue);
00282 double dg = OFstatic_cast(double, y) - 0.3441 * OFstatic_cast(double, cb) - 0.7141 * OFstatic_cast(double, cr) + 0.5291 * OFstatic_cast(double, maxvalue);
00283 double db = OFstatic_cast(double, y) + 1.7720 * OFstatic_cast(double, cb) - 0.8859 * OFstatic_cast(double, maxvalue);
00284 red = (dr < 0.0) ? 0 : (dr > OFstatic_cast(double, maxvalue)) ? maxvalue : OFstatic_cast(T2, dr);
00285 green = (dg < 0.0) ? 0 : (dg > OFstatic_cast(double, maxvalue)) ? maxvalue : OFstatic_cast(T2, dg);
00286 blue = (db < 0.0) ? 0 : (db > OFstatic_cast(double, maxvalue)) ? maxvalue : OFstatic_cast(T2, db);
00287 }
00288 };
00289
00290
00291 #endif
00292
00293
00294
00295
00296
00297
00298
00299
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315
00316
00317
00318
00319
00320
00321
00322
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332
00333
00334
00335
00336
00337
00338
00339
00340
00341
00342
00343
00344
00345
00346
00347
00348
00349
00350
00351
00352
00353
00354
00355
00356
00357
00358
00359
00360
00361