128 |
|
static void getpicture(void); |
129 |
|
static void getcolors(void); |
130 |
|
static void bresp(COLOR y, COLOR x); |
131 |
+ |
static void ibresp(COLOR y, COLOR x); |
132 |
|
static void compute(void); |
133 |
|
static void putmapping(void); |
134 |
|
static void compsoln(COLOR cin[], COLOR cout[], int n); |
435 |
|
|
436 |
|
static void |
437 |
|
bresp( /* piecewise linear interpolation of primaries */ |
438 |
< |
COLOR y, |
439 |
< |
COLOR x |
438 |
> |
COLOR y, /* y is linear output */ |
439 |
> |
COLOR x /* x is non-linear input */ |
440 |
|
) |
441 |
|
{ |
442 |
|
int i, n; |
455 |
|
|
456 |
|
|
457 |
|
static void |
458 |
+ |
ibresp( /* inverse mapping (delinearization) */ |
459 |
+ |
COLOR x, /* x is non-linear output */ |
460 |
+ |
COLOR y /* y is linear input */ |
461 |
+ |
) |
462 |
+ |
{ |
463 |
+ |
int i, n; |
464 |
+ |
|
465 |
+ |
for (i = 0; i < 3; i++) { |
466 |
+ |
for (n = 0; n < NMBNEU-2; n++) |
467 |
+ |
if (colval(y,i) < colval(bramp[n+1][1],i)) |
468 |
+ |
break; |
469 |
+ |
colval(x,i) = ((colval(bramp[n+1][1],i) - colval(y,i)) * |
470 |
+ |
colval(bramp[n][0],i) + |
471 |
+ |
(colval(y,i) - colval(bramp[n][1],i)) * |
472 |
+ |
colval(bramp[n+1][0],i)) / |
473 |
+ |
(colval(bramp[n+1][1],i) - colval(bramp[n][1],i)); |
474 |
+ |
} |
475 |
+ |
} |
476 |
+ |
|
477 |
+ |
|
478 |
+ |
static void |
479 |
|
compute(void) /* compute color mapping */ |
480 |
|
{ |
481 |
|
COLOR clrin[24], clrout[24]; |
501 |
|
copycolor(bramp[i][1], mbRGB[mbneu[i]]); |
502 |
|
} |
503 |
|
/* compute color space gamut */ |
504 |
< |
if (scanning) { |
505 |
< |
copycolor(colmin, cblack); |
506 |
< |
copycolor(colmax, cwhite); |
507 |
< |
scalecolor(colmax, irrad); |
508 |
< |
} else |
509 |
< |
for (i = 0; i < 3; i++) { |
510 |
< |
colval(colmin,i) = colval(bramp[0][0],i) - |
489 |
< |
colval(bramp[0][1],i) * |
490 |
< |
(colval(bramp[1][0],i)-colval(bramp[0][0],i)) / |
491 |
< |
(colval(bramp[1][1],i)-colval(bramp[0][1],i)); |
492 |
< |
colval(colmax,i) = colval(bramp[NMBNEU-2][0],i) + |
493 |
< |
(1.-colval(bramp[NMBNEU-2][1],i)) * |
494 |
< |
(colval(bramp[NMBNEU-1][0],i) - |
495 |
< |
colval(bramp[NMBNEU-2][0],i)) / |
496 |
< |
(colval(bramp[NMBNEU-1][1],i) - |
497 |
< |
colval(bramp[NMBNEU-2][1],i)); |
498 |
< |
} |
504 |
> |
copycolor(colmin, cblack); |
505 |
> |
copycolor(colmax, cwhite); |
506 |
> |
scalecolor(colmax, irrad); |
507 |
> |
if (!scanning) { |
508 |
> |
ibresp(colmin, colmin); |
509 |
> |
ibresp(colmax, colmax); |
510 |
> |
} |
511 |
|
/* compute color mapping */ |
512 |
|
do { |
513 |
|
cflags = inpflags & ~gmtflags; |
514 |
|
n = 0; /* compute transform matrix */ |
515 |
< |
for (i = 0; i < 24; i++) |
516 |
< |
if (cflags & 1L<<i) { |
515 |
> |
for (i = 0; i < 24; i++) { |
516 |
> |
if (!(cflags & 1L<<i)) |
517 |
> |
continue; |
518 |
> |
if (scanning) { |
519 |
|
bresp(clrin[n], inpRGB[i]); |
520 |
|
copycolor(clrout[n], mbRGB[i]); |
521 |
< |
n++; |
521 |
> |
} else { |
522 |
> |
copycolor(clrin[n], inpRGB[i]); |
523 |
> |
ibresp(clrout[n], mbRGB[i]); |
524 |
|
} |
525 |
+ |
n++; |
526 |
+ |
} |
527 |
|
compsoln(clrin, clrout, n); |
528 |
|
for (i = 0; i < 24; i++) /* check gamut */ |
529 |
|
if (cflags & 1L<<i && cvtcolor(ctmp, inpRGB[i])) |
530 |
|
gmtflags |= 1L<<i; |
531 |
|
} while (cflags & gmtflags); |
532 |
+ |
|
533 |
|
if (gmtflags & MODFLGS) |
534 |
|
fprintf(stderr, |
535 |
|
"%s: warning - some moderate colors are out of gamut\n", |