src/px/ra_ps.c

/* Copyright (c) 1992 Regents of the University of California */

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/*
 *  Radiance picture to PostScript file translator -- one way!
 */

#include  <stdio.h>
#ifdef MSDOS
#include  <fcntl.h>
#endif
#include  "color.h"

#define HMARGIN         (.5*72)                 /* horizontal margin */
#define VMARGIN         (.5*72)                 /* vertical margin */
#define PWIDTH          (8.5*72-2*HMARGIN)      /* width of device */
#define PHEIGHT         (11*72-2*VMARGIN)       /* height of device */

#define RUNCHR          '*'                     /* character to start rle */
#define MINRUN          4                       /* minimum run-length */
#define RSTRT           '!'                     /* character for MINRUN */
#define MAXRUN          (MINRUN+'~'-RSTRT)      /* maximum run-length */

char  code[] =                  /* 6-bit code lookup table */
        "0123456789abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ@+";

int  wrongformat = 0;                   /* input in wrong format? */
double  pixaspect = 1.0;                /* pixel aspect ratio */

int  bradj = 0;                         /* brightness adjustment */
int  ncopies = 1;                       /* number of copies */

char  *progname;

int  xmax, ymax;

extern char  *malloc();


headline(s)             /* check header line */
char  *s;
{
        char  fmt[32];

        if (isformat(s)) {
                formatval(fmt, s);
                wrongformat = strcmp(fmt, COLRFMT);
        } else if (isaspect(s))
                pixaspect *= aspectval(s);
}


main(argc, argv)
int  argc;
char  *argv[];
{
        int  i;
        
        progname = argv[0];

        for (i = 1; i < argc; i++)
                if (argv[i][0] == '-')
                        switch (argv[i][1]) {
                        case 'e':               /* exposure adjustment */
                                if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
                                        goto userr;
                                bradj = atoi(argv[++i]);
                                break;
                        case 'n':               /* number of copies */
                                ncopies = atoi(argv[++i]);
                                break;
                        default:
                                goto userr;
                        }
                else
                        break;

        if (i < argc-2)
                goto userr;
        if (i <= argc-1 && freopen(argv[i], "r", stdin) == NULL) {
                fprintf(stderr, "%s: can't open input \"%s\"\n",
                                progname, argv[i]);
                exit(1);
        }
        if (i == argc-2 && freopen(argv[i+1], "w", stdout) == NULL) {
                fprintf(stderr, "%s: can't open output \"%s\"\n",
                                progname, argv[i+1]);
                exit(1);
        }
#ifdef MSDOS
        setmode(fileno(stdin), O_BINARY);
#endif
                                /* get our header */
        getheader(stdin, headline, NULL);
        if (wrongformat || fgetresolu(&xmax, &ymax, stdin) < 0)
                quiterr("bad picture format");
                                /* write header */
        PSheader(i <= argc-1 ? argv[i] : "<stdin>");
                                /* convert file */
        ra2ps();
                                /* write trailer */
        PStrailer();
        exit(0);
userr:
        fprintf(stderr, "Usage: %s [-e +/-stops] [input [output]]\n", progname);
        exit(1);
}


quiterr(err)            /* print message and exit */
char  *err;
{
        if (err != NULL) {
                fprintf(stderr, "%s: %s\n", progname, err);
                exit(1);
        }
        exit(0);
}


