src/px/pcomb.c

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

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

/*
 *  Combine picture files according to calcomp functions.
 *
 *      1/4/89
 */

#include <stdio.h>

#include <errno.h>

#include "color.h"

#include "calcomp.h"

#define MAXINP          16              /* maximum number of input files */

struct {
        char    *name;          /* file name */
        FILE    *fp;            /* stream pointer */
        COLOR   *scan;          /* input scanline */
}       input[MAXINP];                  /* input pictures */

int     nfiles;                         /* number of input files */

char    *vcolin[3] = {"ri", "gi", "bi"};
char    *vcolout[3] = {"ro", "go", "bo"};

#define vxpos           "x"
#define vypos           "y"

int     nowarn = 0;                     /* no warning messages? */

int     xres=0, yres=0;                 /* picture resolution */

int     xpos, ypos;                     /* picture position */


tputs(s)                        /* put out string preceded by a tab */
char    *s;
{
        putchar('\t');
        fputs(s, stdout);
}


main(argc, argv)
int     argc;
char    *argv[];
{
        extern double   l_redin(), l_grnin(), l_bluin();
        int     a;
        
        funset(vcolin[RED], 1, l_redin);
        funset(vcolin[GRN], 1, l_grnin);
        funset(vcolin[BLU], 1, l_bluin);
        
        for (a = 1; a < argc; a++)
                if (argv[a][0] == '-')
                        switch (argv[a][1]) {
                        case '\0':
                                goto getfiles;
                        case 'x':
                                xres = atoi(argv[++a]);
                                break;
                        case 'y':
                                yres = atoi(argv[++a]);
                                break;
                        case 'w':
                                nowarn = !nowarn;
                                break;
                        case 'f':
                                fcompile(argv[++a]);
                                break;
                        case 'e':
                                scompile(NULL, argv[++a]);
                                break;
                        default:
                                eputs("Usage: ");
                                eputs(argv[0]);
        eputs(" [-w][-x xres][-y yres][-e expr][-f file] [picture ..]\n");
                                quit(1);
                        }
                else
                        break;
getfiles:
        nfiles = 0;
        for ( ; a < argc; a++) {
                if (nfiles >= MAXINP) {
                        eputs(argv[0]);
                        eputs(": too many picture files\n");
                        quit(1);
                }
                if (argv[a][0] == '-') {
                        input[nfiles].name = "<stdin>";
                        input[nfiles].fp = stdin;
                } else {
                        input[nfiles].name = argv[a];
                        input[nfiles].fp = fopen(argv[a], "r");
                        if (input[nfiles].fp == NULL) {
                                eputs(argv[a]);
                                eputs(": cannot open\n");
                                quit(1);
                        }
                }
                fputs(input[nfiles].name, stdout);
                fputs(":\n", stdout);
                getheader(input[nfiles].fp, tputs);
                if (fscanf(input[nfiles].fp, "-Y %d +X %d\n", &ypos, &xpos) != 2) {
                        eputs(input[nfiles].name);
                        eputs(": bad picture size\n");
                        quit(1);
                }
                if (xres == 0 && yres == 0) {
                        xres = xpos;
                        yres = ypos;
                } else if (xpos != xres || ypos != yres) {
                        eputs(argv[0]);
                        eputs(": resolution mismatch\n");
                        quit(1);
                }
                input[nfiles].scan = (COLOR *)emalloc(xres*sizeof(COLOR));
                nfiles++;
        }
        printargs(argc, argv, stdout);
        putchar('\n');
        printf("-Y %d +X %d\n", yres, xres);
        combine();
        quit(0);
}


