ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/color.h
Revision: 2.21
Committed: Fri Jun 6 16:38:47 2003 UTC (20 years, 10 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.20: +9 -40 lines
Log Message:
*** empty log message ***

File Contents

# User Rev Content
1 schorsch 2.21 /* RCSid $Id: color.h,v 2.20 2003/05/13 17:58:32 greg Exp $ */
2 greg 1.1 /*
3     * color.h - header for routines using pixel color values.
4     *
5 gwlarson 2.17 * Must be included after X11 headers, since they declare a BYTE type.
6 greg 1.1 *
7     * Two color representations are used, one for calculation and
8     * another for storage. Calculation is done with three floats
9     * for speed. Stored color values use 4 bytes which contain
10     * three single byte mantissas and a common exponent.
11     */
12 schorsch 2.21 #ifndef _RAD_COLOR_H_
13     #define _RAD_COLOR_H_
14     #ifdef __cplusplus
15     extern "C" {
16     #endif
17    
18 greg 1.1
19 greg 2.19 #include "copyright.h"
20 greg 2.18
21     #include <stdio.h>
22     #include <stdlib.h>
23    
24 greg 1.1 #define RED 0
25     #define GRN 1
26     #define BLU 2
27 greg 2.10 #define CIEX 0 /* or, if input is XYZ... */
28     #define CIEY 1
29     #define CIEZ 2
30     #define EXP 3 /* exponent same for either format */
31 greg 1.1 #define COLXS 128 /* excess used for exponent */
32 greg 2.10 #define WHT 3 /* used for RGBPRIMS type */
33 greg 1.1
34 gwlarson 2.17 #undef BYTE
35     #define BYTE unsigned char /* 8-bit unsigned integer */
36 greg 1.1
37 greg 2.10 typedef BYTE COLR[4]; /* red, green, blue (or X,Y,Z), exponent */
38 greg 1.1
39 greg 2.10 typedef float COLOR[3]; /* red, green, blue (or X,Y,Z) */
40    
41     typedef float RGBPRIMS[4][2]; /* (x,y) chromaticities for RGBW */
42     typedef float (*RGBPRIMP)[2]; /* pointer to RGBPRIMS array */
43    
44     typedef float COLORMAT[3][3]; /* color coordinate conversion matrix */
45    
46 greg 1.1 #define copycolr(c1,c2) (c1[0]=c2[0],c1[1]=c2[1], \
47     c1[2]=c2[2],c1[3]=c2[3])
48    
49     #define colval(col,pri) ((col)[pri])
50    
51     #define setcolor(col,r,g,b) ((col)[RED]=(r),(col)[GRN]=(g),(col)[BLU]=(b))
52    
53     #define copycolor(c1,c2) ((c1)[0]=(c2)[0],(c1)[1]=(c2)[1],(c1)[2]=(c2)[2])
54    
55     #define scalecolor(col,sf) ((col)[0]*=(sf),(col)[1]*=(sf),(col)[2]*=(sf))
56    
57     #define addcolor(c1,c2) ((c1)[0]+=(c2)[0],(c1)[1]+=(c2)[1],(c1)[2]+=(c2)[2])
58    
59     #define multcolor(c1,c2) ((c1)[0]*=(c2)[0],(c1)[1]*=(c2)[1],(c1)[2]*=(c2)[2])
60    
61 greg 1.7 #ifdef NTSC
62 greg 2.7 #define CIE_x_r 0.670 /* standard NTSC primaries */
63     #define CIE_y_r 0.330
64     #define CIE_x_g 0.210
65     #define CIE_y_g 0.710
66     #define CIE_x_b 0.140
67     #define CIE_y_b 0.080
68 greg 2.8 #define CIE_x_w 0.3333 /* use true white */
69     #define CIE_y_w 0.3333
70 greg 1.7 #else
71 greg 2.7 #define CIE_x_r 0.640 /* nominal CRT primaries */
72     #define CIE_y_r 0.330
73     #define CIE_x_g 0.290
74     #define CIE_y_g 0.600
75     #define CIE_x_b 0.150
76     #define CIE_y_b 0.060
77 greg 2.8 #define CIE_x_w 0.3333 /* use true white */
78     #define CIE_y_w 0.3333
79 greg 1.7 #endif
80 greg 2.7
81 greg 2.14 #define STDPRIMS {{CIE_x_r,CIE_y_r},{CIE_x_g,CIE_y_g}, \
82     {CIE_x_b,CIE_y_b},{CIE_x_w,CIE_y_w}}
83 greg 2.10
84 greg 2.7 #define CIE_D ( CIE_x_r*(CIE_y_g - CIE_y_b) + \
85     CIE_x_g*(CIE_y_b - CIE_y_r) + \
86     CIE_x_b*(CIE_y_r - CIE_y_g) )
87     #define CIE_C_rD ( (1./CIE_y_w) * \
88     ( CIE_x_w*(CIE_y_g - CIE_y_b) - \
89     CIE_y_w*(CIE_x_g - CIE_x_b) + \
90     CIE_x_g*CIE_y_b - CIE_x_b*CIE_y_g ) )
91     #define CIE_C_gD ( (1./CIE_y_w) * \
92     ( CIE_x_w*(CIE_y_b - CIE_y_r) - \
93     CIE_y_w*(CIE_x_b - CIE_x_r) - \
94     CIE_x_r*CIE_y_b + CIE_x_b*CIE_y_r ) )
95     #define CIE_C_bD ( (1./CIE_y_w) * \
96     ( CIE_x_w*(CIE_y_r - CIE_y_g) - \
97     CIE_y_w*(CIE_x_r - CIE_x_g) + \
98     CIE_x_r*CIE_y_g - CIE_x_g*CIE_y_r ) )
99    
100     #define CIE_rf (CIE_y_r*CIE_C_rD/CIE_D)
101     #define CIE_gf (CIE_y_g*CIE_C_gD/CIE_D)
102     #define CIE_bf (CIE_y_b*CIE_C_bD/CIE_D)
103    
104 greg 2.9 /* As of 9-94, CIE_rf=.265074126, CIE_gf=.670114631 and CIE_bf=.064811243 */
105    
106 greg 2.10 /***** The following definitions are valid for RGB colors only... *****/
107    
108 greg 2.7 #define bright(col) (CIE_rf*(col)[RED]+CIE_gf*(col)[GRN]+CIE_bf*(col)[BLU])
109     #define normbright(c) ( ( (long)(CIE_rf*256.+.5)*(c)[RED] + \
110     (long)(CIE_gf*256.+.5)*(c)[GRN] + \
111     (long)(CIE_bf*256.+.5)*(c)[BLU] ) >> 8 )
112 greg 1.1
113 greg 1.14 /* luminous efficacies over visible spectrum */
114     #define MAXEFFICACY 683. /* defined maximum at 550 nm */
115     #define WHTEFFICACY 179. /* uniform white light */
116     #define D65EFFICACY 203. /* standard illuminant D65 */
117     #define INCEFFICACY 160. /* illuminant A (incand.) */
118     #define SUNEFFICACY 208. /* illuminant B (solar dir.) */
119     #define SKYEFFICACY D65EFFICACY /* skylight */
120     #define DAYEFFICACY D65EFFICACY /* combined sky and solar */
121    
122     #define luminance(col) (WHTEFFICACY * bright(col))
123 greg 1.9
124 greg 2.10 /***** ...end of stuff specific to RGB colors *****/
125    
126 greg 1.1 #define intens(col) ( (col)[0] > (col)[1] \
127     ? (col)[0] > (col)[2] ? (col)[0] : (col)[2] \
128     : (col)[1] > (col)[2] ? (col)[1] : (col)[2] )
129    
130 greg 1.3 #define colrval(c,p) ( (c)[EXP] ? \
131     ldexp((c)[p]+.5,(int)(c)[EXP]-(COLXS+8)) : \
132     0. )
133    
134 greg 1.1 #define WHTCOLOR {1.0,1.0,1.0}
135     #define BLKCOLOR {0.0,0.0,0.0}
136     #define WHTCOLR {128,128,128,COLXS+1}
137     #define BLKCOLR {0,0,0,0}
138 greg 1.3
139 greg 1.11 /* picture format identifier */
140     #define COLRFMT "32-bit_rle_rgbe"
141 greg 2.10 #define CIEFMT "32-bit_rle_xyze"
142     #define PICFMT "32-bit_rle_???e" /* matches either */
143     #define LPICFMT 15 /* max format id len */
144 greg 1.11
145 greg 1.6 /* macros for exposures */
146 greg 1.5 #define EXPOSSTR "EXPOSURE="
147 greg 1.12 #define LEXPOSSTR 9
148     #define isexpos(hl) (!strncmp(hl,EXPOSSTR,LEXPOSSTR))
149     #define exposval(hl) atof((hl)+LEXPOSSTR)
150 greg 1.5 #define fputexpos(ex,fp) fprintf(fp,"%s%e\n",EXPOSSTR,ex)
151 greg 1.6
152     /* macros for pixel aspect ratios */
153     #define ASPECTSTR "PIXASPECT="
154 greg 1.12 #define LASPECTSTR 10
155     #define isaspect(hl) (!strncmp(hl,ASPECTSTR,LASPECTSTR))
156     #define aspectval(hl) atof((hl)+LASPECTSTR)
157 greg 1.6 #define fputaspect(pa,fp) fprintf(fp,"%s%f\n",ASPECTSTR,pa)
158 greg 1.5
159 greg 2.10 /* macros for primary specifications */
160     #define PRIMARYSTR "PRIMARIES="
161     #define LPRIMARYSTR 10
162     #define isprims(hl) (!strncmp(hl,PRIMARYSTR,LPRIMARYSTR))
163     #define primsval(p,hl) sscanf(hl+LPRIMARYSTR, \
164     "%f %f %f %f %f %f %f %f", \
165     &(p)[RED][CIEX],&(p)[RED][CIEY], \
166     &(p)[GRN][CIEX],&(p)[GRN][CIEY], \
167     &(p)[BLU][CIEX],&(p)[BLU][CIEY], \
168     &(p)[WHT][CIEX],&(p)[WHT][CIEY])
169     #define fputprims(p,fp) fprintf(fp, \
170     "%s %.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f\n",\
171     PRIMARYSTR, \
172     (p)[RED][CIEX],(p)[RED][CIEY], \
173     (p)[GRN][CIEX],(p)[GRN][CIEY], \
174     (p)[BLU][CIEX],(p)[BLU][CIEY], \
175     (p)[WHT][CIEX],(p)[WHT][CIEY])
176    
177 greg 1.8 /* macros for color correction */
178     #define COLCORSTR "COLORCORR="
179 greg 1.12 #define LCOLCORSTR 10
180     #define iscolcor(hl) (!strncmp(hl,COLCORSTR,LCOLCORSTR))
181     #define colcorval(cc,hl) sscanf(hl+LCOLCORSTR,"%f %f %f", \
182 greg 1.8 &(cc)[RED],&(cc)[GRN],&(cc)[BLU])
183     #define fputcolcor(cc,fp) fprintf(fp,"%s %f %f %f\n",COLCORSTR, \
184     (cc)[RED],(cc)[GRN],(cc)[BLU])
185 greg 2.6
186 greg 2.10 /*
187     * Conversions to and from XYZ space generally don't apply WHTEFFICACY.
188     * If you need Y to be luminance (cd/m^2), this must be applied when
189     * converting from radiance (watts/sr/m^2).
190     */
191    
192 greg 2.18 extern RGBPRIMS stdprims; /* standard primary chromaticities */
193     extern COLORMAT rgb2xyzmat; /* RGB to XYZ conversion matrix */
194     extern COLORMAT xyz2rgbmat; /* XYZ to RGB conversion matrix */
195     extern COLOR cblack, cwhite; /* black (0,0,0) and white (1,1,1) */
196 greg 2.10
197 greg 2.13 #define CGAMUT_LOWER 01
198     #define CGAMUT_UPPER 02
199     #define CGAMUT (CGAMUT_LOWER|CGAMUT_UPPER)
200    
201     #define rgb_cie(xyz,rgb) colortrans(xyz,rgb2xyzmat,rgb)
202 greg 2.10
203     #ifdef BSD
204 greg 2.20 #define cpcolormat(md,ms) bcopy((void *)ms,(void *)md,sizeof(COLORMAT))
205 greg 2.10 #else
206 greg 2.20 #define cpcolormat(md,ms) memcpy((void *)md,(void *)ms,sizeof(COLORMAT))
207 greg 2.18 #endif
208    
209     /* defined in color.c */
210     extern char *tempbuffer(unsigned int len);
211     extern int fwritecolrs(COLR *scanline, int len, FILE *fp);
212     extern int freadcolrs(COLR *scanline, int len, FILE *fp);
213     extern int fwritescan(COLOR *scanline, int len, FILE *fp);
214     extern int freadscan(COLOR *scanline, int len, FILE *fp);
215     extern void setcolr(COLR clr, double r, double g, double b);
216     extern void colr_color(COLOR col, COLR clr);
217     extern int bigdiff(COLOR c1, COLOR c2, double md);
218     /* defined in spec_rgb.c */
219     extern void spec_rgb(COLOR col, int s, int e);
220     extern void spec_cie(COLOR col, int s, int e);
221     extern void cie_rgb(COLOR rgb, COLOR xyz);
222     extern int clipgamut(COLOR col, double brt, int gamut,
223     COLOR lower, COLOR upper);
224     extern void colortrans(COLOR c2, COLORMAT mat, COLOR c1);
225     extern void multcolormat(COLORMAT m3, COLORMAT m2,
226     COLORMAT m1);
227     extern void compxyz2rgbmat(COLORMAT mat, RGBPRIMS pr);
228     extern void comprgb2xyzmat(COLORMAT mat, RGBPRIMS pr);
229     extern void comprgb2rgbmat(COLORMAT mat, RGBPRIMS pr1, RGBPRIMS pr2);
230     extern void compxyzWBmat(COLORMAT mat, float wht1[2],
231     float wht2[2]);
232     extern void compxyz2rgbWBmat(COLORMAT mat, RGBPRIMS pr);
233     extern void comprgb2xyzWBmat(COLORMAT mat, RGBPRIMS pr);
234     extern void comprgb2rgbWBmat(COLORMAT mat, RGBPRIMS pr1, RGBPRIMS pr2);
235     /* defined in colrops.c */
236     extern int setcolrcor(double (*f)(), double a2);
237     extern int setcolrinv(double (*f)(), double a2);
238     extern int setcolrgam(double g);
239     extern int colrs_gambs(COLR *scan, int len);
240     extern int gambs_colrs(COLR *scan, int len);
241     extern void shiftcolrs(COLR *scan, int len, int adjust);
242     extern void normcolrs(COLR *scan, int len, int adjust);
243    
244    
245     #ifdef __cplusplus
246     }
247 greg 2.6 #endif
248 schorsch 2.21 #endif /* _RAD_COLOR_H_ */
249