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 2.4 by greg, Wed Feb 19 17:08:00 1992 UTC vs.
Revision 2.32 by greg, Fri Sep 8 21:38:25 2006 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1991 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   *
# Line 11 | Line 8 | static char SCCSid[] = "$SunId$ LBL";
8   */
9  
10   #include  <stdio.h>
11 + #include  <math.h>
12 + #include  <time.h>
13 + #include  <string.h>
14  
15 < #include  "color.h"
15 > #include "copyright.h"
16  
17 + #include  "platform.h"
18 + #include  "rtprocess.h"
19 + #include  "rterror.h"
20 + #include  "color.h"
21   #include  "resolu.h"
22  
23 < #define  MAXFILE        32
23 > #define  MAXFILE        512
24  
25 + #define  HASMIN         1
26 + #define  HASMAX         2
27 +
28                                          /* output picture size */
29   int  xsiz = 0;
30   int  ysiz = 0;
# Line 33 | Line 40 | 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 {
# Line 40 | Line 51 | struct {
51          FILE  *fp;                      /* stream pointer */
52          int  xres, yres;                /* picture size */
53          int  xloc, yloc;                /* anchor point */
54 <        int  hasmin, hasmax;            /* has threshold values */
55 <        COLR  thmin, thmax;             /* thresholds */
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  
51 – FILE  *popen(), *lblopen();
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 < tabputs(s)                      /* print line preceded by a tab */
71 < char  *s;
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 (isformat(s)) {
81 <                formatval(fmt, s);
82 <                wrongformat = strcmp(fmt, COLRFMT);
83 <        } else {
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 < main(argc, argv)
97 < int  argc;
98 < char  *argv[];
96 > int
97 > main(
98 >        int  argc,
99 >        char  *argv[]
100 > )
101   {
102          int  ncolumns = 0;
103          int  autolabel = 0;
104 <        int  curcol = 0, curx = 0, cury = 0;
104 >        int  curcol = 0, x0 = 0, curx = 0, cury = 0, spacing = 0;
105 >        int  xsgn, ysgn;
106          char  *thislabel;
107          int  an;
108 <
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 <                        xmax = xsiz = atoi(argv[++an]);
119 >                        xsiz = atoi(argv[++an]);
120                          break;
121                  case 'y':
122 <                        ymax = ysiz = atoi(argv[++an]);
122 >                        ysiz = atoi(argv[++an]);
123                          break;
124                  case 'b':
125                          setcolr(bgcolr, atof(argv[an+1]),
# Line 95 | Line 130 | char  *argv[];
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':
# Line 116 | Line 158 | char  *argv[];
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 <                if (autolabel)
167 <                        thislabel = argv[an];
168 <                else
169 <                        thislabel = NULL;
170 <                input[nfile].hasmin = input[nfile].hasmax = 0;
127 <                while (an < argc && (argv[an][0] == '-' || argv[an][0] == '+'))
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].hasmin = 1;
176 <                                        setcolr(input[nfile].thmin,
177 <                                                        atof(argv[an+1]),
178 <                                                        atof(argv[an+1]),
179 <                                                        atof(argv[an+1]));
180 <                                } else {
181 <                                        input[nfile].hasmax = 1;
182 <                                        setcolr(input[nfile].thmax,
140 <                                                        atof(argv[an+1]),
141 <                                                        atof(argv[an+1]),
142 <                                                        atof(argv[an+1]));
143 <                                }
144 <                                an += 2;
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+1];
150 <                                an += 2;
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;
# Line 156 | Line 203 | dofiles:
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 = "<stdin>";
214 >                        input[nfile].name = StandardInput;
215                          input[nfile].fp = stdin;
216                  } else {
217                          if (argv[an][0] == '!') {
218 <                                input[nfile].name = "<Command>";
218 >                                input[nfile].name = Command;
219                                  input[nfile].fp = popen(argv[an]+1, "r");
220                          } else {
221                                  input[nfile].name = argv[an];
# Line 177 | Line 228 | getfile:
228                  }
229                  an++;
230                                                  /* get header */
231 <                printf("%s:\n", input[nfile].name);
232 <                getheader(input[nfile].fp, tabputs, NULL);
231 >                if (echoheader)
232 >                        printf("%s:\n", input[nfile].name);
233 >                getheader(input[nfile].fp, headline, NULL);
234                  if (wrongformat) {
235 <                        fprintf(stderr, "%s: not a Radiance picture\n",
235 >                        fprintf(stderr, "%s: incompatible input format\n",
236                                          input[nfile].name);
237                          quit(1);
238                  }
# Line 193 | Line 245 | getfile:
245                  }
246                  if (ncolumns > 0) {
247                          if (curcol >= ncolumns) {
248 <                                cury = ymax;
249 <                                curx = 0;
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;
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;
# Line 216 | Line 276 | getfile:
276                  if (thislabel != NULL) {
277                          if (++nfile >= MAXFILE)
278                                  goto toomany;
279 <                        input[nfile].name = "<Label>";
280 <                        input[nfile].hasmin = input[nfile].hasmax = 0;
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)
# Line 229 | Line 291 | getfile:
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                                          /* add new header info. */
301          printargs(argc, argv, stdout);
302 <        fputformat(COLRFMT, stdout);
302 >        if (strcmp(ourfmt, PICFMT))
303 >                fputformat(ourfmt, stdout);     /* print format if known */
304          putchar('\n');
305          fprtresolu(xsiz, ysiz, stdout);
306  
# Line 241 | Line 308 | getfile:
308          
309          quit(0);
310   userr:
311 <        fprintf(stderr, "Usage: %s [-x xr][-y yr][-b r g b][-a n][-la][-lh h] ",
311 >        fprintf(stderr,
312 >        "Usage: %s [-h][-x xr][-y yr][-b r g b][-a n][-s p][-o x0 y0][-la][-lh h] ",
313                          progname);
314 <        fprintf(stderr, "[-t min1][+t max1][-l lab] pic1 x1 y1 ..\n");
314 >        fprintf(stderr, "[-t min1][+t max1][-l lab][=SS] pic1 x1 y1 ..\n");
315          quit(1);
316   toomany:
317          fprintf(stderr, "%s: only %d files and labels allowed\n",
# Line 252 | Line 320 | toomany:
320   labelerr:
321          fprintf(stderr, "%s: error opening label\n", progname);
322          quit(1);
323 +        return 1; /* pro forma return */
324   }
325  
326  
327 < compos()                                /* composite pictures */
327 > static void
328 > compos(void)                            /* composite pictures */
329   {
330          COLR  *scanin, *scanout;
331          int  y;
# Line 292 | Line 362 | compos()                               /* composite pictures */
362                                  if (x > xsiz)
363                                          x = xsiz;
364                                  for (x--; x >= 0 && x >= input[i].xloc; x--) {
365 <                                        if (input[i].hasmin &&
365 >                                        if (input[i].flags & HASMIN &&
366                                          cmpcolr(scanin[x], input[i].thmin) <= 0)
367                                                  continue;
368 <                                        if (input[i].hasmax &&
368 >                                        if (input[i].flags & HASMAX &&
369                                          cmpcolr(scanin[x], input[i].thmax) >= 0)
370                                                  continue;
371                                          copycolr(scanout[x], scanin[x]);
# Line 309 | Line 379 | compos()                               /* composite pictures */
379                          quit(1);
380                  }
381          }
382 +                                        /* read remainders from streams */
383 +        for (i = 0; i < nfile; i++)
384 +                if (input[i].name[0] == '<')
385 +                        while (getc(input[i].fp) != EOF)
386 +                                ;
387          return;
388   memerr:
389          perror(progname);
# Line 316 | Line 391 | memerr:
391   }
392  
393  
394 < int
395 < cmpcolr(c1, c2)                 /* compare two colr's (improvisation) */
396 < register COLR  c1, c2;
394 > static int
395 > cmpcolr(                        /* compare COLR to luminance */
396 >        register COLR  c1,
397 >        double lv2
398 > )
399   {
400 <        register int  i, j;
401 <
402 <        j = 4;                          /* check exponents first! */
403 <        while (j--)
404 <                if (i = c1[j] - c2[j])
405 <                        return(i);
400 >        double  lv1 = .0;
401 >        
402 >        if (c1[EXP])
403 >                lv1 = ldexp((double)normbright(c1), (int)c1[EXP]-(COLXS+8));
404 >        if (lv1 < lv2) return(-1);
405 >        if (lv1 > lv2) return(1);
406          return(0);
407   }
408  
409  
410 < FILE *
411 < lblopen(s, xp, yp)              /* open pipe to label generator */
412 < char  *s;
413 < int  *xp, *yp;
410 > static FILE *
411 > lblopen(                /* open pipe to label generator */
412 >        char  *s,
413 >        int  *xp,
414 >        int  *yp
415 > )
416   {
417 <        char  com[128];
417 >        char  com[PATH_MAX];
418          FILE  *fp;
419  
420 <        sprintf(com, "psign -h %d '%.30s' | pfilt -1 -x /2 -y /2",
342 <                        2*labelht, s);
420 >        sprintf(com, "psign -s -.15 -a 2 -x %d -y %d '%.90s'", *xp, *yp, s);
421          if ((fp = popen(com, "r")) == NULL)
422                  return(NULL);
423          if (checkheader(fp, COLRFMT, NULL) < 0)
# Line 353 | Line 431 | err:
431   }
432  
433  
434 + void
435   quit(code)              /* exit gracefully */
436   int  code;
437   {
438 +        register int  i;
439 +                                /* close input files */
440 +        for (i = 0; i < nfile; i++)
441 +                if (input[i].name == Command || input[i].name == Label)
442 +                        pclose(input[i].fp);
443 +                else
444 +                        fclose(input[i].fp);
445          exit(code);
446   }

Diff Legend

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