combine()                       /* combine pictures */
{
        int     coldef[3];
        COLOR   *scanout;
        register int    i, j;
                                                /* check defined variables */
        for (j = 0; j < 3; j++)
                coldef[j] = vardefined(vcolout[j]);
                                                /* allocate scanline */
        scanout = (COLOR *)emalloc(xres*sizeof(COLOR));
                                                /* combine files */
        for (ypos = yres-1; ypos >= 0; ypos--) {
                for (i = 0; i < nfiles; i++)
                        if (freadscan(input[i].scan, xres, input[i].fp) < 0) {
                                eputs(input[i].name);
                                eputs(": read error\n");
                                quit(1);
                        }
                varset(vypos, (double)ypos);
                for (xpos = 0; xpos < xres; xpos++) {
                        varset(vxpos, (double)xpos);
                        eclock++;
                        for (j = 0; j < 3; j++)
                                if (coldef[j]) {
                                        colval(scanout[xpos],j) = varvalue(vcolout[j]);
                                        if (colval(scanout[xpos],j) < 0.0)
                                                colval(scanout[xpos],j) = 0.0;
                                } else {
                                        colval(scanout[xpos],j) = 0.0;
                                        for (i = 0; i < nfiles; i++)
                                                colval(scanout[xpos],j) += colval(input[i].scan[xpos],j);
                                        colval(scanout[xpos],j) /= (double)nfiles;
                                }
                }
                if (fwritescan(scanout, xres, stdout) < 0) {
                        eputs("write error\n");
                        quit(1);
                }
        }
        efree(scanout);
}


double
colin(fn, ci)                   /* return color value for picture */
register int    fn, ci;
{
        if (fn == 0)
                return((double)nfiles);
        if (fn < 1 || fn > nfiles) {
                errno = EDOM;
                return(0.0);
        }
        return(colval(input[fn-1].scan[xpos],ci));
}


double
l_redin()                       /* get red color */
{
        return(colin((int)(argument(1)+.5), RED));
}


double
l_grnin()                       /* get green color */
{
        return(colin((int)(argument(1)+.5), GRN));
}


double
l_bluin()                       /* get blue color */
{
        return(colin((int)(argument(1)+.5), BLU));
}


wputs(msg)
char    *msg;
{
        if (!nowarn)
                eputs(msg);
}


eputs(msg)
char    *msg;
{
        fputs(msg, stderr);
}


