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.11 by greg, Wed Mar 19 13:04:09 1997 UTC vs.
Revision 3.32 by greg, Fri Oct 27 16:35:02 2023 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 "pcond.h"
12  
13  
14   #define LDMAX           100             /* default max. display luminance */
15 < #define LDDYN           32              /* default dynamic range */
15 > #define LDDYN           100             /* default dynamic range */
16  
17   int     what2do = 0;                    /* desired adjustments */
18  
# 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 */
43  
44   double  scalef = 0.;                    /* linear scaling factor */
45  
46 + static gethfunc headline;
47 + static void getahead(void);
48 + static void mapimage(void);
49 + static void getfovimg(void);
50 + static void check2do(void);
51  
52 < main(argc, argv)
53 < int     argc;
54 < char    *argv[];
52 >
53 >
54 > int
55 > main(
56 >        int     argc,
57 >        char    *argv[]
58 > )
59   {
60          static RGBPRIMS outprimS;
61          int     i;
62 < #define bool(flg)               switch (argv[i][2]) { \
62 > #define check_bool(flg)         switch (argv[i][2]) { \
63                                  case '\0': what2do ^= flg; break; \
64                                  case 'y': case 'Y': case 't': case 'T': \
65                                  case '+': case '1': what2do |= flg; break; \
# Line 62 | Line 72 | char   *argv[];
72          for (i = 1; i < argc && argv[i][0] == '-'; i++)
73                  switch (argv[i][1]) {
74                  case 'h':
75 <                        bool(DO_HUMAN);
75 >                        check_bool(DO_HUMAN);
76                          break;
77                  case 'a':
78 <                        bool(DO_ACUITY);
78 >                        check_bool(DO_ACUITY);
79                          break;
80                  case 'v':
81 <                        bool(DO_VEIL);
81 >                        check_bool(DO_VEIL);
82                          break;
83                  case 's':
84 <                        bool(DO_HSENS);
84 >                        check_bool(DO_HSENS);
85                          break;
86                  case 'c':
87 <                        bool(DO_COLOR);
87 >                        check_bool(DO_COLOR);
88                          break;
89                  case 'w':
90 <                        bool(DO_CWEIGHT);
90 >                        check_bool(DO_CWEIGHT);
91                          break;
92                  case 'i':
93                          if (i+1 >= argc) goto userr;
# Line 85 | Line 95 | char   *argv[];
95                          if (fixfrac > FTINY) what2do |= DO_FIXHIST;
96                          else what2do &= ~DO_FIXHIST;
97                          break;
98 +                case 'I':
99 +                        check_bool(DO_PREHIST);
100 +                        break;
101                  case 'l':
102 <                        bool(DO_LINEAR);
102 >                        check_bool(DO_LINEAR);
103                          break;
104                  case 'p':
105                          if (i+8 >= argc) goto userr;
# Line 103 | Line 116 | char   *argv[];
116                  case 'e':
117                          if (i+1 >= argc) goto userr;
118                          scalef = atof(argv[++i]);
119 <                        if (argv[i][0] == '+' | argv[i][0] == '-')
119 >                        if ((argv[i][0] == '+') | (argv[i][0] == '-'))
120                                  scalef = pow(2.0, scalef);
121                          what2do |= DO_LINEAR;
122                          break;
# Line 137 | Line 150 | char   *argv[];
150                  default:
151                          goto userr;
152                  }
153 +        if ((what2do & (DO_FIXHIST|DO_PREHIST)) == (DO_FIXHIST|DO_PREHIST)) {
154 +                fprintf(stderr, "%s: only one of -i or -I option\n", progname);
155 +                exit(1);
156 +        }
157          if ((mbcalfile != NULL) + (cwarpfile != NULL) +
158                          (outprims != stdprims) > 1) {
159                  fprintf(stderr,
# Line 144 | Line 161 | char   *argv[];
161                                  progname);
162                  exit(1);
163          }
164 <        if (outprims == stdprims & inprims != stdprims)
164 >        if ((outprims == stdprims) & (inprims != stdprims))
165                  outprims = inprims;
166          Bldmin = Bl(ldmax/lddyn);
167          Bldmax = Bl(ldmax);
# Line 156 | Line 173 | char   *argv[];
173                                          /* open output file */
174          if (i+2 == argc && freopen(argv[i+1], "w", stdout) == NULL)
175                  syserror(argv[i+1]);
176 < #ifdef MSDOS
177 <        setmode(fileno(infp), O_BINARY);
161 <        setmode(fileno(stdout), O_BINARY);
162 < #endif
176 >        SET_FILE_BINARY(infp);
177 >        SET_FILE_BINARY(stdout);
178          getahead();                     /* load input header */
179          printargs(argc, argv, stdout);  /* add to output header */
180 <        if (mbcalfile == NULL & outprims != stdprims)
180 >        if ((mbcalfile == NULL) & (outprims != stdprims))
181                  fputprims(outprims, stdout);
182 <        getfovimg();                    /* get foveal sample image */
183 <        if (what2do&DO_FIXHIST)         /* get fixation history? */
182 >        if ((what2do & (DO_PREHIST|DO_VEIL|DO_ACUITY)) != DO_PREHIST)
183 >                getfovimg();            /* get foveal sample image? */
184 >        if (what2do&DO_PREHIST)         /* get histogram? */
185 >                gethisto(stdin);
186 >        else if (what2do&DO_FIXHIST)    /* get fixation history? */
187                  getfixations(stdin);
188          mapimage();                     /* map the picture */
189          if (mapfp != NULL)              /* write out basic mapping */
190                  putmapping(mapfp);
191          exit(0);
192   userr:
193 <        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",
193 >        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",
194                          progname);
195          exit(1);
196 < #undef bool
196 >        return 1; /* pro forma return */
197 > #undef check_bool
198   }
199  
200  
201 < syserror(s)                             /* report system error and exit */
202 < char    *s;
201 > void
202 > syserror(                               /* report system error and exit */
203 >        char    *s
204 > )
205   {
206          fprintf(stderr, "%s: ", progname);
207          perror(s);
# Line 188 | Line 209 | char   *s;
209   }
210  
211  
212 < headline(s)                             /* process header line */
213 < char    *s;
212 > static int
213 > headline(                               /* process header line */
214 >        char    *s,
215 >        void    *p
216 > )
217   {
218          static RGBPRIMS inprimS;
219 <        char    fmt[32];
219 >        char    fmt[MAXFMTLEN];
220  
221          if (formatval(fmt, s)) {        /* check if format string */
222                  if (!strcmp(fmt,COLRFMT)) lumf = rgblum;
223                  else if (!strcmp(fmt,CIEFMT)) lumf = cielum;
224                  else lumf = NULL;
225 <                return;                 /* don't echo */
225 >                return(0);              /* don't echo */
226          }
227 <        if (isprims(s)) {               /* get input primaries */
228 <                primsval(inprimS, s);
229 <                inprims= inprimS;
230 <                return;                 /* don't echo */
227 >                                        /* get input primaries */
228 >        if (isprims(s) && primsval(inprimS, s)) {
229 >                inprims = inprimS;
230 >                return(0);              /* don't echo */
231          }
232          if (isexpos(s)) {               /* picture exposure */
233                  inpexp *= exposval(s);
234 <                return;                 /* don't echo */
234 >                return(0);              /* don't echo */
235          }
236          if (isaspect(s))                /* pixel aspect ratio */
237                  pixaspect *= aspectval(s);
238          if (isview(s))                  /* image view */
239                  gotview += sscanview(&ourview, s);
240 <        fputs(s, stdout);
240 >        return(fputs(s, stdout));
241   }
242  
243  
244 < getahead()                      /* load picture header */
244 > static void
245 > getahead(void)                  /* load picture header */
246   {
247          char    *err;
248  
# Line 228 | Line 253 | getahead()                     /* load picture header */
253                  exit(1);
254          }
255          if (lumf == rgblum)
256 <                comprgb2xyzmat(inrgb2xyz, inprims);
256 >                comprgb2xyzWBmat(inrgb2xyz, inprims);
257          else if (mbcalfile != NULL) {
258                  fprintf(stderr, "%s: macbethcal only works with RGB pictures\n",
259                                  progname);
260                  exit(1);
261          }
262 <        if (!gotview || ourview.type == VT_PAR) {
263 <                copystruct(&ourview, &stdview);
262 >        if (!gotview || ourview.type == VT_PAR ||
263 >                        (ourview.horiz <= 5.) | (ourview.vert <= 5.)) {
264 >                ourview = stdview;
265                  ourview.type = VT_PER;
266                  if (pixaspect*inpres.yr < inpres.xr) {
267                          ourview.horiz = 40.0;
# Line 255 | Line 281 | getahead()                     /* load picture header */
281   }
282  
283  
284 < mapimage()                              /* map picture and send to stdout */
284 > static void
285 > mapimage(void)                          /* map picture and send to stdout */
286   {
287          COLOR   *scan;
288  
289          comphist();                     /* generate adaptation histogram */
290          check2do();                     /* modify what2do flags */
291          if (what2do&DO_VEIL)
292 <                compveil();
293 <        if (!(what2do&DO_LINEAR) && mkbrmap() < 0)      /* make tone map */
294 <                what2do |= DO_LINEAR;   /* failed! -- use linear scaling */
292 >                compveil();             /* compute veil image */
293 >        if (!(what2do&DO_LINEAR))
294 >                if (mkbrmap() < 0)      /* make tone map */
295 >                        what2do |= DO_LINEAR;   /* failed! -- use linear */
296 > #if ADJ_VEIL
297 >                else if (what2do&DO_VEIL)
298 >                        adjveil();      /* else adjust veil image */
299 > #endif
300          if (what2do&DO_LINEAR) {
301                  if (scalef <= FTINY) {
302                          if (what2do&DO_HSENS)
# Line 290 | Line 322 | mapimage()                             /* map picture and send to stdout */
322   }
323  
324  
325 < getfovimg()                     /* load foveal sampled image */
325 > static void
326 > getfovimg(void)                 /* load foveal sampled image */
327   {
328 <        extern FILE     *popen();
296 <        char    combuf[128];
328 >        char    combuf[PATH_MAX];
329          FILE    *fp;
330          int     x, y;
331                                                  /* compute image size */
# Line 301 | Line 333 | getfovimg()                    /* load foveal sampled image */
333          if (fvxr < 2) fvxr = 2;
334          fvyr = sqrt(ourview.vn2)/FOVDIA + 0.5;
335          if (fvyr < 2) fvyr = 2;
336 <        if (!(inpres.or & YMAJOR)) {            /* picture is rotated? */
336 >        if (!(inpres.rt & YMAJOR)) {            /* picture is rotated? */
337                  y = fvyr;
338                  fvyr = fvxr;
339                  fvxr = y;
340          }
341          if ((fovimg = (COLOR *)malloc(fvxr*fvyr*sizeof(COLOR))) == NULL)
342                  syserror("malloc");
343 <        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d %s", fvxr, fvyr, infn);
343 >        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d \"%s\"", fvxr, fvyr, infn);
344          if ((fp = popen(combuf, "r")) == NULL)
345                  syserror("popen");
346 +        SET_FILE_BINARY(fp);
347          getheader(fp, NULL, NULL);      /* skip header */
348 <        if (fgetresolu(&x, &y, fp) < 0 || x != fvxr | y != fvyr)
348 >        if (fgetresolu(&x, &y, fp) < 0 || (x != fvxr) | (y != fvyr))
349                  goto readerr;
350          for (y = 0; y < fvyr; y++)
351                  if (freadscan(fovscan(y), fvxr, fp) < 0)
# Line 326 | Line 359 | readerr:
359   }
360  
361  
362 < check2do()              /* check histogram to see what isn't worth doing */
362 > static void
363 > check2do(void)          /* check histogram to see what isn't worth doing */
364   {
365          double  sum;
366          double  b, l;
367 <        register int    i;
367 >        int     i;
368  
369                                          /* check for within display range */
370          if (bwmax - bwmin <= Bldmax - Bldmin)

Diff Legend

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