ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/lampcolor.c
Revision: 2.3
Committed: Fri Jan 10 11:17:13 1992 UTC (32 years, 10 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.2: +2 -2 lines
Log Message:
change white to D65WHITE to match lamp.tab

File Contents

# User Rev Content
1 greg 1.1 /* Copyright (c) 1991 Regents of the University of California */
2    
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8 greg 1.2 * Program to convert lamp color from table and compute radiance.
9 greg 1.1 */
10    
11     #include <stdio.h>
12    
13 greg 1.8 #include "color.h"
14    
15 greg 1.2 #define PI 3.14159265358979323846
16    
17     extern char *gets(), *strcpy();
18 greg 2.2 #ifndef atof
19 greg 1.2 extern double atof();
20 greg 2.2 #endif
21 greg 1.1 extern float *matchlamp();
22    
23 greg 1.2 /* lamp parameters */
24     #define LTYPE 0
25     #define LUNIT 1
26     #define LGEOM 2
27     #define LOUTP 3
28     #define NPARAMS 4
29 greg 1.1
30 greg 1.2 int typecheck(), unitcheck(), geomcheck(), outpcheck();
31    
32     float *lampcolor; /* the lamp color (RGB) */
33     double unit2meter; /* conversion from units to meters */
34 greg 1.7 double area; /* radiating area for this geometry */
35 greg 1.2 double lumens; /* total lamp lumens */
36    
37     struct {
38     char *name;
39     char value[64];
40     int (*check)();
41     char *help;
42     } param[NPARAMS] = {
43 greg 2.3 { "lamp type", "D65WHITE", typecheck,
44 greg 1.2 "The lamp type is a string which corresponds to one of the types registered\n\
45 greg 2.3 in the lamp table file. A value of \"D65WHITE\" means an uncolored source,\n\
46 greg 1.2 which may be preferable because it results in a color balanced image." },
47     { "length unit", "meter", unitcheck,
48     "Unit must be one of: \"meter\", \"centimeter\", \"foot\", or \"inch\".\n\
49     These may be abbreviated as a single letter." },
50     { "lamp geometry", "polygon", geomcheck,
51     "The lamp geometry must be one of: \"polygon\", \"sphere\", \"cylinder\"\n\
52     or \"ring\". These may be abbreviated as a single letter." },
53     { "total lamp lumens", "0", outpcheck,
54     "This is the overall light output of the lamp and its fixture. If you do\n\
55     not know this value explicitly, you can compute the approximate lumens\n\
56 greg 1.3 by multiplying the input wattage by 14 for incandescent fixtures or 70\n\
57 greg 1.2 for fluorescent fixtures." },
58     };
59    
60    
61 greg 1.1 main(argc, argv)
62     int argc;
63     char *argv[];
64     {
65     char *lamptab = "lamp.tab";
66 greg 1.2 char buf[64];
67     int i;
68 greg 1.1
69     if (argc > 1) lamptab = argv[1];
70     if (loadlamps(lamptab) == 0) {
71 greg 1.2 fprintf(stderr, "%s: no such lamp table\n", lamptab);
72 greg 1.1 exit(1);
73     }
74 greg 1.2 printf("Program to compute lamp radiance. Enter '?' for help.\n");
75 greg 1.1 for ( ; ; ) {
76 greg 1.2 i = 0;
77     while (i < NPARAMS) {
78     printf("Enter %s [%s]: ", param[i].name,
79     param[i].value);
80     if (gets(buf) == NULL)
81     exit(0);
82     if (buf[0] == '?') {
83     puts(param[i].help);
84     continue;
85     }
86     if (buf[0])
87     strcpy(param[i].value, buf);
88     if (!(*param[i].check)(param[i].value)) {
89     fprintf(stderr, "%s: bad value for %s\n",
90     param[i].value, param[i].name);
91     continue;
92     }
93     i++;
94     }
95     compute();
96 greg 1.1 }
97 greg 1.2 }
98    
99    
100     typecheck(s) /* check lamp type */
101     char *s;
102     {
103     lampcolor = matchlamp(s);
104     return(lampcolor != NULL);
105     }
106    
107    
108     unitcheck(s) /* compute conversion to meters */
109     char *s;
110     {
111     int len = strlen(s);
112    
113     switch (*s) {
114     case 'm':
115     if (strncmp(s, "meters", len))
116     return(0);
117     unit2meter = 1.0;
118     return(1);
119     case 'c':
120     if (strncmp(s, "centimeters", len) && strncmp(s, "cms", len))
121     return(0);
122     unit2meter = 0.01;
123     return(1);
124     case 'f':
125     if (strncmp(s, "foot", len) && strncmp(s, "feet", len))
126     return(0);
127     unit2meter = 0.3048;
128     return(1);
129     case 'i':
130     if (strncmp(s, "inches", len))
131     return(0);
132     unit2meter = 0.0254;
133     return(1);
134     }
135     return(0);
136     }
137    
138    
139     geomcheck(s) /* check/set lamp geometry */
140     char *s;
141     {
142     int len = strlen(s);
143    
144     switch (*s) {
145     case 'p':
146     if (strncmp(s, "polygonal", len))
147     return(0);
148     return(getpolygon());
149     case 's':
150     if (strncmp(s, "sphere", len) && strncmp(s, "spherical", len))
151     return(0);
152     return(getsphere());
153     case 'c':
154     if (strncmp(s,"cylinder",len) && strncmp(s,"cylindrical",len))
155     return(0);
156     return(getcylinder());
157     case 'r':
158     if (strncmp(s, "ring", len) && strncmp(s, "disk", len))
159     return(0);
160     return(getring());
161     }
162     return(0);
163     }
164    
165    
166     outpcheck(s) /* check lumen output value */
167     register char *s;
168     {
169     if ((*s < '0' || *s > '9') && *s != '.')
170     return(0);
171     lumens = atof(s);
172     return(1);
173     }
174    
175    
176     compute() /* compute lamp radiance */
177     {
178     double whiteval;
179    
180 greg 1.8 whiteval = lumens/area/(WHTEFFICACY*PI);
181 greg 1.2
182     printf("Lamp color (RGB) = %f %f %f\n",
183     lampcolor[0]*whiteval,
184     lampcolor[1]*whiteval,
185     lampcolor[2]*whiteval);
186     }
187    
188    
189 greg 1.6 getd(name, dp, help) /* get a positive double from stdin */
190     char *name;
191 greg 1.2 double *dp;
192 greg 1.6 char *help;
193 greg 1.2 {
194     char buf[32];
195 greg 1.6 again:
196     printf("%s [%g]: ", name, *dp);
197 greg 1.2 if (gets(buf) == NULL)
198     return(0);
199 greg 1.6 if (buf[0] == '?') {
200     puts(help);
201     goto again;
202     }
203 greg 1.2 if ((buf[0] < '0' || buf[0] > '9') && buf[0] != '.')
204     return(0);
205     *dp = atof(buf);
206     return(1);
207     }
208    
209    
210     getpolygon() /* get projected area for a polygon */
211     {
212 greg 1.9 static double parea = 1.0;
213 greg 1.5
214 greg 1.9 getd("Polygon area", &parea,
215 greg 1.6 "Enter the total radiating area of the polygon.");
216 greg 1.9 area = unit2meter*unit2meter * parea;
217 greg 1.2 return(1);
218     }
219    
220    
221     getsphere() /* get projected area for a sphere */
222     {
223 greg 1.5 static double radius = 1.0;
224 greg 1.2
225 greg 1.6 getd("Sphere radius", &radius,
226     "Enter the distance from the sphere's center to its surface.");
227 greg 1.7 area = 4.*PI*unit2meter*unit2meter * radius*radius;
228 greg 1.2 return(1);
229     }
230    
231    
232     getcylinder() /* get projected area for a cylinder */
233     {
234 greg 1.5 static double length = 1.0, radius = 0.1;
235 greg 1.2
236 greg 1.6 getd("Cylinder length", &length,
237     "Enter the length of the cylinder.");
238     getd("Cylinder radius", &radius,
239     "Enter the distance from the cylinder's axis to its surface.");
240 greg 1.7 area = 2.*PI*unit2meter*unit2meter * radius*length;
241 greg 1.2 return(1);
242     }
243    
244    
245     getring() /* get projected area for a ring */
246     {
247 greg 1.5 static double radius = 1.0;
248 greg 1.2
249 greg 1.6 getd("Disk radius", &radius,
250     "Enter the distance from the ring's center to its outer edge.\n\
251     The inner radius must be zero.");
252 greg 1.7 area = PI*unit2meter*unit2meter * radius*radius;
253 greg 1.2 return(1);
254 greg 1.1 }