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.13 by gwlarson, Tue Oct 27 09:08:26 1998 UTC vs.
Revision 3.29 by greg, Thu Aug 2 18:33:44 2018 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 "paths.h"
12   #include "pcond.h"
13  
14  
# 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 + static void getahead(void);
49 + static void mapimage(void);
50 + static void getfovimg(void);
51 + static void check2do(void);
52  
53 < main(argc, argv)
54 < int     argc;
55 < char    *argv[];
53 >
54 >
55 > int
56 > main(
57 >        int     argc,
58 >        char    *argv[]
59 > )
60   {
61          static RGBPRIMS outprimS;
62          int     i;
63 < #define bool(flg)               switch (argv[i][2]) { \
63 > #define check_bool(flg)         switch (argv[i][2]) { \
64                                  case '\0': what2do ^= flg; break; \
65                                  case 'y': case 'Y': case 't': case 'T': \
66                                  case '+': case '1': what2do |= flg; break; \
# Line 62 | Line 73 | char   *argv[];
73          for (i = 1; i < argc && argv[i][0] == '-'; i++)
74                  switch (argv[i][1]) {
75                  case 'h':
76 <                        bool(DO_HUMAN);
76 >                        check_bool(DO_HUMAN);
77                          break;
78                  case 'a':
79 <                        bool(DO_ACUITY);
79 >                        check_bool(DO_ACUITY);
80                          break;
81                  case 'v':
82 <                        bool(DO_VEIL);
82 >                        check_bool(DO_VEIL);
83                          break;
84                  case 's':
85 <                        bool(DO_HSENS);
85 >                        check_bool(DO_HSENS);
86                          break;
87                  case 'c':
88 <                        bool(DO_COLOR);
88 >                        check_bool(DO_COLOR);
89                          break;
90                  case 'w':
91 <                        bool(DO_CWEIGHT);
91 >                        check_bool(DO_CWEIGHT);
92                          break;
93                  case 'i':
94                          if (i+1 >= argc) goto userr;
# Line 86 | Line 97 | char   *argv[];
97                          else what2do &= ~DO_FIXHIST;
98                          break;
99                  case 'I':
100 <                        bool(DO_PREHIST);
100 >                        check_bool(DO_PREHIST);
101                          break;
102                  case 'l':
103 <                        bool(DO_LINEAR);
103 >                        check_bool(DO_LINEAR);
104                          break;
105                  case 'p':
106                          if (i+8 >= argc) goto userr;
# Line 106 | Line 117 | char   *argv[];
117                  case 'e':
118                          if (i+1 >= argc) goto userr;
119                          scalef = atof(argv[++i]);
120 <                        if (argv[i][0] == '+' | argv[i][0] == '-')
120 >                        if ((argv[i][0] == '+') | (argv[i][0] == '-'))
121                                  scalef = pow(2.0, scalef);
122                          what2do |= DO_LINEAR;
123                          break;
# Line 151 | Line 162 | char   *argv[];
162                                  progname);
163                  exit(1);
164          }
165 <        if (outprims == stdprims & inprims != stdprims)
165 >        if ((outprims == stdprims) & (inprims != stdprims))
166                  outprims = inprims;
167          Bldmin = Bl(ldmax/lddyn);
168          Bldmax = Bl(ldmax);
# Line 163 | Line 174 | char   *argv[];
174                                          /* open output file */
175          if (i+2 == argc && freopen(argv[i+1], "w", stdout) == NULL)
176                  syserror(argv[i+1]);
177 < #ifdef MSDOS
178 <        setmode(fileno(infp), O_BINARY);
168 <        setmode(fileno(stdout), O_BINARY);
169 < #endif
177 >        SET_FILE_BINARY(infp);
178 >        SET_FILE_BINARY(stdout);
179          getahead();                     /* load input header */
180          printargs(argc, argv, stdout);  /* add to output header */
181 <        if (mbcalfile == NULL & outprims != stdprims)
181 >        if ((mbcalfile == NULL) & (outprims != stdprims))
182                  fputprims(outprims, stdout);
183          if ((what2do & (DO_PREHIST|DO_VEIL|DO_ACUITY)) != DO_PREHIST)
184                  getfovimg();            /* get foveal sample image? */
# Line 185 | Line 194 | userr:
194          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",
195                          progname);
196          exit(1);
197 < #undef bool
197 >        return 1; /* pro forma return */
198 > #undef check_bool
199   }
200  
201  
202 < syserror(s)                             /* report system error and exit */
203 < char    *s;
202 > void
203 > syserror(                               /* report system error and exit */
204 >        char    *s
205 > )
206   {
207          fprintf(stderr, "%s: ", progname);
208          perror(s);
# Line 198 | Line 210 | char   *s;
210   }
211  
212  
213 < headline(s)                             /* process header line */
214 < char    *s;
213 > static int
214 > headline(                               /* process header line */
215 >        char    *s,
216 >        void    *p
217 > )
218   {
219          static RGBPRIMS inprimS;
220 <        char    fmt[32];
220 >        char    fmt[MAXFMTLEN];
221  
222          if (formatval(fmt, s)) {        /* check if format string */
223                  if (!strcmp(fmt,COLRFMT)) lumf = rgblum;
# Line 227 | Line 242 | char   *s;
242   }
243  
244  
245 < getahead()                      /* load picture header */
245 > static void
246 > getahead(void)                  /* load picture header */
247   {
248          char    *err;
249  
# Line 238 | Line 254 | getahead()                     /* load picture header */
254                  exit(1);
255          }
256          if (lumf == rgblum)
257 <                comprgb2xyzmat(inrgb2xyz, inprims);
257 >                comprgb2xyzWBmat(inrgb2xyz, inprims);
258          else if (mbcalfile != NULL) {
259                  fprintf(stderr, "%s: macbethcal only works with RGB pictures\n",
260                                  progname);
261                  exit(1);
262          }
263 <        if (!gotview || ourview.type == VT_PAR) {
264 <                copystruct(&ourview, &stdview);
263 >        if (!gotview || ourview.type == VT_PAR ||
264 >                        (ourview.horiz <= 5.) | (ourview.vert <= 5.)) {
265 >                ourview = stdview;
266                  ourview.type = VT_PER;
267                  if (pixaspect*inpres.yr < inpres.xr) {
268                          ourview.horiz = 40.0;
# Line 265 | Line 282 | getahead()                     /* load picture header */
282   }
283  
284  
285 < mapimage()                              /* map picture and send to stdout */
285 > static void
286 > mapimage(void)                          /* map picture and send to stdout */
287   {
288          COLOR   *scan;
289  
290          comphist();                     /* generate adaptation histogram */
291          check2do();                     /* modify what2do flags */
292          if (what2do&DO_VEIL)
293 <                compveil();
294 <        if (!(what2do&DO_LINEAR) && mkbrmap() < 0)      /* make tone map */
295 <                what2do |= DO_LINEAR;   /* failed! -- use linear scaling */
293 >                compveil();             /* compute veil image */
294 >        if (!(what2do&DO_LINEAR))
295 >                if (mkbrmap() < 0)      /* make tone map */
296 >                        what2do |= DO_LINEAR;   /* failed! -- use linear */
297 > #if ADJ_VEIL
298 >                else if (what2do&DO_VEIL)
299 >                        adjveil();      /* else adjust veil image */
300 > #endif
301          if (what2do&DO_LINEAR) {
302                  if (scalef <= FTINY) {
303                          if (what2do&DO_HSENS)
# Line 300 | Line 323 | mapimage()                             /* map picture and send to stdout */
323   }
324  
325  
326 < getfovimg()                     /* load foveal sampled image */
326 > static void
327 > getfovimg(void)                 /* load foveal sampled image */
328   {
329 <        extern FILE     *popen();
306 <        char    combuf[128];
329 >        char    combuf[PATH_MAX];
330          FILE    *fp;
331          int     x, y;
332                                                  /* compute image size */
# Line 311 | Line 334 | getfovimg()                    /* load foveal sampled image */
334          if (fvxr < 2) fvxr = 2;
335          fvyr = sqrt(ourview.vn2)/FOVDIA + 0.5;
336          if (fvyr < 2) fvyr = 2;
337 <        if (!(inpres.or & YMAJOR)) {            /* picture is rotated? */
337 >        if (!(inpres.rt & YMAJOR)) {            /* picture is rotated? */
338                  y = fvyr;
339                  fvyr = fvxr;
340                  fvxr = y;
341          }
342          if ((fovimg = (COLOR *)malloc(fvxr*fvyr*sizeof(COLOR))) == NULL)
343                  syserror("malloc");
344 <        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d %s", fvxr, fvyr, infn);
344 >        sprintf(combuf, "pfilt -1 -b -pa 0 -x %d -y %d \"%s\"", fvxr, fvyr, infn);
345          if ((fp = popen(combuf, "r")) == NULL)
346                  syserror("popen");
347 +        SET_FILE_BINARY(fp);
348          getheader(fp, NULL, NULL);      /* skip header */
349 <        if (fgetresolu(&x, &y, fp) < 0 || x != fvxr | y != fvyr)
349 >        if (fgetresolu(&x, &y, fp) < 0 || (x != fvxr) | (y != fvyr))
350                  goto readerr;
351          for (y = 0; y < fvyr; y++)
352                  if (freadscan(fovscan(y), fvxr, fp) < 0)
# Line 336 | Line 360 | readerr:
360   }
361  
362  
363 < check2do()              /* check histogram to see what isn't worth doing */
363 > static void
364 > check2do(void)          /* check histogram to see what isn't worth doing */
365   {
366          double  sum;
367          double  b, l;
368 <        register int    i;
368 >        int     i;
369  
370                                          /* check for within display range */
371          if (bwmax - bwmin <= Bldmax - Bldmin)

Diff Legend

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