ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pvalue.c
Revision: 2.38
Committed: Thu Feb 6 19:33:32 2020 UTC (4 years, 3 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.37: +15 -10 lines
Log Message:
Fixed issue with BigEndian handling in header

File Contents

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