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

Comparing ray/src/px/pcomb.c (file contents):
Revision 1.1 by greg, Thu Feb 2 10:49:22 1989 UTC vs.
Revision 1.10 by greg, Thu Apr 18 14:35:17 1991 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1989 Regents of the University of California */
1 > /* Copyright (c) 1991 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 18 | Line 18 | static char SCCSid[] = "$SunId$ LBL";
18  
19   #include "calcomp.h"
20  
21 < #define MAXINP          16              /* maximum number of input files */
21 > #define MAXINP          32              /* maximum number of input files */
22  
23   struct {
24          char    *name;          /* file name */
25          FILE    *fp;            /* stream pointer */
26          COLOR   *scan;          /* input scanline */
27 +        COLOR   coef;           /* coefficient */
28   }       input[MAXINP];                  /* input pictures */
29  
30   int     nfiles;                         /* number of input files */
31  
32   char    *vcolin[3] = {"ri", "gi", "bi"};
33   char    *vcolout[3] = {"ro", "go", "bo"};
34 + #define vbrtin          "li"
35 + #define vbrtout         "lo"
36  
37 + #define vnfiles         "nfiles"
38 + #define vxres           "xres"
39 + #define vyres           "yres"
40   #define vxpos           "x"
41   #define vypos           "y"
42  
43   int     nowarn = 0;                     /* no warning messages? */
44  
45 + int     original = 0;                   /* origninal values? */
46 +
47   int     xres=0, yres=0;                 /* picture resolution */
48  
49   int     xpos, ypos;                     /* picture position */
50  
51 + int     wrongformat = 0;
52  
53 +
54   tputs(s)                        /* put out string preceded by a tab */
55   char    *s;
56   {
57 <        putchar('\t');
58 <        fputs(s, stdout);
57 >        char    fmt[32];
58 >        double  d;
59 >        COLOR   ctmp;
60 >
61 >        if (isformat(s)) {                      /* check format */
62 >                formatval(fmt, s);
63 >                wrongformat = strcmp(fmt, COLRFMT);
64 >        } else if (original && isexpos(s)) {    /* exposure */
65 >                d = 1.0/exposval(s);
66 >                scalecolor(input[nfiles].coef, d);
67 >        } else if (original && iscolcor(s)) {   /* color correction */
68 >                colcorval(ctmp, s);
69 >                colval(input[nfiles].coef,RED) /= colval(ctmp,RED);
70 >                colval(input[nfiles].coef,GRN) /= colval(ctmp,GRN);
71 >                colval(input[nfiles].coef,BLU) /= colval(ctmp,BLU);
72 >        } else {                                /* echo unaffected line */
73 >                putchar('\t');
74 >                fputs(s, stdout);
75 >        }
76 >
77   }
78  
79  
# Line 53 | Line 81 | main(argc, argv)
81   int     argc;
82   char    *argv[];
83   {
84 <        extern double   l_redin(), l_grnin(), l_bluin();
84 >        extern double   l_redin(), l_grnin(), l_bluin(), l_brtin(), atof();
85 >        double  f;
86          int     a;
87          
88          funset(vcolin[RED], 1, l_redin);
89          funset(vcolin[GRN], 1, l_grnin);
90          funset(vcolin[BLU], 1, l_bluin);
91 +        funset(vbrtin, 1, l_brtin);
92          
93          for (a = 1; a < argc; a++)
94                  if (argv[a][0] == '-')
95                          switch (argv[a][1]) {
96                          case '\0':
97 +                        case 's':
98 +                        case 'c':
99                                  goto getfiles;
100                          case 'x':
101                                  xres = atoi(argv[++a]);
# Line 74 | Line 106 | char   *argv[];
106                          case 'w':
107                                  nowarn = !nowarn;
108                                  break;
109 +                        case 'o':
110 +                                original = !original;
111 +                                break;
112                          case 'f':
113                                  fcompile(argv[++a]);
114                                  break;
115                          case 'e':
116 <                                scompile(NULL, argv[++a]);
116 >                                scompile(argv[++a], NULL, 0);
117                                  break;
118                          default:
119 <                                eputs("Usage: ");
85 <                                eputs(argv[0]);
86 <        eputs(" [-w][-x xres][-y yres][-e expr][-f file] [picture ..]\n");
87 <                                quit(1);
119 >                                goto usage;
120                          }
121                  else
122                          break;
123   getfiles:
124 +        for (nfiles = 0; nfiles < MAXINP; nfiles++)
125 +                setcolor(input[nfiles].coef, 1.0, 1.0, 1.0);
126          nfiles = 0;
127          for ( ; a < argc; a++) {
128                  if (nfiles >= MAXINP) {
# Line 96 | Line 130 | getfiles:
130                          eputs(": too many picture files\n");
131                          quit(1);
132                  }
133 <                if (argv[a][0] == '-') {
134 <                        input[nfiles].name = "<stdin>";
135 <                        input[nfiles].fp = stdin;
136 <                } else {
133 >                if (argv[a][0] == '-')
134 >                        switch (argv[a][1]) {
135 >                        case '\0':
136 >                                input[nfiles].name = "<stdin>";
137 >                                input[nfiles].fp = stdin;
138 >                                break;
139 >                        case 's':
140 >                                f = atof(argv[++a]);
141 >                                scalecolor(input[nfiles].coef, f);
142 >                                continue;
143 >                        case 'c':
144 >                                colval(input[nfiles].coef,RED)*=atof(argv[++a]);
145 >                                colval(input[nfiles].coef,GRN)*=atof(argv[++a]);
146 >                                colval(input[nfiles].coef,BLU)*=atof(argv[++a]);
147 >                                continue;
148 >                        default:
149 >                                goto usage;
150 >                        }
151 >                else {
152                          input[nfiles].name = argv[a];
153                          input[nfiles].fp = fopen(argv[a], "r");
154                          if (input[nfiles].fp == NULL) {
155 <                                eputs(argv[a]);
107 <                                eputs(": cannot open\n");
155 >                                perror(argv[a]);
156                                  quit(1);
157                          }
158                  }
159                  fputs(input[nfiles].name, stdout);
160                  fputs(":\n", stdout);
161 <                getheader(input[nfiles].fp, tputs);
162 <                if (fscanf(input[nfiles].fp, "-Y %d +X %d\n", &ypos, &xpos) != 2) {
161 >                getheader(input[nfiles].fp, tputs, NULL);
162 >                if (wrongformat) {
163                          eputs(input[nfiles].name);
164 +                        eputs(": not in Radiance picture format\n");
165 +                        quit(1);
166 +                }
167 +                if (fgetresolu(&xpos, &ypos, input[nfiles].fp) !=
168 +                                (YMAJOR|YDECR)) {
169 +                        eputs(input[nfiles].name);
170                          eputs(": bad picture size\n");
171                          quit(1);
172                  }
# Line 128 | Line 182 | getfiles:
182                  nfiles++;
183          }
184          printargs(argc, argv, stdout);
185 +        fputformat(COLRFMT, stdout);
186          putchar('\n');
187 <        printf("-Y %d +X %d\n", yres, xres);
187 >        fputresolu(YMAJOR|YDECR, xres, yres, stdout);
188          combine();
189          quit(0);
190 + usage:
191 +        eputs("Usage: ");
192 +        eputs(argv[0]);
193 +        eputs(
194 + " [-w][-h][-x xr][-y yr][-e expr][-f file] [ [-s f][-c r g b] pic ..]\n");
195 +        quit(1);
196   }
197  
198  
199   combine()                       /* combine pictures */
200   {
201 <        int     coldef[3];
201 >        EPNODE  *coldef[3], *brtdef;
202          COLOR   *scanout;
203 +        double  d;
204          register int    i, j;
205                                                  /* check defined variables */
206 <        for (j = 0; j < 3; j++)
207 <                coldef[j] = vardefined(vcolout[j]);
206 >        for (j = 0; j < 3; j++) {
207 >                if (vardefined(vcolout[j]))
208 >                        coldef[j] = eparse(vcolout[j]);
209 >                else
210 >                        coldef[j] = NULL;
211 >        }
212 >        if (vardefined(vbrtout))
213 >                brtdef = eparse(vbrtout);
214 >        else
215 >                brtdef = NULL;
216 >                                                /* predefine variables */
217 >        varset(vnfiles, (double)nfiles);
218 >        varset(vxres, (double)xres);
219 >        varset(vyres, (double)yres);
220                                                  /* allocate scanline */
221          scanout = (COLOR *)emalloc(xres*sizeof(COLOR));
222                                                  /* combine files */
223          for (ypos = yres-1; ypos >= 0; ypos--) {
224 +            for (i = 0; i < nfiles; i++)
225 +                    if (freadscan(input[i].scan, xres, input[i].fp) < 0) {
226 +                            eputs(input[i].name);
227 +                            eputs(": read error\n");
228 +                            quit(1);
229 +                    }
230 +            varset(vypos, (double)ypos);
231 +            for (xpos = 0; xpos < xres; xpos++) {
232                  for (i = 0; i < nfiles; i++)
233 <                        if (freadscan(input[i].scan, xres, input[i].fp) < 0) {
234 <                                eputs(input[i].name);
235 <                                eputs(": read error\n");
236 <                                quit(1);
233 >                        multcolor(input[i].scan[xpos],input[i].coef);
234 >                varset(vxpos, (double)xpos);
235 >                eclock++;
236 >                if (brtdef != NULL) {
237 >                    d = evalue(brtdef);
238 >                    if (d < 0.0)
239 >                        d = 0.0;
240 >                    setcolor(scanout[xpos], d, d, d);
241 >                } else {
242 >                    for (j = 0; j < 3; j++) {
243 >                        if (coldef[j] != NULL) {
244 >                                colval(scanout[xpos],j) = evalue(coldef[j]);
245 >                        } else {
246 >                            colval(scanout[xpos],j) = 0.0;
247 >                            for (i = 0; i < nfiles; i++)
248 >                                colval(scanout[xpos],j) +=
249 >                                        colval(input[i].scan[xpos],j);
250                          }
251 <                varset(vypos, (double)ypos);
252 <                for (xpos = 0; xpos < xres; xpos++) {
253 <                        varset(vxpos, (double)xpos);
159 <                        eclock++;
160 <                        for (j = 0; j < 3; j++)
161 <                                if (coldef[j]) {
162 <                                        colval(scanout[xpos],j) = varvalue(vcolout[j]);
163 <                                        if (colval(scanout[xpos],j) < 0.0)
164 <                                                colval(scanout[xpos],j) = 0.0;
165 <                                } else {
166 <                                        colval(scanout[xpos],j) = 0.0;
167 <                                        for (i = 0; i < nfiles; i++)
168 <                                                colval(scanout[xpos],j) += colval(input[i].scan[xpos],j);
169 <                                        colval(scanout[xpos],j) /= (double)nfiles;
170 <                                }
251 >                        if (colval(scanout[xpos],j) < 0.0)
252 >                            colval(scanout[xpos],j) = 0.0;
253 >                    }
254                  }
255 <                if (fwritescan(scanout, xres, stdout) < 0) {
256 <                        eputs("write error\n");
257 <                        quit(1);
258 <                }
255 >            }
256 >            if (fwritescan(scanout, xres, stdout) < 0) {
257 >                    perror("write error");
258 >                    quit(1);
259 >            }
260          }
261          efree(scanout);
262   }
# Line 210 | Line 294 | double
294   l_bluin()                       /* get blue color */
295   {
296          return(colin((int)(argument(1)+.5), BLU));
297 + }
298 +
299 +
300 + double
301 + l_brtin()                       /* get brightness value */
302 + {
303 +        register int    fn;
304 +
305 +        fn = argument(1)+.5;
306 +        if (fn == 0)
307 +                return((double)nfiles);
308 +        if (fn < 1 || fn > nfiles) {
309 +                errno = EDOM;
310 +                return(0.0);
311 +        }
312 +        return(bright(input[fn-1].scan[xpos]));
313   }
314  
315  

Diff Legend

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