ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pcompos.c
Revision: 2.35
Committed: Tue Feb 2 18:02:32 2016 UTC (8 years, 3 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad5R1
Changes since 2.34: +2 -2 lines
Log Message:
Moved declaration of popen to paths.h and put convert_command() into module

File Contents

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