ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pcompos.c
(Generate patch)

Comparing ray/src/px/pcompos.c (file contents):
Revision 1.5 by greg, Thu Sep 14 09:20:36 1989 UTC vs.
Revision 2.38 by greg, Fri Jul 19 17:37:56 2019 UTC

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

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines