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

Comparing ray/src/px/macbethcal.c (file contents):
Revision 2.10 by greg, Wed Nov 1 09:04:16 1995 UTC vs.
Revision 2.13 by greg, Fri Jan 31 15:56:17 1997 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1995 Regents of the University of California */
1 > /* Copyright (c) 1997 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 95 | Line 95 | short  mbneu[NMBNEU] = {Black,Neutral35,Neutral5,Neutra
95   #define  RG_ORIG        02      /* original color region */
96   #define  RG_CORR        04      /* corrected color region */
97  
98 + int     scanning = 1;           /* scanned input (or recorded output)? */
99 +
100   int     xmax, ymax;             /* input image dimensions */
101   int     bounds[4][2];           /* image coordinates of chart corners */
102   double  imgxfm[3][3];           /* coordinate transformation matrix */
# Line 104 | Line 106 | long   inpflags = 0;           /* flags of which colors were inpu
106   long    gmtflags = 0;           /* flags of out-of-gamut colors */
107  
108   COLOR   bramp[NMBNEU][2];       /* brightness ramp (per primary) */
109 < double  solmat[3][3];           /* color mapping matrix */
109 > COLORMAT        solmat;         /* color mapping matrix */
110 > COLOR   colmin, colmax;         /* gamut limits */
111  
112   FILE    *debugfp = NULL;        /* debug output picture */
113   char    *progname;
# Line 116 | Line 119 | main(argc, argv)
119   int     argc;
120   char    **argv;
121   {
119 –        int     inpispic = 1;
122          int     i;
123  
124          progname = argv[0];
# Line 147 | Line 149 | char   **argv;
149                          bounds[2][1] = atoi(argv[++i]);
150                          bounds[3][0] = atoi(argv[++i]);
151                          bounds[3][1] = atoi(argv[++i]);
152 <                        inpispic = 2;
152 >                        scanning = 2;
153                          break;
154                  case 'c':                               /* color input */
155 <                        inpispic = 0;
155 >                        scanning = 0;
156                          break;
157                  default:
158                          goto userr;
159                  }
160                                                          /* open files */
161          if (i < argc && freopen(argv[i], "r", stdin) == NULL) {
162 <                perror(argv[1]);
162 >                perror(argv[i]);
163                  exit(1);
164          }
165          if (i+1 < argc && freopen(argv[i+1], "w", stdout) == NULL) {
166 <                perror(argv[2]);
166 >                perror(argv[i+1]);
167                  exit(1);
168          }
169 <        if (inpispic) {                 /* load input picture header */
169 >        if (scanning) {                 /* load input picture header */
170   #ifdef MSDOS
171                  setmode(fileno(stdin), O_BINARY);
172   #endif
# Line 177 | Line 179 | char   **argv;
179                  xmax = 512;
180                  ymax = 2*512/3;
181          }
182 <        if (inpispic != 2) {            /* use default boundaries */
182 >        if (scanning != 2) {            /* use default boundaries */
183                  bounds[0][0] = bounds[2][0] = .029*xmax + .5;
184                  bounds[0][1] = bounds[1][1] = .956*ymax + .5;
185                  bounds[1][0] = bounds[3][0] = .971*xmax + .5;
186                  bounds[2][1] = bounds[3][1] = .056*ymax + .5;
187          }
188          init();                         /* initialize */
189 <        if (inpispic)                   /* get picture colors */
189 >        if (scanning)                   /* get picture colors */
190                  getpicture();
191          else
192                  getcolors();
# Line 194 | Line 196 | char   **argv;
196          printargs(argc, argv, stdout);
197          printf("\n\tUsage: pcomb -f %s uncorrected.pic > corrected.pic\n",
198                          i+1 < argc ? argv[i+1] : "{this_file}");
199 +        if (!scanning)
200 +                printf("\t   Or: pcond [options] -f %s orig.pic > output.pic\n",
201 +                                i+1 < argc ? argv[i+1] : "{this_file}");
202          printf("}\n");
203          putmapping();                   /* put out color mapping */
204          if (debugfp != NULL)            /* put out debug picture */
205 <                if (inpispic)
205 >                if (scanning)
206                          picdebug();
207                  else
208                          clrdebug();
# Line 322 | Line 327 | getcolors()                    /* get xyY colors from standard input */
327                                  fgetval(stdin, 'f', &xyYin[0]) != 1 ||
328                                  fgetval(stdin, 'f', &xyYin[1]) != 1 ||
329                                  fgetval(stdin, 'f', &xyYin[2]) != 1 ||
330 <                                xyYin[0] < 0. | xyYin[0] > 1. |
331 <                                xyYin[1] < 0. | xyYin[1] > 1.) {
330 >                                xyYin[0] < 0. | xyYin[1] < 0. ||
331 >                                xyYin[0] + xyYin[1] > 1.) {
332                          fprintf(stderr, "%s: bad color input data\n",
333                                          progname);
334                          exit(1);
# Line 381 | Line 386 | compute()                      /* compute color mapping */
386          COLOR   clrin[24], clrout[24];
387          long    cflags;
388          COLOR   ctmp;
389 <        register int    i, j, n;
389 >        register int    i, n;
390                                          /* did we get what we need? */
391          if ((inpflags & REQFLGS) != REQFLGS) {
392                  fprintf(stderr, "%s: missing required input colors\n",
# Line 393 | Line 398 | compute()                      /* compute color mapping */
398                  copycolor(bramp[i][0], inpRGB[mbneu[i]]);
399                  copycolor(bramp[i][1], mbRGB[mbneu[i]]);
400          }
401 +                                        /* compute color space gamut */
402 +        if (scanning) {
403 +                copycolor(colmin, cblack);
404 +                copycolor(colmax, cwhite);
405 +        } else
406 +                for (i = 0; i < 3; i++) {
407 +                        colval(colmin,i) = colval(bramp[0][0],i) -
408 +                                colval(bramp[0][1],i) *
409 +                                (colval(bramp[1][0],i)-colval(bramp[0][0],i)) /
410 +                                (colval(bramp[1][1],i)-colval(bramp[1][0],i));
411 +                        colval(colmax,i) = colval(bramp[NMBNEU-2][0],i) +
412 +                                (1.-colval(bramp[NMBNEU-2][1],i)) *
413 +                                (colval(bramp[NMBNEU-1][0],i) -
414 +                                        colval(bramp[NMBNEU-2][0],i)) /
415 +                                (colval(bramp[NMBNEU-1][1],i) -
416 +                                        colval(bramp[NMBNEU-2][1],i));
417 +                }
418                                          /* compute color mapping */
419          do {
420                  cflags = inpflags & ~gmtflags;
# Line 406 | Line 428 | compute()                      /* compute color mapping */
428                  compsoln(clrin, clrout, n);
429                                                  /* check out-of-gamut colors */
430                  for (i = 0; i < 24; i++)
431 <                        if (cflags & 1L<<i) {
432 <                                cresp(ctmp, mbRGB[i]);
411 <                                for (j = 0; j < 3; j++)
412 <                                        if (colval(ctmp,j) <= 0. ||
413 <                                                colval(ctmp,j) >= 1.) {
414 <                                                gmtflags |= 1L<<i;
415 <                                                break;
416 <                                        }
417 <                        }
431 >                        if (cflags & 1L<<i && cvtcolor(ctmp, mbRGB[i]))
432 >                                gmtflags |= 1L<<i;
433          } while (cflags & gmtflags);
434          if (gmtflags & MODFLGS)
435                  fprintf(stderr,
# Line 437 | Line 452 | putmapping()                   /* put out color mapping for pcomb -f *
452                  for (i = 0; i < NMBNEU; i++)
453                          printf(",%g", colval(bramp[i][1],j));
454                  printf(");\n");
440 –                printf("%c = %ci(1);\n", cchar[j], cchar[j]);
455                  printf("%cfi(n) = if(n-%g, %d, if(%cxa(n+1)-%c, n, %cfi(n+1)));\n",
456                                  cchar[j], NMBNEU-1.5, NMBNEU-1, cchar[j],
457                                  cchar[j], cchar[j]);
458                  printf("%cndx = %cfi(1);\n", cchar[j], cchar[j]);
459 <                printf("%cn = ((%cxa(%cndx+1)-%c)*%cya(%cndx) + ",
460 <                                cchar[j], cchar[j], cchar[j],
461 <                                cchar[j], cchar[j], cchar[j]);
459 >                printf("%c%c = ((%cxa(%cndx+1)-%c)*%cya(%cndx) + ",
460 >                                cchar[j], scanning?'n':'o', cchar[j],
461 >                                cchar[j], cchar[j], cchar[j], cchar[j]);
462                  printf("(%c-%cxa(%cndx))*%cya(%cndx+1)) /\n",
463                                  cchar[j], cchar[j], cchar[j],
464                                  cchar[j], cchar[j]);
# Line 452 | Line 466 | putmapping()                   /* put out color mapping for pcomb -f *
466                                  cchar[j], cchar[j], cchar[j], cchar[j]);
467          }
468                                          /* print color mapping */
469 <        printf("ro = %g*rn + %g*gn + %g*bn ;\n",
470 <                        solmat[0][0], solmat[0][1], solmat[0][2]);
471 <        printf("go = %g*rn + %g*gn + %g*bn ;\n",
472 <                        solmat[1][0], solmat[1][1], solmat[1][2]);
473 <        printf("bo = %g*rn + %g*gn + %g*bn ;\n",
474 <                        solmat[2][0], solmat[2][1], solmat[2][2]);
469 >        if (scanning) {
470 >                printf("r = ri(1); g = gi(1); b = bi(1);\n");
471 >                printf("ro = %g*rn + %g*gn + %g*bn ;\n",
472 >                                solmat[0][0], solmat[0][1], solmat[0][2]);
473 >                printf("go = %g*rn + %g*gn + %g*bn ;\n",
474 >                                solmat[1][0], solmat[1][1], solmat[1][2]);
475 >                printf("bo = %g*rn + %g*gn + %g*bn ;\n",
476 >                                solmat[2][0], solmat[2][1], solmat[2][2]);
477 >        } else {
478 >                printf("r1 = ri(1); g1 = gi(1); b1 = bi(1);\n");
479 >                printf("r = %g*r1 + %g*g1 + %g*b1 ;\n",
480 >                                solmat[0][0], solmat[0][1], solmat[0][2]);
481 >                printf("g = %g*r1 + %g*g1 + %g*b1 ;\n",
482 >                                solmat[1][0], solmat[1][1], solmat[1][2]);
483 >                printf("b = %g*r1 + %g*g1 + %g*b1 ;\n",
484 >                                solmat[2][0], solmat[2][1], solmat[2][2]);
485 >        }
486   }
487  
488  
# Line 521 | Line 546 | int    n;
546   }
547  
548  
549 + int
550   cvtcolor(cout, cin)             /* convert color according to our mapping */
551   COLOR   cout, cin;
552   {
553          COLOR   ctmp;
554 +        int     clipped;
555  
556 <        bresp(ctmp, cin);
557 <        cresp(cout, ctmp);
558 <        if (colval(cout,RED) < 0.)
559 <                colval(cout,RED) = 0.;
560 <        if (colval(cout,GRN) < 0.)
561 <                colval(cout,GRN) = 0.;
562 <        if (colval(cout,BLU) < 0.)
563 <                colval(cout,BLU) = 0.;
556 >        if (scanning) {
557 >                bresp(ctmp, cin);
558 >                clipped = cresp(cout, ctmp);
559 >        } else {
560 >                clipped = cresp(ctmp, cin);
561 >                bresp(cout, ctmp);
562 >        }
563 >        return(clipped);
564   }
565  
566  
567 + int
568   cresp(cout, cin)                /* transform color according to matrix */
569   COLOR   cout, cin;
570   {
571 <        double  r, g, b;
572 <
545 <        r = colval(cin,0)*solmat[0][0] + colval(cin,1)*solmat[0][1]
546 <                        + colval(cin,2)*solmat[0][2];
547 <        g = colval(cin,0)*solmat[1][0] + colval(cin,1)*solmat[1][1]
548 <                        + colval(cin,2)*solmat[1][2];
549 <        b = colval(cin,0)*solmat[2][0] + colval(cin,1)*solmat[2][1]
550 <                        + colval(cin,2)*solmat[2][2];
551 <        setcolor(cout, r, g, b);
571 >        colortrans(cout, solmat, cin);
572 >        return(clipgamut(cout, bright(cout), CGAMUT, colmin, colmax));
573   }
574  
575  
# Line 563 | Line 584 | register float xyYin[3];
584          ctmp[0] = xyYin[0] * d;
585          ctmp[1] = xyYin[2];
586          ctmp[2] = (1. - xyYin[0] - xyYin[1]) * d;
587 <        cie_rgb(rgbout, ctmp);
587 >                                /* allow negative values */
588 >        colortrans(rgbout, xyz2rgbmat, ctmp);
589   }
590  
591  
# Line 600 | Line 622 | picdebug()                     /* put out debugging picture */
622                  for (x = 0; x < xmax; x++) {
623                          rg = chartndx(x, y, &i);
624                          if (rg == RG_CENT) {
625 <                                if (!(1L<<i & gmtflags) || (x+y)&07)
625 >                                if (!(1L<<i & gmtflags) || (x+y)&07) {
626                                          copycolor(scan[x], mbRGB[i]);
627 <                                else
627 >                                        clipgamut(scan[x], bright(scan[x]),
628 >                                                CGAMUT, cblack, cwhite);
629 >                                } else
630                                          copycolor(scan[x], blkcol);
631                          } else if (rg == RG_CORR)
632                                  cvtcolor(scan[x], scan[x]);
# Line 626 | Line 650 | clrdebug()                     /* put out debug picture from color input
650          static COLR     blkclr = BLKCOLR;
651          COLR    mbclr[24], cvclr[24], orclr[24];
652          COLR    *scan;
653 <        COLOR   ctmp;
653 >        COLOR   ctmp, ct2;
654          int     y, i;
655          register int    x, rg;
656                                                  /* convert colors */
657          for (i = 0; i < 24; i++) {
658 <                setcolr(mbclr[i], colval(mbRGB[i],RED),
659 <                                colval(mbRGB[i],GRN), colval(mbRGB[i],BLU));
658 >                copycolor(ctmp, mbRGB[i]);
659 >                clipgamut(ctmp, bright(ctmp), CGAMUT, cblack, cwhite);
660 >                setcolr(mbclr[i], colval(ctmp,RED),
661 >                                colval(ctmp,GRN), colval(ctmp,BLU));
662                  if (inpflags & 1L<<i) {
663 <                        setcolr(orclr[i], colval(inpRGB[i],RED),
664 <                                        colval(inpRGB[i],GRN),
665 <                                        colval(inpRGB[i],BLU));
640 <                        cvtcolor(ctmp, inpRGB[i]);
641 <                        setcolr(cvclr[i], colval(ctmp,RED),
663 >                        copycolor(ctmp, inpRGB[i]);
664 >                        clipgamut(ctmp, bright(ctmp), CGAMUT, cblack, cwhite);
665 >                        setcolr(orclr[i], colval(ctmp,RED),
666                                          colval(ctmp,GRN), colval(ctmp,BLU));
667 +                        bresp(ctmp, inpRGB[i]);
668 +                        colortrans(ct2, solmat, ctmp);
669 +                        clipgamut(ct2, bright(ct2), CGAMUT, cblack, cwhite);
670 +                        setcolr(cvclr[i], colval(ct2,RED),
671 +                                        colval(ct2,GRN), colval(ct2,BLU));
672                  }
673          }
674                                                  /* allocate scanline */

Diff Legend

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