PSheader(name)                  /* print PostScript header */
char  *name;
{
        int  landscape = 0;
        double  pwidth, pheight;
        double  iwidth, iheight;

        printf("%%!\n");
        printf("%%%%Title: %s\n", name);
        printf("%%%%Creator: %s\n", progname);
        printf("%%%%Pages: %d\n", ncopies);
        if (landscape = xmax > pixaspect*ymax)
                printf("%%%%Landscape\n");
        else
                printf("%%%%Portrait\n");
        printf("%%%%EndComments\n");
        printf("64 dict begin\n");
                                        /* define image reader */
        PSprocdef("read6bitRLE");
                                        /* set up transformation matrix */
        printf("%f %f translate\n", HMARGIN, VMARGIN);
        if (PWIDTH > PHEIGHT ^ landscape) {
                printf("0 %f translate\n", PHEIGHT);
                printf("-90 rotate\n");
                pwidth = PHEIGHT;
                pheight = PWIDTH;
        } else {
                pwidth = PWIDTH;
                pheight = PHEIGHT;
        }
        if (pheight/pwidth > pixaspect*ymax/xmax) {
                iwidth = pwidth;
                iheight = pwidth*pixaspect*ymax/xmax;
        } else {
                iheight = pheight;
                iwidth = pheight*xmax/(pixaspect*ymax);
        }
        printf("%f %f translate\n", (pwidth-iwidth)*.5, (pheight-iheight)*.5);
        printf("%f %f scale\n", iwidth, iheight);
        printf("%%%%EndProlog\n");
                                        /* start image procedure */
        printf("%d %d 8 [%d 0 0 %d 0 %d] {read6bitRLE} image\n", xmax, ymax,
                        xmax, -ymax, ymax);
}


PStrailer()                     /* print PostScript trailer */
{
        puts("%%Trailer");
        if (ncopies > 1)
                printf("/#copies %d def\n", ncopies);
        puts("showpage");
        puts("end");
        puts("%%EOF");
}


PSprocdef(nam)                  /* define PS procedure to read image */
char  *nam;
{
        short  itab[128];
        register int  i;
                                /* assign code values */
        for (i = 0; i < 128; i++)       /* clear */
                itab[i] = -1;
        for (i = 1; i < 63; i++)        /* assign greys */
                itab[code[i]] = i<<2 | 2;
        itab[code[0]] = 0;              /* black is black */
        itab[code[63]] = 255;           /* and white is white */
        printf("/codetab [");
        for (i = 0; i < 128; i++) {
                if (!(i & 0xf))
                        putchar('\n');
                printf(" %3d", itab[i]);
        }
        printf("\n] def\n");
        printf("/scanline %d string def\n", xmax);
        printf("/nrept 0 def\n");
        printf("/readbyte { currentfile read not {stop} if } def\n");
        printf("/decode { codetab exch get } def\n");
        printf("/%s {\n", nam);
        printf("\t{ 0 1 %d { scanline exch\n", xmax-1);
        printf("\t\tnrept 0 le\n");
        printf("\t\t\t{ { readbyte dup %d eq\n", RUNCHR);
        printf("\t\t\t\t\t{ pop /nrept readbyte %d sub def\n", RSTRT-MINRUN+1);
        printf("\t\t\t\t\t\t/reptv readbyte decode def\n");
        printf("\t\t\t\t\t\treptv exit }\n");
        printf("\t\t\t\t\t{ decode dup 0 lt {pop} {exit} ifelse }\n");
        printf("\t\t\t\tifelse } loop }\n");
        printf("\t\t\t{ /nrept nrept 1 sub def reptv }\n");
        printf("\t\tifelse put\n");
        printf("\t\t} for\n");
        printf("\t} stopped {pop pop 0 string} {scanline} ifelse\n");
        printf("} bind def\n");
}


ra2ps()                         /* convert Radiance scanlines to 6-bit */
{
        register COLR   *scanin;
        register int    c;
        int     lastc, cnt;
        register int    x;
        int     y;
                                                /* allocate scanline */
        scanin = (COLR *)malloc(xmax*sizeof(COLR));
        if (scanin == NULL)
                quiterr("out of memory in ra2ps");
                                                /* convert image */
        for (y = ymax-1; y >= 0; y--) {
                if (freadcolrs(scanin, xmax, stdin) < 0)
                        quiterr("error reading Radiance picture");
                normcolrs(scanin, xmax, bradj); /* normalize */
                lastc = -1; cnt = 0;
                for (x = 0; x < xmax; x++) {
                        c = normbright(scanin[x]) + 2;
                        if (c > 255) c = 255;
                        c = code[c>>2];
                        if (c == lastc && cnt < MAXRUN)
                                cnt++;
                        else {
                                putrle(cnt, lastc);
                                lastc = c;
                                cnt = 1;
                        }
                }
                putrle(cnt, lastc);
                if (ferror(stdout))
                        quiterr("error writing PostScript file");
        }
        putchar('\n');
                                                /* free scanline */
        free((char *)scanin);
}


