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.9 by greg, Wed Jan 29 13:22:08 1997 UTC vs.
Revision 3.20 by schorsch, Fri Jan 2 12:47:01 2004 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 "paths.h"
11 + #include "rtprocess.h"
12   #include "pcond.h"
13  
14  
15   #define LDMAX           100             /* default max. display luminance */
16 < #define LDMINF          0.01            /* default min. display lum. factor */
16 > #define LDDYN           32              /* default dynamic range */
17  
18   int     what2do = 0;                    /* desired adjustments */
19  
20   double  ldmax = LDMAX;                  /* maximum output luminance */
21 < double  ldmin = 0.;                     /* minimum output luminance */
22 < double  Bldmin, Bldmax;                 /* Bl(ldmin) and Bl(ldmax) */
21 > double  lddyn = LDDYN;                  /* display dynamic range */
22 > double  Bldmin, Bldmax;                 /* Bl(ldmax/lddyn) and Bl(ldmax) */
23  
24   char    *progname;                      /* global argv[0] */
25  
# Line 32 | Line 33 | double fixfrac = 0.;                   /* histogram share due to fixat
33   RESOLU  inpres;                         /* input picture resolution */
34  
35   COLOR   *fovimg;                        /* foveal (1 degree) averaged image */
36 < short   fvxr, fvyr;                     /* foveal image resolution */
36 > int     fvxr, fvyr;                     /* foveal image resolution */
37 > float   *crfimg;                        /* contrast reduction factors */
38   short   (*fixlst)[2];                   /* fixation history list */
39   int     nfixations;                     /* number of fixation points */
40 < float   bwhist[HISTRES];                /* luminance histogram */
40 > double  bwhist[HISTRES];                /* luminance histogram */
41   double  histot;                         /* total count of histogram */
42   double  bwmin, bwmax;                   /* histogram limits */
43   double  bwavg;                          /* mean brightness */
44  
45   double  scalef = 0.;                    /* linear scaling factor */
46  
47 + static gethfunc headline;
48  
49 +
50   main(argc, argv)
51   int     argc;
52   char    *argv[];
# Line 85 | Line 89 | char   *argv[];
89                          if (fixfrac > FTINY) what2do |= DO_FIXHIST;
90                          else what2do &= ~DO_FIXHIST;
91                          break;
92 +                case 'I':
93 +                        bool(DO_PREHIST);
94 +                        break;
95                  case 'l':
96                          bool(DO_LINEAR);
97                          break;
# Line 103 | Line 110 | char   *argv[];
110                  case 'e':
111                          if (i+1 >= argc) goto userr;
112                          scalef = atof(argv[++i]);
113 <                        if (argv[i][0] == '+' | argv[i][0] == '-')
113 >                        if ((argv[i][0] == '+') | (argv[i][0] == '-'))
114                                  scalef = pow(2.0, scalef);
115                          what2do |= DO_LINEAR;
116                          break;
# Line 111 | Line 118 | char   *argv[];
118                          if (i+1 >= argc) goto userr;
119                          mbcalfile = argv[++i];
120                          break;
121 <                case 't':
121 >                case 'm':
122                          if (i+1 >= argc) goto userr;
123 +                        cwarpfile = argv[++i];
124 +                        break;
125 +                case 'u':
126 +                        if (i+1 >= argc) goto userr;
127                          ldmax = atof(argv[++i]);
128                          if (ldmax <= FTINY)
129                                  goto userr;
130                          break;
131 <                case 'b':
131 >                case 'd':
132                          if (i+1 >= argc) goto userr;
133 <                        ldmin = atof(argv[++i]);
133 >                        lddyn = atof(argv[++i]);
134                          break;
135 <                case 'm':
135 >                case 'x':
136                          if (i+1 >= argc) goto userr;
137                          if ((mapfp = fopen(argv[++i], "w")) == NULL) {
138                                  fprintf(stderr,
# Line 133 | Line 144 | char   *argv[];
144                  default:
145                          goto userr;
146                  }
147 <        if (mbcalfile != NULL & outprims != stdprims) {
148 <                fprintf(stderr, "%s: only one of -p or -f option supported\n",
138 <                                progname);
147 >        if ((what2do & (DO_FIXHIST|DO_PREHIST)) == (DO_FIXHIST|DO_PREHIST)) {
148 >                fprintf(stderr, "%s: only one of -i or -I option\n", progname);
149                  exit(1);
150          }
151 <        if (outprims == stdprims & inprims != stdprims)
152 <                outprims = inprims;
153 <        if (ldmin <= FTINY)
154 <                ldmin = ldmax*LDMINF;
155 <        else if (ldmin >= ldmax) {
146 <                fprintf(stderr, "%s: Ldmin (%f) >= Ldmax (%f)!\n", progname,
147 <                                ldmin, ldmax);
151 >        if ((mbcalfile != NULL) + (cwarpfile != NULL) +
152 >                        (outprims != stdprims) > 1) {
153 >                fprintf(stderr,
154 >                        "%s: only one of -p, -m or -f option supported\n",
155 >                                progname);
156                  exit(1);
157          }
158 <        Bldmin = Bl(ldmin);
158 >        if ((outprims == stdprims) & (inprims != stdprims))
159 >                outprims = inprims;
160 >        Bldmin = Bl(ldmax/lddyn);
161          Bldmax = Bl(ldmax);
162          if (i >= argc || i+2 < argc)
163                  goto userr;
# Line 157 | Line 167 | char   *argv[];
167                                          /* open output file */
168          if (i+2 == argc && freopen(argv[i+1], "w", stdout) == NULL)
169                  syserror(argv[i+1]);
170 < #ifdef MSDOS
171 <        setmode(fileno(infp), O_BINARY);
162 <        setmode(fileno(stdout), O_BINARY);
163 < #endif
170 >        SET_FILE_BINARY(infp);
171 >        SET_FILE_BINARY(stdout);
172          getahead();                     /* load input header */
173          printargs(argc, argv, stdout);  /* add to output header */
174 <        if (mbcalfile == NULL & outprims != stdprims)
174 >        if ((mbcalfile == NULL) & (outprims != stdprims))
175                  fputprims(outprims, stdout);
176 <        getfovimg();                    /* get foveal sample image */
177 <        if (what2do&DO_FIXHIST)         /* get fixation history? */
176 >        if ((what2do & (DO_PREHIST|DO_VEIL|DO_ACUITY)) != DO_PREHIST)
177 >                getfovimg();            /* get foveal sample image? */
178 >        if (what2do&DO_PREHIST)         /* get histogram? */
179 >                gethisto(stdin);
180 >        else if (what2do&DO_FIXHIST)    /* get fixation history? */
181                  getfixations(stdin);
182          mapimage();                     /* map the picture */
183          if (mapfp != NULL)              /* write out basic mapping */
184                  putmapping(mapfp);
185          exit(0);
186   userr:
187 <        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][-t Ldmax][-b Ldmin][-m mapfile] inpic [outpic]\n",
187 >        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",
188                          progname);
189          exit(1);
190   #undef bool
# Line 189 | Line 200 | char   *s;
200   }
201  
202  
203 < headline(s)                             /* process header line */
204 < char    *s;
203 > static int
204 > headline(                               /* process header line */
205 >        char    *s,
206 >        void    *p
207 > )
208   {
209          static RGBPRIMS inprimS;
210          char    fmt[32];
# Line 199 | Line 213 | char   *s;
213                  if (!strcmp(fmt,COLRFMT)) lumf = rgblum;
214                  else if (!strcmp(fmt,CIEFMT)) lumf = cielum;
215                  else lumf = NULL;
216 <                return;                 /* don't echo */
216 >                return(0);              /* don't echo */
217          }
218          if (isprims(s)) {               /* get input primaries */
219                  primsval(inprimS, s);
220                  inprims= inprimS;
221 <                return;                 /* don't echo */
221 >                return(0);              /* don't echo */
222          }
223          if (isexpos(s)) {               /* picture exposure */
224                  inpexp *= exposval(s);
225 <                return;                 /* don't echo */
225 >                return(0);              /* don't echo */
226          }
227          if (isaspect(s))                /* pixel aspect ratio */
228                  pixaspect *= aspectval(s);
229          if (isview(s))                  /* image view */
230                  gotview += sscanview(&ourview, s);
231 <        fputs(s, stdout);
231 >        return(fputs(s, stdout));
232   }
233  
234  
# Line 229 | Line 243 | getahead()                     /* load picture header */
243                  exit(1);
244          }
245          if (lumf == rgblum)
246 <                comprgb2xyzmat(inrgb2xyz, inprims);
246 >                comprgb2xyzWBmat(inrgb2xyz, inprims);
247          else if (mbcalfile != NULL) {
248                  fprintf(stderr, "%s: macbethcal only works with RGB pictures\n",
249                                  progname);
250                  exit(1);
251          }
252          if (!gotview || ourview.type == VT_PAR) {
253 <                copystruct(&ourview, &stdview);
253 >                ourview = stdview;
254                  ourview.type = VT_PER;
255                  if (pixaspect*inpres.yr < inpres.xr) {
256                          ourview.horiz = 40.0;
# Line 263 | Line 277 | mapimage()                             /* map picture and send to stdout */
277          comphist();                     /* generate adaptation histogram */
278          check2do();                     /* modify what2do flags */
279          if (what2do&DO_VEIL)
280 <                compveil();
281 <        if (!(what2do&DO_LINEAR) && mkbrmap() < 0)      /* make tone map */
282 <                what2do |= DO_LINEAR;   /* failed! -- use linear scaling */
280 >                compveil();             /* compute veil image */
281 >        if (!(what2do&DO_LINEAR))
282 >                if (mkbrmap() < 0)      /* make tone map */
283 >                        what2do |= DO_LINEAR;   /* failed! -- use linear */
284 > #if ADJ_VEIL
285 >                else if (what2do&DO_VEIL)
286 >                        adjveil();      /* else adjust veil image */
287 > #endif
288          if (what2do&DO_LINEAR) {
289                  if (scalef <= FTINY) {
290                          if (what2do&DO_HSENS)
# Line 293 | Line 312 | mapimage()                             /* map picture and send to stdout */
312  
313   getfovimg()                     /* load foveal sampled image */
314   {
315 <        extern FILE     *popen();
297 <        char    combuf[128];
315 >        char    combuf[PATH_MAX];
316          FILE    *fp;
317          int     x, y;
318                                                  /* compute image size */
# Line 302 | Line 320 | getfovimg()                    /* load foveal sampled image */
320          if (fvxr < 2) fvxr = 2;
321          fvyr = sqrt(ourview.vn2)/FOVDIA + 0.5;
322          if (fvyr < 2) fvyr = 2;
323 <        if (!(inpres.or & YMAJOR)) {            /* picture is rotated? */
323 >        if (!(inpres.rt & YMAJOR)) {            /* picture is rotated? */
324                  y = fvyr;
325                  fvyr = fvxr;
326                  fvxr = y;
327          }
328          if ((fovimg = (COLOR *)malloc(fvxr*fvyr*sizeof(COLOR))) == NULL)
329                  syserror("malloc");
330 <        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d %s", fvxr, fvyr, infn);
330 >        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d \"%s\"", fvxr, fvyr, infn);
331          if ((fp = popen(combuf, "r")) == NULL)
332                  syserror("popen");
333          getheader(fp, NULL, NULL);      /* skip header */
334 <        if (fgetresolu(&x, &y, fp) < 0 || x != fvxr | y != fvyr)
334 >        if (fgetresolu(&x, &y, fp) < 0 || (x != fvxr) | (y != fvyr))
335                  goto readerr;
336          for (y = 0; y < fvyr; y++)
337                  if (freadscan(fovscan(y), fvxr, fp) < 0)

Diff Legend

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