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

Comparing ray/src/px/ra_t8.c (file contents):
Revision 1.4 by greg, Thu Jan 18 23:58:24 1990 UTC vs.
Revision 2.4 by greg, Tue Oct 13 11:36:13 1992 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 15 | Line 15 | static char SCCSid[] = "$SunId$ LBL";
15  
16   #include  "color.h"
17  
18 < #include  "pic.h"
18 > #include  "resolu.h"
19  
20   #include  "targa.h"
21  
22 < #ifndef  BSD
23 < #define  bcopy(s,d,n)           (void)memcpy(d,s,n)
22 > #ifdef MSDOS
23 > #include  <fcntl.h>
24 > #endif
25 >
26 > #include  <math.h>
27 >
28 > #ifndef  BSD
29 > #define  bcopy(s,d,n)           (void)memcpy(d,s,n)
30   extern char  *memcpy();
31   #endif
26 –                        /* descriptor for a picture file or frame buffer */
27 – typedef struct {
28 –        char    *name;                  /* file name */
29 –        FILE    *fp;                    /* file pointer */
30 –        int     nexty;                  /* file positioning */
31 –        int     bytes_line;             /* 0 == variable length lines */
32 –        union {
33 –                long    b;                      /* initial scanline */
34 –                long    *y;                     /* individual scanline */
35 –        } pos;                          /* position(s) */
36 – } pic;
32  
33 < #define  goodpic(h)     (my_imType(h) && my_mapType(h))
34 < #define  my_imType(h)   (((h)->dataType==IM_CMAP || (h)->dataType==IM_CCMAP) \
33 > #define  goodpic(h)     (my_imType(h) && my_mapType(h))
34 > #define  my_imType(h)   (((h)->dataType==IM_CMAP || (h)->dataType==IM_CCMAP) \
35                                  && (h)->dataBits==8 && (h)->imType==0)
36 < #define  my_mapType(h)  ((h)->mapType==CM_HASMAP && \
36 > #define  my_mapType(h)  ((h)->mapType==CM_HASMAP && \
37                                  ((h)->CMapBits==24 || (h)->CMapBits==32))
38  
39 < #define  taralloc(h)    (pixel *)emalloc((h)->x*(h)->y*sizeof(pixel))
39 > #define  taralloc(h)    (BYTE *)emalloc((h)->x*(h)->y)
40  
41 < extern pic      *openinput();
41 > extern BYTE  clrtab[][3];
42  
43   extern char     *ecalloc(), *emalloc();
44  
45   extern long  ftell();
46  
47 < extern double  atof(), pow();
47 > double  gamv = 2.2;                     /* gamv correction */
48  
49 < double  gamma = 2.0;                    /* gamma correction */
49 > int  bradj = 0;                         /* brightness adjustment */
50  
56 – pic     *inpic;
57 –
51   char  *progname;
52  
53   char  errmsg[128];
54  
55 < COLR    *inline;
55 > COLR    *inl;
56  
57 < pixel   *tarData;
57 > BYTE    *tarData;
58  
59   int  xmax, ymax;
60  
# Line 70 | Line 63 | main(argc, argv)
63   int  argc;
64   char  *argv[];
65   {
66 <        colormap  rasmap;
74 <        struct hdStruct  head;
66 >        struct hdStruct  head;
67          int  dither = 1;
68          int  reverse = 0;
69          int  ncolors = 256;
70          int  greyscale = 0;
71          int  i;
72 <        
72 > #ifdef MSDOS
73 >        extern int  _fmode;
74 >        _fmode = O_BINARY;
75 >        setmode(fileno(stdin), O_BINARY);
76 >        setmode(fileno(stdout), O_BINARY);
77 > #endif
78          progname = argv[0];
79  
80          for (i = 1; i < argc; i++)
# Line 87 | Line 84 | char  *argv[];
84                                  dither = !dither;
85                                  break;
86                          case 'g':
87 <                                gamma = atof(argv[++i]);
87 >                                gamv = atof(argv[++i]);
88                                  break;
89                          case 'r':
90                                  reverse = !reverse;
# Line 95 | Line 92 | char  *argv[];
92                          case 'b':
93                                  greyscale = 1;
94                                  break;
95 +                        case 'e':
96 +                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
97 +                                        goto userr;
98 +                                bradj = atoi(argv[++i]);
99 +                                break;
100                          case 'c':
101                                  ncolors = atoi(argv[++i]);
102                                  break;
# Line 104 | Line 106 | char  *argv[];
106                  else
107                          break;
108  
109 +        if (i < argc-2)
110 +                goto userr;
111 +        if (i <= argc-1 && freopen(argv[i], "r", stdin) == NULL) {
112 +                sprintf(errmsg, "cannot open input \"%s\"", argv[i]);
113 +                quiterr(errmsg);
114 +        }
115 +        if (i == argc-2 && freopen(argv[i+1], "w", stdout) == NULL) {
116 +                sprintf(errmsg, "cannot open output \"%s\"", argv[i+1]);
117 +                quiterr(errmsg);
118 +        }
119          if (reverse) {
108 –                if (i < argc-2)
109 –                        goto userr;
110 –                if (i <= argc-1 && freopen(argv[i], "r", stdin) == NULL) {
111 –                        sprintf(errmsg, "can't open input \"%s\"", argv[i]);
112 –                        quiterr(errmsg);
113 –                }
120                                          /* get header */
121                  if (getthead(&head, NULL, stdin) < 0)
122                          quiterr("bad targa file");
# Line 118 | Line 124 | char  *argv[];
124                          quiterr("incompatible format");
125                  xmax = head.x;
126                  ymax = head.y;
121 –                                        /* open output file */
122 –                if (i == argc-2 && freopen(argv[i+1], "w", stdout) == NULL) {
123 –                        sprintf(errmsg, "can't open output \"%s\"", argv[i+1]);
124 –                        quiterr(errmsg);
125 –                }
127                                          /* put header */
128 <                printargs(argc, argv, stdout);
128 >                printargs(i, argv, stdout);
129 >                fputformat(COLRFMT, stdout);
130                  putchar('\n');
131 <                fputresolu(YMAJOR|YDECR, xmax, ymax, stdout);
131 >                fprtresolu(xmax, ymax, stdout);
132                                          /* convert file */
133                  tg2ra(&head);
134          } else {
135 <                if (i > argc-1 || i < argc-2)
136 <                        goto userr;
135 <                if ((inpic = openinput(argv[i], &head)) == NULL) {
136 <                        sprintf(errmsg, "can't open input \"%s\"", argv[i]);
137 <                        quiterr(errmsg);
138 <                }
139 <                if (i == argc-2 && freopen(argv[i+1], "w", stdout) == NULL) {
140 <                        sprintf(errmsg, "can't open output \"%s\"", argv[i+1]);
141 <                        quiterr(errmsg);
142 <                }
135 >                if (getrhead(&head, stdin) < 0)
136 >                        quiterr("bad Radiance input");
137                                          /* write header */
138                  putthead(&head, NULL, stdout);
139                                          /* convert file */
140                  if (greyscale)
141 <                        biq(dither,ncolors,1,rasmap);
141 >                        getgrey(ncolors);
142                  else
143 <                        ciq(dither,ncolors,1,rasmap);
143 >                        getmapped(ncolors, dither);
144                                          /* write data */
145                  writetarga(&head, tarData, stdout);
146          }
147          quiterr(NULL);
148   userr:
149          fprintf(stderr,
150 <        "Usage: %s [-d][-c ncolors][-b][-g gamma] input [output]\n",
150 >        "Usage: %s [-d][-c ncolors][-b][-g gamv][-e +/-stops] input [output]\n",
151                          progname);
152 <        fprintf(stderr, "   Or: %s -r [-g gamma] [input [output]]\n",
152 >        fprintf(stderr, "   Or: %s -r [-g gamv][-e +/-stops] [input [output]]\n",
153                          progname);
154          exit(1);
155   }
# Line 209 | Line 203 | int code;
203  
204  
205   getthead(hp, ip, fp)            /* read header from input */
206 < struct hdStruct  *hp;
206 > struct hdStruct  *hp;
207   char  *ip;
208   register FILE  *fp;
209   {
# Line 242 | Line 236 | register FILE  *fp;
236  
237  
238   putthead(hp, ip, fp)            /* write header to output */
239 < struct hdStruct  *hp;
239 > struct hdStruct  *hp;
240   char  *ip;
241   register FILE  *fp;
242   {
# Line 269 | Line 263 | register FILE  *fp;
263   }
264  
265  
266 < pic *
273 < openinput(fname, h)             /* open RADIANCE input file */
274 < char  *fname;
266 > getrhead(h, fp)                 /* load RADIANCE input file header */
267   register struct hdStruct  *h;
268 + FILE  *fp;
269   {
270 <        register pic  *p;
271 <
272 <        p = (pic *)emalloc(sizeof(pic));
273 <        p->name = fname;
281 <        if (fname == NULL)
282 <                p->fp = stdin;
283 <        else if ((p->fp = fopen(fname, "r")) == NULL)
284 <                return(NULL);
285 <                                        /* discard header */
286 <        getheader(p->fp, NULL);
287 <        if (fgetresolu(&xmax, &ymax, p->fp) != (YMAJOR|YDECR))
288 <                quiterr("bad picture size");
289 <        p->nexty = 0;
290 <        p->bytes_line = 0;              /* variable length lines */
291 <        p->pos.y = (long *)ecalloc(ymax, sizeof(long));
292 <        p->pos.y[0] = ftell(p->fp);
270 >                                        /* get header info. */
271 >        if (checkheader(fp, COLRFMT, NULL) < 0 ||
272 >                        fgetresolu(&xmax, &ymax, fp) < 0)
273 >                return(-1);
274                                          /* assign header */
275          h->textSize = 0;
276          h->mapType = CM_HASMAP;
# Line 304 | Line 285 | register struct hdStruct  *h;
285          h->dataBits = 8;
286          h->imType = 0;
287                                          /* allocate scanline */
288 <        inline = (COLR *)emalloc(xmax*sizeof(COLR));
288 >        inl = (COLR *)emalloc(xmax*sizeof(COLR));
289                                          /* allocate targa data */
290          tarData = taralloc(h);
291  
292 <        return(p);
292 >        return(0);
293   }
294  
295  
296   tg2ra(hp)                       /* targa file to RADIANCE file */
297 < struct hdStruct  *hp;
297 > struct hdStruct  *hp;
298   {
299          union {
300                  BYTE  c3[256][3];
# Line 324 | Line 305 | struct hdStruct  *hp;
305          register int  i, j;
306  
307                                          /* get color table */
308 <        if ((hp->CMapBits==24 ? fread((char *)map.c3,sizeof(map.c3),1,stdin) :
309 <                        fread((char *)map.c4,sizeof(map.c4),1,stdin)) != 1)
308 >        if ((hp->CMapBits==24 ? fread((char *)(map.c3+hp->mapOrig),
309 >                                3*hp->mapLength,1,stdin) :
310 >                        fread((char *)(map.c4+hp->mapOrig),
311 >                                4*hp->mapLength,1,stdin)) != 1)
312                  quiterr("error reading color table");
313                                          /* convert table */
314          for (i = hp->mapOrig; i < hp->mapOrig+hp->mapLength; i++)
315                  if (hp->CMapBits == 24)
316                          setcolr(ctab[i],
317 <                                        pow((map.c3[i][2]+.5)/256.,gamma),
318 <                                        pow((map.c3[i][1]+.5)/256.,gamma),
319 <                                        pow((map.c3[i][0]+.5)/256.,gamma));
317 >                                        pow((map.c3[i][2]+.5)/256.,gamv),
318 >                                        pow((map.c3[i][1]+.5)/256.,gamv),
319 >                                        pow((map.c3[i][0]+.5)/256.,gamv));
320                  else
321                          setcolr(ctab[i],
322 <                                        pow((map.c4[i][3]+.5)/256.,gamma),
323 <                                        pow((map.c4[i][2]+.5)/256.,gamma),
324 <                                        pow((map.c4[i][1]+.5)/256.,gamma));
325 <
322 >                                        pow((map.c4[i][3]+.5)/256.,gamv),
323 >                                        pow((map.c4[i][2]+.5)/256.,gamv),
324 >                                        pow((map.c4[i][1]+.5)/256.,gamv));
325 >        if (bradj)
326 >                shiftcolrs(ctab, 256, bradj);
327                                          /* allocate targa data */
328          tarData = taralloc(hp);
329                                          /* get data */
# Line 358 | Line 342 | struct hdStruct  *hp;
342   }
343  
344  
345 < picreadline3(y, l3)                     /* read in 3-byte scanline */
346 < int  y;
347 < register rgbpixel  *l3;
345 > getmapped(nc, dith)             /* read in and quantize image */
346 > int  nc;                /* number of colors to use */
347 > int  dith;              /* use dithering? */
348   {
349 <        register int    i;
349 >        long  fpos;
350 >        register int  y;
351  
352 <        if (inpic->nexty != y) {                        /* find scanline */
353 <                if (inpic->bytes_line == 0) {
354 <                        if (inpic->pos.y[y] == 0) {
355 <                                while (inpic->nexty < y) {
356 <                                        if (freadcolrs(inline, xmax, inpic->fp) < 0)
357 <                                                quiterr("read error in picreadline3");
358 <                                        inpic->pos.y[++inpic->nexty] = ftell(inpic->fp);
359 <                                }
360 <                        } else if (fseek(inpic->fp, inpic->pos.y[y], 0) == EOF)
361 <                                quiterr("seek error in picreadline3");
377 <                } else if (fseek(inpic->fp, y*inpic->bytes_line+inpic->pos.b, 0) == EOF)
378 <                        quiterr("seek error in picreadline3");
379 <        } else if (inpic->bytes_line == 0 && inpic->pos.y[inpic->nexty] == 0)
380 <                inpic->pos.y[inpic->nexty] = ftell(inpic->fp);
381 <        if (freadcolrs(inline, xmax, inpic->fp) < 0)    /* read scanline */
382 <                quiterr("read error in picreadline3");
383 <        inpic->nexty = y+1;
384 <                                                        /* convert scanline */
385 <        normcolrs(inline, xmax);
386 <        for (i = 0; i < xmax; i++) {
387 <                l3[i].r = inline[i][RED];
388 <                l3[i].g = inline[i][GRN];
389 <                l3[i].b = inline[i][BLU];
352 >        setcolrgam(gamv);
353 >        fpos = ftell(stdin);
354 >        new_histo();                    /* build histogram */
355 >        for (y = ymax-1; y >= 0; y--) {
356 >                if (freadcolrs(inl, xmax, stdin) < 0)
357 >                        quiterr("error reading Radiance input");
358 >                if (bradj)
359 >                        shiftcolrs(inl, xmax, bradj);
360 >                colrs_gambs(inl, xmax);
361 >                cnt_colrs(inl, xmax);
362          }
363 +        if (fseek(stdin, fpos, 0) == EOF)
364 +                quiterr("Radiance input must be from a file");
365 +        new_clrtab(nc);                 /* map colors */
366 +        for (y = ymax-1; y >= 0; y--) {
367 +                if (freadcolrs(inl, xmax, stdin) < 0)
368 +                        quiterr("error reading Radiance input");
369 +                if (bradj)
370 +                        shiftcolrs(inl, xmax, bradj);
371 +                colrs_gambs(inl, xmax);
372 +                if (dith)
373 +                        dith_colrs(tarData+y*xmax, inl, xmax);
374 +                else
375 +                        map_colrs(tarData+y*xmax, inl, xmax);
376 +        }
377   }
378  
379  
380 < picwriteline(y, l)                      /* save output scanline */
381 < int  y;
396 < pixel  *l;
380 > getgrey(nc)                     /* read in and convert to greyscale image */
381 > int  nc;                /* number of colors to use */
382   {
383 <        bcopy((char *)l, (char *)&tarData[(ymax-1-y)*xmax], xmax*sizeof(pixel));
383 >        int  y;
384 >        register BYTE  *dp;
385 >        register int  x;
386 >
387 >        setcolrgam(gamv);
388 >        dp = tarData+xmax*ymax;;
389 >        for (y = ymax-1; y >= 0; y--) {
390 >                if (freadcolrs(inl, xmax, stdin) < 0)
391 >                        quiterr("error reading Radiance input");
392 >                if (bradj)
393 >                        shiftcolrs(inl, xmax, bradj);
394 >                colrs_gambs(inl, xmax);
395 >                x = xmax;
396 >                if (nc < 256)
397 >                        while (x--)
398 >                                *--dp = ((long)normbright(inl[x])*nc+128)>>8;
399 >                else
400 >                        while (x--)
401 >                                *--dp = normbright(inl[x]);
402 >        }
403 >        for (x = 0; x < nc; x++)
404 >                clrtab[x][RED] = clrtab[x][GRN] =
405 >                        clrtab[x][BLU] = ((long)x*256+nc/2)/nc;
406   }
407  
408  
409   writetarga(h, d, fp)            /* write out targa data */
410 < struct hdStruct  *h;
411 < pixel  *d;
410 > struct hdStruct  *h;
411 > BYTE  *d;
412   FILE  *fp;
413   {
414 +        register int  i, j;
415 +
416 +        for (i = 0; i < h->mapLength; i++)      /* write color map */
417 +                for (j = 2; j >= 0; j--)
418 +                        putc(clrtab[i][j], fp);
419          if (h->dataType == IM_CMAP) {           /* uncompressed */
420 <                if (fwrite((char *)d,h->x*sizeof(pixel),h->y,fp) != h->y)
420 >                if (fwrite((char *)d,h->x*sizeof(BYTE),h->y,fp) != h->y)
421                          quiterr("error writing targa file");
422                  return;
423          }
# Line 414 | Line 426 | FILE  *fp;
426  
427  
428   readtarga(h, data, fp)          /* read in targa data */
429 < struct hdStruct  *h;
430 < pixel  *data;
429 > struct hdStruct  *h;
430 > BYTE  *data;
431   FILE  *fp;
432   {
433          register int  cnt, c;
434 <        register pixel  *dp;
434 >        register BYTE   *dp;
435  
436          if (h->dataType == IM_CMAP) {           /* uncompressed */
437 <                if (fread((char *)data,h->x*sizeof(pixel),h->y,fp) != h->y)
437 >                if (fread((char *)data,h->x*sizeof(BYTE),h->y,fp) != h->y)
438                          goto readerr;
439                  return;
440          }
# Line 445 | Line 457 | FILE  *fp;
457          return;
458   readerr:
459          quiterr("error reading targa file");
448 – }
449 –
450 –
451 – picwritecm(cm)                  /* write out color map */
452 – colormap  cm;
453 – {
454 –        register int  i, j;
455 –
456 –        for (j = 0; j < 256; j++)
457 –                for (i = 2; i >= 0; i--)
458 –                        putc(cm[i][j], stdout);
459 – }
460 –
461 –
462 – picreadcm(map)                  /* do gamma correction if requested */
463 – colormap  map;
464 – {
465 –        register int  i, val;
466 –
467 –        for (i = 0; i < 256; i++) {
468 –                val = pow(i/256.0, 1.0/gamma) * 256.0;
469 –                map[0][i] = map[1][i] = map[2][i] = val;
470 –        }
460   }

Diff Legend

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