src/cv/lampcolor.c

#ifndef lint
static const char       RCSid[] = "$Id: lampcolor.c,v 2.7 2003/06/08 12:03:09 schorsch Exp $";
#endif
/*
 * Program to convert lamp color from table and compute radiance.
 */

#include <stdio.h>
#include <string.h>
#include <math.h>

#include "standard.h"
#include "color.h"

                                /* lamp parameters */
#define LTYPE           0
#define LUNIT           1
#define LGEOM           2
#define LOUTP           3
#define NPARAMS         4

static int typecheck(char *s);
static int unitcheck(char *s);
static int geomcheck(char *s);
static int outpcheck(char *s);
static void compute(void);
static int getpolygon(void), getsphere(void), getcylinder(void), getring(void);


float   *lampcolor;             /* the lamp color (RGB) */
double  unit2meter;             /* conversion from units to meters */
double  area;                   /* radiating area for this geometry */
double  lumens;                 /* total lamp lumens */

struct {
        char    *name;
        char    value[64];
        int     (*check)();
        char    *help;
} param[NPARAMS] = {
        { "lamp type", "WHITE", typecheck,
"The lamp type is a string which corresponds to one of the types registered\n\
in the lamp table file.  A value of \"WHITE\" means an uncolored source,\n\
which may be preferable because it results in a color balanced image." },
        { "length unit", "meter", unitcheck,
"Unit must be one of:  \"meter\", \"centimeter\", \"foot\", or \"inch\".\n\
These may be abbreviated as a single letter." },
        { "lamp geometry", "polygon", geomcheck,
"The lamp geometry must be one of:  \"polygon\", \"sphere\", \"cylinder\"\n\
or \"ring\".  These may be abbreviated as a single letter." },
        { "total lamp lumens", "0", outpcheck,
"This is the overall light output of the lamp and its fixture.  If you do\n\
not know this value explicitly, you can compute the approximate lumens\n\
by multiplying the input wattage by 14 for incandescent fixtures or 70\n\
for fluorescent fixtures." },
};


main(argc, argv)
int     argc;
char    *argv[];
{
        char    *lamptab = "lamp.tab";
        char    buf[64];
        int     i;

        if (argc > 1) lamptab = argv[1];
        if (loadlamps(lamptab) == 0) {
                fprintf(stderr, "%s: no such lamp table\n", lamptab);
                exit(1);
        }
        printf("Program to compute lamp radiance.  Enter '?' for help.\n");
        for ( ; ; ) {
                i = 0;
                while (i < NPARAMS) {
                        printf("Enter %s [%s]: ", param[i].name,
                                        param[i].value);
                        if (gets(buf) == NULL)
                                exit(0);
                        if (buf[0] == '?') {
                                puts(param[i].help);
                                continue;
                        }
                        if (buf[0])
                                strcpy(param[i].value, buf);
                        if (!(*param[i].check)(param[i].value)) {
                                fprintf(stderr, "%s: bad value for %s\n",
                                                param[i].value, param[i].name);
                                continue;
                        }
                        i++;
                }
                compute();
        }
}


static int
typecheck(                      /* check lamp type */
char    *s
)
{
        lampcolor = matchlamp(s);
        return(lampcolor != NULL);
}


static int
unitcheck(                      /* compute conversion to meters */
char    *s
)
{
        int     len = strlen(s);

        switch (*s) {
        case 'm':
                if (strncmp(s, "meters", len))
                        return(0);
                unit2meter = 1.0;
                return(1);
        case 'c':
                if (strncmp(s, "centimeters", len) && strncmp(s, "cms", len))
                        return(0);
                unit2meter = 0.01;
                return(1);
        case 'f':
                if (strncmp(s, "foot", len) && strncmp(s, "feet", len))
                        return(0);
                unit2meter = 0.3048;
                return(1);
        case 'i':
                if (strncmp(s, "inches", len))
                        return(0);
                unit2meter = 0.0254;
                return(1);
        }
        return(0);
}


static int
geomcheck(                      /* check/set lamp geometry */
char    *s
)
{
        int     len = strlen(s);

        switch (*s) {
        case 'p':
                if (strncmp(s, "polygonal", len))
                        return(0);
                return(getpolygon());
        case 's':
                if (strncmp(s, "sphere", len) && strncmp(s, "spherical", len))
                        return(0);
                return(getsphere());
        case 'c':
                if (strncmp(s,"cylinder",len) && strncmp(s,"cylindrical",len))
                        return(0);
                return(getcylinder());
        case 'r':
                if (strncmp(s, "ring", len) && strncmp(s, "disk", len))
                        return(0);
                return(getring());
        }
        return(0);
}


static int
outpcheck(                      /* check lumen output value */
register char   *s
)
{
        if ((*s < '0' || *s > '9') && *s != '.')
                return(0);
        lumens = atof(s);
        return(1);
}


static void
compute(void)                   /* compute lamp radiance */
{
        double  whiteval;

        whiteval = lumens/area/(WHTEFFICACY*PI);

        printf("Lamp color (RGB) = %f %f %f\n",
                        lampcolor[0]*whiteval,
                        lampcolor[1]*whiteval,
                        lampcolor[2]*whiteval);
}


getd(name, dp, help)            /* get a positive double from stdin */
char    *name;
double  *dp;
char    *help;
{
        char    buf[32];
again:
        printf("%s [%g]: ", name, *dp);
        if (gets(buf) == NULL)
                return(0);
        if (buf[0] == '?') {
                puts(help);
                goto again;
        }
        if ((buf[0] < '0' || buf[0] > '9') && buf[0] != '.')
                return(0);
        *dp = atof(buf);
        return(1);
}


static int
getpolygon(void)                        /* get projected area for a polygon */
{
        static double   parea = 1.0;

        getd("Polygon area", &parea,
                "Enter the total radiating area of the polygon.");
        area = unit2meter*unit2meter * parea;
        return(1);
}


static int
getsphere(void)                 /* get projected area for a sphere */
{
        static double   radius = 1.0;

        getd("Sphere radius", &radius,
                "Enter the distance from the sphere's center to its surface.");
        area = 4.*PI*unit2meter*unit2meter * radius*radius;
        return(1);
}


static int
getcylinder(void)                       /* get projected area for a cylinder */
{
        static double   length = 1.0, radius = 0.1;

        getd("Cylinder length", &length,
                "Enter the length of the cylinder.");
        getd("Cylinder radius", &radius,
                "Enter the distance from the cylinder's axis to its surface.");
        area = 2.*PI*unit2meter*unit2meter * radius*length;
        return(1);
}


static int
getring(void)                   /* get projected area for a ring */
{
        static double   radius = 1.0;

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