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

Comparing ray/src/px/macbethcal.c (file contents):
Revision 2.16 by gwlarson, Fri Jul 16 16:12:20 1999 UTC vs.
Revision 2.25 by greg, Mon Nov 10 19:08:19 2008 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1999 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Calibrate a scanned MacBeth Color Checker Chart
6   *
# Line 15 | Line 12 | static char SCCSid[] = "$SunId$ SGI";
12  
13   #include <stdio.h>
14   #include <math.h>
15 < #ifdef MSDOS
16 < #include <fcntl.h>
17 < #endif
15 > #include <time.h>
16 >
17 > #include "platform.h"
18 > #include "rtprocess.h"
19 > #include "rtio.h"
20   #include "color.h"
21   #include "resolu.h"
22   #include "pmap.h"
23   #include "warp3d.h"
24 + #include "mx3.h"
25  
26                                  /* MacBeth colors */
27   #define DarkSkin        0
# Line 101 | Line 101 | short  mbneu[NMBNEU] = {Black,Neutral35,Neutral5,Neutra
101   #define  RG_CORR        04      /* corrected color region */
102  
103   #ifndef  DISPCOM
104 < #define  DISPCOM        "ximage -op %s"
104 > #define  DISPCOM        "ximage -op \"%s\""
105   #endif
106  
107   int     scanning = 1;           /* scanned input (or recorded output)? */
# Line 125 | Line 125 | WARP3D *wcor = NULL;           /* color space warp */
125   FILE    *debugfp = NULL;        /* debug output picture */
126   char    *progname;
127  
128 < extern char     *malloc();
128 > static void init(void);
129 > static int chartndx(int x, int y, int   *np);
130 > static void getpicture(void);
131 > static void getcolors(void);
132 > static void bresp(COLOR y, COLOR        x);
133 > static void compute(void);
134 > static void putmapping(void);
135 > static void compsoln(COLOR      cin[], COLOR    cout[], int     n);
136 > static void cwarp(void);
137 > static int cvtcolor(COLOR       cout, COLOR     cin);
138 > static int cresp(COLOR  cout, COLOR     cin);
139 > static void xyY2RGB(COLOR       rgbout, float   xyYin[3]);
140 > static void picdebug(void);
141 > static void clrdebug(void);
142 > static void getpos(char *name, int      bnds[2], FILE   *fp);
143 > static void pickchartpos(char   *pfn);
144  
145  
146 < main(argc, argv)
147 < int     argc;
148 < char    **argv;
146 > int
147 > main(
148 >        int     argc,
149 >        char    **argv
150 > )
151   {
152          int     i;
153  
# Line 145 | Line 162 | char   **argv;
162                                  perror(argv[i]);
163                                  exit(1);
164                          }
165 < #ifdef MSDOS
149 <                        setmode(fileno(debugfp), O_BINARY);
150 < #endif
165 >                        SET_FILE_BINARY(debugfp);
166                          newheader("RADIANCE", debugfp);         /* start */
167                          printargs(argc, argv, debugfp);         /* header */
168                          break;
# Line 192 | Line 207 | char   **argv;
207                  exit(1);
208          }
209          if (scanning) {                 /* load input picture header */
210 < #ifdef MSDOS
196 <                setmode(fileno(stdin), O_BINARY);
197 < #endif
210 >                SET_FILE_BINARY(stdin);
211                  if (checkheader(stdin, COLRFMT, NULL) < 0 ||
212                                  fgetresolu(&xmax, &ymax, stdin) < 0) {
213                          fprintf(stderr, "%s: bad input picture\n", progname);
# Line 227 | Line 240 | char   **argv;
240  
241                  printf("# Color correspondence produced by:\n#\t\t");
242                  printargs(argc, argv, stdout);
243 <                printf("#\tUsage: pcwarp %s uncorrected.pic > corrected.pic\n",
243 >                printf("#\tUsage: pcwarp %s uncorrected.hdr > corrected.hdr\n",
244                                  i+1 < argc ? argv[i+1] : "{this_file}");
245 <                printf("#\t   Or: pcond [options] -m %s orig.pic > output.pic\n",
245 >                printf("#\t   Or: pcond [options] -m %s orig.hdr > output.hdr\n",
246                                  i+1 < argc ? argv[i+1] : "{this_file}");
247                  for (j = 0; j < 24; j++)
248                          printf("%f %f %f    %f %f %f\n",
# Line 242 | Line 255 | char   **argv;
255                                                  /* print header */
256                  printf("{\n\tColor correction file computed by:\n\t\t");
257                  printargs(argc, argv, stdout);
258 <                printf("\n\tUsage: pcomb -f %s uncorrected.pic > corrected.pic\n",
258 >                printf("\n\tUsage: pcomb -f %s uncorrected.hdr > corrected.hdr\n",
259                                  i+1 < argc ? argv[i+1] : "{this_file}");
260                  if (!scanning)
261 <                        printf("\t   Or: pcond [options] -f %s orig.pic > output.pic\n",
261 >                        printf("\t   Or: pcond [options] -f %s orig.hdr > output.hdr\n",
262                                          i+1 < argc ? argv[i+1] : "{this_file}");
263                  printf("}\n");
264                  putmapping();                   /* put out color mapping */
265          }
266 <        if (debugfp != NULL)            /* put out debug picture */
266 >        if (debugfp != NULL) {          /* put out debug picture */
267                  if (scanning)
268                          picdebug();
269                  else
270                          clrdebug();
271 +        }
272          exit(0);
273   userr:
274          fprintf(stderr,
275 < "Usage: %s [-d dbg.pic][-P | -p xul yul xur yur xll yll xlr ylr][-i irrad][-m] input.pic [output.{cal|cwp}]\n",
275 > "Usage: %s [-d dbg.hdr][-P | -p xul yul xur yur xll yll xlr ylr][-i irrad][-m] input.hdr [output.{cal|cwp}]\n",
276                          progname);
277 <        fprintf(stderr, "   or: %s [-d dbg.pic][-i irrad][-m] -c [xyY.dat [output.{cal|cwp}]]\n",
277 >        fprintf(stderr, "   or: %s [-d dbg.hdr][-i irrad][-m] -c [xyY.dat [output.{cal|cwp}]]\n",
278                          progname);
279          exit(1);
280 +        return 1; /* pro forma return */
281   }
282  
283  
284 < init()                          /* initialize */
284 > static void
285 > init(void)                              /* initialize */
286   {
287          double  quad[4][2];
288          register int    i;
# Line 292 | Line 308 | init()                         /* initialize */
308   }
309  
310  
311 < int
312 < chartndx(x, y, np)                      /* find color number for position */
313 < int     x, y;
314 < int     *np;
311 > static int
312 > chartndx(                       /* find color number for position */
313 >        int     x,
314 >        int y,
315 >        int     *np
316 > )
317   {
318          double  ipos[3], cpos[3];
319          int     ix, iy;
# Line 324 | Line 342 | int    *np;
342   }
343  
344  
345 < getpicture()                            /* load in picture colors */
345 > static void
346 > getpicture(void)                                /* load in picture colors */
347   {
348          COLR    *scanln;
349          COLOR   pval;
# Line 362 | Line 381 | getpicture()                           /* load in picture colors */
381                  scalecolor(inpRGB[i], d);
382                  inpflags |= 1L<<i;
383          }
384 <        free((char *)scanln);
384 >        free((void *)scanln);
385   }
386  
387  
388 < getcolors()                     /* get xyY colors from standard input */
388 > static void
389 > getcolors(void)                 /* get xyY colors from standard input */
390   {
391          int     gotwhite = 0;
392          COLOR   whiteclr;
# Line 374 | Line 394 | getcolors()                    /* get xyY colors from standard input */
394          float   xyYin[3];
395  
396          while (fgetval(stdin, 'i', &n) == 1) {          /* read colors */
397 <                if (n < 0 | n > 24 ||
397 >                if ((n < 0) | (n > 24) ||
398                                  fgetval(stdin, 'f', &xyYin[0]) != 1 ||
399                                  fgetval(stdin, 'f', &xyYin[1]) != 1 ||
400                                  fgetval(stdin, 'f', &xyYin[2]) != 1 ||
401 <                                xyYin[0] < 0. | xyYin[1] < 0. ||
401 >                                (xyYin[0] < 0.) | (xyYin[1] < 0.) ||
402                                  xyYin[0] + xyYin[1] > 1.) {
403                          fprintf(stderr, "%s: bad color input data\n",
404                                          progname);
# Line 414 | Line 434 | getcolors()                    /* get xyY colors from standard input */
434   }
435  
436  
437 < bresp(y, x)             /* piecewise linear interpolation of primaries */
438 < COLOR   y, x;
437 > static void
438 > bresp(          /* piecewise linear interpolation of primaries */
439 >        COLOR   y,
440 >        COLOR   x
441 > )
442   {
443          register int    i, n;
444  
# Line 432 | Line 455 | COLOR  y, x;
455   }
456  
457  
458 < compute()                       /* compute color mapping */
458 > static void
459 > compute(void)                   /* compute color mapping */
460   {
461          COLOR   clrin[24], clrout[24];
462          long    cflags;
# Line 447 | Line 471 | compute()                      /* compute color mapping */
471                                          /* compute piecewise luminance curve */
472          for (i = 0; i < NMBNEU; i++) {
473                  copycolor(bramp[i][0], inpRGB[mbneu[i]]);
474 +                for (n = i ? 3 : 0; n--; )
475 +                        if (colval(bramp[i][0],n) <=
476 +                                        colval(bramp[i-1][0],n)+1e-7) {
477 +                                fprintf(stderr,
478 +                "%s: non-increasing neutral patch\n", progname);
479 +                                exit(1);
480 +                        }
481                  copycolor(bramp[i][1], mbRGB[mbneu[i]]);
482          }
483                                          /* compute color space gamut */
# Line 490 | Line 521 | compute()                      /* compute color mapping */
521   }
522  
523  
524 < putmapping()                    /* put out color mapping */
524 > static void
525 > putmapping(void)                        /* put out color mapping */
526   {
527          static char     cchar[3] = {'r', 'g', 'b'};
528          register int    i, j;
# Line 538 | Line 570 | putmapping()                   /* put out color mapping */
570   }
571  
572  
573 < compsoln(cin, cout, n)          /* solve 3xN system using least-squares */
574 < COLOR   cin[], cout[];
575 < int     n;
573 > static void
574 > compsoln(               /* solve 3xN system using least-squares */
575 >        COLOR   cin[],
576 >        COLOR   cout[],
577 >        int     n
578 > )
579   {
580          extern double   mx3d_adjoint(), fabs();
581          double  mat[3][3], invmat[3][3];
# Line 598 | Line 633 | int    n;
633   }
634  
635  
636 < cwarp()                         /* compute color warp map */
636 > static void
637 > cwarp(void)                             /* compute color warp map */
638   {
639          register int    i;
640  
# Line 614 | Line 650 | memerr:
650   }
651  
652  
653 < int
654 < cvtcolor(cout, cin)             /* convert color according to our mapping */
655 < COLOR   cout, cin;
653 > static int
654 > cvtcolor(               /* convert color according to our mapping */
655 >        COLOR   cout,
656 >        COLOR   cin
657 > )
658   {
659          COLOR   ctmp;
660          int     clipped;
# Line 635 | Line 673 | COLOR  cout, cin;
673   }
674  
675  
676 < int
677 < cresp(cout, cin)                /* transform color according to matrix */
678 < COLOR   cout, cin;
676 > static int
677 > cresp(          /* transform color according to matrix */
678 >        COLOR   cout,
679 >        COLOR   cin
680 > )
681   {
682          colortrans(cout, solmat, cin);
683          return(clipgamut(cout, bright(cout), CGAMUT, colmin, colmax));
684   }
685  
686  
687 < xyY2RGB(rgbout, xyYin)          /* convert xyY to RGB */
688 < COLOR   rgbout;
689 < register float  xyYin[3];
687 > static void
688 > xyY2RGB(                /* convert xyY to RGB */
689 >        COLOR   rgbout,
690 >        register float  xyYin[3]
691 > )
692   {
693          COLOR   ctmp;
694          double  d;
# Line 660 | Line 702 | register float xyYin[3];
702   }
703  
704  
705 < picdebug()                      /* put out debugging picture */
705 > static void
706 > picdebug(void)                  /* put out debugging picture */
707   {
708          static COLOR    blkcol = BLKCOLOR;
709          COLOR   *scan;
# Line 712 | Line 755 | picdebug()                     /* put out debugging picture */
755          }
756                                                  /* clean up */
757          fclose(debugfp);
758 <        free((char *)scan);
758 >        free((void *)scan);
759   }
760  
761  
762 < clrdebug()                      /* put out debug picture from color input */
762 > static void
763 > clrdebug(void)                  /* put out debug picture from color input */
764   {
765          static COLR     blkclr = BLKCOLR;
766          COLR    mbclr[24], cvclr[24], orclr[24];
# Line 782 | Line 826 | clrdebug()                     /* put out debug picture from color input
826          }
827                                                  /* clean up */
828          fclose(debugfp);
829 <        free((char *)scan);
829 >        free((void *)scan);
830   }
831  
832  
833 < getpos(name, bnds, fp)          /* get boundary position */
834 < char    *name;
835 < int     bnds[2];
836 < FILE    *fp;
833 > static void
834 > getpos(         /* get boundary position */
835 >        char    *name,
836 >        int     bnds[2],
837 >        FILE    *fp
838 > )
839   {
840          char    buf[64];
841  
# Line 803 | Line 849 | FILE   *fp;
849   }
850  
851  
852 < pickchartpos(pfn)               /* display picture and pick chart location */
853 < char    *pfn;
852 > static void
853 > pickchartpos(           /* display picture and pick chart location */
854 >        char    *pfn
855 > )
856   {
857 <        char    combuf[512];
857 >        char    combuf[PATH_MAX];
858          FILE    *pfp;
859  
860          sprintf(combuf, DISPCOM, pfn);

Diff Legend

– Removed lines
+ Added lines
< Changed lines (old)
> Changed lines (new)