ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pvalue.c
Revision: 2.7
Committed: Tue Nov 9 15:21:00 1993 UTC (30 years, 5 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.6: +26 -3 lines
Log Message:
added +/-e option to correct/adjust exposure

File Contents

# Content
1 /* Copyright (c) 1992 Regents of the University of California */
2
3 #ifndef lint
4 static char SCCSid[] = "$SunId$ LBL";
5 #endif
6
7 /*
8 * pvalue.c - program to print pixel values.
9 *
10 * 4/23/86
11 */
12
13 #include "standard.h"
14
15 #include "color.h"
16
17 #include "resolu.h"
18
19 #define min(a,b) ((a)<(b)?(a):(b))
20
21 RESOLU picres; /* resolution of picture */
22
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? */
30
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? */
37
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;
54
55 int (*getval)(), (*putval)();
56
57
58 main(argc, argv)
59 int argc;
60 char **argv;
61 {
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] == '-' || argv[i][0] == '+')
70 switch (argv[i][1]) {
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 = argv[i][0] == '-';
79 break;
80 case 'o': /* original values */
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 = argv[i][0] == '-';
109 break;
110 case 'b': /* brightness values */
111 brightonly = argv[i][0] == '-';
112 break;
113 case 'd': /* data only (no indices) */
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 dataonly = 1;
137 format = 'd';
138 fmtid = "double";
139 break;
140 default:
141 goto unkopt;
142 }
143 break;
144 case 'x': /* x resolution */
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 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:
162 fprintf(stderr, "%s: unknown option: %s\n",
163 progname, argv[i]);
164 quit(1);
165 break;
166 }
167 else
168 break;
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) {
183 if ((fin = fopen(argv[i], "r")) == NULL) {
184 fprintf(stderr, "%s: can't open file \"%s\"\n",
185 progname, argv[i]);
186 quit(1);
187 }
188 } else {
189 fprintf(stderr, "%s: bad # file arguments\n", progname);
190 quit(1);
191 }
192
193 set_io();
194
195 if (reverse) {
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 /* 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, 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 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
249 quit(0);
250 }
251
252
253 checkhead(line) /* deal with line from header */
254 char *line;
255 {
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);
275 }
276
277
278 pixtoval() /* convert picture to values */
279 {
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(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 = 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 < scanlen(&picres); x++) {
300 if (uniq)
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 (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);
325 }
326 }
327 }
328 free((char *)scanln);
329 }
330
331
332 valtopix() /* convert values to a pixel file */
333 {
334 int dogamma;
335 register COLOR *scanln;
336 int y;
337 register int x;
338
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 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, scanlen(&picres), stdout) < 0) {
362 fprintf(stderr, "%s: write error\n", progname);
363 quit(1);
364 }
365 }
366 free((char *)scanln);
367 }
368
369
370 quit(code)
371 int code;
372 {
373 exit(code);
374 }
375
376
377 getcascii(col, fp) /* get an ascii color value from fp */
378 COLOR col;
379 FILE *fp;
380 {
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[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
386 return(0);
387 }
388
389
390 getcdouble(col, fp) /* get a double color value from fp */
391 COLOR col;
392 FILE *fp;
393 {
394 double vd[3];
395
396 if (fread((char *)vd, sizeof(double), 3, fp) != 3)
397 return(-1);
398 setcolor(col, vd[rord[RED]], vd[rord[GRN]], vd[rord[BLU]]);
399 return(0);
400 }
401
402
403 getcfloat(col, fp) /* get a float color value from fp */
404 COLOR col;
405 FILE *fp;
406 {
407 float vf[3];
408
409 if (fread((char *)vf, sizeof(float), 3, fp) != 3)
410 return(-1);
411 setcolor(col, vf[rord[RED]], vf[rord[GRN]], vf[rord[BLU]]);
412 return(0);
413 }
414
415
416 getcint(col, fp) /* get an int color value from fp */
417 COLOR col;
418 FILE *fp;
419 {
420 int vi[3];
421
422 if (fscanf(fp, "%d %d %d", &vi[0], &vi[1], &vi[2]) != 3)
423 return(-1);
424 setcolor(col, (vi[rord[RED]]+.5)/256.,
425 (vi[rord[GRN]]+.5)/256., (vi[rord[BLU]]+.5)/256.);
426 return(0);
427 }
428
429
430 getcbyte(col, fp) /* get a byte color value from fp */
431 COLOR col;
432 FILE *fp;
433 {
434 BYTE vb[3];
435
436 if (fread((char *)vb, sizeof(BYTE), 3, fp) != 3)
437 return(-1);
438 setcolor(col, (vb[rord[RED]]+.5)/256.,
439 (vb[rord[GRN]]+.5)/256., (vb[rord[BLU]]+.5)/256.);
440 return(0);
441 }
442
443
444 getbascii(col, fp) /* get an ascii brightness value from fp */
445 COLOR col;
446 FILE *fp;
447 {
448 double vd;
449
450 if (fscanf(fp, "%lf", &vd) != 1)
451 return(-1);
452 setcolor(col, vd, vd, vd);
453 return(0);
454 }
455
456
457 getbdouble(col, fp) /* get a double brightness value from fp */
458 COLOR col;
459 FILE *fp;
460 {
461 double vd;
462
463 if (fread((char *)&vd, sizeof(double), 1, fp) != 1)
464 return(-1);
465 setcolor(col, vd, vd, vd);
466 return(0);
467 }
468
469
470 getbfloat(col, fp) /* get a float brightness value from fp */
471 COLOR col;
472 FILE *fp;
473 {
474 float vf;
475
476 if (fread((char *)&vf, sizeof(float), 1, fp) != 1)
477 return(-1);
478 setcolor(col, vf, vf, vf);
479 return(0);
480 }
481
482
483 getbint(col, fp) /* get an int brightness value from fp */
484 COLOR col;
485 FILE *fp;
486 {
487 int vi;
488 double d;
489
490 if (fscanf(fp, "%d", &vi) != 1)
491 return(-1);
492 d = (vi+.5)/256.;
493 setcolor(col, d, d, d);
494 return(0);
495 }
496
497
498 getbbyte(col, fp) /* get a byte brightness value from fp */
499 COLOR col;
500 FILE *fp;
501 {
502 BYTE vb;
503 double d;
504
505 if (fread((char *)&vb, sizeof(BYTE), 1, fp) != 1)
506 return(-1);
507 d = (vb+.5)/256.;
508 setcolor(col, d, d, d);
509 return(0);
510 }
511
512
513 putcascii(col, fp) /* put an ascii color to fp */
514 COLOR col;
515 FILE *fp;
516 {
517 fprintf(fp, "%15.3e %15.3e %15.3e\n",
518 colval(col,ord[0]),
519 colval(col,ord[1]),
520 colval(col,ord[2]));
521
522 return(ferror(fp) ? -1 : 0);
523 }
524
525
526 putcfloat(col, fp) /* put a float color to fp */
527 COLOR col;
528 FILE *fp;
529 {
530 float vf[3];
531
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 }
539
540
541 putcdouble(col, fp) /* put a double color to fp */
542 COLOR col;
543 FILE *fp;
544 {
545 double vd[3];
546
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 }
554
555
556 putcint(col, fp) /* put an int color to fp */
557 COLOR col;
558 FILE *fp;
559 {
560 fprintf(fp, "%d %d %d\n",
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 }
567
568
569 putcbyte(col, fp) /* put a byte color to fp */
570 COLOR col;
571 FILE *fp;
572 {
573 register int i;
574 BYTE vb[3];
575
576 i = colval(col,ord[0])*256.;
577 vb[0] = min(i,255);
578 i = colval(col,ord[1])*256.;
579 vb[1] = min(i,255);
580 i = colval(col,ord[2])*256.;
581 vb[2] = min(i,255);
582 fwrite((char *)vb, sizeof(BYTE), 3, fp);
583
584 return(ferror(fp) ? -1 : 0);
585 }
586
587
588 putbascii(col, fp) /* put an ascii brightness to fp */
589 COLOR col;
590 FILE *fp;
591 {
592 fprintf(fp, "%15.3e\n", bright(col));
593
594 return(ferror(fp) ? -1 : 0);
595 }
596
597
598 putbfloat(col, fp) /* put a float brightness to fp */
599 COLOR col;
600 FILE *fp;
601 {
602 float vf;
603
604 vf = bright(col);
605 fwrite((char *)&vf, sizeof(float), 1, fp);
606
607 return(ferror(fp) ? -1 : 0);
608 }
609
610
611 putbdouble(col, fp) /* put a double brightness to fp */
612 COLOR col;
613 FILE *fp;
614 {
615 double vd;
616
617 vd = bright(col);
618 fwrite((char *)&vd, sizeof(double), 1, fp);
619
620 return(ferror(fp) ? -1 : 0);
621 }
622
623
624 putbint(col, fp) /* put an int brightness to fp */
625 COLOR col;
626 FILE *fp;
627 {
628 fprintf(fp, "%d\n", (int)(bright(col)*256.));
629
630 return(ferror(fp) ? -1 : 0);
631 }
632
633
634 putbbyte(col, fp) /* put a byte brightness to fp */
635 COLOR col;
636 FILE *fp;
637 {
638 register int i;
639 BYTE vb;
640
641 i = bright(col)*256.;
642 vb = min(i,255);
643 fwrite((char *)&vb, sizeof(BYTE), 1, fp);
644
645 return(ferror(fp) ? -1 : 0);
646 }
647
648
649 set_io() /* set put and get functions */
650 {
651 switch (format) {
652 case 'a': /* ascii */
653 if (brightonly) {
654 getval = getbascii;
655 putval = putbascii;
656 } else {
657 getval = getcascii;
658 putval = putcascii;
659 }
660 return;
661 case 'f': /* binary float */
662 if (brightonly) {
663 getval = getbfloat;
664 putval = putbfloat;
665 } else {
666 getval = getcfloat;
667 putval = putcfloat;
668 }
669 return;
670 case 'd': /* binary double */
671 if (brightonly) {
672 getval = getbdouble;
673 putval = putbdouble;
674 } else {
675 getval = getcdouble;
676 putval = putcdouble;
677 }
678 return;
679 case 'i': /* integer */
680 if (brightonly) {
681 getval = getbint;
682 putval = putbint;
683 } else {
684 getval = getcint;
685 putval = putcint;
686 }
687 return;
688 case 'b': /* byte */
689 if (brightonly) {
690 getval = getbbyte;
691 putval = putbbyte;
692 } else {
693 getval = getcbyte;
694 putval = putcbyte;
695 }
696 return;
697 }
698 }