ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pvalue.c
Revision: 2.21
Committed: Thu Jun 26 00:58:10 2003 UTC (20 years, 10 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.20: +2 -2 lines
Log Message:
Abstracted process and path handling for Windows.
Renamed FLOAT to RREAL because of conflict on Windows.
Added conditional compiles for some signal handlers.

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2 schorsch 2.21 static const char RCSid[] = "$Id: pvalue.c,v 2.20 2003/06/20 00:25:50 greg Exp $";
3 greg 1.1 #endif
4     /*
5     * pvalue.c - program to print pixel values.
6     *
7     * 4/23/86
8     */
9    
10 schorsch 2.18 #include <time.h>
11    
12 greg 1.8 #include "standard.h"
13 schorsch 2.19 #include "platform.h"
14 greg 1.1 #include "color.h"
15 greg 1.8 #include "resolu.h"
16 greg 2.15
17 greg 2.5 #define min(a,b) ((a)<(b)?(a):(b))
18 greg 1.1
19 greg 2.10 /* what to put out (also RED, GRN, BLU) */
20     #define ALL 3
21     #define BRIGHT 4
22    
23 greg 2.5 RESOLU picres; /* resolution of picture */
24 greg 1.1
25 greg 1.6 int uniq = 0; /* print only unique values? */
26 greg 1.1
27 greg 2.7 int doexposure = 0; /* exposure change? (>100 to print) */
28 greg 1.1
29     int dataonly = 0; /* data only format? */
30    
31 greg 2.10 int putprim = ALL; /* what to put out */
32 greg 1.1
33     int reverse = 0; /* reverse conversion? */
34    
35     int format = 'a'; /* input/output format */
36 greg 1.6 char *fmtid = "ascii"; /* format identifier for header */
37 greg 1.1
38 greg 2.3 int header = 1; /* do header? */
39 greg 1.1
40 greg 2.10 long skipbytes = 0; /* skip bytes in input? */
41    
42 greg 2.16 int swapbytes = 0; /* swap bytes in 16-bit words? */
43    
44 greg 2.12 int interleave = 1; /* file is interleaved? */
45    
46 greg 2.3 int resolution = 1; /* put/get resolution string? */
47    
48 greg 2.7 int original = 0; /* convert to original values? */
49    
50 greg 1.6 int wrongformat = 0; /* wrong input format? */
51    
52 greg 2.5 double gamcor = 1.0; /* gamma correction */
53 greg 1.7
54 greg 2.4 int ord[3] = {RED, GRN, BLU}; /* RGB ordering */
55     int rord[4]; /* reverse ordering */
56    
57 greg 1.5 COLOR exposure = WHTCOLOR;
58 greg 1.1
59     char *progname;
60    
61     FILE *fin;
62 greg 2.10 FILE *fin2 = NULL, *fin3 = NULL; /* for other color channels */
63 greg 1.1
64     int (*getval)(), (*putval)();
65    
66 greg 2.14 double
67     rgb_bright(clr)
68     COLOR clr;
69     {
70     return(bright(clr));
71     }
72    
73     double
74     xyz_bright(clr)
75     COLOR clr;
76     {
77     return(clr[CIEY]);
78     }
79    
80     double (*mybright)() = &rgb_bright;
81    
82 greg 1.1
83     main(argc, argv)
84     int argc;
85     char **argv;
86     {
87     extern int checkhead();
88 greg 2.10 extern long atol();
89 greg 2.7 double d, expval = 1.0;
90 greg 1.1 int i;
91    
92     progname = argv[0];
93    
94     for (i = 1; i < argc; i++)
95 greg 1.8 if (argv[i][0] == '-' || argv[i][0] == '+')
96 greg 1.1 switch (argv[i][1]) {
97 greg 1.8 case 'h': /* header */
98     header = argv[i][0] == '+';
99 greg 1.1 break;
100 greg 2.3 case 'H': /* resolution string */
101     resolution = argv[i][0] == '+';
102     break;
103 greg 2.10 case 's': /* skip bytes in header */
104     skipbytes = atol(argv[++i]);
105     break;
106 greg 1.1 case 'u': /* unique values */
107 greg 1.8 uniq = argv[i][0] == '-';
108 greg 1.1 break;
109     case 'o': /* original values */
110 greg 1.8 original = argv[i][0] == '-';
111 greg 1.1 break;
112 greg 1.7 case 'g': /* gamma correction */
113 greg 2.7 gamcor = atof(argv[i+1]);
114 greg 1.8 if (argv[i][0] == '+')
115     gamcor = 1.0/gamcor;
116 greg 2.7 i++;
117 greg 1.7 break;
118 greg 2.7 case 'e': /* exposure correction */
119     d = atof(argv[i+1]);
120     if (argv[i+1][0] == '-' || argv[i+1][0] == '+')
121     d = pow(2.0, d);
122     if (argv[i][0] == '-')
123 greg 2.9 expval *= d;
124 greg 2.7 scalecolor(exposure, d);
125     doexposure++;
126     i++;
127     break;
128 greg 2.4 case 'R': /* reverse byte sequence */
129     if (argv[i][0] == '-') {
130     ord[0]=BLU; ord[1]=GRN; ord[2]=RED;
131     } else {
132     ord[0]=RED; ord[1]=GRN; ord[2]=BLU;
133     }
134     break;
135 greg 1.1 case 'r': /* reverse conversion */
136 greg 1.8 reverse = argv[i][0] == '-';
137 greg 1.1 break;
138 greg 2.12 case 'n': /* non-interleaved RGB */
139     interleave = argv[i][0] == '+';
140     break;
141 greg 1.1 case 'b': /* brightness values */
142 greg 2.10 putprim = argv[i][0] == '-' ? BRIGHT : ALL;
143 greg 1.1 break;
144 greg 2.10 case 'p': /* put primary */
145     switch (argv[i][2]) {
146     case 'r': case 'R': putprim = RED; break;
147     case 'g': case 'G': putprim = GRN; break;
148     case 'b': case 'B': putprim = BLU; break;
149     default: goto unkopt;
150     }
151     break;
152 greg 1.1 case 'd': /* data only (no indices) */
153 greg 1.8 dataonly = argv[i][0] == '-';
154 greg 1.1 switch (argv[i][2]) {
155     case '\0':
156     case 'a': /* ascii */
157     format = 'a';
158 greg 1.6 fmtid = "ascii";
159 greg 1.1 break;
160     case 'i': /* integer */
161 greg 1.6 format = 'i';
162     fmtid = "ascii";
163     break;
164 greg 1.1 case 'b': /* byte */
165 greg 1.8 dataonly = 1;
166 greg 1.6 format = 'b';
167     fmtid = "byte";
168     break;
169 greg 2.16 case 'W': /* 16-bit swapped */
170     swapbytes = 1;
171 greg 2.15 case 'w': /* 16-bit */
172     dataonly = 1;
173     format = 'w';
174     fmtid = "16-bit";
175     break;
176 greg 1.1 case 'f': /* float */
177 greg 1.8 dataonly = 1;
178 greg 1.6 format = 'f';
179     fmtid = "float";
180     break;
181 greg 1.1 case 'd': /* double */
182 greg 1.8 dataonly = 1;
183 greg 1.6 format = 'd';
184     fmtid = "double";
185 greg 1.1 break;
186     default:
187     goto unkopt;
188     }
189     break;
190     case 'x': /* x resolution */
191 greg 2.6 case 'X': /* x resolution */
192 greg 2.3 resolution = 0;
193 greg 1.8 if (argv[i][0] == '-')
194 greg 2.14 picres.rt |= XDECR;
195 greg 1.8 picres.xr = atoi(argv[++i]);
196 greg 1.1 break;
197     case 'y': /* y resolution */
198 greg 2.6 case 'Y': /* y resolution */
199 greg 2.3 resolution = 0;
200 greg 1.8 if (argv[i][0] == '-')
201 greg 2.14 picres.rt |= YDECR;
202 greg 1.8 if (picres.xr == 0)
203 greg 2.14 picres.rt |= YMAJOR;
204 greg 1.8 picres.yr = atoi(argv[++i]);
205 greg 1.1 break;
206     default:
207     unkopt:
208     fprintf(stderr, "%s: unknown option: %s\n",
209     progname, argv[i]);
210     quit(1);
211     break;
212     }
213     else
214     break;
215 greg 1.6 /* recognize special formats */
216     if (dataonly && format == 'b')
217 greg 2.17 if (putprim == ALL)
218     fmtid = "24-bit_rgb";
219     else
220 greg 1.6 fmtid = "8-bit_grey";
221 greg 2.17 if (dataonly && format == 'w')
222     if (putprim == ALL)
223     fmtid = "48-bit_rgb";
224 greg 1.6 else
225 greg 2.15 fmtid = "16-bit_grey";
226 greg 2.4 /* assign reverse ordering */
227     rord[ord[0]] = 0;
228     rord[ord[1]] = 1;
229     rord[ord[2]] = 2;
230     /* get input */
231 greg 1.1 if (i == argc) {
232     fin = stdin;
233 greg 2.10 } else if (i < argc) {
234 greg 1.1 if ((fin = fopen(argv[i], "r")) == NULL) {
235     fprintf(stderr, "%s: can't open file \"%s\"\n",
236     progname, argv[i]);
237     quit(1);
238     }
239 greg 2.17 if (reverse && putprim != BRIGHT && i == argc-3) {
240 greg 2.10 if ((fin2 = fopen(argv[i+1], "r")) == NULL) {
241     fprintf(stderr, "%s: can't open file \"%s\"\n",
242     progname, argv[i+1]);
243     quit(1);
244     }
245     if ((fin3 = fopen(argv[i+2], "r")) == NULL) {
246     fprintf(stderr, "%s: can't open file \"%s\"\n",
247     progname, argv[i+2]);
248     quit(1);
249     }
250 greg 2.12 interleave = -1;
251 greg 2.10 } else if (i != argc-1)
252     fin = NULL;
253 greg 2.17 if (reverse && putprim != BRIGHT && !interleave) {
254 greg 2.12 fin2 = fopen(argv[i], "r");
255     fin3 = fopen(argv[i], "r");
256     }
257     if (skipbytes && (fseek(fin, skipbytes, 0) || (fin2 != NULL &&
258     (fseek(fin2, skipbytes, 0) ||
259     fseek(fin3, skipbytes, 0))))) {
260     fprintf(stderr, "%s: cannot skip %ld bytes on input\n",
261     progname, skipbytes);
262     quit(1);
263     }
264 greg 2.10 }
265     if (fin == NULL) {
266 greg 1.1 fprintf(stderr, "%s: bad # file arguments\n", progname);
267     quit(1);
268     }
269    
270     if (reverse) {
271 schorsch 2.18 #ifdef _WIN32
272     SET_FILE_BINARY(stdout);
273 greg 2.5 if (format != 'a' && format != 'i')
274 schorsch 2.18 SET_FILE_BINARY(fin);
275 greg 2.5 #endif
276 greg 1.6 /* get header */
277 greg 2.8 if (header) {
278     if (checkheader(fin, fmtid, stdout) < 0) {
279     fprintf(stderr, "%s: wrong input format\n",
280     progname);
281     quit(1);
282     }
283 greg 2.10 if (fin2 != NULL) {
284     getheader(fin2, NULL, NULL);
285     getheader(fin3, NULL, NULL);
286     }
287 greg 2.8 } else
288     newheader("RADIANCE", stdout);
289 greg 2.3 /* get resolution */
290     if ((resolution && !fgetsresolu(&picres, fin)) ||
291     picres.xr <= 0 || picres.yr <= 0) {
292     fprintf(stderr, "%s: missing resolution\n", progname);
293     quit(1);
294     }
295 greg 2.10 if (resolution && fin2 != NULL) {
296     RESOLU pres2;
297     if (!fgetsresolu(&pres2, fin2) ||
298 greg 2.14 pres2.rt != picres.rt ||
299 greg 2.10 pres2.xr != picres.xr ||
300     pres2.yr != picres.yr ||
301     !fgetsresolu(&pres2, fin3) ||
302 greg 2.14 pres2.rt != picres.rt ||
303 greg 2.10 pres2.xr != picres.xr ||
304     pres2.yr != picres.yr) {
305     fprintf(stderr, "%s: resolution mismatch\n",
306     progname);
307     quit(1);
308     }
309     }
310 greg 1.6 /* add to header */
311     printargs(i, argv, stdout);
312 greg 2.9 if (expval < .99 || expval > 1.01)
313 greg 2.7 fputexpos(expval, stdout);
314 greg 1.6 fputformat(COLRFMT, stdout);
315 greg 1.8 putchar('\n');
316 greg 2.3 fputsresolu(&picres, stdout); /* always put resolution */
317 greg 1.1 valtopix();
318     } else {
319 schorsch 2.18 #ifdef _WIN32
320     SET_FILE_BINARY(fin);
321 greg 2.5 if (format != 'a' && format != 'i')
322 schorsch 2.18 SET_FILE_BINARY(stdout);
323 greg 2.5 #endif
324 greg 1.1 /* get header */
325 greg 1.6 getheader(fin, checkhead, NULL);
326     if (wrongformat) {
327 greg 2.14 fprintf(stderr,
328     "%s: input not a Radiance RGBE picture\n",
329 greg 1.6 progname);
330     quit(1);
331     }
332 greg 2.3 if (!fgetsresolu(&picres, fin)) {
333     fprintf(stderr, "%s: missing resolution\n", progname);
334     quit(1);
335     }
336 greg 1.1 if (header) {
337     printargs(i, argv, stdout);
338 greg 2.9 if (expval < .99 || expval > 1.01)
339 greg 2.7 fputexpos(expval, stdout);
340 greg 1.6 fputformat(fmtid, stdout);
341 greg 1.8 putchar('\n');
342 greg 1.1 }
343 greg 2.3 if (resolution) /* put resolution */
344     fputsresolu(&picres, stdout);
345 greg 1.1 pixtoval();
346     }
347    
348     quit(0);
349     }
350    
351    
352 gwlarson 2.13 int
353 greg 1.1 checkhead(line) /* deal with line from header */
354     char *line;
355     {
356 greg 1.6 char fmt[32];
357 greg 1.5 double d;
358     COLOR ctmp;
359    
360 greg 2.14 if (formatval(fmt, line)) {
361     if (!strcmp(fmt, CIEFMT)) {
362     mybright = &xyz_bright;
363     if (original) {
364     scalecolor(exposure, 1./WHTEFFICACY);
365     doexposure++;
366     }
367     } else if (!strcmp(fmt, COLRFMT))
368     mybright = &rgb_bright;
369     else
370     wrongformat++;
371     } else if (original && isexpos(line)) {
372 greg 1.5 d = 1.0/exposval(line);
373     scalecolor(exposure, d);
374 greg 2.7 doexposure++;
375 greg 1.6 } else if (original && iscolcor(line)) {
376 greg 1.5 colcorval(ctmp, line);
377 greg 1.7 setcolor(exposure, colval(exposure,RED)/colval(ctmp,RED),
378     colval(exposure,GRN)/colval(ctmp,GRN),
379     colval(exposure,BLU)/colval(ctmp,BLU));
380 greg 2.7 doexposure++;
381 greg 1.6 } else if (header)
382     fputs(line, stdout);
383 gwlarson 2.13 return(0);
384 greg 1.1 }
385    
386    
387     pixtoval() /* convert picture to values */
388     {
389 greg 2.5 register COLOR *scanln;
390 greg 1.7 int dogamma;
391 greg 1.1 COLOR lastc;
392 schorsch 2.21 RREAL hv[2];
393 greg 2.12 int startprim, endprim;
394     long startpos;
395 greg 1.1 int y;
396     register int x;
397    
398 greg 1.8 scanln = (COLOR *)malloc(scanlen(&picres)*sizeof(COLOR));
399 greg 1.1 if (scanln == NULL) {
400     fprintf(stderr, "%s: out of memory\n", progname);
401     quit(1);
402     }
403 greg 1.8 dogamma = gamcor < .95 || gamcor > 1.05;
404 greg 2.12 if (putprim == ALL && !interleave) {
405     startprim = RED; endprim = BLU;
406     startpos = ftell(fin);
407     } else {
408     startprim = putprim; endprim = putprim;
409     }
410     for (putprim = startprim; putprim <= endprim; putprim++) {
411     if (putprim != startprim && fseek(fin, startpos, 0)) {
412     fprintf(stderr, "%s: seek error on input file\n",
413     progname);
414 greg 1.1 quit(1);
415     }
416 greg 2.12 set_io();
417     setcolor(lastc, 0.0, 0.0, 0.0);
418     for (y = 0; y < numscans(&picres); y++) {
419     if (freadscan(scanln, scanlen(&picres), fin) < 0) {
420     fprintf(stderr, "%s: read error\n", progname);
421 greg 1.1 quit(1);
422     }
423 greg 2.12 for (x = 0; x < scanlen(&picres); x++) {
424     if (uniq)
425     if ( colval(scanln[x],RED) ==
426     colval(lastc,RED) &&
427     colval(scanln[x],GRN) ==
428     colval(lastc,GRN) &&
429     colval(scanln[x],BLU) ==
430     colval(lastc,BLU) )
431     continue;
432     else
433     copycolor(lastc, scanln[x]);
434     if (doexposure)
435     multcolor(scanln[x], exposure);
436     if (dogamma)
437     setcolor(scanln[x],
438     pow(colval(scanln[x],RED), 1.0/gamcor),
439     pow(colval(scanln[x],GRN), 1.0/gamcor),
440     pow(colval(scanln[x],BLU), 1.0/gamcor));
441     if (!dataonly) {
442     pix2loc(hv, &picres, x, y);
443     printf("%7d %7d ",
444     (int)(hv[0]*picres.xr),
445     (int)(hv[1]*picres.yr));
446     }
447     if ((*putval)(scanln[x]) < 0) {
448     fprintf(stderr, "%s: write error\n",
449     progname);
450     quit(1);
451     }
452     }
453 greg 1.1 }
454     }
455 greg 2.14 free((void *)scanln);
456 greg 1.1 }
457    
458    
459     valtopix() /* convert values to a pixel file */
460     {
461 greg 1.7 int dogamma;
462 greg 2.5 register COLOR *scanln;
463 greg 1.1 int y;
464     register int x;
465    
466 greg 1.8 scanln = (COLOR *)malloc(scanlen(&picres)*sizeof(COLOR));
467 greg 1.1 if (scanln == NULL) {
468     fprintf(stderr, "%s: out of memory\n", progname);
469     quit(1);
470     }
471 greg 1.8 dogamma = gamcor < .95 || gamcor > 1.05;
472 greg 2.12 set_io();
473 greg 1.8 for (y = 0; y < numscans(&picres); y++) {
474     for (x = 0; x < scanlen(&picres); x++) {
475 greg 2.10 if (!dataonly) {
476 greg 1.1 fscanf(fin, "%*d %*d");
477 greg 2.10 if (fin2 != NULL) {
478     fscanf(fin2, "%*d %*d");
479     fscanf(fin3, "%*d %*d");
480     }
481     }
482 greg 2.11 if ((*getval)(scanln[x]) < 0) {
483 greg 1.1 fprintf(stderr, "%s: read error\n", progname);
484     quit(1);
485     }
486 greg 1.7 if (dogamma)
487     setcolor(scanln[x],
488 greg 1.8 pow(colval(scanln[x],RED), gamcor),
489     pow(colval(scanln[x],GRN), gamcor),
490     pow(colval(scanln[x],BLU), gamcor));
491 greg 2.7 if (doexposure)
492     multcolor(scanln[x], exposure);
493 greg 1.1 }
494 greg 1.8 if (fwritescan(scanln, scanlen(&picres), stdout) < 0) {
495 greg 1.1 fprintf(stderr, "%s: write error\n", progname);
496     quit(1);
497     }
498     }
499 greg 2.14 free((void *)scanln);
500 greg 1.1 }
501    
502    
503 greg 2.14 void
504 greg 1.1 quit(code)
505     int code;
506     {
507     exit(code);
508     }
509    
510    
511 greg 2.16 swap16(wp, n) /* swap n 16-bit words */
512     register uint16 *wp;
513     int n;
514     {
515     while (n-- > 0) {
516     *wp = *wp << 8 | ((*wp >> 8) & 0xff);
517     wp++;
518     }
519     }
520    
521 greg 2.11 getcascii(col) /* get an ascii color value from stream(s) */
522 greg 1.1 COLOR col;
523     {
524 greg 2.5 double vd[3];
525 greg 1.1
526 greg 2.11 if (fin2 == NULL) {
527     if (fscanf(fin, "%lf %lf %lf", &vd[0], &vd[1], &vd[2]) != 3)
528 greg 2.10 return(-1);
529     } else {
530 greg 2.11 if (fscanf(fin, "%lf", &vd[0]) != 1 ||
531     fscanf(fin2, "%lf", &vd[1]) != 1 ||
532     fscanf(fin3, "%lf", &vd[2]) != 1)
533 greg 2.10 return(-1);
534     }
535 greg 2.4 setcolor(col, vd[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
536 greg 1.1 return(0);
537     }
538    
539    
540 greg 2.11 getcdouble(col) /* get a double color value from stream(s) */
541 greg 1.1 COLOR col;
542     {
543 greg 2.5 double vd[3];
544 greg 1.1
545 greg 2.11 if (fin2 == NULL) {
546     if (fread((char *)vd, sizeof(double), 3, fin) != 3)
547 greg 2.10 return(-1);
548     } else {
549 greg 2.11 if (fread((char *)vd, sizeof(double), 1, fin) != 1 ||
550     fread((char *)(vd+1), sizeof(double), 1, fin2) != 1 ||
551     fread((char *)(vd+2), sizeof(double), 1, fin3) != 1)
552 greg 2.10 return(-1);
553     }
554 greg 2.4 setcolor(col, vd[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
555 greg 1.1 return(0);
556     }
557    
558    
559 greg 2.11 getcfloat(col) /* get a float color value from stream(s) */
560 greg 1.1 COLOR col;
561     {
562     float vf[3];
563    
564 greg 2.11 if (fin2 == NULL) {
565     if (fread((char *)vf, sizeof(float), 3, fin) != 3)
566 greg 2.10 return(-1);
567     } else {
568 greg 2.11 if (fread((char *)vf, sizeof(float), 1, fin) != 1 ||
569     fread((char *)(vf+1), sizeof(float), 1, fin2) != 1 ||
570     fread((char *)(vf+2), sizeof(float), 1, fin3) != 1)
571 greg 2.10 return(-1);
572     }
573 greg 2.4 setcolor(col, vf[rord[RED]], vf[rord[GRN]], vf[rord[BLU]]);
574 greg 1.1 return(0);
575     }
576    
577    
578 greg 2.11 getcint(col) /* get an int color value from stream(s) */
579 greg 1.1 COLOR col;
580     {
581     int vi[3];
582    
583 greg 2.11 if (fin2 == NULL) {
584     if (fscanf(fin, "%d %d %d", &vi[0], &vi[1], &vi[2]) != 3)
585 greg 2.10 return(-1);
586     } else {
587 greg 2.11 if (fscanf(fin, "%d", &vi[0]) != 1 ||
588     fscanf(fin2, "%d", &vi[1]) != 1 ||
589     fscanf(fin3, "%d", &vi[2]) != 1)
590 greg 2.10 return(-1);
591     }
592 greg 2.4 setcolor(col, (vi[rord[RED]]+.5)/256.,
593     (vi[rord[GRN]]+.5)/256., (vi[rord[BLU]]+.5)/256.);
594 greg 1.1 return(0);
595     }
596    
597    
598 greg 2.11 getcbyte(col) /* get a byte color value from stream(s) */
599 greg 1.1 COLOR col;
600     {
601     BYTE vb[3];
602    
603 greg 2.11 if (fin2 == NULL) {
604     if (fread((char *)vb, sizeof(BYTE), 3, fin) != 3)
605 greg 2.10 return(-1);
606     } else {
607 greg 2.11 if (fread((char *)vb, sizeof(BYTE), 1, fin) != 1 ||
608     fread((char *)(vb+1), sizeof(BYTE), 1, fin2) != 1 ||
609     fread((char *)(vb+2), sizeof(BYTE), 1, fin3) != 1)
610 greg 2.10 return(-1);
611     }
612 greg 2.4 setcolor(col, (vb[rord[RED]]+.5)/256.,
613     (vb[rord[GRN]]+.5)/256., (vb[rord[BLU]]+.5)/256.);
614 greg 1.1 return(0);
615     }
616    
617    
618 greg 2.15 getcword(col) /* get a 16-bit color value from stream(s) */
619     COLOR col;
620     {
621     uint16 vw[3];
622    
623     if (fin2 == NULL) {
624     if (fread((char *)vw, sizeof(uint16), 3, fin) != 3)
625     return(-1);
626     } else {
627     if (fread((char *)vw, sizeof(uint16), 1, fin) != 1 ||
628     fread((char *)(vw+1), sizeof(uint16), 1, fin2) != 1 ||
629     fread((char *)(vw+2), sizeof(uint16), 1, fin3) != 1)
630     return(-1);
631     }
632 greg 2.16 if (swapbytes)
633     swap16(vw, 3);
634 greg 2.15 setcolor(col, (vw[rord[RED]]+.5)/65536.,
635     (vw[rord[GRN]]+.5)/65536., (vw[rord[BLU]]+.5)/65536.);
636     return(0);
637     }
638    
639    
640 greg 2.11 getbascii(col) /* get an ascii brightness value from fin */
641 greg 1.1 COLOR col;
642     {
643 greg 2.5 double vd;
644 greg 1.1
645 greg 2.11 if (fscanf(fin, "%lf", &vd) != 1)
646 greg 1.1 return(-1);
647     setcolor(col, vd, vd, vd);
648     return(0);
649     }
650    
651    
652 greg 2.11 getbdouble(col) /* get a double brightness value from fin */
653 greg 1.1 COLOR col;
654     {
655 greg 2.5 double vd;
656 greg 1.1
657 greg 2.11 if (fread((char *)&vd, sizeof(double), 1, fin) != 1)
658 greg 1.1 return(-1);
659     setcolor(col, vd, vd, vd);
660     return(0);
661     }
662    
663    
664 greg 2.11 getbfloat(col) /* get a float brightness value from fin */
665 greg 1.1 COLOR col;
666     {
667     float vf;
668    
669 greg 2.11 if (fread((char *)&vf, sizeof(float), 1, fin) != 1)
670 greg 1.1 return(-1);
671     setcolor(col, vf, vf, vf);
672     return(0);
673     }
674    
675    
676 greg 2.11 getbint(col) /* get an int brightness value from fin */
677 greg 1.1 COLOR col;
678     {
679     int vi;
680 greg 2.5 double d;
681 greg 1.1
682 greg 2.11 if (fscanf(fin, "%d", &vi) != 1)
683 greg 1.1 return(-1);
684     d = (vi+.5)/256.;
685     setcolor(col, d, d, d);
686     return(0);
687     }
688    
689    
690 greg 2.11 getbbyte(col) /* get a byte brightness value from fin */
691 greg 1.1 COLOR col;
692     {
693     BYTE vb;
694 greg 2.5 double d;
695 greg 1.1
696 greg 2.11 if (fread((char *)&vb, sizeof(BYTE), 1, fin) != 1)
697 greg 1.1 return(-1);
698     d = (vb+.5)/256.;
699     setcolor(col, d, d, d);
700     return(0);
701     }
702    
703    
704 greg 2.15 getbword(col) /* get a 16-bit brightness value from fin */
705     COLOR col;
706     {
707     uint16 vw;
708     double d;
709    
710     if (fread((char *)&vw, sizeof(uint16), 1, fin) != 1)
711     return(-1);
712 greg 2.16 if (swapbytes)
713     swap16(&vw, 1);
714 greg 2.15 d = (vw+.5)/65536.;
715     setcolor(col, d, d, d);
716     return(0);
717     }
718    
719    
720 greg 2.11 putcascii(col) /* put an ascii color to stdout */
721 greg 1.1 COLOR col;
722     {
723 greg 2.11 fprintf(stdout, "%15.3e %15.3e %15.3e\n",
724 greg 2.4 colval(col,ord[0]),
725     colval(col,ord[1]),
726     colval(col,ord[2]));
727 greg 1.1
728 greg 2.11 return(ferror(stdout) ? -1 : 0);
729 greg 1.1 }
730    
731    
732 greg 2.11 putcfloat(col) /* put a float color to stdout */
733 greg 1.1 COLOR col;
734     {
735     float vf[3];
736    
737 greg 2.4 vf[0] = colval(col,ord[0]);
738     vf[1] = colval(col,ord[1]);
739     vf[2] = colval(col,ord[2]);
740 greg 2.11 fwrite((char *)vf, sizeof(float), 3, stdout);
741 greg 1.1
742 greg 2.11 return(ferror(stdout) ? -1 : 0);
743 greg 1.1 }
744    
745    
746 greg 2.11 putcdouble(col) /* put a double color to stdout */
747 greg 1.1 COLOR col;
748     {
749 greg 2.5 double vd[3];
750 greg 1.1
751 greg 2.4 vd[0] = colval(col,ord[0]);
752     vd[1] = colval(col,ord[1]);
753     vd[2] = colval(col,ord[2]);
754 greg 2.11 fwrite((char *)vd, sizeof(double), 3, stdout);
755 greg 1.1
756 greg 2.11 return(ferror(stdout) ? -1 : 0);
757 greg 1.1 }
758    
759    
760 greg 2.11 putcint(col) /* put an int color to stdout */
761 greg 1.1 COLOR col;
762     {
763 greg 2.11 fprintf(stdout, "%d %d %d\n",
764 greg 2.4 (int)(colval(col,ord[0])*256.),
765     (int)(colval(col,ord[1])*256.),
766     (int)(colval(col,ord[2])*256.));
767 greg 1.1
768 greg 2.11 return(ferror(stdout) ? -1 : 0);
769 greg 1.1 }
770    
771    
772 greg 2.11 putcbyte(col) /* put a byte color to stdout */
773 greg 1.1 COLOR col;
774     {
775 greg 2.15 long i;
776 greg 1.1 BYTE vb[3];
777    
778 greg 2.4 i = colval(col,ord[0])*256.;
779 greg 1.1 vb[0] = min(i,255);
780 greg 2.4 i = colval(col,ord[1])*256.;
781 greg 1.1 vb[1] = min(i,255);
782 greg 2.4 i = colval(col,ord[2])*256.;
783 greg 1.1 vb[2] = min(i,255);
784 greg 2.11 fwrite((char *)vb, sizeof(BYTE), 3, stdout);
785 greg 1.1
786 greg 2.11 return(ferror(stdout) ? -1 : 0);
787 greg 1.1 }
788    
789    
790 greg 2.15 putcword(col) /* put a 16-bit color to stdout */
791     COLOR col;
792     {
793     long i;
794     uint16 vw[3];
795    
796     i = colval(col,ord[0])*65536.;
797     vw[0] = min(i,65535);
798     i = colval(col,ord[1])*65536.;
799     vw[1] = min(i,65535);
800     i = colval(col,ord[2])*65536.;
801     vw[2] = min(i,65535);
802 greg 2.16 if (swapbytes)
803     swap16(vw, 3);
804 greg 2.15 fwrite((char *)vw, sizeof(uint16), 3, stdout);
805    
806     return(ferror(stdout) ? -1 : 0);
807     }
808    
809    
810 greg 2.11 putbascii(col) /* put an ascii brightness to stdout */
811 greg 1.1 COLOR col;
812     {
813 greg 2.14 fprintf(stdout, "%15.3e\n", (*mybright)(col));
814 greg 1.1
815 greg 2.11 return(ferror(stdout) ? -1 : 0);
816 greg 1.1 }
817    
818    
819 greg 2.11 putbfloat(col) /* put a float brightness to stdout */
820 greg 1.1 COLOR col;
821     {
822     float vf;
823    
824 greg 2.14 vf = (*mybright)(col);
825 greg 2.11 fwrite((char *)&vf, sizeof(float), 1, stdout);
826 greg 1.1
827 greg 2.11 return(ferror(stdout) ? -1 : 0);
828 greg 1.1 }
829    
830    
831 greg 2.11 putbdouble(col) /* put a double brightness to stdout */
832 greg 1.1 COLOR col;
833     {
834 greg 2.5 double vd;
835 greg 1.1
836 greg 2.14 vd = (*mybright)(col);
837 greg 2.11 fwrite((char *)&vd, sizeof(double), 1, stdout);
838 greg 1.1
839 greg 2.11 return(ferror(stdout) ? -1 : 0);
840 greg 1.1 }
841    
842    
843 greg 2.11 putbint(col) /* put an int brightness to stdout */
844 greg 1.1 COLOR col;
845     {
846 greg 2.14 fprintf(stdout, "%d\n", (int)((*mybright)(col)*256.));
847 greg 1.1
848 greg 2.11 return(ferror(stdout) ? -1 : 0);
849 greg 1.1 }
850    
851    
852 greg 2.11 putbbyte(col) /* put a byte brightness to stdout */
853 greg 1.1 COLOR col;
854     {
855     register int i;
856     BYTE vb;
857    
858 greg 2.14 i = (*mybright)(col)*256.;
859 greg 1.1 vb = min(i,255);
860 greg 2.11 fwrite((char *)&vb, sizeof(BYTE), 1, stdout);
861 greg 1.1
862 greg 2.11 return(ferror(stdout) ? -1 : 0);
863 greg 1.1 }
864    
865    
866 greg 2.15 putbword(col) /* put a 16-bit brightness to stdout */
867     COLOR col;
868     {
869     long i;
870     uint16 vw;
871    
872     i = (*mybright)(col)*65536.;
873     vw = min(i,65535);
874 greg 2.16 if (swapbytes)
875     swap16(&vw, 1);
876 greg 2.15 fwrite((char *)&vw, sizeof(uint16), 1, stdout);
877    
878     return(ferror(stdout) ? -1 : 0);
879     }
880    
881    
882 greg 2.11 putpascii(col) /* put an ascii primary to stdout */
883 greg 2.10 COLOR col;
884     {
885 greg 2.11 fprintf(stdout, "%15.3e\n", colval(col,putprim));
886 greg 2.10
887 greg 2.11 return(ferror(stdout) ? -1 : 0);
888 greg 2.10 }
889    
890    
891 greg 2.11 putpfloat(col) /* put a float primary to stdout */
892 greg 2.10 COLOR col;
893     {
894     float vf;
895    
896     vf = colval(col,putprim);
897 greg 2.11 fwrite((char *)&vf, sizeof(float), 1, stdout);
898 greg 2.10
899 greg 2.11 return(ferror(stdout) ? -1 : 0);
900 greg 2.10 }
901    
902    
903 greg 2.11 putpdouble(col) /* put a double primary to stdout */
904 greg 2.10 COLOR col;
905     {
906     double vd;
907    
908     vd = colval(col,putprim);
909 greg 2.11 fwrite((char *)&vd, sizeof(double), 1, stdout);
910 greg 2.10
911 greg 2.11 return(ferror(stdout) ? -1 : 0);
912 greg 2.10 }
913    
914    
915 greg 2.11 putpint(col) /* put an int primary to stdout */
916 greg 2.10 COLOR col;
917     {
918 greg 2.11 fprintf(stdout, "%d\n", (int)(colval(col,putprim)*256.));
919 greg 2.10
920 greg 2.11 return(ferror(stdout) ? -1 : 0);
921 greg 2.10 }
922    
923    
924 greg 2.11 putpbyte(col) /* put a byte primary to stdout */
925 greg 2.10 COLOR col;
926     {
927 greg 2.15 long i;
928 greg 2.10 BYTE vb;
929    
930     i = colval(col,putprim)*256.;
931     vb = min(i,255);
932 greg 2.11 fwrite((char *)&vb, sizeof(BYTE), 1, stdout);
933 greg 2.10
934 greg 2.11 return(ferror(stdout) ? -1 : 0);
935 greg 2.10 }
936    
937    
938 greg 2.15 putpword(col) /* put a 16-bit primary to stdout */
939     COLOR col;
940     {
941     long i;
942     uint16 vw;
943    
944     i = colval(col,putprim)*65536.;
945     vw = min(i,65535);
946 greg 2.16 if (swapbytes)
947     swap16(&vw, 1);
948 greg 2.15 fwrite((char *)&vw, sizeof(uint16), 1, stdout);
949    
950     return(ferror(stdout) ? -1 : 0);
951     }
952    
953    
954 greg 1.1 set_io() /* set put and get functions */
955     {
956     switch (format) {
957     case 'a': /* ascii */
958 greg 2.10 if (putprim == BRIGHT) {
959 greg 1.1 getval = getbascii;
960     putval = putbascii;
961 greg 2.10 } else if (putprim != ALL) {
962     getval = getbascii;
963     putval = putpascii;
964 greg 1.1 } else {
965     getval = getcascii;
966     putval = putcascii;
967 greg 2.12 if (reverse && !interleave) {
968     fprintf(stderr,
969     "%s: ASCII input files must be interleaved\n",
970     progname);
971     quit(1);
972     }
973 greg 1.1 }
974     return;
975     case 'f': /* binary float */
976 greg 2.10 if (putprim == BRIGHT) {
977 greg 1.1 getval = getbfloat;
978     putval = putbfloat;
979 greg 2.10 } else if (putprim != ALL) {
980     getval = getbfloat;
981     putval = putpfloat;
982 greg 1.1 } else {
983     getval = getcfloat;
984     putval = putcfloat;
985 greg 2.12 if (reverse && !interleave) {
986     if (fin2 == NULL)
987     goto namerr;
988     if (fseek(fin2,
989     (long)sizeof(float)*picres.xr*picres.yr, 1))
990     goto seekerr;
991     if (fseek(fin3,
992     (long)sizeof(float)*2*picres.xr*picres.yr, 1))
993     goto seekerr;
994     }
995 greg 1.1 }
996     return;
997     case 'd': /* binary double */
998 greg 2.10 if (putprim == BRIGHT) {
999 greg 1.1 getval = getbdouble;
1000     putval = putbdouble;
1001 greg 2.10 } else if (putprim != ALL) {
1002     getval = getbdouble;
1003     putval = putpdouble;
1004 greg 1.1 } else {
1005     getval = getcdouble;
1006     putval = putcdouble;
1007 greg 2.12 if (reverse && !interleave) {
1008     if (fin2 == NULL)
1009     goto namerr;
1010     if (fseek(fin2,
1011     (long)sizeof(double)*picres.xr*picres.yr, 1))
1012     goto seekerr;
1013     if (fseek(fin3,
1014     (long)sizeof(double)*2*picres.xr*picres.yr, 1))
1015     goto seekerr;
1016     }
1017 greg 1.1 }
1018     return;
1019     case 'i': /* integer */
1020 greg 2.10 if (putprim == BRIGHT) {
1021 greg 1.1 getval = getbint;
1022     putval = putbint;
1023 greg 2.10 } else if (putprim != ALL) {
1024     getval = getbint;
1025     putval = putpint;
1026 greg 1.1 } else {
1027     getval = getcint;
1028     putval = putcint;
1029 greg 2.12 if (reverse && !interleave) {
1030     fprintf(stderr,
1031     "%s: integer input files must be interleaved\n",
1032     progname);
1033     quit(1);
1034     }
1035 greg 1.1 }
1036     return;
1037     case 'b': /* byte */
1038 greg 2.10 if (putprim == BRIGHT) {
1039 greg 1.1 getval = getbbyte;
1040     putval = putbbyte;
1041 greg 2.10 } else if (putprim != ALL) {
1042     getval = getbbyte;
1043     putval = putpbyte;
1044 greg 1.1 } else {
1045     getval = getcbyte;
1046     putval = putcbyte;
1047 greg 2.12 if (reverse && !interleave) {
1048     if (fin2 == NULL)
1049     goto namerr;
1050     if (fseek(fin2,
1051     (long)sizeof(BYTE)*picres.xr*picres.yr, 1))
1052     goto seekerr;
1053     if (fseek(fin3,
1054     (long)sizeof(BYTE)*2*picres.xr*picres.yr, 1))
1055 greg 2.15 goto seekerr;
1056     }
1057     }
1058     return;
1059     case 'w': /* 16-bit */
1060     if (putprim == BRIGHT) {
1061     getval = getbword;
1062     putval = putbword;
1063     } else if (putprim != ALL) {
1064     getval = getbword;
1065     putval = putpword;
1066     } else {
1067     getval = getcword;
1068     putval = putcword;
1069     if (reverse && !interleave) {
1070     if (fin2 == NULL)
1071     goto namerr;
1072     if (fseek(fin2,
1073     (long)sizeof(uint16)*picres.xr*picres.yr, 1))
1074     goto seekerr;
1075     if (fseek(fin3,
1076     (long)sizeof(uint16)*2*picres.xr*picres.yr, 1))
1077 greg 2.12 goto seekerr;
1078     }
1079 greg 1.1 }
1080     return;
1081     }
1082 greg 2.12 badopt:
1083     fprintf(stderr, "%s: botched file type\n", progname);
1084     quit(1);
1085     namerr:
1086     fprintf(stderr, "%s: non-interleaved file(s) must be named\n",
1087     progname);
1088     quit(1);
1089     seekerr:
1090     fprintf(stderr, "%s: cannot seek on interleaved input file\n",
1091     progname);
1092     quit(1);
1093 greg 1.1 }