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 1.2 by greg, Tue Sep 12 13:04:30 1989 UTC vs.
Revision 2.7 by greg, Tue Nov 9 15:21:00 1993 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1986 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 10 | Line 10 | static char SCCSid[] = "$SunId$ LBL";
10   *     4/23/86
11   */
12  
13 < #include  <stdio.h>
13 > #include  "standard.h"
14  
15   #include  "color.h"
16  
17 < #define  min(a,b)               ((a)<(b)?(a):(b))
17 > #include  "resolu.h"
18  
19 < int  xres = 0;                  /* resolution of input */
20 < int  yres = 0;
19 > #define  min(a,b)               ((a)<(b)?(a):(b))
20  
21 < int  uniq = 0;                  /* unique values? */
21 > RESOLU  picres;                 /* resolution of picture */
22  
23 < int  original = 0;              /* original values? */
23 > int  uniq = 0;                  /* print only unique values? */
24  
25 + int  doexposure = 0;            /* exposure change? (>100 to print) */
26 +
27   int  dataonly = 0;              /* data only format? */
28  
29   int  brightonly = 0;            /* only brightness values? */
# Line 30 | Line 31 | int  brightonly = 0;           /* only brightness values? */
31   int  reverse = 0;               /* reverse conversion? */
32  
33   int  format = 'a';              /* input/output format */
34 + char  *fmtid = "ascii";         /* format identifier for header */
35  
36 < int  header = 1;                /* do header */
36 > int  header = 1;                /* do header? */
37  
38 < double  exposure = 1.0;
38 > int  resolution = 1;            /* put/get resolution string? */
39  
40 + int  original = 0;              /* convert to original values? */
41 +
42 + int  wrongformat = 0;           /* wrong input format? */
43 +
44 + double  gamcor = 1.0;           /* gamma correction */
45 +
46 + int  ord[3] = {RED, GRN, BLU};  /* RGB ordering */
47 + int  rord[4];                   /* reverse ordering */
48 +
49 + COLOR  exposure = WHTCOLOR;
50 +
51   char  *progname;
52  
53   FILE  *fin;
# Line 46 | Line 59 | main(argc, argv)
59   int  argc;
60   char  **argv;
61   {
49 –        extern double  atof();
62          extern int  checkhead();
63 +        double  d, expval = 1.0;
64          int  i;
65  
66          progname = argv[0];
67  
68          for (i = 1; i < argc; i++)
69 <                if (argv[i][0] == '-')
69 >                if (argv[i][0] == '-' || argv[i][0] == '+')
70                          switch (argv[i][1]) {
71 <                        case 'h':               /* no header */
72 <                                header = 0;
71 >                        case 'h':               /* header */
72 >                                header = argv[i][0] == '+';
73                                  break;
74 +                        case 'H':               /* resolution string */
75 +                                resolution = argv[i][0] == '+';
76 +                                break;
77                          case 'u':               /* unique values */
78 <                                uniq = 1;
78 >                                uniq = argv[i][0] == '-';
79                                  break;
80                          case 'o':               /* original values */
81 <                                original = 1;
81 >                                original = argv[i][0] == '-';
82                                  break;
83 +                        case 'g':               /* gamma correction */
84 +                                gamcor = atof(argv[i+1]);
85 +                                if (argv[i][0] == '+')
86 +                                        gamcor = 1.0/gamcor;
87 +                                i++;
88 +                                break;
89 +                        case 'e':               /* exposure correction */
90 +                                d = atof(argv[i+1]);
91 +                                if (argv[i+1][0] == '-' || argv[i+1][0] == '+')
92 +                                        d = pow(2.0, d);
93 +                                if (argv[i][0] == '-')
94 +                                        doexposure = 100;
95 +                                scalecolor(exposure, d);
96 +                                expval *= d;
97 +                                doexposure++;
98 +                                i++;
99 +                                break;
100 +                        case 'R':               /* reverse byte sequence */
101 +                                if (argv[i][0] == '-') {
102 +                                        ord[0]=BLU; ord[1]=GRN; ord[2]=RED;
103 +                                } else {
104 +                                        ord[0]=RED; ord[1]=GRN; ord[2]=BLU;
105 +                                }
106 +                                break;
107                          case 'r':               /* reverse conversion */
108 <                                reverse = 1;
108 >                                reverse = argv[i][0] == '-';
109                                  break;
110                          case 'b':               /* brightness values */
111 <                                brightonly = 1;
111 >                                brightonly = argv[i][0] == '-';
112                                  break;
113                          case 'd':               /* data only (no indices) */
114 <                                dataonly = 1;
114 >                                dataonly = argv[i][0] == '-';
115                                  switch (argv[i][2]) {
116                                  case '\0':
117                                  case 'a':               /* ascii */
118                                          format = 'a';
119 +                                        fmtid = "ascii";
120                                          break;
121                                  case 'i':               /* integer */
122 +                                        format = 'i';
123 +                                        fmtid = "ascii";
124 +                                        break;
125                                  case 'b':               /* byte */
126 +                                        dataonly = 1;
127 +                                        format = 'b';
128 +                                        fmtid = "byte";
129 +                                        break;
130                                  case 'f':               /* float */
131 +                                        dataonly = 1;
132 +                                        format = 'f';
133 +                                        fmtid = "float";
134 +                                        break;
135                                  case 'd':               /* double */
136 <                                        format = argv[i][2];
136 >                                        dataonly = 1;
137 >                                        format = 'd';
138 >                                        fmtid = "double";
139                                          break;
140                                  default:
141                                          goto unkopt;
142                                  }
143                                  break;
144                          case 'x':               /* x resolution */
145 <                                xres = atoi(argv[++i]);
145 >                        case 'X':               /* x resolution */
146 >                                resolution = 0;
147 >                                if (argv[i][0] == '-')
148 >                                        picres.or |= XDECR;
149 >                                picres.xr = atoi(argv[++i]);
150                                  break;
151                          case 'y':               /* y resolution */
152 <                                yres = atoi(argv[++i]);
152 >                        case 'Y':               /* y resolution */
153 >                                resolution = 0;
154 >                                if (argv[i][0] == '-')
155 >                                        picres.or |= YDECR;
156 >                                if (picres.xr == 0)
157 >                                        picres.or |= YMAJOR;
158 >                                picres.yr = atoi(argv[++i]);
159                                  break;
160                          default:
161   unkopt:
# Line 102 | Line 166 | unkopt:
166                          }
167                  else
168                          break;
169 <                        
169 >                                        /* recognize special formats */
170 >        if (dataonly && format == 'b')
171 >                if (brightonly)
172 >                        fmtid = "8-bit_grey";
173 >                else
174 >                        fmtid = "24-bit_rgb";
175 >                                        /* assign reverse ordering */
176 >        rord[ord[0]] = 0;
177 >        rord[ord[1]] = 1;
178 >        rord[ord[2]] = 2;
179 >                                        /* get input */
180          if (i == argc) {
181                  fin = stdin;
182          } else if (i == argc-1) {
# Line 119 | Line 193 | unkopt:
193          set_io();
194  
195          if (reverse) {
196 <                if (header)                     /* get header */
197 <                        copyheader(fin, stdout);
198 <                                                /* add to header */
199 <                printargs(i, argv, stdout);
200 <                printf("\n");
201 <                if (yres <= 0 || xres <= 0) {
202 <                        fprintf(stderr, "%s: missing x and y resolution\n",
203 <                                        progname);
196 > #ifdef MSDOS
197 >                setmode(fileno(stdout), O_BINARY);
198 >                if (format != 'a' && format != 'i')
199 >                        setmode(fileno(fin), O_BINARY);
200 > #endif
201 >                                        /* get header */
202 >                if (header && checkheader(fin, fmtid, stdout) < 0) {
203 >                        fprintf(stderr, "%s: wrong input format\n", progname);
204                          quit(1);
205                  }
206 <                fputresolu(YMAJOR|YDECR, xres, yres, stdout);
206 >                                        /* get resolution */
207 >                if ((resolution && !fgetsresolu(&picres, fin)) ||
208 >                                picres.xr <= 0 || picres.yr <= 0) {
209 >                        fprintf(stderr, "%s: missing resolution\n", progname);
210 >                        quit(1);
211 >                }
212 >                                                /* add to header */
213 >                printargs(i, argv, stdout);
214 >                if (doexposure > 100)
215 >                        fputexpos(expval, stdout);
216 >                fputformat(COLRFMT, stdout);
217 >                putchar('\n');
218 >                fputsresolu(&picres, stdout);   /* always put resolution */
219                  valtopix();
220          } else {
221 + #ifdef MSDOS
222 +                setmode(fileno(fin), O_BINARY);
223 +                if (format != 'a' && format != 'i')
224 +                        setmode(fileno(stdout), O_BINARY);
225 + #endif
226                                                  /* get header */
227 <                getheader(fin, checkhead);
228 <
229 <                if (xres <= 0 || yres <= 0)             /* get picture size */
230 <                        if (fgetresolu(&xres, &yres, fin) != (YMAJOR|YDECR)) {
231 <                                fprintf(stderr,
232 <                                "%s: missing x and y resolution\n",
233 <                                                progname);
234 <                                quit(1);
235 <                        }
227 >                getheader(fin, checkhead, NULL);
228 >                if (wrongformat) {
229 >                        fprintf(stderr, "%s: input not a Radiance picture\n",
230 >                                        progname);
231 >                        quit(1);
232 >                }
233 >                if (!fgetsresolu(&picres, fin)) {
234 >                        fprintf(stderr, "%s: missing resolution\n", progname);
235 >                        quit(1);
236 >                }
237                  if (header) {
238                          printargs(i, argv, stdout);
239 <                        printf("\n");
239 >                        if (doexposure > 100)
240 >                                fputexpos(expval, stdout);
241 >                        fputformat(fmtid, stdout);
242 >                        putchar('\n');
243                  }
244 +                if (resolution)                 /* put resolution */
245 +                        fputsresolu(&picres, stdout);
246                  pixtoval();
247          }
248  
# Line 156 | Line 253 | unkopt:
253   checkhead(line)                         /* deal with line from header */
254   char  *line;
255   {
256 <        if (header)
256 >        char    fmt[32];
257 >        double  d;
258 >        COLOR   ctmp;
259 >
260 >        if (isformat(line)) {
261 >                formatval(fmt, line);
262 >                wrongformat = strcmp(fmt, COLRFMT);
263 >        } else if (original && isexpos(line)) {
264 >                d = 1.0/exposval(line);
265 >                scalecolor(exposure, d);
266 >                doexposure++;
267 >        } else if (original && iscolcor(line)) {
268 >                colcorval(ctmp, line);
269 >                setcolor(exposure, colval(exposure,RED)/colval(ctmp,RED),
270 >                                colval(exposure,GRN)/colval(ctmp,GRN),
271 >                                colval(exposure,BLU)/colval(ctmp,BLU));
272 >                doexposure++;
273 >        } else if (header)
274                  fputs(line, stdout);
161 –        if (!strncmp(line, "EXPOSURE=", 9))
162 –                exposure *= atof(line+9);
275   }
276  
277  
278   pixtoval()                              /* convert picture to values */
279   {
280 <        COLOR  *scanln;
280 >        register COLOR  *scanln;
281 >        int  dogamma;
282          COLOR  lastc;
283 +        FLOAT  hv[2];
284          int  y;
285          register int  x;
286  
287 <        scanln = (COLOR *)malloc(xres*sizeof(COLOR));
287 >        scanln = (COLOR *)malloc(scanlen(&picres)*sizeof(COLOR));
288          if (scanln == NULL) {
289                  fprintf(stderr, "%s: out of memory\n", progname);
290                  quit(1);
291          }
292 +        dogamma = gamcor < .95 || gamcor > 1.05;
293          setcolor(lastc, 0.0, 0.0, 0.0);
294 <        for (y = yres-1; y >= 0; y--) {
295 <                if (freadscan(scanln, xres, fin) < 0) {
294 >        for (y = 0; y < numscans(&picres); y++) {
295 >                if (freadscan(scanln, scanlen(&picres), fin) < 0) {
296                          fprintf(stderr, "%s: read error\n", progname);
297                          quit(1);
298                  }
299 <                for (x = 0; x < xres; x++) {
299 >                for (x = 0; x < scanlen(&picres); x++) {
300                          if (uniq)
301 <                                if (    scanln[x][RED] == lastc[RED] &&
302 <                                        scanln[x][GRN] == lastc[GRN] &&
303 <                                        scanln[x][BLU] == lastc[BLU]    )
301 >                                if (    colval(scanln[x],RED) ==
302 >                                                colval(lastc,RED) &&
303 >                                        colval(scanln[x],GRN) ==
304 >                                                colval(lastc,GRN) &&
305 >                                        colval(scanln[x],BLU) ==
306 >                                                colval(lastc,BLU)       )
307                                          continue;
308                                  else
309                                          copycolor(lastc, scanln[x]);
310 <                        if (original)
311 <                                scalecolor(scanln[x], 1.0/exposure);
312 <                        if (!dataonly)
313 <                                printf("%7d %7d ", x, y);
310 >                        if (doexposure)
311 >                                multcolor(scanln[x], exposure);
312 >                        if (dogamma)
313 >                                setcolor(scanln[x],
314 >                                pow(colval(scanln[x],RED), 1.0/gamcor),
315 >                                pow(colval(scanln[x],GRN), 1.0/gamcor),
316 >                                pow(colval(scanln[x],BLU), 1.0/gamcor));
317 >                        if (!dataonly) {
318 >                                pix2loc(hv, &picres, x, y);
319 >                                printf("%7d %7d ", (int)(hv[0]*picres.xr),
320 >                                                (int)(hv[1]*picres.yr));
321 >                        }
322                          if ((*putval)(scanln[x], stdout) < 0) {
323                                  fprintf(stderr, "%s: write error\n", progname);
324                                  quit(1);
# Line 205 | Line 331 | pixtoval()                             /* convert picture to values */
331  
332   valtopix()                      /* convert values to a pixel file */
333   {
334 <        COLOR  *scanln;
334 >        int  dogamma;
335 >        register COLOR  *scanln;
336          int  y;
337          register int  x;
338  
339 <        scanln = (COLOR *)malloc(xres*sizeof(COLOR));
339 >        scanln = (COLOR *)malloc(scanlen(&picres)*sizeof(COLOR));
340          if (scanln == NULL) {
341                  fprintf(stderr, "%s: out of memory\n", progname);
342                  quit(1);
343          }
344 <        for (y = yres-1; y >= 0; y--) {
345 <                for (x = 0; x < xres; x++) {
344 >        dogamma = gamcor < .95 || gamcor > 1.05;
345 >        for (y = 0; y < numscans(&picres); y++) {
346 >                for (x = 0; x < scanlen(&picres); x++) {
347                          if (!dataonly)
348                                  fscanf(fin, "%*d %*d");
349                          if ((*getval)(scanln[x], fin) < 0) {
350                                  fprintf(stderr, "%s: read error\n", progname);
351                                  quit(1);
352                          }
353 +                        if (dogamma)
354 +                                setcolor(scanln[x],
355 +                                        pow(colval(scanln[x],RED), gamcor),
356 +                                        pow(colval(scanln[x],GRN), gamcor),
357 +                                        pow(colval(scanln[x],BLU), gamcor));
358 +                        if (doexposure)
359 +                                multcolor(scanln[x], exposure);
360                  }
361 <                if (fwritescan(scanln, xres, stdout) < 0
227 <                                || fflush(stdout) < 0) {
361 >                if (fwritescan(scanln, scanlen(&picres), stdout) < 0) {
362                          fprintf(stderr, "%s: write error\n", progname);
363                          quit(1);
364                  }
# Line 244 | Line 378 | getcascii(col, fp)             /* get an ascii color value from f
378   COLOR  col;
379   FILE  *fp;
380   {
381 <        double  vd[3];
381 >        double  vd[3];
382  
383          if (fscanf(fp, "%lf %lf %lf", &vd[0], &vd[1], &vd[2]) != 3)
384                  return(-1);
385 <        setcolor(col, vd[0], vd[1], vd[2]);
385 >        setcolor(col, vd[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
386          return(0);
387   }
388  
# Line 257 | Line 391 | getcdouble(col, fp)            /* get a double color value from
391   COLOR  col;
392   FILE  *fp;
393   {
394 <        double  vd[3];
394 >        double  vd[3];
395  
396 <        if (fread(vd, sizeof(double), 3, fp) != 3)
396 >        if (fread((char *)vd, sizeof(double), 3, fp) != 3)
397                  return(-1);
398 <        setcolor(col, vd[0], vd[1], vd[2]);
398 >        setcolor(col, vd[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
399          return(0);
400   }
401  
# Line 272 | Line 406 | FILE  *fp;
406   {
407          float  vf[3];
408  
409 <        if (fread(vf, sizeof(float), 3, fp) != 3)
409 >        if (fread((char *)vf, sizeof(float), 3, fp) != 3)
410                  return(-1);
411 <        setcolor(col, vf[0], vf[1], vf[2]);
411 >        setcolor(col, vf[rord[RED]], vf[rord[GRN]], vf[rord[BLU]]);
412          return(0);
413   }
414  
# Line 287 | Line 421 | FILE  *fp;
421  
422          if (fscanf(fp, "%d %d %d", &vi[0], &vi[1], &vi[2]) != 3)
423                  return(-1);
424 <        setcolor(col,(vi[0]+.5)/256.,(vi[1]+.5)/256.,(vi[2]+.5)/256.);
424 >        setcolor(col, (vi[rord[RED]]+.5)/256.,
425 >                        (vi[rord[GRN]]+.5)/256., (vi[rord[BLU]]+.5)/256.);
426          return(0);
427   }
428  
# Line 298 | Line 433 | FILE  *fp;
433   {
434          BYTE  vb[3];
435  
436 <        if (fread(vb, sizeof(BYTE), 3, fp) != 3)
436 >        if (fread((char *)vb, sizeof(BYTE), 3, fp) != 3)
437                  return(-1);
438 <        setcolor(col,(vb[0]+.5)/256.,(vb[1]+.5)/256.,(vb[2]+.5)/256.);
438 >        setcolor(col, (vb[rord[RED]]+.5)/256.,
439 >                        (vb[rord[GRN]]+.5)/256., (vb[rord[BLU]]+.5)/256.);
440          return(0);
441   }
442  
# Line 309 | Line 445 | getbascii(col, fp)             /* get an ascii brightness value f
445   COLOR  col;
446   FILE  *fp;
447   {
448 <        double  vd;
448 >        double  vd;
449  
450          if (fscanf(fp, "%lf", &vd) != 1)
451                  return(-1);
# Line 322 | Line 458 | getbdouble(col, fp)            /* get a double brightness value
458   COLOR  col;
459   FILE  *fp;
460   {
461 <        double  vd;
461 >        double  vd;
462  
463 <        if (fread(&vd, sizeof(double), 1, fp) != 1)
463 >        if (fread((char *)&vd, sizeof(double), 1, fp) != 1)
464                  return(-1);
465          setcolor(col, vd, vd, vd);
466          return(0);
# Line 337 | Line 473 | FILE  *fp;
473   {
474          float  vf;
475  
476 <        if (fread(&vf, sizeof(float), 1, fp) != 1)
476 >        if (fread((char *)&vf, sizeof(float), 1, fp) != 1)
477                  return(-1);
478          setcolor(col, vf, vf, vf);
479          return(0);
# Line 349 | Line 485 | COLOR  col;
485   FILE  *fp;
486   {
487          int  vi;
488 <        double  d;
488 >        double  d;
489  
490          if (fscanf(fp, "%d", &vi) != 1)
491                  return(-1);
# Line 364 | Line 500 | COLOR  col;
500   FILE  *fp;
501   {
502          BYTE  vb;
503 <        double  d;
503 >        double  d;
504  
505 <        if (fread(&vb, sizeof(BYTE), 1, fp) != 1)
505 >        if (fread((char *)&vb, sizeof(BYTE), 1, fp) != 1)
506                  return(-1);
507          d = (vb+.5)/256.;
508          setcolor(col, d, d, d);
# Line 379 | Line 515 | COLOR  col;
515   FILE  *fp;
516   {
517          fprintf(fp, "%15.3e %15.3e %15.3e\n",
518 <                        colval(col,RED),
519 <                        colval(col,GRN),
520 <                        colval(col,BLU));
518 >                        colval(col,ord[0]),
519 >                        colval(col,ord[1]),
520 >                        colval(col,ord[2]));
521  
522          return(ferror(fp) ? -1 : 0);
523   }
# Line 393 | Line 529 | FILE  *fp;
529   {
530          float  vf[3];
531  
532 <        vf[0] = colval(col,RED);
533 <        vf[1] = colval(col,GRN);
534 <        vf[2] = colval(col,BLU);
535 <        fwrite(vf, sizeof(float), 3, fp);
532 >        vf[0] = colval(col,ord[0]);
533 >        vf[1] = colval(col,ord[1]);
534 >        vf[2] = colval(col,ord[2]);
535 >        fwrite((char *)vf, sizeof(float), 3, fp);
536  
537          return(ferror(fp) ? -1 : 0);
538   }
# Line 406 | Line 542 | putcdouble(col, fp)                    /* put a double color to fp */
542   COLOR  col;
543   FILE  *fp;
544   {
545 <        double  vd[3];
545 >        double  vd[3];
546  
547 <        vd[0] = colval(col,RED);
548 <        vd[1] = colval(col,GRN);
549 <        vd[2] = colval(col,BLU);
550 <        fwrite(vd, sizeof(double), 3, fp);
547 >        vd[0] = colval(col,ord[0]);
548 >        vd[1] = colval(col,ord[1]);
549 >        vd[2] = colval(col,ord[2]);
550 >        fwrite((char *)vd, sizeof(double), 3, fp);
551  
552          return(ferror(fp) ? -1 : 0);
553   }
# Line 422 | Line 558 | COLOR  col;
558   FILE  *fp;
559   {
560          fprintf(fp, "%d %d %d\n",
561 <                        (int)(colval(col,RED)*256.),
562 <                        (int)(colval(col,GRN)*256.),
563 <                        (int)(colval(col,BLU)*256.));
561 >                        (int)(colval(col,ord[0])*256.),
562 >                        (int)(colval(col,ord[1])*256.),
563 >                        (int)(colval(col,ord[2])*256.));
564  
565          return(ferror(fp) ? -1 : 0);
566   }
# Line 437 | Line 573 | FILE  *fp;
573          register int  i;
574          BYTE  vb[3];
575  
576 <        i = colval(col,RED)*256.;
576 >        i = colval(col,ord[0])*256.;
577          vb[0] = min(i,255);
578 <        i = colval(col,GRN)*256.;
578 >        i = colval(col,ord[1])*256.;
579          vb[1] = min(i,255);
580 <        i = colval(col,BLU)*256.;
580 >        i = colval(col,ord[2])*256.;
581          vb[2] = min(i,255);
582 <        fwrite(vb, sizeof(BYTE), 3, fp);
582 >        fwrite((char *)vb, sizeof(BYTE), 3, fp);
583  
584          return(ferror(fp) ? -1 : 0);
585   }
# Line 466 | Line 602 | FILE  *fp;
602          float  vf;
603  
604          vf = bright(col);
605 <        fwrite(&vf, sizeof(float), 1, fp);
605 >        fwrite((char *)&vf, sizeof(float), 1, fp);
606  
607          return(ferror(fp) ? -1 : 0);
608   }
# Line 476 | Line 612 | putbdouble(col, fp)                    /* put a double brightness to fp
612   COLOR  col;
613   FILE  *fp;
614   {
615 <        double  vd;
615 >        double  vd;
616  
617          vd = bright(col);
618 <        fwrite(&vd, sizeof(double), 1, fp);
618 >        fwrite((char *)&vd, sizeof(double), 1, fp);
619  
620          return(ferror(fp) ? -1 : 0);
621   }
# Line 504 | Line 640 | FILE  *fp;
640  
641          i = bright(col)*256.;
642          vb = min(i,255);
643 <        fwrite(&vb, sizeof(BYTE), 1, fp);
643 >        fwrite((char *)&vb, sizeof(BYTE), 1, fp);
644  
645          return(ferror(fp) ? -1 : 0);
646   }

Diff Legend

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