ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pvalue.c
(Generate patch)

Comparing ray/src/px/pvalue.c (file contents):
Revision 2.8 by greg, Sun Feb 27 10:16:48 1994 UTC vs.
Revision 2.16 by greg, Tue Feb 25 16:47:23 2003 UTC

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

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines