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

Diff Legend

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