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

Comparing ray/src/px/pcond.c (file contents):
Revision 3.8 by greg, Thu Jan 9 13:56:23 1997 UTC vs.
Revision 3.11 by greg, Wed Mar 19 13:04:09 1997 UTC

# Line 12 | Line 12 | static char SCCSid[] = "$SunId$ LBL";
12  
13  
14   #define LDMAX           100             /* default max. display luminance */
15 < #define LDMINF          0.01            /* default min. display lum. factor */
15 > #define LDDYN           32              /* default dynamic range */
16  
17   int     what2do = 0;                    /* desired adjustments */
18  
19   double  ldmax = LDMAX;                  /* maximum output luminance */
20 < double  ldmin = 0.;                     /* minimum output luminance */
21 < double  Bldmin, Bldmax;                 /* Bl(ldmin) and Bl(ldmax) */
20 > double  lddyn = LDDYN;                  /* display dynamic range */
21 > double  Bldmin, Bldmax;                 /* Bl(ldmax/lddyn) and Bl(ldmax) */
22  
23   char    *progname;                      /* global argv[0] */
24  
# Line 28 | Line 28 | FILE   *mapfp = NULL;                  /* tone-mapping function stream
28   VIEW    ourview = STDVIEW;              /* picture view */
29   int     gotview = 0;                    /* picture has view */
30   double  pixaspect = 1.0;                /* pixel aspect ratio */
31 + double  fixfrac = 0.;                   /* histogram share due to fixations */
32   RESOLU  inpres;                         /* input picture resolution */
33  
34   COLOR   *fovimg;                        /* foveal (1 degree) averaged image */
35   short   fvxr, fvyr;                     /* foveal image resolution */
36 + short   (*fixlst)[2];                   /* fixation history list */
37 + int     nfixations;                     /* number of fixation points */
38   float   bwhist[HISTRES];                /* luminance histogram */
39   double  histot;                         /* total count of histogram */
40   double  bwmin, bwmax;                   /* histogram limits */
# Line 76 | Line 79 | char   *argv[];
79                  case 'w':
80                          bool(DO_CWEIGHT);
81                          break;
82 +                case 'i':
83 +                        if (i+1 >= argc) goto userr;
84 +                        fixfrac = atof(argv[++i]);
85 +                        if (fixfrac > FTINY) what2do |= DO_FIXHIST;
86 +                        else what2do &= ~DO_FIXHIST;
87 +                        break;
88                  case 'l':
89                          bool(DO_LINEAR);
90                          break;
# Line 102 | Line 111 | char   *argv[];
111                          if (i+1 >= argc) goto userr;
112                          mbcalfile = argv[++i];
113                          break;
114 <                case 't':
114 >                case 'm':
115                          if (i+1 >= argc) goto userr;
116 +                        cwarpfile = argv[++i];
117 +                        break;
118 +                case 'u':
119 +                        if (i+1 >= argc) goto userr;
120                          ldmax = atof(argv[++i]);
121                          if (ldmax <= FTINY)
122                                  goto userr;
123                          break;
124 <                case 'b':
124 >                case 'd':
125                          if (i+1 >= argc) goto userr;
126 <                        ldmin = atof(argv[++i]);
126 >                        lddyn = atof(argv[++i]);
127                          break;
128 <                case 'm':
128 >                case 'x':
129                          if (i+1 >= argc) goto userr;
130                          if ((mapfp = fopen(argv[++i], "w")) == NULL) {
131                                  fprintf(stderr,
# Line 124 | Line 137 | char   *argv[];
137                  default:
138                          goto userr;
139                  }
140 <        if (mbcalfile != NULL & outprims != stdprims) {
141 <                fprintf(stderr, "%s: only one of -p or -f option supported\n",
140 >        if ((mbcalfile != NULL) + (cwarpfile != NULL) +
141 >                        (outprims != stdprims) > 1) {
142 >                fprintf(stderr,
143 >                        "%s: only one of -p, -m or -f option supported\n",
144                                  progname);
145                  exit(1);
146          }
147          if (outprims == stdprims & inprims != stdprims)
148                  outprims = inprims;
149 <        if (ldmin <= FTINY)
135 <                ldmin = ldmax*LDMINF;
136 <        else if (ldmin >= ldmax) {
137 <                fprintf(stderr, "%s: Ldmin (%f) >= Ldmax (%f)!\n", progname,
138 <                                ldmin, ldmax);
139 <                exit(1);
140 <        }
141 <        Bldmin = Bl(ldmin);
149 >        Bldmin = Bl(ldmax/lddyn);
150          Bldmax = Bl(ldmax);
151          if (i >= argc || i+2 < argc)
152                  goto userr;
153 +                                        /* open input file */
154          if ((infp = fopen(infn=argv[i], "r")) == NULL)
155                  syserror(infn);
156 +                                        /* open output file */
157          if (i+2 == argc && freopen(argv[i+1], "w", stdout) == NULL)
158                  syserror(argv[i+1]);
159   #ifdef MSDOS
# Line 152 | Line 162 | char   *argv[];
162   #endif
163          getahead();                     /* load input header */
164          printargs(argc, argv, stdout);  /* add to output header */
165 <        if (outprims != inprims)
165 >        if (mbcalfile == NULL & outprims != stdprims)
166                  fputprims(outprims, stdout);
167 +        getfovimg();                    /* get foveal sample image */
168 +        if (what2do&DO_FIXHIST)         /* get fixation history? */
169 +                getfixations(stdin);
170          mapimage();                     /* map the picture */
171          if (mapfp != NULL)              /* write out basic mapping */
172                  putmapping(mapfp);
173          exit(0);
174   userr:
175 <        fprintf(stderr, "Usage: %s [-{h|a|v|s|c|l|w}[+-]][-e ev][-p xr yr xg yg xb yb xw yw|-f mbf.cal][-t Ldmax][-b Ldmin][-m mapfile] inpic [outpic]\n",
175 >        fprintf(stderr, "Usage: %s [-{h|a|v|s|c|l|w}[+-]][-i ffrac][-e ev][-p xr yr xg yg xb yb xw yw|-f mbf.cal|-m rgb.cwp][-u Ldmax][-d Lddyn][-x mapfile] inpic [outpic]\n",
176                          progname);
177          exit(1);
178   #undef bool
# Line 246 | Line 259 | mapimage()                             /* map picture and send to stdout */
259   {
260          COLOR   *scan;
261  
249 – #ifdef DEBUG
250 –        fprintf(stderr, "%s: generating histogram...", progname);
251 – #endif
252 –        getfovimg();                    /* get foveal sample image */
262          comphist();                     /* generate adaptation histogram */
254 – #ifdef DEBUG
255 –        fputs("done\n", stderr);
256 – #endif
263          check2do();                     /* modify what2do flags */
264 <        if (what2do&DO_VEIL) {
259 < #ifdef DEBUG
260 <                fprintf(stderr, "%s: computing veiling...", progname);
261 < #endif
264 >        if (what2do&DO_VEIL)
265                  compveil();
266 < #ifdef DEBUG
267 <                fputs("done\n", stderr);
265 < #endif
266 <        }
267 < #ifdef DEBUG
268 <        fprintf(stderr, "%s: computing brightness mapping...", progname);
269 < #endif
270 <        if (!(what2do&DO_LINEAR) && mkbrmap() < 0) {    /* make tone map */
271 <                what2do |= DO_LINEAR;           /* use linear scaling */
272 < #ifdef DEBUG
273 <                fputs("failed!\n", stderr);
274 <        } else
275 <                fputs("done\n", stderr);
276 < #else
277 <        }
278 < #endif
266 >        if (!(what2do&DO_LINEAR) && mkbrmap() < 0)      /* make tone map */
267 >                what2do |= DO_LINEAR;   /* failed! -- use linear scaling */
268          if (what2do&DO_LINEAR) {
269                  if (scalef <= FTINY) {
270                          if (what2do&DO_HSENS)
# Line 285 | Line 274 | mapimage()                             /* map picture and send to stdout */
274                                  scalef = Lb(0.5*(Bldmax+Bldmin)) / Lb(bwavg);
275                          scalef *= WHTEFFICACY/(inpexp*ldmax);
276                  }
288 – #ifdef DEBUG
289 –                fprintf(stderr, "%s: linear scaling factor = %f\n",
290 –                                progname, scalef);
291 – #endif
277                  fputexpos(inpexp*scalef, stdout);       /* record exposure */
278                  if (lumf == cielum) scalef /= WHTEFFICACY;
279          }
280 <        putchar('\n');                  /* complete header */
280 >        fputformat(COLRFMT, stdout);    /* complete header */
281 >        putchar('\n');
282          fputsresolu(&inpres, stdout);   /* resolution doesn't change */
283 <
283 >                                        /* condition our image */
284          for (scan = firstscan(); scan != NULL; scan = nextscan())
285                  if (fwritescan(scan, scanlen(&inpres), stdout) < 0) {
286                          fprintf(stderr, "%s: scanline write error\n",
# Line 304 | Line 290 | mapimage()                             /* map picture and send to stdout */
290   }
291  
292  
307 – double
308 – centprob(x, y)                  /* center-weighting probability function */
309 – int     x, y;
310 – {
311 –        double  xr, yr, p;
312 –                                /* paraboloid, 0 at 90 degrees from center */
313 –        xr = (x - .5*(fvxr-1))/90.;     /* 180 degree fisheye has fv?r == 90 */
314 –        yr = (y - .5*(fvyr-1))/90.;
315 –        p = 1. - xr*xr - yr*yr;
316 –        return(p < 0. ? 0. : p);
317 – }
318 –
319 –
293   getfovimg()                     /* load foveal sampled image */
294   {
295          extern FILE     *popen();
# Line 335 | Line 308 | getfovimg()                    /* load foveal sampled image */
308          }
309          if ((fovimg = (COLOR *)malloc(fvxr*fvyr*sizeof(COLOR))) == NULL)
310                  syserror("malloc");
311 <        sprintf(combuf, "pfilt -1 -b -x %d -y %d %s", fvxr, fvyr, infn);
311 >        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d %s", fvxr, fvyr, infn);
312          if ((fp = popen(combuf, "r")) == NULL)
313                  syserror("popen");
314          getheader(fp, NULL, NULL);      /* skip header */
# Line 353 | Line 326 | readerr:
326   }
327  
328  
356 – comphist()                      /* create foveal sampled image and histogram */
357 – {
358 –        double  l, b, lwmin, lwmax;
359 –        register int    x, y;
360 –
361 –        lwmin = 1e10;                   /* find extrema */
362 –        lwmax = 0.;
363 –        for (y = 0; y < fvyr; y++)
364 –                for (x = 0; x < fvxr; x++) {
365 –                        l = plum(fovscan(y)[x]);
366 –                        if (l < lwmin) lwmin = l;
367 –                        if (l > lwmax) lwmax = l;
368 –                }
369 –        lwmin -= FTINY;
370 –        lwmax += FTINY;
371 –        if (lwmin < LMIN) lwmin = LMIN;
372 –        if (lwmax > LMAX) lwmax = LMAX;
373 –                                        /* compute histogram */
374 –        bwmin = Bl(lwmin);
375 –        bwmax = Bl(lwmax);
376 –        bwavg = 0.;
377 –        for (y = 0; y < fvyr; y++)
378 –                for (x = 0; x < fvxr; x++) {
379 –                        l = plum(fovscan(y)[x]);
380 –                        if (l < lwmin) continue;
381 –                        if (l > lwmax) continue;
382 –                        b = Bl(l);
383 –                        bwavg += b;
384 –                        l = what2do&DO_CWEIGHT ? centprob(x,y) : 1.;
385 –                        bwhist[bwhi(b)] += l;
386 –                        histot += l;
387 –                }
388 –        bwavg /= histot;
389 – }
390 –
391 –
329   check2do()              /* check histogram to see what isn't worth doing */
330   {
331          double  sum;
# Line 396 | Line 333 | check2do()             /* check histogram to see what isn't worth
333          register int    i;
334  
335                                          /* check for within display range */
336 <        l = Lb(bwmax)/Lb(bwmin);
400 <        if (l <= ldmax/ldmin)
336 >        if (bwmax - bwmin <= Bldmax - Bldmin)
337                  what2do |= DO_LINEAR;
338                                          /* determine if veiling significant */
339 <        if (l < 100.)                   /* heuristic */
339 >        if (bwmax - bwmin < 4.5)                /* heuristic */
340                  what2do &= ~DO_VEIL;
341  
342          if (!(what2do & (DO_ACUITY|DO_COLOR)))

Diff Legend

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