ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pcompos.c
Revision: 2.30
Committed: Wed May 25 04:44:26 2005 UTC (18 years, 11 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad3R7P2, rad3R7P1
Changes since 2.29: +2 -2 lines
Log Message:
Created rtcontrib program for computing ray contributions and coefficients

File Contents

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