ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/tmapluv.c
(Generate patch)

Comparing ray/src/common/tmapluv.c (file contents):
Revision 3.1 by gwlarson, Thu Oct 22 13:53:06 1998 UTC vs.
Revision 3.5 by greg, Tue Feb 25 02:47:22 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Routines for tone-mapping LogLuv encoded pixels.
6 + *
7 + * Externals declared in tmaptiff.h
8   */
9  
10 + #include "copyright.h"
11 +
12 + #define LOGLUV_PUBLIC           1
13 +
14   #include <stdio.h>
15 + #include <string.h>
16   #include <math.h>
17   #include "tiffio.h"
18   #include "tmprivat.h"
19 < #include "uvcode.h"
19 > #include "tmaptiff.h"
20  
21   #ifndef BSD
22   #define bzero(d,n)              (void)memset(d,0,n)
# Line 23 | Line 27 | static char SCCSid[] = "$SunId$ SGI";
27   #define clruv(p,uv)             uvflgop(p,uv,&=~)
28   #define clruvall(p)             bzero((MEM_PTR)(p)->rgbflg,sizeof((p)->rgbflg))
29  
30 < #define U_NEU                   0.210526316
27 < #define V_NEU                   0.473684211
28 <
30 > #ifdef NOPROTO
31   static MEM_PTR  luv32Init();
32   static void     luv32NewSpace();
33   static MEM_PTR  luv24Init();
34   static void     luv24NewSpace();
35 < extern void     free();
35 > #else
36 > static MEM_PTR  luv32Init(struct tmStruct *);
37 > static void     luv32NewSpace(struct tmStruct *);
38 > static MEM_PTR  luv24Init(struct tmStruct *);
39 > static void     luv24NewSpace(struct tmStruct *);
40 > #endif
41  
42   typedef struct {
43          int     offset;                 /* computed luminance offset */
# Line 38 | Line 45 | typedef struct {
45          uint32  rgbflg[1<<(16-5)];      /* flags for computed values */
46   } LUV32DATA;                    /* LogLuv 32-bit conversion data */
47  
41 – #define UVSCALE                 410.
48   #define UVNEU                   ((int)(UVSCALE*U_NEU)<<8 \
49                                          | (int)(UVSCALE*V_NEU))
50  
# Line 61 | Line 67 | static int     luv24Reg = -1;          /* 24-bit package reg. numb
67   static int      uv14neu = -1;           /* neutral index for 14-bit (u',v') */
68  
69  
70 < static
70 > static void
71   uv2rgb(rgb, tm, uvp)                    /* compute RGB from uv coordinate */
72   BYTE    rgb[3];
73   register struct tmStruct        *tm;
74   double  uvp[2];
75   {       /* Should check that tm->inppri==TM_XYZPRIM beforehand... */
76          double  d, x, y;
77 <        COLR    XYZ, RGB;
77 >        COLOR   XYZ, RGB;
78                                          /* convert to XYZ */
79          d = 1./(6.*uvp[0] - 16.*uvp[1] + 12.);
80          x = 9.*uvp[0] * d;
# Line 78 | Line 84 | double uvp[2];
84          XYZ[CIEZ] = (1.-x-y)/y * XYZ[CIEY];
85                                          /* convert to RGB and clip */
86          colortrans(RGB, tm->cmat, XYZ);
87 <        clipgamut(RGB, XYZ[CIEY], CGAMUT_LOWER, cblack, cwhite);
87 >        clipgamut(RGB, 1., CGAMUT_LOWER, cblack, cwhite);
88                                          /* perform final scaling & gamma */
89          d = tm->clf[RED] * RGB[RED];
90          rgb[RED] = d>=.999 ? 255 : (int)(256.*pow(d, 1./tm->mongam));
# Line 94 | Line 100 | compmeshift(li, uvp)                   /* compute mesopic color shift
100   TMbright        li;             /* encoded world luminance */
101   double  uvp[2];                 /* world (u',v') -> returned desaturated */
102   {
103 <        /* UNIMPLEMENTED */
104 <        return(li);
99 < }
103 >        double  scotrat;
104 >        register double d;
105  
106 <
107 < static int
108 < uvpencode(uvp)                  /* encode (u',v') coordinates */
109 < double  uvp[2];
110 < {
111 <        register int    vi, ui;
112 <
113 <        if (uvp[1] < UV_VSTART)
114 <                return(-1);
115 <        vi = (uvp[1] - UV_VSTART)*(1./UV_SQSIZ);
116 <        if (vi >= UV_NVS)
117 <                return(-1);
118 <        if (uvp[0] < uv_row[vi].ustart)
119 <                return(-1);
120 <        ui = (uvp[0] - uv_row[vi].ustart)*(1./UV_SQSIZ);
121 <        if (ui >= uv_row[vi].nus)
122 <                return(-1);
123 <        return(uv_row[vi].ncum + ui);
106 >        if (li >= BMESUPPER)
107 >                return(li);
108 >        scotrat = (.767676768 - 1.02356902*uvp[1])/uvp[0] - .343434343;
109 >        if (li <= BMESLOWER) {
110 >                d = 0.;
111 >                uvp[0] = U_NEU; uvp[1] = V_NEU;
112 >        } else {
113 >                d = (tmMesofact[li-BMESLOWER] + .5) * (1./256.);
114 >                uvp[0] = d*uvp[0] + (1.-d)*U_NEU;
115 >                uvp[1] = d*uvp[1] + (1.-d)*V_NEU;
116 >        }
117 >        /*
118 >        d = li + (double)TM_BRTSCALE*log(d + (1.-d)*scotrat);
119 >        */
120 >        d = d + (1.-d)*scotrat;
121 >        d -= 1.;                        /* Taylor expansion of log(x) about 1 */
122 >        d = d*(1. + d*(-.5 + d*(1./3. + d*-.125)));
123 >        d = li + (double)TM_BRTSCALE*d;
124 >        return((TMbright)(d>0. ? d+.5 : d-.5));
125   }
126  
127  
122 – static int
123 – uvpdecode(uvp, c)               /* decode (u',v') index */
124 – double  uvp[2];
125 – int     c;
126 – {
127 –        int     upper, lower;
128 –        register int    ui, vi;
129 –
130 –        if (c < 0 || c >= UV_NDIVS)
131 –                return(-1);
132 –        lower = 0;                      /* binary search */
133 –        upper = UV_NVS;
134 –        do {
135 –                vi = (lower + upper) >> 1;
136 –                ui = c - uv_row[vi].ncum;
137 –                if (ui > 0)
138 –                        lower = vi;
139 –                else if (ui < 0)
140 –                        upper = vi;
141 –                else
142 –                        break;
143 –        } while (upper - lower > 1);
144 –        vi = lower;
145 –        ui = c - uv_row[vi].ncum;
146 –        uvp[0] = uv_row[vi].ustart + (ui+.5)*UV_SQSIZ;
147 –        uvp[1] = UV_VSTART + (vi+.5)*UV_SQSIZ;
148 –        return(0);
149 – }
150 –
151 –
128   int
129   tmCvLuv32(ls, cs, luvs, len)            /* convert raw 32-bit LogLuv values */
130   TMbright        *ls;
# Line 166 | Line 142 | int    len;
142          if (ls == NULL | luvs == NULL | len < 0)
143                  returnErr(TM_E_ILLEGAL);
144                                          /* check package registration */
145 <        if (luv32Reg < 0 && (luv32Reg = tmRegPkg(&luv32Pkg)) < 0)
146 <                returnErr(TM_E_CODERR1);
145 >        if (luv32Reg < 0) {
146 >                if ((luv32Reg = tmRegPkg(&luv32Pkg)) < 0)
147 >                        returnErr(TM_E_CODERR1);
148 >                tmMkMesofact();
149 >        }
150                                          /* get package data */
151          if ((ld = (LUV32DATA *)tmPkgData(tmTop,luv32Reg)) == NULL)
152                  returnErr(TM_E_NOMEM);
# Line 175 | Line 154 | int    len;
154          for (i = len; i--; ) {
155                  j = luvs[i] >> 16;              /* get luminance */
156                  if (j & 0x8000)                 /* negative luminance */
157 <                        ls[i] = MINBRT-1;       /* assign bogus value */
157 >                        ls[i] = TM_NOBRT;       /* assign bogus value */
158                  else                            /* else convert to lnL */
159 <                        ls[i] = (BRT2SCALE*j >> 8) - ld->offset;
159 >                        ls[i] = (BRT2SCALE(j) >> 8) - ld->offset;
160                  if (cs == TM_NOCHROM)           /* no color? */
161                          continue;
162                                                  /* get chrominance */
# Line 185 | Line 164 | int    len;
164                          uvp[0] = 1./UVSCALE*((luvs[i]>>8 & 0xff) + .5);
165                          uvp[1] = 1./UVSCALE*((luvs[i] & 0xff) + .5);
166                          ls[i] = compmeshift(ls[i], uvp);
167 <                        j = tmTop->flags & TM_F_BW || ls[i] < BMESLOWER
167 >                        j = tmTop->flags&TM_F_BW || ls[i]<BMESLOWER
168                                          ? UVNEU
169                                          : (int)(uvp[0]*UVSCALE)<<8
170                                                  | (int)(uvp[1]*UVSCALE);
# Line 193 | Line 172 | int    len;
172                          j = tmTop->flags&TM_F_BW ? UVNEU : luvs[i]&0xffff;
173                  }
174                  if (!isuvset(ld, j)) {
175 <                        uvp[0] = 1./UVSCALE*((luvs[i]>>8 & 0xff) + .5);
176 <                        uvp[1] = 1./UVSCALE*((luvs[i] & 0xff) + .5);
175 >                        uvp[0] = 1./UVSCALE*((j>>8) + .5);
176 >                        uvp[1] = 1./UVSCALE*((j & 0xff) + .5);
177                          uv2rgb(ld->rgbval[j], tmTop, uvp);
178                          setuv(ld, j);
179                  }
# Line 223 | Line 202 | int    len;
202          if (ls == NULL | luvs == NULL | len < 0)
203                  returnErr(TM_E_ILLEGAL);
204                                          /* check package registration */
205 <        if (luv24Reg < 0 && (luv24Reg = tmRegPkg(&luv24Pkg)) < 0)
206 <                returnErr(TM_E_CODERR1);
205 >        if (luv24Reg < 0) {
206 >                if ((luv24Reg = tmRegPkg(&luv24Pkg)) < 0)
207 >                        returnErr(TM_E_CODERR1);
208 >                tmMkMesofact();
209 >        }
210                                          /* get package data */
211          if ((ld = (LUV24DATA *)tmPkgData(tmTop,luv24Reg)) == NULL)
212                  returnErr(TM_E_NOMEM);
213                                          /* convert each pixel */
214          for (i = len; i--; ) {
215                  j = luvs[i] >> 14;              /* get luminance */
216 <                ls[i] = (BRT2SCALE*j >> 6) - ld->offset;
216 >                ls[i] = (BRT2SCALE(j) >> 6) - ld->offset;
217                  if (cs == TM_NOCHROM)           /* no color? */
218                          continue;
219                                                  /* get chrominance */
220                  if (tmTop->flags & TM_F_MESOPIC && ls[i] < BMESUPPER) {
221 <                        if (uvpdecode(uvp, luvs[i]&0x3fff) < 0) {
221 >                        if (uv_decode(uvp, uvp+1, luvs[i]&0x3fff) < 0) {
222                                  uvp[0] = U_NEU;         /* should barf? */
223                                  uvp[1] = V_NEU;
224                          }
225                          ls[i] = compmeshift(ls[i], uvp);
226 <                        if (tmTop->flags & TM_F_BW || ls[i] < BMESLOWER
227 <                                        || (j = uvpencode(uvp)) < 0)
226 >                        if (tmTop->flags&TM_F_BW || ls[i]<BMESLOWER
227 >                                        || (j = uv_encode(uvp[0], uvp[1],
228 >                                                SGILOGENCODE_NODITHER)) < 0)
229                                  j = uv14neu;
230                  } else {
231                          j = tmTop->flags&TM_F_BW ? uv14neu : luvs[i]&0x3fff;
232                  }
233                  if (!isuvset(ld, j)) {
234 <                        if (uvpdecode(uvp, j) < 0) {
235 <                                uvp[0] = U_NEU;
253 <                                uvp[1] = V_NEU;
234 >                        if (uv_decode(&uvp[0], &uvp[1], j) < 0) {
235 >                                uvp[0] = U_NEU; uvp[1] = V_NEU;
236                          }
237                          uv2rgb(ld->rgbval[j], tmTop, uvp);
238                          setuv(ld, j);
# Line 280 | Line 262 | int    len;
262                  returnErr(TM_E_ILLEGAL);
263                                          /* check scaling offset */
264          if (!FEQ(tmTop->inpsf, lastsf)) {
265 <                offset = BRT2SCALE*64;
265 >                offset = BRT2SCALE(64);
266                  if (tmTop->inpsf > 1.0001)
267                          offset -= (int)(TM_BRTSCALE*log(tmTop->inpsf)+.5);
268                  else if (tmTop->inpsf < 0.9999)
# Line 290 | Line 272 | int    len;
272                                          /* convert each pixel */
273          for (i = len; i--; ) {
274                  if (l16s[i] & 0x8000)           /* negative luminance */
275 <                        ls[i] = MINBRT-1;       /* assign bogus value */
275 >                        ls[i] = TM_NOBRT;       /* assign bogus value */
276                  else                            /* else convert to lnL */
277 <                        ls[i] = (BRT2SCALE*l16s[i] >> 8) - offset;
277 >                        ls[i] = (BRT2SCALE(l16s[i]) >> 8) - offset;
278          }
279          returnOK;
280   }
# Line 309 | Line 291 | struct tmStruct        *tm;
291                  exit(1);
292          }
293          ld = (LUV32DATA *)tm->pd[luv32Reg];
294 <        ld->offset = BRT2SCALE*64;
294 >        ld->offset = BRT2SCALE(64);
295          if (tm->inpsf > 1.0001)
296                  ld->offset -= (int)(TM_BRTSCALE*log(tmTop->inpsf)+.5);
297          else if (tm->inpsf < 0.9999)
# Line 338 | Line 320 | luv24NewSpace(tm)              /* initialize 24-bit LogLuv color s
320   struct tmStruct *tm;
321   {
322          register LUV24DATA      *ld;
341 –        double  uvp[2];
323  
324          if (tm->inppri != TM_XYZPRIM) {         /* panic time! */
325                  fputs("Improper input color space in luv24NewSpace!\n", stderr);
326                  exit(1);
327          }
328          ld = (LUV24DATA *)tm->pd[luv24Reg];
329 <        ld->offset = BRT2SCALE*12;
329 >        ld->offset = BRT2SCALE(12);
330          if (tm->inpsf > 1.0001)
331                  ld->offset -= (int)(TM_BRTSCALE*log(tmTop->inpsf)+.5);
332          else if (tm->inpsf < 0.9999)
# Line 364 | Line 345 | struct tmStruct        *tm;
345          if (ld == NULL)
346                  return(NULL);
347          tm->pd[luv24Reg] = (MEM_PTR)ld;
348 <        if (uv14neu < 0) {              /* initialize neutral color index */
349 <                double  uvp[2];
369 <                uvp[0] = U_NEU; uvp[1] = V_NEU;
370 <                uv14neu = uvpencode(uvp);
371 <        }
348 >        if (uv14neu < 0)                /* initialize neutral color index */
349 >                uv14neu = uv_encode(U_NEU, V_NEU, SGILOGENCODE_NODITHER);
350          luv24NewSpace(tm);
351          return((MEM_PTR)ld);
352   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines