ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pcompos.c
Revision: 2.18
Committed: Sun Feb 27 10:17:04 1994 UTC (30 years, 2 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.17: +5 -3 lines
Log Message:
Added new ID to first line of header and changed use of formatval

File Contents

# User Rev Content
1 greg 2.17 /* Copyright (c) 1993 Regents of the University of California */
2 greg 1.1
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * pcompos.c - program to composite pictures.
9     *
10     * 6/30/87
11     */
12    
13     #include <stdio.h>
14    
15 greg 2.16 #include <math.h>
16    
17 greg 2.15 #ifdef MSDOS
18     #include <fcntl.h>
19     #endif
20    
21 greg 1.1 #include "color.h"
22    
23 greg 1.10 #include "resolu.h"
24 greg 1.1
25 greg 2.15 #define MAXFILE 64
26 greg 1.1
27 greg 2.17 #define HASMIN 1
28     #define HASMAX 2
29    
30 greg 1.1 /* output picture size */
31     int xsiz = 0;
32     int ysiz = 0;
33     /* input dimensions */
34     int xmin = 0;
35     int ymin = 0;
36     int xmax = 0;
37     int ymax = 0;
38    
39     COLR bgcolr = BLKCOLR; /* background color */
40    
41 greg 1.9 int labelht = 24; /* label height */
42    
43 greg 1.1 int checkthresh = 0; /* check threshold value */
44    
45 greg 2.14 char Command[] = "<Command>";
46     char Label[] = "<Label>";
47    
48 greg 1.1 char *progname;
49    
50     struct {
51 greg 1.7 char *name; /* file or command name */
52 greg 1.1 FILE *fp; /* stream pointer */
53     int xres, yres; /* picture size */
54     int xloc, yloc; /* anchor point */
55 greg 2.17 int flags; /* HASMIN, HASMAX */
56 greg 1.1 COLR thmin, thmax; /* thresholds */
57     } input[MAXFILE]; /* our input files */
58    
59     int nfile; /* number of files */
60    
61 greg 1.6 int wrongformat = 0;
62 greg 1.1
63 greg 1.9 FILE *popen(), *lblopen();
64 greg 1.6
65 greg 2.15 extern char *malloc();
66 greg 1.7
67 greg 2.15
68 greg 1.1 tabputs(s) /* print line preceded by a tab */
69     char *s;
70     {
71 greg 1.6 char fmt[32];
72    
73 greg 2.18 if (isheadid(s))
74     return;
75     if (formatval(fmt, s))
76 greg 1.6 wrongformat = strcmp(fmt, COLRFMT);
77 greg 2.18 else {
78 greg 1.6 putc('\t', stdout);
79     fputs(s, stdout);
80     }
81 greg 1.1 }
82    
83    
84     main(argc, argv)
85     int argc;
86     char *argv[];
87     {
88 greg 1.7 int ncolumns = 0;
89 greg 1.9 int autolabel = 0;
90 greg 2.11 int curcol = 0, x0 = 0, curx = 0, cury = 0, spacing = 0;
91 greg 2.17 int xsgn, ysgn;
92 greg 1.9 char *thislabel;
93 greg 1.1 int an;
94 greg 2.15 #ifdef MSDOS
95     extern int _fmode;
96     _fmode = O_BINARY;
97     setmode(fileno(stdin), O_BINARY);
98     setmode(fileno(stdout), O_BINARY);
99     #endif
100 greg 1.1 progname = argv[0];
101    
102     for (an = 1; an < argc && argv[an][0] == '-'; an++)
103     switch (argv[an][1]) {
104     case 'x':
105 greg 2.11 xsiz = atoi(argv[++an]);
106 greg 1.1 break;
107     case 'y':
108 greg 2.11 ysiz = atoi(argv[++an]);
109 greg 1.1 break;
110     case 'b':
111     setcolr(bgcolr, atof(argv[an+1]),
112     atof(argv[an+2]),
113     atof(argv[an+3]));
114     an += 3;
115     break;
116 greg 1.7 case 'a':
117     ncolumns = atoi(argv[++an]);
118     break;
119 greg 2.11 case 's':
120     spacing = atoi(argv[++an]);
121     break;
122     case 'o':
123     curx = x0 = atoi(argv[++an]);
124     cury = atoi(argv[++an]);
125     break;
126 greg 1.7 case 'l':
127 greg 1.9 switch (argv[an][2]) {
128     case 'a':
129     autolabel++;
130     break;
131     case 'h':
132     labelht = atoi(argv[++an]);
133     break;
134     case '\0':
135     goto dofiles;
136     default:
137     goto userr;
138     }
139 greg 1.7 break;
140 greg 1.1 case '\0':
141     case 't':
142     goto dofiles;
143     default:
144     goto userr;
145     }
146     dofiles:
147 greg 2.18 newheader("RADIANCE", stdout);
148 greg 1.1 for (nfile = 0; an < argc; nfile++) {
149 greg 1.9 if (nfile >= MAXFILE)
150     goto toomany;
151 greg 2.6 thislabel = NULL;
152 greg 2.17 input[nfile].flags = 0;
153     xsgn = ysgn = '-';
154     while (an < argc && (argv[an][0] == '-' || argv[an][0] == '+'
155     || argv[an][0] == '=')) {
156 greg 1.1 switch (argv[an][1]) {
157     case 't':
158     checkthresh = 1;
159     if (argv[an][0] == '-') {
160 greg 2.17 input[nfile].flags |= HASMIN;
161 greg 1.1 setcolr(input[nfile].thmin,
162     atof(argv[an+1]),
163     atof(argv[an+1]),
164     atof(argv[an+1]));
165 greg 2.17 } else if (argv[an][0] == '+') {
166     input[nfile].flags |= HASMAX;
167 greg 1.1 setcolr(input[nfile].thmax,
168     atof(argv[an+1]),
169     atof(argv[an+1]),
170     atof(argv[an+1]));
171 greg 2.17 } else
172     goto userr;
173     an++;
174 greg 1.1 break;
175 greg 1.9 case 'l':
176     if (strcmp(argv[an], "-l"))
177     goto userr;
178 greg 2.17 thislabel = argv[++an];
179 greg 1.9 break;
180 greg 2.17 case '+':
181     case '-':
182     case '0':
183     if (argv[an][0] != '=')
184     goto userr;
185     xsgn = argv[an][1];
186     ysgn = argv[an][2];
187     if (ysgn != '+' && ysgn != '-' && ysgn != '0')
188     goto userr;
189     break;
190 greg 1.1 case '\0':
191     if (argv[an][0] == '-')
192     goto getfile;
193 greg 1.9 goto userr;
194 greg 1.1 default:
195     goto userr;
196     }
197 greg 2.17 an++;
198     }
199 greg 1.1 getfile:
200 greg 1.7 if (argc-an < (ncolumns ? 1 : 3))
201 greg 1.1 goto userr;
202 greg 2.6 if (autolabel && thislabel == NULL)
203     thislabel = argv[an];
204 greg 1.1 if (!strcmp(argv[an], "-")) {
205     input[nfile].name = "<stdin>";
206     input[nfile].fp = stdin;
207     } else {
208 greg 2.4 if (argv[an][0] == '!') {
209 greg 2.14 input[nfile].name = Command;
210 greg 2.4 input[nfile].fp = popen(argv[an]+1, "r");
211     } else {
212     input[nfile].name = argv[an];
213     input[nfile].fp = fopen(argv[an], "r");
214     }
215     if (input[nfile].fp == NULL) {
216 greg 1.5 perror(argv[an]);
217 greg 1.1 quit(1);
218     }
219     }
220     an++;
221     /* get header */
222     printf("%s:\n", input[nfile].name);
223 greg 1.6 getheader(input[nfile].fp, tabputs, NULL);
224     if (wrongformat) {
225     fprintf(stderr, "%s: not a Radiance picture\n",
226     input[nfile].name);
227     quit(1);
228     }
229 greg 1.1 /* get picture size */
230 greg 1.4 if (fgetresolu(&input[nfile].xres, &input[nfile].yres,
231 greg 1.10 input[nfile].fp) < 0) {
232 greg 1.1 fprintf(stderr, "%s: bad picture size\n",
233     input[nfile].name);
234     quit(1);
235     }
236 greg 1.7 if (ncolumns > 0) {
237     if (curcol >= ncolumns) {
238 greg 2.11 cury = ymax + spacing;
239     curx = x0;
240 greg 1.7 curcol = 0;
241     }
242     input[nfile].xloc = curx;
243     input[nfile].yloc = cury;
244 greg 2.11 curx += input[nfile].xres + spacing;
245 greg 1.7 curcol++;
246     } else {
247     input[nfile].xloc = atoi(argv[an++]);
248 greg 2.17 if (xsgn == '+')
249     input[nfile].xloc -= input[nfile].xres;
250     else if (xsgn == '0')
251     input[nfile].xloc -= input[nfile].xres/2;
252 greg 1.7 input[nfile].yloc = atoi(argv[an++]);
253 greg 2.17 if (ysgn == '+')
254     input[nfile].yloc -= input[nfile].yres;
255     else if (ysgn == '0')
256     input[nfile].yloc -= input[nfile].yres/2;
257 greg 1.7 }
258 greg 1.1 if (input[nfile].xloc < xmin)
259     xmin = input[nfile].xloc;
260     if (input[nfile].yloc < ymin)
261     ymin = input[nfile].yloc;
262     if (input[nfile].xloc+input[nfile].xres > xmax)
263     xmax = input[nfile].xloc+input[nfile].xres;
264     if (input[nfile].yloc+input[nfile].yres > ymax)
265     ymax = input[nfile].yloc+input[nfile].yres;
266 greg 1.9 if (thislabel != NULL) {
267     if (++nfile >= MAXFILE)
268     goto toomany;
269 greg 2.14 input[nfile].name = Label;
270 greg 2.17 input[nfile].flags = 0;
271 greg 2.7 input[nfile].xres = input[nfile-1].xres;
272     input[nfile].yres = labelht;
273 greg 1.9 if ((input[nfile].fp = lblopen(thislabel,
274     &input[nfile].xres,
275     &input[nfile].yres)) == NULL)
276 greg 1.7 goto labelerr;
277     input[nfile].xloc = input[nfile-1].xloc;
278     input[nfile].yloc = input[nfile-1].yloc +
279 greg 1.9 input[nfile-1].yres-input[nfile].yres;
280 greg 1.7 }
281 greg 1.1 }
282     if (xsiz <= 0)
283     xsiz = xmax;
284 greg 2.11 else if (xsiz > xmax)
285     xmax = xsiz;
286 greg 1.1 if (ysiz <= 0)
287     ysiz = ymax;
288 greg 2.11 else if (ysiz > ymax)
289     ymax = ysiz;
290 greg 1.1 /* add new header info. */
291     printargs(argc, argv, stdout);
292 greg 1.6 fputformat(COLRFMT, stdout);
293 greg 1.10 putchar('\n');
294     fprtresolu(xsiz, ysiz, stdout);
295 greg 1.1
296     compos();
297    
298     quit(0);
299     userr:
300 greg 2.11 fprintf(stderr,
301     "Usage: %s [-x xr][-y yr][-b r g b][-a n][-s p][-o x0 y0][-la][-lh h] ",
302 greg 1.8 progname);
303 greg 2.17 fprintf(stderr, "[-t min1][+t max1][-l lab][=SS] pic1 x1 y1 ..\n");
304 greg 1.7 quit(1);
305 greg 1.9 toomany:
306     fprintf(stderr, "%s: only %d files and labels allowed\n",
307     progname, MAXFILE);
308     quit(1);
309 greg 1.7 labelerr:
310     fprintf(stderr, "%s: error opening label\n", progname);
311 greg 1.1 quit(1);
312     }
313    
314    
315     compos() /* composite pictures */
316     {
317     COLR *scanin, *scanout;
318     int y;
319     register int x, i;
320    
321 greg 1.3 scanin = (COLR *)malloc((xmax-xmin)*sizeof(COLR));
322 greg 1.1 if (scanin == NULL)
323     goto memerr;
324     scanin -= xmin;
325     if (checkthresh) {
326     scanout = (COLR *)malloc(xsiz*sizeof(COLR));
327     if (scanout == NULL)
328     goto memerr;
329     } else
330     scanout = scanin;
331     for (y = ymax-1; y >= 0; y--) {
332     for (x = 0; x < xsiz; x++)
333     copycolr(scanout[x], bgcolr);
334     for (i = 0; i < nfile; i++) {
335     if (input[i].yloc > y ||
336     input[i].yloc+input[i].yres <= y)
337     continue;
338     if (freadcolrs(scanin+input[i].xloc,
339     input[i].xres, input[i].fp) < 0) {
340 greg 1.2 fprintf(stderr, "%s: read error (y==%d)\n",
341     input[i].name,
342     y-input[i].yloc);
343 greg 1.1 quit(1);
344     }
345     if (y >= ysiz)
346     continue;
347     if (checkthresh) {
348     x = input[i].xloc+input[i].xres;
349     if (x > xsiz)
350     x = xsiz;
351     for (x--; x >= 0 && x >= input[i].xloc; x--) {
352 greg 2.17 if (input[i].flags & HASMIN &&
353 greg 1.1 cmpcolr(scanin[x], input[i].thmin) <= 0)
354     continue;
355 greg 2.17 if (input[i].flags & HASMAX &&
356 greg 1.1 cmpcolr(scanin[x], input[i].thmax) >= 0)
357     continue;
358     copycolr(scanout[x], scanin[x]);
359     }
360     }
361     }
362     if (y >= ysiz)
363     continue;
364     if (fwritecolrs(scanout, xsiz, stdout) < 0) {
365 greg 1.5 perror(progname);
366 greg 1.1 quit(1);
367     }
368     }
369     return;
370     memerr:
371 greg 1.5 perror(progname);
372 greg 1.1 quit(1);
373     }
374    
375    
376     int
377     cmpcolr(c1, c2) /* compare two colr's (improvisation) */
378     register COLR c1, c2;
379     {
380     register int i, j;
381    
382     j = 4; /* check exponents first! */
383     while (j--)
384     if (i = c1[j] - c2[j])
385     return(i);
386     return(0);
387 greg 1.9 }
388    
389    
390     FILE *
391     lblopen(s, xp, yp) /* open pipe to label generator */
392     char *s;
393     int *xp, *yp;
394     {
395     char com[128];
396     FILE *fp;
397    
398 greg 2.7 sprintf(com, "psign -s -.15 -a 2 -x %d -y %d '%.90s'", *xp, *yp, s);
399 greg 1.9 if ((fp = popen(com, "r")) == NULL)
400     return(NULL);
401     if (checkheader(fp, COLRFMT, NULL) < 0)
402     goto err;
403 greg 1.10 if (fgetresolu(xp, yp, fp) < 0)
404 greg 1.9 goto err;
405     return(fp);
406     err:
407     pclose(fp);
408     return(NULL);
409 greg 1.1 }
410    
411    
412     quit(code) /* exit gracefully */
413     int code;
414     {
415 greg 2.13 register int i;
416     /* close input files */
417     for (i = 0; i < nfile; i++)
418 greg 2.14 if (input[i].name == Command || input[i].name == Label)
419     pclose(input[i].fp);
420     else
421     fclose(input[i].fp);
422 greg 1.1 exit(code);
423     }