quit(code)
int     code;
{
        exit(code);
}
Revision:	1.1
Committed:	Thu Feb 2 10:49:22 1989 UTC (36 years, 9 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Log Message:	Initial revision
#	Content
1	/* Copyright (c) 1989 Regents of the University of California */
2
3	#ifndef lint
4	static char SCCSid[] = "$SunId$ LBL";
5	#endif
6
7	/*
8	* Combine picture files according to calcomp functions.
9	*
10	* 1/4/89
11	*/
12
13	#include <stdio.h>
14
15	#include <errno.h>
16
17	#include "color.h"
18
19	#include "calcomp.h"
20
21	#define MAXINP 16 /* maximum number of input files */
22
23	struct {
24	char name; / file name */
25	FILE fp; / stream pointer */
26	COLOR scan; / input scanline */
27	} input[MAXINP]; /* input pictures */
28
29	int nfiles; /* number of input files */
30
31	char *vcolin[3] = {"ri", "gi", "bi"};
32	char *vcolout[3] = {"ro", "go", "bo"};
33
34	#define vxpos "x"
35	#define vypos "y"
36
37	int nowarn = 0; /* no warning messages? */
38
39	int xres=0, yres=0; /* picture resolution */
40
41	int xpos, ypos; /* picture position */
42
43
44	tputs(s) /* put out string preceded by a tab */
45	char *s;
46	{
47	putchar('\t');
48	fputs(s, stdout);
49	}
50
51
52	main(argc, argv)
53	int argc;
54	char *argv[];
55	{
56	extern double l_redin(), l_grnin(), l_bluin();
57	int a;
58
59	funset(vcolin[RED], 1, l_redin);
60	funset(vcolin[GRN], 1, l_grnin);
61	funset(vcolin[BLU], 1, l_bluin);
62
63	for (a = 1; a < argc; a++)
64	if (argv[a][0] == '-')
65	switch (argv[a][1]) {
66	case '\0':
67	goto getfiles;
68	case 'x':
69	xres = atoi(argv[++a]);
70	break;
71	case 'y':
72	yres = atoi(argv[++a]);
73	break;
74	case 'w':
75	nowarn = !nowarn;
76	break;
77	case 'f':
78	fcompile(argv[++a]);
79	break;
80	case 'e':
81	scompile(NULL, argv[++a]);
82	break;
83	default:
84	eputs("Usage: ");
85	eputs(argv[0]);
86	eputs(" [-w][-x xres][-y yres][-e expr][-f file] [picture ..]\n");
87	quit(1);
88	}
89	else
90	break;
91	getfiles:
92	nfiles = 0;
93	for ( ; a < argc; a++) {
94	if (nfiles >= MAXINP) {
95	eputs(argv[0]);
96	eputs(": too many picture files\n");
97	quit(1);
98	}
99	if (argv[a][0] == '-') {
100	input[nfiles].name = "<stdin>";
101	input[nfiles].fp = stdin;
102	} else {
103	input[nfiles].name = argv[a];
104	input[nfiles].fp = fopen(argv[a], "r");
105	if (input[nfiles].fp == NULL) {
106	eputs(argv[a]);
107	eputs(": cannot open\n");
108	quit(1);
109	}
110	}
111	fputs(input[nfiles].name, stdout);
112	fputs(":\n", stdout);
113	getheader(input[nfiles].fp, tputs);
114	if (fscanf(input[nfiles].fp, "-Y %d +X %d\n", &ypos, &xpos) != 2) {
115	eputs(input[nfiles].name);
116	eputs(": bad picture size\n");
117	quit(1);
118	}
119	if (xres == 0 && yres == 0) {
120	xres = xpos;
121	yres = ypos;
122	} else if (xpos != xres \|\| ypos != yres) {
123	eputs(argv[0]);
124	eputs(": resolution mismatch\n");
125	quit(1);
126	}
127	input[nfiles].scan = (COLOR )emalloc(xressizeof(COLOR));
128	nfiles++;
129	}
130	printargs(argc, argv, stdout);
131	putchar('\n');
132	printf("-Y %d +X %d\n", yres, xres);
133	combine();
134	quit(0);
135	}
136
137
138	combine() /* combine pictures */
139	{
140	int coldef[3];
141	COLOR *scanout;
142	register int i, j;
143	/* check defined variables */
144	for (j = 0; j < 3; j++)
145	coldef[j] = vardefined(vcolout[j]);
146	/* allocate scanline */
147	scanout = (COLOR )emalloc(xressizeof(COLOR));
148	/* combine files */
149	for (ypos = yres-1; ypos >= 0; ypos--) {
150	for (i = 0; i < nfiles; i++)
151	if (freadscan(input[i].scan, xres, input[i].fp) < 0) {
152	eputs(input[i].name);
153	eputs(": read error\n");
154	quit(1);
155	}
156	varset(vypos, (double)ypos);
157	for (xpos = 0; xpos < xres; xpos++) {
158	varset(vxpos, (double)xpos);
159	eclock++;
160	for (j = 0; j < 3; j++)
161	if (coldef[j]) {
162	colval(scanout[xpos],j) = varvalue(vcolout[j]);
163	if (colval(scanout[xpos],j) < 0.0)
164	colval(scanout[xpos],j) = 0.0;
165	} else {
166	colval(scanout[xpos],j) = 0.0;
167	for (i = 0; i < nfiles; i++)
168	colval(scanout[xpos],j) += colval(input[i].scan[xpos],j);
169	colval(scanout[xpos],j) /= (double)nfiles;
170	}
171	}
172	if (fwritescan(scanout, xres, stdout) < 0) {
173	eputs("write error\n");
174	quit(1);
175	}
176	}
177	efree(scanout);
178	}
179
180
181	double
182	colin(fn, ci) /* return color value for picture */
183	register int fn, ci;
184	{
185	if (fn == 0)
186	return((double)nfiles);
187	if (fn < 1 \|\| fn > nfiles) {
188	errno = EDOM;
189	return(0.0);
190	}
191	return(colval(input[fn-1].scan[xpos],ci));
192	}
193
194
195	double
196	l_redin() /* get red color */
197	{
198	return(colin((int)(argument(1)+.5), RED));
199	}
200
201
202	double
203	l_grnin() /* get green color */
204	{
205	return(colin((int)(argument(1)+.5), GRN));
206	}
207
208
209	double
210	l_bluin() /* get blue color */
211	{
212	return(colin((int)(argument(1)+.5), BLU));
213	}
214
215
216	wputs(msg)
217	char *msg;
218	{
219	if (!nowarn)
220	eputs(msg);
221	}
222
223
224	eputs(msg)
225	char *msg;
226	{
227	fputs(msg, stderr);
228	}
229
230
231	quit(code)
232	int code;
233	{
234	exit(code);
235	}