ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pvalue.c
Revision: 2.18
Committed: Thu Jun 5 19:29:34 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.17: +9 -11 lines
Log Message:
Macros for setting binary file mode. Replacing MSDOS by _WIN32.

File Contents

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