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

Comparing ray/src/common/color.c (file contents):
Revision 1.4 by greg, Thu May 11 22:15:57 1989 UTC vs.
Revision 1.15 by greg, Wed Aug 28 08:39:47 1991 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1986 Regents of the University of California */
1 > /* Copyright (c) 1991 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 14 | Line 14 | static char SCCSid[] = "$SunId$ LBL";
14  
15   #include  "color.h"
16  
17 < #ifdef  SPEC_RGB
18 < /*
19 < *      The following table contains the CIE tristimulus integrals
20 < *  for X, Y, and Z.  The table is cumulative, so that
21 < *  each color coordinate integrates to 1.
22 < */
17 > #define  MINELEN        8       /* minimum scanline length for encoding */
18 > #define  MINRUN         4       /* minimum run length */
19  
24 – #define  STARTWL        380             /* starting wavelength (nanometers) */
25 – #define  INCWL          10              /* wavelength increment */
26 – #define  NINC           40              /* # of values */
20  
21 < static BYTE  chroma[3][NINC] = {
22 <        {                                                       /* X */
23 <                0,   0,   0,   2,   6,   13,  22,  30,  36,  41,
31 <                42,  43,  43,  44,  46,  52,  60,  71,  87,  106,
32 <                128, 153, 178, 200, 219, 233, 243, 249, 252, 254,
33 <                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
34 <        }, {                                                    /* Y */
35 <                0,   0,   0,   0,   0,   1,   2,   4,   7,   11,
36 <                17,  24,  34,  48,  64,  84,  105, 127, 148, 169,
37 <                188, 205, 220, 232, 240, 246, 250, 253, 254, 255,
38 <                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
39 <        }, {                                                    /* Z */
40 <                0,   0,   2,   10,  32,  66,  118, 153, 191, 220,
41 <                237, 246, 251, 253, 254, 255, 255, 255, 255, 255,
42 <                255, 255, 255, 255, 255, 255, 255, 255, 255, 255,
43 <                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
44 <        }
45 < };
46 <
47 <
48 < spec_rgb(col, s, e)             /* comput RGB color from spectral range */
49 < COLOR  col;
50 < int  s, e;
21 > char *
22 > tempbuffer(len)                 /* get a temporary buffer */
23 > unsigned  len;
24   {
25 <        COLOR  ciecolor;
25 >        extern char  *malloc(), *realloc();
26 >        static char  *tempbuf = NULL;
27 >        static int  tempbuflen = 0;
28  
29 <        spec_cie(ciecolor, s, e);
30 <        cie_rgb(col, ciecolor);
29 >        if (len > tempbuflen) {
30 >                if (tempbuflen > 0)
31 >                        tempbuf = realloc(tempbuf, len);
32 >                else
33 >                        tempbuf = malloc(len);
34 >                tempbuflen = tempbuf==NULL ? 0 : len;
35 >        }
36 >        return(tempbuf);
37   }
38  
39  
40 < spec_cie(col, s, e)             /* compute a color from a spectral range */
41 < COLOR  col;             /* returned color */
42 < int  s, e;              /* starting and ending wavelengths */
40 > fwritecolrs(scanline, len, fp)          /* write out a colr scanline */
41 > register COLR  *scanline;
42 > int  len;
43 > register FILE  *fp;
44   {
45 <        register int  i, d, r;
45 >        register int  i, j, beg, cnt;
46 >        int  c2;
47          
48 <        s -= STARTWL;
49 <        if (s < 0)
50 <                s = 0;
51 <
52 <        e -= STARTWL;
53 <        if (e >= INCWL*(NINC - 1))
54 <                e = INCWL*(NINC - 1) - 1;
55 <
56 <        d = e / INCWL;                  /* interpolate values */
57 <        r = e % INCWL;
58 <        for (i = 0; i < 3; i++)
59 <                col[i] = chroma[i][d]*(INCWL - r) + chroma[i][d + 1]*r;
60 <
61 <        d = s / INCWL;
62 <        r = s % INCWL;
63 <        for (i = 0; i < 3; i++)
64 <                col[i] -= chroma[i][d]*(INCWL - r) - chroma[i][d + 1]*r;
65 <
66 <        col[RED] = (col[RED] + 0.5) / (256*INCWL);
67 <        col[GRN] = (col[GRN] + 0.5) / (256*INCWL);
68 <        col[BLU] = (col[BLU] + 0.5) / (256*INCWL);
69 < }
70 <
71 <
72 < cie_rgb(rgbcolor, ciecolor)             /* convert CIE to RGB (NTSC) */
73 < register COLOR  rgbcolor, ciecolor;
74 < {
75 <        static float  cmat[3][3] = {
76 <                1.73, -.48, -.26,
77 <                -.81, 1.65, -.02,
78 <                 .08, -.17, 1.28,
79 <        };
80 <        register int  i;
81 <
82 <        for (i = 0; i < 3; i++) {
83 <                rgbcolor[i] =   cmat[i][0]*ciecolor[0] +
84 <                                cmat[i][1]*ciecolor[1] +
85 <                                cmat[i][2]*ciecolor[2] ;
86 <                if (rgbcolor[i] < 0.0)
87 <                        rgbcolor[i] = 0.0;
48 >        if (len < MINELEN)              /* too small to encode */
49 >                return(fwrite((char *)scanline,sizeof(COLR),len,fp) - len);
50 >        if (len > 32767)                /* too big! */
51 >                return(-1);
52 >        putc(2, fp);                    /* put magic header */
53 >        putc(2, fp);
54 >        putc(len>>8, fp);
55 >        putc(len&255, fp);
56 >                                        /* put components seperately */
57 >        for (i = 0; i < 4; i++) {
58 >            for (j = 0; j < len; j += cnt) {    /* find next run */
59 >                for (beg = j; beg < len; beg += cnt) {
60 >                    for (cnt = 1; cnt < 127 && beg+cnt < len &&
61 >                            scanline[beg+cnt][i] == scanline[beg][i]; cnt++)
62 >                        ;
63 >                    if (cnt >= MINRUN)
64 >                        break;                  /* long enough */
65 >                }
66 >                if (beg-j > 1 && beg-j < MINRUN) {
67 >                    c2 = j+1;
68 >                    while (scanline[c2++][i] == scanline[j][i])
69 >                        if (c2 == beg) {        /* short run */
70 >                            putc(128+beg-j, fp);
71 >                            putc(scanline[j][i], fp);
72 >                            j = beg;
73 >                            break;
74 >                        }
75 >                }
76 >                while (j < beg) {               /* write out non-run */
77 >                    if ((c2 = beg-j) > 128) c2 = 128;
78 >                    putc(c2, fp);
79 >                    while (c2--)
80 >                        putc(scanline[j++][i], fp);
81 >                }
82 >                if (cnt >= MINRUN) {            /* write out run */
83 >                    putc(128+cnt, fp);
84 >                    putc(scanline[beg][i], fp);
85 >                } else
86 >                    cnt = 0;
87 >            }
88          }
89 +        return(ferror(fp) ? -1 : 0);
90   }
107 – #endif
91  
92  
93 < fwritecolrs(scanline, len, fp)          /* write out a colr scanline */
93 > freadcolrs(scanline, len, fp)           /* read in an encoded colr scanline */
94   register COLR  *scanline;
95   int  len;
96   register FILE  *fp;
97   {
98 <        COLR  lastcolr;
99 <        int  rept;
100 <        
101 <        lastcolr[RED] = lastcolr[GRN] = lastcolr[BLU] = 1;
102 <        lastcolr[EXP] = 0;
103 <        rept = 0;
104 <        
105 <        while (len > 0) {
106 <                if (scanline[0][EXP] == lastcolr[EXP] &&
107 <                                scanline[0][RED] == lastcolr[RED] &&
125 <                                scanline[0][GRN] == lastcolr[GRN] &&
126 <                                scanline[0][BLU] == lastcolr[BLU])
127 <                        rept++;
128 <                else {
129 <                        while (rept) {          /* write out count */
130 <                                putc(1, fp);
131 <                                putc(1, fp);
132 <                                putc(1, fp);
133 <                                putc(rept & 255, fp);
134 <                                rept >>= 8;
135 <                        }
136 <                        putc(scanline[0][RED], fp);     /* new color */
137 <                        putc(scanline[0][GRN], fp);
138 <                        putc(scanline[0][BLU], fp);
139 <                        putc(scanline[0][EXP], fp);
140 <                        copycolr(lastcolr, scanline[0]);
141 <                        rept = 0;
142 <                }
143 <                scanline++;
144 <                len--;
98 >        register int  i, j;
99 >        int  code;
100 >                                        /* determine scanline type */
101 >        if (len < MINELEN)
102 >                return(oldreadcolrs(scanline, len, fp));
103 >        if ((i = getc(fp)) == EOF)
104 >                return(-1);
105 >        if (i != 2) {
106 >                ungetc(i, fp);
107 >                return(oldreadcolrs(scanline, len, fp));
108          }
109 <        while (rept) {          /* write out count */
110 <                putc(1, fp);
111 <                putc(1, fp);
112 <                putc(1, fp);
113 <                putc(rept & 255, fp);
114 <                rept >>= 8;
109 >        scanline[0][GRN] = getc(fp);
110 >        scanline[0][BLU] = getc(fp);
111 >        if ((i = getc(fp)) == EOF)
112 >                return(-1);
113 >        if (scanline[0][GRN] != 2 || scanline[0][BLU] & 128) {
114 >                scanline[0][RED] = 2;
115 >                scanline[0][EXP] = i;
116 >                return(oldreadcolrs(scanline+1, len-1, fp));
117          }
118 <        return(ferror(fp) ? -1 : 0);
118 >        if ((scanline[0][BLU]<<8 | i) != len)
119 >                return(-1);             /* length mismatch! */
120 >                                        /* read each component */
121 >        for (i = 0; i < 4; i++)
122 >            for (j = 0; j < len; ) {
123 >                if ((code = getc(fp)) == EOF)
124 >                    return(-1);
125 >                if (code > 128) {       /* run */
126 >                    scanline[j++][i] = getc(fp);
127 >                    for (code &= 127; --code; j++)
128 >                        scanline[j][i] = scanline[j-1][i];
129 >                } else                  /* non-run */
130 >                    while (code--)
131 >                        scanline[j++][i] = getc(fp);
132 >            }
133 >        return(feof(fp) ? -1 : 0);
134   }
135  
136  
137 < freadcolrs(scanline, len, fp)           /* read in a colr scanline */
137 > oldreadcolrs(scanline, len, fp)         /* read in an old colr scanline */
138   register COLR  *scanline;
139   int  len;
140   register FILE  *fp;
# Line 193 | Line 173 | register FILE  *fp;
173   fwritescan(scanline, len, fp)           /* write out a scanline */
174   register COLOR  *scanline;
175   int  len;
176 < register FILE  *fp;
176 > FILE  *fp;
177   {
178 <        COLR  lastcolr, thiscolr;
179 <        int  rept;
180 <        
181 <        lastcolr[RED] = lastcolr[GRN] = lastcolr[BLU] = 1;
182 <        lastcolr[EXP] = 0;
183 <        rept = 0;
184 <        
185 <        while (len > 0) {
186 <                setcolr(thiscolr, scanline[0][RED],
178 >        COLR  *clrscan;
179 >        int  n;
180 >        register COLR  *sp;
181 >                                        /* get scanline buffer */
182 >        if ((sp = (COLR *)tempbuffer(len*sizeof(COLR))) == NULL)
183 >                return(-1);
184 >        clrscan = sp;
185 >                                        /* convert scanline */
186 >        n = len;
187 >        while (n-- > 0) {
188 >                setcolr(sp[0], scanline[0][RED],
189                                    scanline[0][GRN],
190                                    scanline[0][BLU]);
209 –                if (thiscolr[EXP] == lastcolr[EXP] &&
210 –                                thiscolr[RED] == lastcolr[RED] &&
211 –                                thiscolr[GRN] == lastcolr[GRN] &&
212 –                                thiscolr[BLU] == lastcolr[BLU])
213 –                        rept++;
214 –                else {
215 –                        while (rept) {          /* write out count */
216 –                                putc(1, fp);
217 –                                putc(1, fp);
218 –                                putc(1, fp);
219 –                                putc(rept & 255, fp);
220 –                                rept >>= 8;
221 –                        }
222 –                        putc(thiscolr[RED], fp);        /* new color */
223 –                        putc(thiscolr[GRN], fp);
224 –                        putc(thiscolr[BLU], fp);
225 –                        putc(thiscolr[EXP], fp);
226 –                        copycolr(lastcolr, thiscolr);
227 –                        rept = 0;
228 –                }
191                  scanline++;
192 <                len--;
192 >                sp++;
193          }
194 <        while (rept) {          /* write out count */
233 <                putc(1, fp);
234 <                putc(1, fp);
235 <                putc(1, fp);
236 <                putc(rept & 255, fp);
237 <                rept >>= 8;
238 <        }
239 <        return(ferror(fp) ? -1 : 0);
194 >        return(fwritecolrs(clrscan, len, fp));
195   }
196  
197  
198   freadscan(scanline, len, fp)            /* read in a scanline */
199   register COLOR  *scanline;
200   int  len;
201 < register FILE  *fp;
201 > FILE  *fp;
202   {
203 <        COLR  thiscolr;
204 <        int  rshift;
205 <        register int  i;
206 <        
207 <        rshift = 0;
208 <        
209 <        while (len > 0) {
210 <                thiscolr[RED] = getc(fp);
211 <                thiscolr[GRN] = getc(fp);
212 <                thiscolr[BLU] = getc(fp);
213 <                thiscolr[EXP] = getc(fp);
214 <                if (feof(fp) || ferror(fp))
215 <                        return(-1);
216 <                if (thiscolr[RED] == 1 &&
217 <                                thiscolr[GRN] == 1 &&
218 <                                thiscolr[BLU] == 1) {
219 <                        for (i = thiscolr[EXP] << rshift; i > 0; i--) {
265 <                                copycolor(scanline[0], scanline[-1]);
266 <                                scanline++;
267 <                                len--;
268 <                        }
269 <                        rshift += 8;
270 <                } else {
271 <                        colr_color(scanline[0], thiscolr);
272 <                        scanline++;
273 <                        len--;
274 <                        rshift = 0;
275 <                }
203 >        register COLR  *clrscan;
204 >
205 >        if ((clrscan = (COLR *)tempbuffer(len*sizeof(COLR))) == NULL)
206 >                return(-1);
207 >        if (freadcolrs(clrscan, len, fp) < 0)
208 >                return(-1);
209 >                                        /* convert scanline */
210 >        colr_color(scanline[0], clrscan[0]);
211 >        while (--len > 0) {
212 >                scanline++; clrscan++;
213 >                if (clrscan[0][RED] == clrscan[-1][RED] &&
214 >                            clrscan[0][GRN] == clrscan[-1][GRN] &&
215 >                            clrscan[0][BLU] == clrscan[-1][BLU] &&
216 >                            clrscan[0][EXP] == clrscan[-1][EXP])
217 >                        copycolor(scanline[0], scanline[-1]);
218 >                else
219 >                        colr_color(scanline[0], clrscan[0]);
220          }
221          return(0);
222   }
# Line 308 | Line 252 | colr_color(col, clr)           /* convert short to float color
252   register COLOR  col;
253   register COLR  clr;
254   {
255 <        double  ldexp(), f;
255 >        double  f;
256          
257          if (clr[EXP] == 0)
258                  col[RED] = col[GRN] = col[BLU] = 0.0;
# Line 318 | Line 262 | register COLR  clr;
262                  col[GRN] = (clr[GRN] + 0.5)*f;
263                  col[BLU] = (clr[BLU] + 0.5)*f;
264          }
265 + }
266 +
267 +
268 + bigdiff(c1, c2, md)                     /* c1 delta c2 > md? */
269 + register COLOR  c1, c2;
270 + double  md;
271 + {
272 +        register int  i;
273 +
274 +        for (i = 0; i < 3; i++)
275 +                if (colval(c1,i)-colval(c2,i) > md*colval(c2,i) ||
276 +                        colval(c2,i)-colval(c1,i) > md*colval(c1,i))
277 +                        return(1);
278 +        return(0);
279   }

Diff Legend

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