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

Comparing ray/src/px/pfilt.c (file contents):
Revision 1.10 by greg, Thu Sep 13 09:46:00 1990 UTC vs.
Revision 2.8 by greg, Fri Jun 18 16:20:16 1993 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1986 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 10 | Line 10 | static char SCCSid[] = "$SunId$ LBL";
10   *     9/26/85
11   */
12  
13 < #include  <stdio.h>
13 > #include  "standard.h"
14  
15   #include  <signal.h>
16  
17   #include  "color.h"
18  
19 < extern char  *malloc();
19 > #include  "resolu.h"
20  
21 < #define  CHECKRAD       1.5     /* radius to check for filtering */
21 > #include  "paths.h"
22  
23 + extern float  *matchlamp();
24 +
25 + #define  FEQ(a,b)       ((a) >= .98*(b) && (a) <= 1.02*(b))
26 +
27 + #define  CHECKRAD       1.5     /* radius to check for filtering */
28 +
29   COLOR  exposure = WHTCOLOR;     /* exposure for the frame */
30  
31 < double  rad = 0.0;              /* output pixel radius for filtering */
31 > double  rad = 0.0;              /* output pixel radius for filtering */
32  
33   int  nrows = 0;                 /* number of rows for output */
34   int  ncols = 0;                 /* number of columns for output */
35  
36 < double  x_c = 1.0;              /* ratio of output x size to input */
37 < double  y_r = 1.0;              /* ratio of output y size to input */
36 > double  x_c = 1.0;              /* ratio of output x size to input */
37 > double  y_r = 1.0;              /* ratio of output y size to input */
38  
39   int  singlepass = 0;            /* true means skip first pass */
40  
41   int  avghot = 0;                /* true means average in bright spots */
42  
43 < double  hotlvl = 1000.0;        /* level considered "hot" */
43 > double  hotlvl = 100.0;         /* level considered "hot" */
44  
45   int  npts = 0;                  /* (half) number of points for stars */
46  
47 < double  spread = 1e-4;          /* spread for star points */
47 > double  spread = 1e-4;          /* spread for star points */
48  
43 – #define  TEMPLATE       "/usr/tmp/pfXXXXXX"
44 –
49   char  *tfname = NULL;
50  
51 + char  *lampdat = "lamp.tab";    /* lamp data file */
52 +
53 + int  order;                     /* scanline ordering of input */
54   int  xres, yres;                /* resolution of input */
55 < double  inpaspect = 1.0;        /* pixel aspect ratio of input */
55 > double  inpaspect = 1.0;        /* pixel aspect ratio of input */
56 > int  correctaspect = 0;         /* aspect ratio correction? */
57  
58 + int  wrongformat = 0;
59 +
60   int  xrad;                      /* x window size */
61   int  yrad;                      /* y window size */
62  
# Line 61 | Line 71 | main(argc, argv)
71   int  argc;
72   char  **argv;
73   {
64 –        extern char  *mktemp();
65 –        extern double  atof(), pow();
74          extern long  ftell();
75          extern int  quit(), headline();
76          FILE  *fin;
77 +        float  *lampcolor;
78 +        char  *lamptype = NULL;
79          long  fpos;
80 <        double  outaspect = 0.0;
81 <        double  d;
82 <        int  i;
83 <
80 >        double  outaspect = 0.0;
81 >        double  d;
82 >        int  i, j;
83 > #ifdef MSDOS
84 >        extern int  _fmode;
85 >        _fmode = O_BINARY;
86 >        setmode(fileno(stdin), O_BINARY);
87 >        setmode(fileno(stdout), O_BINARY);
88 > #endif
89          if (signal(SIGINT, quit) == SIG_IGN)
90                  signal(SIGINT, SIG_IGN);
91          if (signal(SIGHUP, quit) == SIG_IGN)
92                  signal(SIGINT, SIG_IGN);
93          signal(SIGTERM, quit);
94          signal(SIGPIPE, quit);
95 < #ifdef  SIGXCPU
95 > #ifdef  SIGXCPU
96          signal(SIGXCPU, quit);
97          signal(SIGXFSZ, quit);
98   #endif
99  
100 <        progname = argv[0];
100 >        progname = argv[0] = fixargv0(argv[0]);
101  
102          for (i = 1; i < argc; i++)
103                  if (argv[i][0] == '-')
# Line 103 | Line 118 | char  **argv;
118                                  } else
119                                          nrows = atoi(argv[i]);
120                                  break;
121 +                        case 'c':
122 +                                correctaspect = !correctaspect;
123 +                                break;
124                          case 'p':
125                                  i++;
126                                  outaspect = atof(argv[i]);
# Line 112 | Line 130 | char  **argv;
130                                          d = pow(2.0, atof(argv[i+1]));
131                                  else
132                                          d = atof(argv[i+1]);
133 +                                if (d < 1e-20 || d > 1e20) {
134 +                                        fprintf(stderr,
135 +                                                "%s: exposure out of range\n",
136 +                                                        argv[0]);
137 +                                        exit(1);
138 +                                }
139                                  switch (argv[i][2]) {
140                                  case '\0':
141                                          scalecolor(exposure, d);
# Line 130 | Line 154 | char  **argv;
154                                  }
155                                  i++;
156                                  break;
157 +                        case 'f':
158 +                                lampdat = argv[++i];
159 +                                break;
160 +                        case 't':
161 +                                lamptype = argv[++i];
162 +                                break;
163                          case '1':
164                                  singlepass = 1;
165                                  break;
# Line 163 | Line 193 | char  **argv;
193                          }
194                  else
195                          break;
196 <                        
196 >                                        /* get lamp data (if necessary) */
197 >        if (lamptype != NULL) {
198 >                if (loadlamps(lampdat) < 0)
199 >                        quit(1);
200 >                if ((lampcolor = matchlamp(lamptype)) == NULL) {
201 >                        fprintf(stderr, "%s: unknown lamp type\n", lamptype);
202 >                        quit(1);
203 >                }
204 >                for (i = 0; i < 3; i++)
205 >                        if (lampcolor[i] > 1e-4)
206 >                                colval(exposure,i) /= lampcolor[i];
207 >                freelamps();
208 >        }
209 >                                        /* open input file */
210          if (i == argc) {
211                  if (singlepass)
212                          fin = stdin;
# Line 191 | Line 234 | char  **argv;
234                  quit(1);
235          }
236                                          /* get header */
237 <        getheader(fin, headline);
237 >        getheader(fin, headline, NULL);
238 >        if (wrongformat) {
239 >                fprintf(stderr, "%s: input must be a Radiance picture\n",
240 >                                progname);
241 >                quit(1);
242 >        }
243                                          /* add new header info. */
244          printargs(i, argv, stdout);
245                                          /* get picture size */
246 <        if (fgetresolu(&xres, &yres, fin) != (YMAJOR|YDECR)) {
246 >        if ((order = fgetresolu(&xres, &yres, fin)) < 0) {
247                  fprintf(stderr, "%s: bad picture size\n", progname);
248                  quit(1);
249          }
250 +        if (!(order & YMAJOR))
251 +                inpaspect = 1.0/inpaspect;
252                                          /* compute output resolution */
253          if (ncols <= 0)
254                  ncols = x_c*xres + .5;
# Line 237 | Line 287 | char  **argv;
287   headline(s)                             /* process line from header */
288   char  *s;
289   {
290 +        char  fmt[32];
291 +
292          fputs(s, stdout);               /* copy to output */
293          if (isaspect(s))                /* get aspect ratio */
294                  inpaspect *= aspectval(s);
295 +        else if (isexpos(s))
296 +                hotlvl *= exposval(s);
297 +        else if (isformat(s)) {
298 +                formatval(fmt, s);
299 +                wrongformat = strcmp(fmt, COLRFMT);
300 +        }
301   }
302  
303  
# Line 273 | Line 331 | FILE  *in;
331          }
332          for (i = 0; i < yres; i++) {
333                  if (freadscan(scan, xres, in) < 0) {
334 <                        nrows = nrows * i / yres;       /* adjust frame */
334 >                        nrows = (long)nrows * i / yres; /* adjust frame */
335                          if (nrows <= 0) {
336                                  fprintf(stderr, "%s: empty frame\n", progname);
337                                  quit(1);
338                          }
339                          fprintf(stderr, "%s: warning - partial frame (%d%%)\n",
340 <                                        progname, 100*i/yres);
340 >                                        progname, (int)(100L*i/yres));
341                          yres = i;
342 +                        y_r = (double)nrows/yres;
343                          break;
344                  }
345                  pass1scan(scan, i);
# Line 295 | Line 354 | FILE  *in;
354          int  yread;
355          int  ycent, xcent;
356          int  r, c;
357 <        
357 >        
358          pass2init();
359          scan2init();
360          yread = 0;
# Line 337 | Line 396 | FILE  *in;
396  
397   scan2init()                     /* prepare scanline arrays */
398   {
399 <        double  d;
399 >        COLOR   ctmp;
400 >        double  d;
401          register int  i;
402  
403          if (rad <= 0.0) {
# Line 352 | Line 412 | scan2init()                    /* prepare scanline arrays */
412  
413                  initmask();             /* initialize filter table */
414          }
415 <        barsize = 2 * yrad;
415 >        barsize = 2*yrad + 1;
416          scanin = (COLOR **)malloc(barsize*sizeof(COLOR *));
417          for (i = 0; i < barsize; i++) {
418                  scanin[i] = (COLOR *)malloc(xres*sizeof(COLOR));
# Line 366 | Line 426 | scan2init()                    /* prepare scanline arrays */
426                  fprintf(stderr, "%s: out of memory\n", progname);
427                  quit(1);
428          }
429 <                                        /* record pixel aspect and exposure */
430 <        d = x_c / y_r;
431 <        if (d < .99 || d > 1.01)
432 <                fputaspect(d, stdout);
429 >                                        /* record pixel aspect ratio */
430 >        if (!correctaspect) {
431 >                d = order & YMAJOR ? x_c/y_r : y_r/x_c ;
432 >                if (!FEQ(d,1.0))
433 >                        fputaspect(d, stdout);
434 >        }
435 >                                        /* record exposure */
436          d = bright(exposure);
437 <        if (d < .995 || d > 1.005)
437 >        if (!FEQ(d,1.0))
438                  fputexpos(d, stdout);
439 +                                        /* record color correction */
440 +        copycolor(ctmp, exposure);
441 +        scalecolor(ctmp, 1.0/d);
442 +        if (!FEQ(colval(ctmp,RED),colval(ctmp,GRN)) ||
443 +                        !FEQ(colval(ctmp,GRN),colval(ctmp,BLU)))
444 +                fputcolcor(ctmp, stdout);
445          printf("\n");
446 <        fputresolu(YMAJOR|YDECR, ncols, nrows, stdout); /* resolution */
446 >                                        /* write out resolution */
447 >        fputresolu(order, ncols, nrows, stdout);
448   }
449  
450  

Diff Legend

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