putrle(cnt, cod)                /* put out cnt of cod */
register int    cnt, cod;
{
        static int      col = 0;

        if (cnt >= MINRUN) {
                col += 3;
                putchar(RUNCHR);
                putchar(RSTRT-MINRUN+cnt);
                putchar(cod);
        } else {
                col += cnt;
                while (cnt-- > 0)
                        putchar(cod);
        }
        if (col >= 72) {
                putchar('\n');
                col = 0;
        }
}
Revision:	2.9
Committed:	Thu Feb 10 21:59:53 1994 UTC (30 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.8:	+57 -13 lines
Log Message:	added run-length compression
#	User	Rev	Content
1	greg	2.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			* Radiance picture to PostScript file translator -- one way!
9			*/
10
11			#include <stdio.h>
12	greg	2.6	#ifdef MSDOS
13			#include <fcntl.h>
14			#endif
15	greg	2.1	#include "color.h"
16
17			#define HMARGIN (.572) / horizontal margin */
18			#define VMARGIN (.572) / vertical margin */
19			#define PWIDTH (8.572-2HMARGIN) /* width of device */
20			#define PHEIGHT (1172-2VMARGIN) /* height of device */
21
22	greg	2.9	#define RUNCHR '' / character to start rle */
23			#define MINRUN 4 /* minimum run-length */
24			#define RSTRT '!' /* character for MINRUN */
25			#define MAXRUN (MINRUN+'~'-RSTRT) /* maximum run-length */
26
27	greg	2.1	char code[] = /* 6-bit code lookup table */
28			"0123456789abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ@+";
29
30			int wrongformat = 0; /* input in wrong format? */
31	greg	2.6	double pixaspect = 1.0; /* pixel aspect ratio */
32	greg	2.1
33			int bradj = 0; /* brightness adjustment */
34	greg	2.4	int ncopies = 1; /* number of copies */
35	greg	2.1
36			char *progname;
37
38			int xmax, ymax;
39
40	greg	2.6	extern char *malloc();
41	greg	2.1
42	greg	2.6
43	greg	2.1	headline(s) /* check header line */
44			char *s;
45			{
46			char fmt[32];
47
48			if (isformat(s)) {
49			formatval(fmt, s);
50			wrongformat = strcmp(fmt, COLRFMT);
51			} else if (isaspect(s))
52			pixaspect *= aspectval(s);
53			}
54
55
56			main(argc, argv)
57			int argc;
58			char *argv[];
59			{
60			int i;
61
62			progname = argv[0];
63
64			for (i = 1; i < argc; i++)
65			if (argv[i][0] == '-')
66			switch (argv[i][1]) {
67			case 'e': /* exposure adjustment */
68			if (argv[i+1][0] != '+' && argv[i+1][0] != '-')
69			goto userr;
70			bradj = atoi(argv[++i]);
71			break;
72	greg	2.4	case 'n': /* number of copies */
73			ncopies = atoi(argv[++i]);
74			break;
75	greg	2.1	default:
76			goto userr;
77			}
78			else
79			break;
80
81			if (i < argc-2)
82			goto userr;
83			if (i <= argc-1 && freopen(argv[i], "r", stdin) == NULL) {
84			fprintf(stderr, "%s: can't open input \"%s\"\n",
85			progname, argv[i]);
86			exit(1);
87			}
88			if (i == argc-2 && freopen(argv[i+1], "w", stdout) == NULL) {
89	greg	2.7	fprintf(stderr, "%s: can't open output \"%s\"\n",
90	greg	2.1	progname, argv[i+1]);
91			exit(1);
92			}
93	greg	2.6	#ifdef MSDOS
94			setmode(fileno(stdin), O_BINARY);
95			#endif
96	greg	2.1	/* get our header */
97			getheader(stdin, headline, NULL);
98			if (wrongformat \|\| fgetresolu(&xmax, &ymax, stdin) < 0)
99			quiterr("bad picture format");
100			/* write header */
101			PSheader(i <= argc-1 ? argv[i] : "<stdin>");
102			/* convert file */
103			ra2ps();
104			/* write trailer */
105			PStrailer();
106			exit(0);
107			userr:
108			fprintf(stderr, "Usage: %s [-e +/-stops] [input [output]]\n", progname);
109			exit(1);
110			}
111
112
113			quiterr(err) /* print message and exit */
114			char *err;
115			{
116			if (err != NULL) {
117			fprintf(stderr, "%s: %s\n", progname, err);
118			exit(1);
119			}
120			exit(0);
121			}
122
123
124			PSheader(name) /* print PostScript header */
125			char *name;
126			{
127			int landscape = 0;
128	greg	2.6	double pwidth, pheight;
129			double iwidth, iheight;
130	greg	2.1
131			printf("%%!\n");
132			printf("%%%%Title: %s\n", name);
133			printf("%%%%Creator: %s\n", progname);
134	greg	2.4	printf("%%%%Pages: %d\n", ncopies);
135	greg	2.1	if (landscape = xmax > pixaspect*ymax)
136			printf("%%%%Landscape\n");
137			else
138			printf("%%%%Portrait\n");
139			printf("%%%%EndComments\n");
140			printf("64 dict begin\n");
141	greg	2.3	/* define image reader */
142	greg	2.9	PSprocdef("read6bitRLE");
143	greg	2.1	/* set up transformation matrix */
144	greg	2.2	printf("%f %f translate\n", HMARGIN, VMARGIN);
145	greg	2.1	if (PWIDTH > PHEIGHT ^ landscape) {
146			printf("0 %f translate\n", PHEIGHT);
147	greg	2.2	printf("-90 rotate\n");
148	greg	2.1	pwidth = PHEIGHT;
149			pheight = PWIDTH;
150			} else {
151			pwidth = PWIDTH;
152			pheight = PHEIGHT;
153			}
154			if (pheight/pwidth > pixaspect*ymax/xmax) {
155	greg	2.2	iwidth = pwidth;
156			iheight = pwidthpixaspectymax/xmax;
157	greg	2.1	} else {
158	greg	2.2	iheight = pheight;
159			iwidth = pheightxmax/(pixaspectymax);
160	greg	2.1	}
161	greg	2.2	printf("%f %f translate\n", (pwidth-iwidth).5, (pheight-iheight).5);
162			printf("%f %f scale\n", iwidth, iheight);
163	greg	2.1	printf("%%%%EndProlog\n");
164	greg	2.3	/* start image procedure */
165	greg	2.9	printf("%d %d 8 [%d 0 0 %d 0 %d] {read6bitRLE} image\n", xmax, ymax,
166	greg	2.2	xmax, -ymax, ymax);
167	greg	2.1	}
168
169
170			PStrailer() /* print PostScript trailer */
171			{
172			puts("%%Trailer");
173	greg	2.4	if (ncopies > 1)
174			printf("/#copies %d def\n", ncopies);
175			puts("showpage");
176	greg	2.1	puts("end");
177			puts("%%EOF");
178			}
179
180
181			PSprocdef(nam) /* define PS procedure to read image */
182			char *nam;
183			{
184			short itab[128];
185			register int i;
186	greg	2.5	/* assign code values */
187			for (i = 0; i < 128; i++) /* clear */
188	greg	2.1	itab[i] = -1;
189	greg	2.5	for (i = 1; i < 63; i++) /* assign greys */
190	greg	2.1	itab[code[i]] = i<<2 \| 2;
191	greg	2.5	itab[code[0]] = 0; /* black is black */
192			itab[code[63]] = 255; /* and white is white */
193	greg	2.9	printf("/codetab [");
194	greg	2.1	for (i = 0; i < 128; i++) {
195			if (!(i & 0xf))
196			putchar('\n');
197			printf(" %3d", itab[i]);
198			}
199			printf("\n] def\n");
200			printf("/scanline %d string def\n", xmax);
201	greg	2.9	printf("/nrept 0 def\n");
202			printf("/readbyte { currentfile read not {stop} if } def\n");
203			printf("/decode { codetab exch get } def\n");
204	greg	2.1	printf("/%s {\n", nam);
205			printf("\t{ 0 1 %d { scanline exch\n", xmax-1);
206	greg	2.9	printf("\t\tnrept 0 le\n");
207			printf("\t\t\t{ { readbyte dup %d eq\n", RUNCHR);
208			printf("\t\t\t\t\t{ pop /nrept readbyte %d sub def\n", RSTRT-MINRUN+1);
209			printf("\t\t\t\t\t\t/reptv readbyte decode def\n");
210			printf("\t\t\t\t\t\treptv exit }\n");
211			printf("\t\t\t\t\t{ decode dup 0 lt {pop} {exit} ifelse }\n");
212			printf("\t\t\t\tifelse } loop }\n");
213			printf("\t\t\t{ /nrept nrept 1 sub def reptv }\n");
214			printf("\t\tifelse put\n");
215			printf("\t\t} for\n");
216			printf("\t} stopped {pop pop 0 string} {scanline} ifelse\n");
217	greg	2.8	printf("} bind def\n");
218	greg	2.1	}
219
220
221			ra2ps() /* convert Radiance scanlines to 6-bit */
222			{
223	greg	2.9	register COLR *scanin;
224	greg	2.1	register int c;
225	greg	2.9	int lastc, cnt;
226	greg	2.1	register int x;
227			int y;
228			/* allocate scanline */
229			scanin = (COLR )malloc(xmaxsizeof(COLR));
230			if (scanin == NULL)
231			quiterr("out of memory in ra2ps");
232			/* convert image */
233			for (y = ymax-1; y >= 0; y--) {
234			if (freadcolrs(scanin, xmax, stdin) < 0)
235			quiterr("error reading Radiance picture");
236	greg	2.6	normcolrs(scanin, xmax, bradj); /* normalize */
237	greg	2.9	lastc = -1; cnt = 0;
238	greg	2.1	for (x = 0; x < xmax; x++) {
239	greg	2.9	c = normbright(scanin[x]) + 2;
240	greg	2.1	if (c > 255) c = 255;
241	greg	2.9	c = code[c>>2];
242			if (c == lastc && cnt < MAXRUN)
243			cnt++;
244			else {
245			putrle(cnt, lastc);
246			lastc = c;
247			cnt = 1;
248			}
249	greg	2.1	}
250	greg	2.9	putrle(cnt, lastc);
251	greg	2.1	if (ferror(stdout))
252			quiterr("error writing PostScript file");
253			}
254			putchar('\n');
255			/* free scanline */
256			free((char *)scanin);
257	greg	2.9	}
258
259
260			putrle(cnt, cod) /* put out cnt of cod */
261			register int cnt, cod;
262			{
263			static int col = 0;
264
265			if (cnt >= MINRUN) {
266			col += 3;
267			putchar(RUNCHR);
268			putchar(RSTRT-MINRUN+cnt);
269			putchar(cod);
270			} else {
271			col += cnt;
272			while (cnt-- > 0)
273			putchar(cod);
274			}
275			if (col >= 72) {
276			putchar('\n');
277			col = 0;
278			}
279	greg	2.1	}