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.16 by greg, Fri Jun 4 14:47:01 1993 UTC vs.
Revision 2.33 by greg, Mon Aug 10 04:26:38 2009 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1992 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>
14 –
11   #include  <math.h>
12 + #include  <time.h>
13 + #include  <string.h>
14  
15 < #ifdef MSDOS
18 < #include  <fcntl.h>
19 < #endif
15 > #include "copyright.h"
16  
17 + #include  "platform.h"
18 + #include  "rtprocess.h"
19 + #include  "rterror.h"
20   #include  "color.h"
22 –
21   #include  "resolu.h"
22  
23 < #define  MAXFILE        64
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;
# Line 39 | 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  
# Line 49 | 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  
60 – FILE  *popen(), *lblopen();
64  
65 < extern char  *malloc();
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 < tabputs(s)                      /* print line preceded by a tab */
72 < char  *s;
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, x0 = 0, curx = 0, cury = 0, spacing = 0;
105 +        int  xsgn, ysgn;
106          char  *thislabel;
107          int  an;
108 < #ifdef MSDOS
109 <        extern int  _fmode;
110 <        _fmode = O_BINARY;
92 <        setmode(fileno(stdin), O_BINARY);
93 <        setmode(fileno(stdout), O_BINARY);
94 < #endif
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;
# Line 139 | 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                  thislabel = NULL;
167 <                input[nfile].hasmin = input[nfile].hasmax = 0;
168 <                while (an < argc && (argv[an][0] == '-' || argv[an][0] == '+'))
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,
160 <                                                        atof(argv[an+1]),
161 <                                                        atof(argv[an+1]),
162 <                                                        atof(argv[an+1]));
163 <                                }
164 <                                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];
170 <                                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 176 | 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] == '!') {
# Line 199 | 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 225 | Line 255 | getfile:
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 239 | Line 277 | getfile:
277                          if (++nfile >= MAXFILE)
278                                  goto toomany;
279                          input[nfile].name = Label;
280 <                        input[nfile].hasmin = input[nfile].hasmax = 0;
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,
# Line 261 | Line 299 | getfile:
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 270 | Line 309 | getfile:
309          quit(0);
310   userr:
311          fprintf(stderr,
312 <        "Usage: %s [-x xr][-y yr][-b r g b][-a n][-s p][-o x0 y0][-la][-lh h] ",
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 281 | 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 321 | 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 338 | 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 345 | 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 -s -.15 -a 2 -x %d -y %d '%.90s'", *xp, *yp, s);
# Line 381 | Line 431 | err:
431   }
432  
433  
434 + void
435   quit(code)              /* exit gracefully */
436   int  code;
437   {

Diff Legend

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