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.10 by greg, Wed Feb 5 16:08:14 1997 UTC vs.
Revision 3.18 by schorsch, Mon Oct 27 10:24:51 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1997 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   * Condition Radiance picture for display/output
6 + *  Added white-balance adjustment 10/01 (GW).
7   */
8  
9 + #include "platform.h"
10 + #include "rtprocess.h"
11   #include "pcond.h"
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 32 | Line 32 | double fixfrac = 0.;                   /* histogram share due to fixat
32   RESOLU  inpres;                         /* input picture resolution */
33  
34   COLOR   *fovimg;                        /* foveal (1 degree) averaged image */
35 < short   fvxr, fvyr;                     /* foveal image resolution */
35 > int     fvxr, fvyr;                     /* foveal image resolution */
36 > float   *crfimg;                        /* contrast reduction factors */
37   short   (*fixlst)[2];                   /* fixation history list */
38   int     nfixations;                     /* number of fixation points */
39 < float   bwhist[HISTRES];                /* luminance histogram */
39 > double  bwhist[HISTRES];                /* luminance histogram */
40   double  histot;                         /* total count of histogram */
41   double  bwmin, bwmax;                   /* histogram limits */
42   double  bwavg;                          /* mean brightness */
# Line 85 | Line 86 | char   *argv[];
86                          if (fixfrac > FTINY) what2do |= DO_FIXHIST;
87                          else what2do &= ~DO_FIXHIST;
88                          break;
89 +                case 'I':
90 +                        bool(DO_PREHIST);
91 +                        break;
92                  case 'l':
93                          bool(DO_LINEAR);
94                          break;
# Line 103 | Line 107 | char   *argv[];
107                  case 'e':
108                          if (i+1 >= argc) goto userr;
109                          scalef = atof(argv[++i]);
110 <                        if (argv[i][0] == '+' | argv[i][0] == '-')
110 >                        if ((argv[i][0] == '+') | (argv[i][0] == '-'))
111                                  scalef = pow(2.0, scalef);
112                          what2do |= DO_LINEAR;
113                          break;
# Line 115 | Line 119 | char   *argv[];
119                          if (i+1 >= argc) goto userr;
120                          cwarpfile = argv[++i];
121                          break;
122 <                case 't':
122 >                case 'u':
123                          if (i+1 >= argc) goto userr;
124                          ldmax = atof(argv[++i]);
125                          if (ldmax <= FTINY)
126                                  goto userr;
127                          break;
128 <                case 'b':
128 >                case 'd':
129                          if (i+1 >= argc) goto userr;
130 <                        ldmin = atof(argv[++i]);
130 >                        lddyn = atof(argv[++i]);
131                          break;
132 <                case 'd':
132 >                case 'x':
133                          if (i+1 >= argc) goto userr;
134                          if ((mapfp = fopen(argv[++i], "w")) == NULL) {
135                                  fprintf(stderr,
# Line 137 | Line 141 | char   *argv[];
141                  default:
142                          goto userr;
143                  }
144 +        if ((what2do & (DO_FIXHIST|DO_PREHIST)) == (DO_FIXHIST|DO_PREHIST)) {
145 +                fprintf(stderr, "%s: only one of -i or -I option\n", progname);
146 +                exit(1);
147 +        }
148          if ((mbcalfile != NULL) + (cwarpfile != NULL) +
149                          (outprims != stdprims) > 1) {
150                  fprintf(stderr,
# Line 144 | Line 152 | char   *argv[];
152                                  progname);
153                  exit(1);
154          }
155 <        if (outprims == stdprims & inprims != stdprims)
155 >        if ((outprims == stdprims) & (inprims != stdprims))
156                  outprims = inprims;
157 <        if (ldmin <= FTINY)
150 <                ldmin = ldmax*LDMINF;
151 <        else if (ldmin >= ldmax) {
152 <                fprintf(stderr, "%s: Ldmin (%f) >= Ldmax (%f)!\n", progname,
153 <                                ldmin, ldmax);
154 <                exit(1);
155 <        }
156 <        Bldmin = Bl(ldmin);
157 >        Bldmin = Bl(ldmax/lddyn);
158          Bldmax = Bl(ldmax);
159          if (i >= argc || i+2 < argc)
160                  goto userr;
# Line 163 | Line 164 | char   *argv[];
164                                          /* open output file */
165          if (i+2 == argc && freopen(argv[i+1], "w", stdout) == NULL)
166                  syserror(argv[i+1]);
167 < #ifdef MSDOS
168 <        setmode(fileno(infp), O_BINARY);
168 <        setmode(fileno(stdout), O_BINARY);
169 < #endif
167 >        SET_FILE_BINARY(infp);
168 >        SET_FILE_BINARY(stdout);
169          getahead();                     /* load input header */
170          printargs(argc, argv, stdout);  /* add to output header */
171 <        if (mbcalfile == NULL & outprims != stdprims)
171 >        if ((mbcalfile == NULL) & (outprims != stdprims))
172                  fputprims(outprims, stdout);
173 <        getfovimg();                    /* get foveal sample image */
174 <        if (what2do&DO_FIXHIST)         /* get fixation history? */
173 >        if ((what2do & (DO_PREHIST|DO_VEIL|DO_ACUITY)) != DO_PREHIST)
174 >                getfovimg();            /* get foveal sample image? */
175 >        if (what2do&DO_PREHIST)         /* get histogram? */
176 >                gethisto(stdin);
177 >        else if (what2do&DO_FIXHIST)    /* get fixation history? */
178                  getfixations(stdin);
179          mapimage();                     /* map the picture */
180          if (mapfp != NULL)              /* write out basic mapping */
181                  putmapping(mapfp);
182          exit(0);
183   userr:
184 <        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][-t Ldmax][-b Ldmin][-d mapfile] inpic [outpic]\n",
184 >        fprintf(stderr, "Usage: %s [-{h|a|v|s|c|l|w}[+-]][-I|-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",
185                          progname);
186          exit(1);
187   #undef bool
# Line 205 | Line 207 | char   *s;
207                  if (!strcmp(fmt,COLRFMT)) lumf = rgblum;
208                  else if (!strcmp(fmt,CIEFMT)) lumf = cielum;
209                  else lumf = NULL;
210 <                return;                 /* don't echo */
210 >                return(0);              /* don't echo */
211          }
212          if (isprims(s)) {               /* get input primaries */
213                  primsval(inprimS, s);
214                  inprims= inprimS;
215 <                return;                 /* don't echo */
215 >                return(0);              /* don't echo */
216          }
217          if (isexpos(s)) {               /* picture exposure */
218                  inpexp *= exposval(s);
219 <                return;                 /* don't echo */
219 >                return(0);              /* don't echo */
220          }
221          if (isaspect(s))                /* pixel aspect ratio */
222                  pixaspect *= aspectval(s);
223          if (isview(s))                  /* image view */
224                  gotview += sscanview(&ourview, s);
225 <        fputs(s, stdout);
225 >        return(fputs(s, stdout));
226   }
227  
228  
# Line 235 | Line 237 | getahead()                     /* load picture header */
237                  exit(1);
238          }
239          if (lumf == rgblum)
240 <                comprgb2xyzmat(inrgb2xyz, inprims);
240 >                comprgb2xyzWBmat(inrgb2xyz, inprims);
241          else if (mbcalfile != NULL) {
242                  fprintf(stderr, "%s: macbethcal only works with RGB pictures\n",
243                                  progname);
244                  exit(1);
245          }
246          if (!gotview || ourview.type == VT_PAR) {
247 <                copystruct(&ourview, &stdview);
247 >                ourview = stdview;
248                  ourview.type = VT_PER;
249                  if (pixaspect*inpres.yr < inpres.xr) {
250                          ourview.horiz = 40.0;
# Line 269 | Line 271 | mapimage()                             /* map picture and send to stdout */
271          comphist();                     /* generate adaptation histogram */
272          check2do();                     /* modify what2do flags */
273          if (what2do&DO_VEIL)
274 <                compveil();
275 <        if (!(what2do&DO_LINEAR) && mkbrmap() < 0)      /* make tone map */
276 <                what2do |= DO_LINEAR;   /* failed! -- use linear scaling */
274 >                compveil();             /* compute veil image */
275 >        if (!(what2do&DO_LINEAR))
276 >                if (mkbrmap() < 0)      /* make tone map */
277 >                        what2do |= DO_LINEAR;   /* failed! -- use linear */
278 > #if ADJ_VEIL
279 >                else if (what2do&DO_VEIL)
280 >                        adjveil();      /* else adjust veil image */
281 > #endif
282          if (what2do&DO_LINEAR) {
283                  if (scalef <= FTINY) {
284                          if (what2do&DO_HSENS)
# Line 299 | Line 306 | mapimage()                             /* map picture and send to stdout */
306  
307   getfovimg()                     /* load foveal sampled image */
308   {
302 –        extern FILE     *popen();
309          char    combuf[128];
310          FILE    *fp;
311          int     x, y;
# Line 308 | Line 314 | getfovimg()                    /* load foveal sampled image */
314          if (fvxr < 2) fvxr = 2;
315          fvyr = sqrt(ourview.vn2)/FOVDIA + 0.5;
316          if (fvyr < 2) fvyr = 2;
317 <        if (!(inpres.or & YMAJOR)) {            /* picture is rotated? */
317 >        if (!(inpres.rt & YMAJOR)) {            /* picture is rotated? */
318                  y = fvyr;
319                  fvyr = fvxr;
320                  fvxr = y;
# Line 319 | Line 325 | getfovimg()                    /* load foveal sampled image */
325          if ((fp = popen(combuf, "r")) == NULL)
326                  syserror("popen");
327          getheader(fp, NULL, NULL);      /* skip header */
328 <        if (fgetresolu(&x, &y, fp) < 0 || x != fvxr | y != fvyr)
328 >        if (fgetresolu(&x, &y, fp) < 0 || (x != fvxr) | (y != fvyr))
329                  goto readerr;
330          for (y = 0; y < fvyr; y++)
331                  if (freadscan(fovscan(y), fvxr, fp) < 0)

Diff Legend

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