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

Comparing ray/src/px/ra_ps.c (file contents):
Revision 2.14 by greg, Wed Oct 4 16:01:32 1995 UTC vs.
Revision 2.26 by schorsch, Sun Jul 27 22:12:03 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   *  Radiance picture to PostScript file translator -- one way!
6   */
7  
8   #include  <stdio.h>
9   #include  <math.h>
10 < #ifdef MSDOS
11 < #include  <fcntl.h>
12 < #endif
10 > #include  <ctype.h>
11 >
12 > #include  "platform.h"
13   #include  "color.h"
14  
15 < #define GAMVAL          1.5                     /* gamma value for pixels */
15 > #define UPPER(c)        ((c)&~0x20)             /* ASCII trick */
16  
17 + #define CODE6GAM        1.47                    /* gamma for 6-bit codes */
18 + #define DEFGGAM         1.0                     /* greyscale device gamma */
19 + #define DEFCGAM         1.8                     /* color device gamma */
20 +
21   #define GRY             -1                      /* artificial index for grey */
22  
23 < #define HMARGIN         (.5*72)                 /* horizontal margin */
24 < #define VMARGIN         (.5*72)                 /* vertical margin */
25 < #define PWIDTH          (8.5*72-2*HMARGIN)      /* width of device */
26 < #define PHEIGHT         (11*72-2*VMARGIN)       /* height of device */
23 > #define DEFMARG         0.5                     /* default margin (inches) */
24 > #define DEFWIDTH        8.5                     /* default page width */
25 > #define DEFHEIGHT       11                      /* default page height */
26 > #define HMARGIN         (hmarg*72)              /* horizontal margin */
27 > #define VMARGIN         (vmarg*72)              /* vertical margin */
28 > #define PWIDTH          (width*72-2*HMARGIN)    /* width of device */
29 > #define PHEIGHT         (height*72-2*VMARGIN)   /* height of device */
30  
31   #define RUNCHR          '*'                     /* character to start rle */
32   #define MINRUN          4                       /* minimum run-length */
# Line 35 | Line 39 | char  code[] =                 /* 6-bit code lookup table */
39   int  wrongformat = 0;                   /* input in wrong format? */
40   double  pixaspect = 1.0;                /* pixel aspect ratio */
41  
42 < int  docolor = 0;                       /* produce color image? */
42 > double  devgam = 0.;                    /* device gamma response */
43 > double  hmarg = DEFMARG,
44 >        vmarg = DEFMARG;                /* horizontal and vertical margins */
45 > double  width = DEFWIDTH,
46 >        height = DEFHEIGHT;             /* default paper width and height */
47 > double  dpi = 0;                        /* print density (0 if unknown) */
48 > int  docolor = 1;                       /* produce color image? */
49   int  bradj = 0;                         /* brightness adjustment */
50   int  ncopies = 1;                       /* number of copies */
51  
52 + extern int      Aputprim(), Bputprim(), Cputprim();
53 +
54 + int  (*putprim)() = Aputprim;           /* function for writing scanline */
55 +
56   char  *progname;
57  
58 < int  xmax, ymax;
58 > int  xmax, ymax;                        /* input image dimensions */
59  
60 < extern char  *malloc();
60 > extern double   unit2inch();
61  
62  
63 + int
64   headline(s)             /* check header line */
65   char  *s;
66   {
# Line 56 | Line 71 | char  *s;
71                  wrongformat = strcmp(fmt, COLRFMT);
72          } else if (isaspect(s))
73                  pixaspect *= aspectval(s);
74 +        return(0);
75   }
76  
77  
# Line 64 | Line 80 | int  argc;
80   char  *argv[];
81   {
82          int  i;
83 +        double  d;
84          
85          progname = argv[0];
86  
87          for (i = 1; i < argc; i++)
88                  if (argv[i][0] == '-')
89                          switch (argv[i][1]) {
90 +                        case 'b':               /* produce b&w PostScript */
91 +                                docolor = 0;
92 +                                break;
93                          case 'c':               /* produce color PostScript */
94 <                                docolor++;
94 >                                docolor = 1;
95                                  break;
96 +                        case 'A':               /* standard ASCII encoding */
97 +                                putprim = Aputprim;
98 +                                break;
99 +                        case 'B':               /* standard binary encoding */
100 +                                putprim = Bputprim;
101 +                                break;
102 +                        case 'C':               /* compressed ASCII encoding */
103 +                                putprim = Cputprim;
104 +                                break;
105 +                        case 'd':               /* print density */
106 +                                dpi = atof(argv[++i]);
107 +                                break;
108 +                        case 'g':               /* device gamma adjustment */
109 +                                devgam = atof(argv[++i]);
110 +                                break;
111 +                        case 'p':               /* paper size */
112 +                                parsepaper(argv[++i]);
113 +                                break;
114 +                        case 'm':               /* margin */
115 +                                d = atof(argv[i+1]);
116 +                                d *= unit2inch(argv[i+1]);
117 +                                switch (argv[i][2]) {
118 +                                case '\0':
119 +                                        hmarg = vmarg = d;
120 +                                        break;
121 +                                case 'h':
122 +                                        hmarg = d;
123 +                                        break;
124 +                                case 'v':
125 +                                        vmarg = d;
126 +                                        break;
127 +                                default:
128 +                                        goto userr;
129 +                                }
130 +                                i++;
131 +                                break;
132                          case 'e':               /* exposure adjustment */
133                                  if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
134                                          goto userr;
# Line 99 | Line 155 | char  *argv[];
155                                  progname, argv[i+1]);
156                  exit(1);
157          }
158 < #ifdef MSDOS
103 <        setmode(fileno(stdin), O_BINARY);
104 < #endif
158 >        SET_FILE_BINARY(stdin);
159                                  /* get our header */
160          getheader(stdin, headline, NULL);
161          if (wrongformat || fgetresolu(&xmax, &ymax, stdin) < 0)
162                  quiterr("bad picture format");
163                                  /* gamma compression */
164 <        setcolrgam(GAMVAL);
164 >        if (devgam <= 0.05)
165 >                devgam = docolor ? DEFCGAM : DEFGGAM;
166 >        if (putprim == Cputprim)
167 >                setcolrgam(CODE6GAM);
168 >        else if (devgam != 1.)
169 >                setcolrgam(devgam);
170                                  /* write header */
171 <        PSheader(i <= argc-1 ? argv[i] : "<stdin>");
171 >        PSheader(argc, argv);
172                                  /* convert file */
173          ra2ps();
174                                  /* write trailer */
175          PStrailer();
176          exit(0);
177   userr:
178 <        fprintf(stderr, "Usage: %s [-c][-e +/-stops] [input [output]]\n",
178 >        fprintf(stderr,
179 > "Usage: %s [-b|c][-A|B|C][-e +/-stops][-p paper][-m[h|v] margin][-d dpi][-g gamma] [input [output]]\n",
180                          progname);
181          exit(1);
182   }
183  
184  
185 + double
186 + unit2inch(s)            /* determine unit */
187 + register char   *s;
188 + {
189 +        static struct unit {char n; float f;} u[] = {
190 +                'i', 1.,
191 +                'm', 1./25.4,
192 +                'c', 1./2.54,
193 +                '\0' };
194 +        register struct unit    *up;
195 +
196 +        while (*s && !isalpha(*s))
197 +                s++;
198 +        for (up = u; up->n; up++)
199 +                if (up->n == *s)
200 +                        return(up->f);
201 +        return(1.);
202 + }
203 +
204 +
205 + int
206 + matchid(name, id)       /* see if name matches id (case insensitive) */
207 + char    *name;
208 + register char   *id;
209 + {
210 +        register char   *s = name;
211 +
212 +        while (*s) {
213 +                if (isalpha(*s)) {
214 +                        if (!isalpha(*id) || UPPER(*s) != UPPER(*id))
215 +                                return(0);
216 +                } else if (*s != *id)
217 +                        return(0);
218 +                s++; id++;
219 +        }
220 +        return(!*id || s-name >= 3);    /* substrings >= 3 chars OK */
221 + }
222 +
223 +
224 + parsepaper(ps)          /* determine paper size from name */
225 + char    *ps;
226 + {
227 +        static struct psize {char n[12]; float w,h;} p[] = {
228 +                "envelope", 4.12, 9.5,
229 +                "executive", 7.25, 10.5,
230 +                "letter", 8.5, 11.,
231 +                "lettersmall", 7.68, 10.16,
232 +                "legal", 8.5, 14.,
233 +                "monarch", 3.87, 7.5,
234 +                "statement", 5.5, 8.5,
235 +                "tabloid", 11., 17.,
236 +                "A3", 11.69, 16.54,
237 +                "A4", 8.27, 11.69,
238 +                "A4small", 7.47, 10.85,
239 +                "A5", 6.00, 8.27,
240 +                "A6", 4.13, 6.00,
241 +                "B4", 10.12, 14.33,
242 +                "B5", 7.17, 10.12,
243 +                "C5", 6.38, 9.01,
244 +                "C6", 4.49, 6.38,
245 +                "DL", 4.33, 8.66,
246 +                "hagaki", 3.94, 5.83,
247 +                "" };
248 +        register struct psize   *pp;
249 +        register char   *s = ps;
250 +        double  d;
251 +
252 +        if (isdigit(*s)) {              /* check for WWxHH specification */
253 +                width = atof(s);
254 +                while (*s && !isalpha(*s))
255 +                        s++;
256 +                d = unit2inch(s);
257 +                height = atof(++s);
258 +                width *= d;
259 +                height *= d;
260 +                if ((width >= 1.) & (height >= 1.))
261 +                        return;
262 +        } else                          /* check for match to standard size */
263 +                for (pp = p; pp->n[0]; pp++)
264 +                        if (matchid(s, pp->n)) {
265 +                                width = pp->w;
266 +                                height = pp->h;
267 +                                return;
268 +                        }
269 +        fprintf(stderr, "%s: unknown paper size \"%s\" -- known sizes:\n",
270 +                        progname, ps);
271 +        fprintf(stderr, "_Name________Width_Height_(inches)\n");
272 +        for (pp = p; pp->n[0]; pp++)
273 +                fprintf(stderr, "%-11s  %5.2f  %5.2f\n", pp->n, pp->w, pp->h);
274 +        fprintf(stderr, "Or use WWxHH size specification\n");
275 +        exit(1);
276 + }
277 +
278 +
279   quiterr(err)            /* print message and exit */
280   char  *err;
281   {
# Line 133 | Line 287 | char  *err;
287   }
288  
289  
290 < PSheader(name)                  /* print PostScript header */
291 < char  *name;
290 > PSheader(ac, av)                /* print PostScript header */
291 > int  ac;
292 > char  **av;
293   {
294 <        int  landscape = 0;
294 >        char  *rstr;
295 >        int  landscape, rotate, n;
296          double  pwidth, pheight;
297          double  iwidth, iheight;
298                                          /* EPS comments */
299          puts("%!PS-Adobe-2.0 EPSF-2.0");
300 <        printf("%%%%Title: %s\n", name);
301 <        printf("%%%%Creator: %s = %s\n", progname, SCCSid);
300 >        printf("%%%%Title: "); printargs(ac, av, stdout);
301 >        printf("%%%%Creator: %s\n", progname);
302          printf("%%%%Pages: %d\n", ncopies);
303 <        if (landscape = xmax > pixaspect*ymax)
303 >        if ( (landscape = xmax > pixaspect*ymax) )
304                  puts("%%Orientation: Landscape");
305          else
306                  puts("%%Orientation: Portrait");
307 <        if (PWIDTH > PHEIGHT ^ landscape) {
307 >        if ( (rotate = (PWIDTH > PHEIGHT) ^ landscape) ) {
308                  pwidth = PHEIGHT;
309                  pheight = PWIDTH;
310          } else {
311                  pwidth = PWIDTH;
312                  pheight = PHEIGHT;
313          }
314 <        if (pheight/pwidth > pixaspect*ymax/xmax) {
315 <                iwidth = pwidth;
316 <                iheight = pwidth*pixaspect*ymax/xmax;
317 <        } else {
318 <                iheight = pheight;
319 <                iwidth = pheight*xmax/(pixaspect*ymax);
320 <        }
321 <        if (pwidth == PHEIGHT)
314 >        if (dpi > 100 && (pixaspect >= 0.99) & (pixaspect <= 1.01))
315 >                if (pheight/pwidth > ymax/xmax) {
316 >                        n = pwidth*dpi/xmax;    /* floor */
317 >                        iwidth = n > 0 ? (double)(n*xmax)/dpi : pwidth;
318 >                        iheight = iwidth*ymax/xmax;
319 >                } else {
320 >                        n = pheight*dpi/ymax;   /* floor */
321 >                        iheight = n > 0 ? (double)(n*ymax)/dpi : pheight;
322 >                        iwidth = iheight*xmax/ymax;
323 >                }
324 >        else
325 >                if (pheight/pwidth > pixaspect*ymax/xmax) {
326 >                        iwidth = pwidth;
327 >                        iheight = iwidth*pixaspect*ymax/xmax;
328 >                } else {
329 >                        iheight = pheight;
330 >                        iwidth = iheight*xmax/(pixaspect*ymax);
331 >                }
332 >        if (rotate)
333                  printf("%%%%BoundingBox: %.0f %.0f %.0f %.0f\n",
334                                  HMARGIN+(pheight-iheight)*.5,
335                                  VMARGIN+(pwidth-iwidth)*.5,
# Line 176 | Line 343 | char  *name;
343                                  VMARGIN+(pheight-iheight)*.5+iheight);
344          puts("%%EndComments");
345          puts("gsave save");
346 <        puts("64 dict begin");
346 >        puts("17 dict begin");
347                                          /* define image reader */
348          if (docolor) {
349                  printf("/redline %d string def\n", xmax);
350                  printf("/grnline %d string def\n", xmax);
351                  printf("/bluline %d string def\n", xmax);
352 <                printf("/rgbline %d string def\n", 3*xmax);
353 <                PSprocdef("read6bitRLE", "interleave");
354 <        } else {
355 <                printf("/greyline %d string def\n", xmax);
356 <                PSprocdef("read6bitRLE", NULL);
190 <        }
352 >        } else
353 >                printf("/gryline %d string def\n", xmax);
354 >                                        /* use compressed encoding? */
355 >        if (putprim == Cputprim)
356 >                PSprocdef("read6bitRLE");
357                                          /* set up transformation matrix */
358          printf("%f %f translate\n", HMARGIN, VMARGIN);
359 <        if (pwidth == PHEIGHT) {
359 >        if (rotate) {
360                  printf("%f 0 translate\n", PWIDTH);
361                  puts("90 rotate");
362          }
# Line 198 | Line 364 | char  *name;
364          printf("%f %f scale\n", iwidth, iheight);
365          puts("%%EndProlog");
366                                          /* start image procedure */
367 <        printf("%d %d 8 [%d 0 0 %d 0 %d] ", xmax, ymax, xmax, -ymax, ymax);
368 <        if (docolor) {
369 <                puts("{redline read6bitRLE grnline read6bitRLE");
370 <                puts("\tbluline read6bitRLE rgbline interleave} false 3 colorimage");
371 <        } else
372 <                puts("{greyline read6bitRLE} image");
367 >        printf("%d %d 8 [%d 0 0 %d 0 %d]\n", xmax, ymax, xmax, -ymax, ymax);
368 >        if (putprim == Cputprim) {
369 >                if (docolor) {
370 >                        puts("{redline read6bitRLE}");
371 >                        puts("{grnline read6bitRLE}");
372 >                        puts("{bluline read6bitRLE}");
373 >                        puts("true 3 colorimage");
374 >                } else
375 >                        puts("{gryline read6bitRLE} image");
376 >        } else {
377 >                rstr = putprim==Aputprim ? "readhexstring" : "readstring";
378 >                if (docolor) {
379 >                        printf("{currentfile redline %s pop}\n", rstr);
380 >                        printf("{currentfile grnline %s pop}\n", rstr);
381 >                        printf("{currentfile bluline %s pop}\n", rstr);
382 >                        puts("true 3 colorimage");
383 >                } else
384 >                        printf("{currentfile gryline %s pop} image\n", rstr);
385 >        }
386   }
387  
388  
# Line 219 | Line 398 | PStrailer()                    /* print PostScript trailer */
398   }
399  
400  
401 < PSprocdef(nam, inam)            /* define PS procedure to read image */
402 < char  *nam, *inam;
401 > PSprocdef(nam)                  /* define PS procedure to read image */
402 > char  *nam;
403   {
404          short  itab[128];
405          register int  i;
# Line 228 | Line 407 | char  *nam, *inam;
407          for (i = 0; i < 128; i++)       /* clear */
408                  itab[i] = -1;
409          for (i = 1; i < 63; i++)        /* assign greys */
410 <                itab[code[i]] = 256.0*pow(((i<<2)+2.5)/256.0, GAMVAL);
410 >                itab[code[i]] = 256.0*pow((i+.5)/64.0, CODE6GAM/devgam);
411          itab[code[0]] = 0;              /* black is black */
412          itab[code[63]] = 255;           /* and white is white */
413          printf("/codetab [");
# Line 256 | Line 435 | char  *nam, *inam;
435          printf("\t\t} for\n");
436          printf("\t} stopped {pop pop 0 string} {scanline} ifelse\n");
437          printf("} bind def\n");
259 –        if (inam == NULL)
260 –                return;
261 –                                        /* define interleaving procedure */
262 –        printf("/%s {\t%% redscn grnscn bluscn rgbscn\n", inam);
263 –        printf("\t/rgbscan exch def /bscan exch def\n");
264 –        printf("\t/gscan exch def /rscan exch def\n");
265 –        printf("\trscan length %d eq gscan length %d eq and ", xmax, xmax);
266 –        printf("bscan length %d eq and \n", xmax);
267 –        printf("\t{ 0 1 %d { /ndx exch def\n", xmax-1);
268 –        printf("\t\trgbscan ndx 3 mul rscan ndx get put\n");
269 –        printf("\t\trgbscan ndx 3 mul 1 add gscan ndx get put\n");
270 –        printf("\t\trgbscan ndx 3 mul 2 add bscan ndx get put\n");
271 –        printf("\t\t} for rgbscan }\n");
272 –        printf("\t{0 string} ifelse\n");
273 –        printf("} bind def\n");
438   }
439  
440  
# Line 286 | Line 450 | ra2ps()                                /* convert Radiance scanlines to 6-bit */
450          for (y = ymax-1; y >= 0; y--) {
451                  if (freadcolrs(scanin, xmax, stdin) < 0)
452                          quiterr("error reading Radiance picture");
453 <                if (bradj)                      /* adjust exposure */
454 <                        shiftcolrs(scanin, xmax, bradj);
455 <                colrs_gambs(scanin, xmax);      /* gamma compression */
453 >                if (putprim == Cputprim || devgam != 1.) {
454 >                        if (bradj)                      /* adjust exposure */
455 >                                shiftcolrs(scanin, xmax, bradj);
456 >                        colrs_gambs(scanin, xmax);      /* gamma compression */
457 >                } else
458 >                        normcolrs(scanin, xmax, bradj);
459                  if (docolor) {
460 <                        putprim(scanin, RED);
461 <                        putprim(scanin, GRN);
462 <                        putprim(scanin, BLU);
460 >                        (*putprim)(scanin, RED);
461 >                        (*putprim)(scanin, GRN);
462 >                        (*putprim)(scanin, BLU);
463                  } else
464 <                        putprim(scanin, GRY);
464 >                        (*putprim)(scanin, GRY);
465                  if (ferror(stdout))
466                          quiterr("error writing PostScript file");
467          }
468          putchar('\n');
469                                                  /* free scanline */
470 <        free((char *)scanin);
470 >        free((void *)scanin);
471   }
472  
473  
474 < putprim(scn, pri)               /* put out one primary from scanline */
474 > int
475 > Aputprim(scn, pri)              /* put out hex ASCII primary from scanline */
476 > COLR    *scn;
477 > int     pri;
478 > {
479 >        static char     hexdigit[] = "0123456789ABCDEF";
480 >        static int      col = 0;
481 >        register int    x, c;
482 >
483 >        for (x = 0; x < xmax; x++) {
484 >                if (pri == GRY)
485 >                        c = normbright(scn[x]);
486 >                else
487 >                        c = scn[x][pri];
488 >                if (c > 255) c = 255;
489 >                putchar(hexdigit[c>>4]);
490 >                putchar(hexdigit[c&0xf]);
491 >                if ((col += 2) >= 72) {
492 >                        putchar('\n');
493 >                        col = 0;
494 >                }
495 >        }
496 > }
497 >
498 >
499 > int
500 > Bputprim(scn, pri)              /* put out binary primary from scanline */
501 > COLR    *scn;
502 > int     pri;
503 > {
504 >        register int    x, c;
505 >
506 >        for (x = 0; x < xmax; x++) {
507 >                if (pri == GRY)
508 >                        c = normbright(scn[x]);
509 >                else
510 >                        c = scn[x][pri];
511 >                if (c > 255) c = 255;
512 >                putchar(c);
513 >        }
514 > }
515 >
516 >
517 > int
518 > Cputprim(scn, pri)              /* put out compressed primary from scanline */
519   COLR    *scn;
520   int     pri;
521   {

Diff Legend

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