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

Diff Legend

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