ImagerColorImpl.mesa
Copyright © 1984, 1985 by Xerox Corporation. All rights reserved.
Doug Wyatt, May 31, 1985 3:13:20 pm PDT
Stone, June 25, 1985 5:15:17 pm PDT
Michael Plass, July 23, 1985 1:52:36 pm PDT
DIRECTORY
Imager,
ImagerColor,
ImagerColorDefs,
ImagerColorOperator,
ImagerColorPrivate,
ImagerPixelArray USING [MaxSampleValue],
PrincOps,
Real,
RefTab,
RuntimeError;
ImagerColorImpl: CEDAR PROGRAM
IMPORTS Imager, ImagerColorOperator, ImagerPixelArray, Real, RefTab, RuntimeError
EXPORTS ImagerColor, ImagerColorDefs, ImagerColorPrivate
~ BEGIN OPEN ImagerColorPrivate, ImagerColor, ImagerColorDefs;
ConstantColorImpl: TYPE ~ ImagerColorPrivate.ConstantColorImpl;
ConstantColorImplRep: PUBLIC TYPE ~ ImagerColorPrivate.ConstantColorImplRep;
ConstantColorClass: TYPE ~ ImagerColorPrivate.ConstantColorClass;
ConstantColorClassRep: PUBLIC TYPE ~ ImagerColorPrivate.ConstantColorClassRep;
RGBCalibrationImpl: TYPE ~ ImagerColorPrivate.RGBCalibrationImpl;
RGBCalibrationImplRep: PUBLIC TYPE ~ ImagerColorPrivate.RGBCalibrationImplRep;
defaultCalibration: RGBCalibration ← CreateCalibration[xxxx];
defaultCalibration: RGBCalibration ←
NIL;
GetDefaultCalibration:
PUBLIC
PROC
RETURNS [RGBCalibration] ~ {
RETURN[defaultCalibration];
};
CreateCalibration:
PUBLIC
PROC [type:
ATOM, red, green, blue: CIEChromaticity,
white: CIEChromaticity, YMax:
REAL ← 100]
RETURNS [RGBCalibration] ~ {
ERROR;
};
CIEFromRGB:
PUBLIC
PROC [rgb:
RGB, calibration: RGBCalibration ←
NIL]
RETURNS [
CIE] ~ {
cal: RGBCalibration ~ IF calibration=NIL THEN defaultCalibration ELSE calibration;
impl: RGBCalibrationImpl ~ cal.impl;
RETURN[[
X: impl.cXR*rgb.R+impl.cXG*rgb.G+impl.cXB*rgb.B,
Y: impl.cYR*rgb.R+impl.cYG*rgb.G+impl.cYB*rgb.B,
Z: impl.cZR*rgb.R+impl.cZG*rgb.G+impl.cZB*rgb.B
]];
};
RGBFromCIE:
PUBLIC
PROC [cie:
CIE, calibration: RGBCalibration ←
NIL]
RETURNS [
RGB] ~ {
cal: RGBCalibration ~ IF calibration=NIL THEN defaultCalibration ELSE calibration;
impl: RGBCalibrationImpl ~ cal.impl;
RETURN[[
R: impl.cRX*cie.X+impl.cRY*cie.Y+impl.cRZ*cie.Z,
G: impl.cGX*cie.X+impl.cGY*cie.Y+impl.cGZ*cie.Z,
B: impl.cBX*cie.X+impl.cBY*cie.Y+impl.cBZ*cie.Z
]];
};
RGBMaxY:
PUBLIC
PROC [c: CIEChromaticity, calibration: RGBCalibration ←
NIL]
RETURNS [Y:
REAL] ~ {
We want to find the maximum value of Y such that c.x, c.y, Y is inside of the current RGB gamut. (x,y,z)*S=(X,Y,Z). (x,y,z)*S*CIEToRGB=(r,g,b)*S. We want to find the maximum value for S such that R, G and B <=1. Find the MAX[r,g,b]. S=1/max. Y=y*S.
NOTE: there is no chromaticity value for black (r=g=b=0). Black is Y=0.
cie: CIE ← [X: c.x, Y: c.y, Z: 1-(c.x+c.y)];
rgb: RGB ← RGBFromCIE[cie, calibration];
max: REAL ← MAX[MAX[rgb.R, rgb.G],rgb.B];
Y ← c.y/max; --it would be an unusual device that had max=0
};
ChromaticityFromCIE:
PUBLIC
PROC [c:
CIE]
RETURNS [CIEChromaticity] ~ {
sum: REAL ~ c.X+c.Y+c.Z;
RETURN[[x: c.X/sum, y: c.Y/sum]];
};
CIEFromChromaticity:
PUBLIC
PROC [c: CIEChromaticity, Y:
REAL]
RETURNS [
CIE] ~ {
scale: REAL ~ Y/c.y;
RETURN[[X: c.x*scale, Y: Y, Z: (1-c.x-c.y)*scale]];
};
ToRange:
PROC[v:
REAL]
RETURNS[
REAL] =
INLINE {
IF v IN[0..1] THEN RETURN[v] ELSE ERROR RuntimeError.BoundsFault;
};
ensures that v is in [0..1]; raises BoundsFault if not
These algorithms use the hexacone model described in
"Color Gamut Transform Pairs" by Alvy Ray Smith
Siggraph 1978, p. 12.
Algorithms from Foley and van Dam
HSLFromRGB:
PUBLIC
PROC [val:
RGB]
RETURNS [
HSL] ~ {
max,min,rc,gc,bc,del, h, s, l: REAL;
red: REAL ← ToRange[val.R];
green: REAL ← ToRange[val.G];
blue: REAL ← ToRange[val.B];
max ← MAX[red,MAX[green,blue]];
min ← MIN[red,MIN[green,blue]];
l ← (max+min)/2;
IF max=min THEN RETURN[[0,0,l]]; --gray
del ← max-min;
s ← IF l <= 0.5 THEN del/(max+min) ELSE del/(2-max-min);
rc ← (max-red)/del;
gc ← (max-green)/del;
bc ← (max-blue)/del;
IF max = red THEN h ← bc-gc --between yellow and magenta
ELSE IF max = green THEN h ← 2+rc-bc --between cyan and yellow
ELSE IF max = blue THEN h ← 4+gc-rc --between magenta and cyan
ELSE ERROR Imager.Error[[$invalidColor, "Invalid RGB color"]];
h ← h/6.0;
IF h < 0 THEN h ← h+1;
RETURN[[h, s, l]];
};
RGBFromHSL:
PUBLIC
PROC [val:
HSL]
RETURNS [
RGB] ~ {
m1,m2,hue,saturation, lightness, r, g, b: REAL;
Value:
PROC[n1,n2,h1:
REAL]
RETURNS[v:
REAL] = {
IF h1 > 360 THEN h1 ← h1-360;
IF h1 < 0 THEN h1 ← h1+360;
v ←
SELECT
TRUE
FROM
h1 IN [0..60) => n1+(n2-n1)*h1/60,
h1 IN [60..180) => n2,
h1 IN [180..240) => n1+(n2-n1)*(240-h1)/60,
ENDCASE => n1;
};
IF val.S=0 THEN RETURN[[val.L, val.L, val.L]];
saturation ← ToRange[val.S];
lightness ← ToRange[val.L];
hue ← 360*ToRange[val.H];
m2 ←
IF lightness <= 0.5
THEN lightness*(1+saturation)
ELSE lightness+saturation-lightness*saturation;
m1 ← 2*lightness-m2;
r ← Value[m1,m2,hue+120];
g ← Value[m1,m2,hue];
b ← Value[m1,m2,hue-120];
RETURN[[r, g, b]];
};
HSVFromRGB:
PUBLIC
PROC [val:
RGB]
RETURNS [
HSV] ~ {
r, g, b, h, s, v: REAL ← 0;
max,min,rc,gc,bc: REAL;
r ← ToRange[val.R]; g ← ToRange[val.G]; b ← ToRange[val.B];
min ← MIN[MIN[r,g],b]; --amount of white
v ← max ← MAX[MAX[r,g],b];--maximum "brightness"
IF max#0 THEN s ← (max-min)/max
ELSE s ← 0;
IF s=0 THEN RETURN[[0,0,v]]; --gray
rc ← (max - r)/(max - min);
gc ← (max - g)/(max - min);
bc ← (max - b)/(max - min);
IF r=max THEN h�-gc
ELSE IF g=max THEN h𡤂+rc-bc
ELSE IF b=max THEN h𡤄+gc-rc;
h ← h / 6.0;
IF h<0 THEN h←h+1;
RETURN[[h, s, v]];
};
RGBFromHSV:
PUBLIC
PROC [val:
HSV]
RETURNS [
RGB] ~ {
hue, saturation, value: REAL;
ihue: INTEGER;
fhue,m,n,k: REAL;
IF val.V=0 OR val.S=0 THEN RETURN[[val.V, val.V, val.V]];
hue ← ToRange[val.H];
saturation ← ToRange[val.S];
value ← ToRange[val.V];
hue ← hue*6;
ihue ← Real.FixI[hue]; --integer hue
fhue ← hue-ihue; --fractional hue
IF ihue=6 THEN ihue ← 0;
m ← value*(1-saturation);
n ← value*(1-(saturation*fhue));
k ← value*(1-(saturation*(1-fhue)));
SELECT ihue
FROM
0 => RETURN[[value,k,m]];
1 => RETURN[[n,value,m]];
2 => RETURN[[m,value,k]];
3 => RETURN[[m,n,value]];
4 => RETURN[[k,m,value]];
5 => RETURN[[value,m,n]];
ENDCASE => RETURN[[0,0,0]];
};
YIQFromRGB:
PUBLIC
PROC [val:
RGB]
RETURNS [
YIQ] ~ {
0.3, 0.59, 0.11
0.6, -0.28, -0.32
0.21, -0.52, 0.31
RETURN[[
Y: 0.30*val.R+0.59*val.G+0.11*val.B,
I: 0.60*val.R-0.28*val.G-0.32*val.B,
Q: 0.21*val.R-0.52*val.G+0.31*val.B
]];
};
RGBFromYIQ:
PUBLIC
PROC [val:
YIQ]
RETURNS [
RGB] ~ {
1.0, 0.9482623, 0.6240127
1.0, -0.2760664, -0.6398104
1.0, -1.10545, 1.729858
RETURN[[
R: 1.0*val.Y+0.9482623*val.I+0.6240127*val.Q,
G: 1.0*val.Y-0.2760664*val.I-0.6398104*val.Q,
B: 1.0*val.Y-1.10545*val.I+1.729858*val.Q
]];
};
colorTable: RefTab.Ref ~ RefTab.Create[mod: 31];
PutColor:
PROC [atom:
ATOM, color: ConstantColor] ~ {
impl: ConstantColorImpl ~ color.impl;
impl.atom ← atom;
[] ← RefTab.Store[x: colorTable, key: atom, val: color];
};
GetColor:
PROC [atom:
ATOM]
RETURNS [color: ConstantColor] ~ {
found: BOOL; val: REF;
[found: found, val: val] ← RefTab.Fetch[x: colorTable, key: atom];
IF found
THEN
WITH val
SELECT
FROM
color: ConstantColor => RETURN[color];
ENDCASE => ERROR;
RETURN[NIL];
};
IntensityFromGray:
PROC [f:
REAL]
RETURNS [
REAL] ~ {
IF f>=1 THEN RETURN[0];
IF f<=0 THEN RETURN[1];
RETURN[1-f];
};
NewGray:
PROC [f:
REAL]
RETURNS [ConstantColor] ~ {
impl: ConstantColorImpl ~ NEW[ConstantColorImplRep.gray ← [
Y: IntensityFromGray[f], variant: gray[f]]];
RETURN[NEW[ColorRep.constant ← [variant: constant[impl: impl]]]];
};
IntensityFromRGB:
PROC [val:
RGB]
RETURNS [
REAL] ~ {
Y: REAL ~ 0.30*val.R+0.59*val.G+0.11*val.B;
IF Y<=0 THEN RETURN[0];
IF Y>=1 THEN RETURN[1];
RETURN[Y];
};
NewRGB:
PROC [rgb:
RGB]
RETURNS [ConstantColor] ~ {
impl: ConstantColorImpl ~ NEW[ConstantColorImplRep.rgb ← [
Y: IntensityFromRGB[rgb], variant: rgb[val: rgb]]];
RETURN[NEW[ColorRep.constant ← [variant: constant[impl: impl]]]];
};
NewHSV:
PROC [hsv:
HSV]
RETURNS [ConstantColor] ~ {
RETURN[NewRGB[RGBFromHSV[hsv]]];
};
NewHSL:
PROC [hsl:
HSL]
RETURNS [ConstantColor] ~ {
RETURN[NewRGB[RGBFromHSL[hsl]]];
};
ColorFromAtom:
PUBLIC
PROC [atom:
ATOM]
RETURNS [ConstantColor] ~ {
RETURN[GetColor[atom]];
};
ColorFromGray:
PUBLIC
PROC [f:
REAL]
RETURNS [ConstantColor] ~ {
IF f<=0.0 THEN RETURN[GetColor[$White]];
IF f>=1.0 THEN RETURN[GetColor[$Black]];
IF f=0.5 THEN RETURN[GetColor[$Gray]];
RETURN[NewGray[f]];
};
ColorFromRGB:
PUBLIC
PROC [rgb:
RGB, calibration: RGBCalibration ←
NIL]
RETURNS [ConstantColor] ~ {
IF rgb.R=0 AND rgb.G=0 AND rgb.B=0 THEN RETURN[GetColor[$Black]];
IF rgb.R=1 AND rgb.G=1 AND rgb.B=1 THEN RETURN[GetColor[$White]];
IF rgb.R=1 AND rgb.G=0 AND rgb.B=0 THEN RETURN[GetColor[$Red]];
IF rgb.R=0 AND rgb.G=1 AND rgb.B=0 THEN RETURN[GetColor[$Green]];
IF rgb.R=0 AND rgb.G=0 AND rgb.B=1 THEN RETURN[GetColor[$Blue]];
IF rgb.R=0 AND rgb.G=1 AND rgb.B=1 THEN RETURN[GetColor[$Cyan]];
IF rgb.R=1 AND rgb.G=0 AND rgb.B=1 THEN RETURN[GetColor[$Magenta]];
IF rgb.R=1 AND rgb.G=1 AND rgb.B=0 THEN RETURN[GetColor[$Yellow]];
RETURN[NewRGB[rgb]];
};
ColorFromHSV:
PUBLIC
PROC [hsv:
HSV]
RETURNS [ConstantColor] ~ {
IF hsv.V=0 THEN RETURN[GetColor[$Black]];
IF hsv.V=1 AND hsv.S=0 THEN RETURN[GetColor[$White]];
RETURN[NewHSV[hsv]];
};
GrayFromColor:
PUBLIC
PROC [color: ConstantColor]
RETURNS [
REAL] ~ {
impl: ConstantColorImpl ~ color.impl;
RETURN[1-impl.Y];
};
AtomFromColor: PUBLIC PROC [color: ConstantColor] RETURNS [ATOM] ~ {
impl: ConstantColorImpl ~ color.impl;
RETURN[impl.atom];
};
ComponentFromColor:
PUBLIC
PROC [color: ConstantColor, name:
ATOM]
RETURNS [
REAL] ~ {
class: ConstantColorClass ~ color.class;
impl: REF ConstantColorImplRep ~ color.impl;
IF class#
NIL
AND class.ComponentFromColor #
NIL
THEN RETURN [class.ComponentFromColor[color, name]];
WITH impl
SELECT
FROM
stipple: REF ConstantColorImplRep.stipple => RETURN [impl.Y];
gray: REF ConstantColorImplRep.gray => RETURN [impl.Y];
rgb:
REF ConstantColorImplRep.rgb => {
SELECT name
FROM
$Intensity => RETURN [impl.Y];
$Red => RETURN [rgb.val.R];
$Green => RETURN [rgb.val.G];
$Blue => RETURN [rgb.val.B];
ENDCASE => RETURN [impl.Y];
};
ENDCASE => RETURN [impl.Y];
};
IntensityFromStipple:
PROC [word:
WORD]
RETURNS [
REAL] ~ {
nBits: NAT ~ 16;
bits: PACKED ARRAY [0..nBits) OF [0..1] ~ LOOPHOLE[word];
count: NAT ← 0; -- count the number of 1 bits
FOR i: NAT IN[0..nBits) DO count ← count+bits[i] ENDLOOP;
RETURN[REAL[nBits-count]/nBits];
};
ColorFromStipple:
PUBLIC
PROC [word:
WORD, function: StippleFunction]
RETURNS [ConstantColor] ~ {
impl: ConstantColorImpl ~ NEW[ConstantColorImplRep ← [
Y: IntensityFromStipple[word], variant: stipple[word: word, function: function]]];
RETURN[NEW[ColorRep.constant ← [variant: constant[impl: impl]]]];
};
MakeSampledBlack:
PUBLIC
PROC[
pa: PixelArray, um: Transformation, clear:
BOOL ←
FALSE
]
RETURNS[SampledColor] ~ {
IF pa.samplesPerPixel#1 THEN ERROR;
IF ImagerPixelArray.MaxSampleValue[pa, 0]#1 THEN ERROR;
RETURN[NEW[ColorRep.sampled ← [variant: sampled[
pa: pa, um: um, colorOperator: ImagerColorOperator.BlackColorModel[clear]]]]];
};
MakeSampledColor:
PUBLIC
PROC[
pa: PixelArray, um: Transformation, colorOperator: ColorOperator
]
RETURNS[SampledColor] ~ {
RETURN[NEW[ColorRep[sampled] ← [variant: sampled[
pa: pa, um: um, colorOperator: colorOperator]]]];
};
InitColorTable:
PROC ~ {
PutColor[$White, NewGray[0.0]];
PutColor[$Black, NewGray[1.0]];
PutColor[$Invert, ColorFromStipple[word: WORD.LAST, function: invert]];
PutColor[$Clear, ???];
PutColor[$Gray, NewGray[0.5]];
PutColor[$Red, NewRGB[[R: 1, G: 0, B: 0]]];
PutColor[$Green, NewRGB[[R: 0, G: 1, B: 0]]];
PutColor[$Blue, NewRGB[[R: 0, G: 0, B: 1]]];
PutColor[$Cyan, NewRGB[[R: 0, G: 1, B: 1]]];
PutColor[$Magenta, NewRGB[[R: 1, G: 0, B: 1]]];
PutColor[$Yellow, NewRGB[[R: 1, G: 1, B: 0]]];
PutColor[$Pink, NewHSL[[H: 0.0, S: 0.5, L: 0.7]]]; -- ???
PutColor[$Orange, NewHSL[[H: 0.04, S: 0.6, L: 0.4]]];
PutColor[$Brown, NewHSL[[H: 0.08, S: 0.6, L: 0.2]]];
PutColor[$Olive, NewHSL[[H: 0.25, S: 0.6, L: 0.2]]]; -- ???
PutColor[$YellowGreen, NewHSL[[H: 0.25, S: 0.6, L: 0.5]]]; -- ???
PutColor[$Purple, NewHSL[[H: 0.73, S: 0.6, L: 0.4]]];
};
InitColorTable[];
END.
CubeRootFn: PROC [r: REAL] RETURNS [REAL] ~ {
IF r>0.008856 THEN RETURN[RealFns.Root[index: 3, arg: r]]
ELSE RETURN[(903.29*r+16.0)/116.0];
};
LStar: PUBLIC PROC [Y: REAL] RETURNS [REAL] ~ {
RETURN[116.0*CubeRootFn[Y/100]-16.0];
};
LABFromCIE: PUBLIC PROC [val, illum: CIE] RETURNS [CIELAB] ~ {
fX: REAL ~ CubeRootFn[val.X/illum.X];
fY: REAL ~ CubeRootFn[val.Y/illum.Y];
fZ: REAL ~ CubeRootFn[val.Z/illum.Z];
RETURN[[L: 116*fY-16, a: 500*(fX-fY), b: 200*(fY-fZ)]];
};
CIEFromLAB: PUBLIC PROC [val: CIELAB, illum: CIE] RETURNS [CIE] ~ {
fY: REAL ~ (val.L+16)/116;
fX: REAL ~ fY+val.a/500;
fZ: REAL ~ fY-val.b/200;
RETURN[[X: illum.X*(fX*fX*fX), Y: illum.Y*(fY*fY*fY), Z: illum.Z*(fZ*fZ*fZ)]];
};
LUVFromCIE: PUBLIC PROC [val, illum: CIE] RETURNS [CIELUV] ~ {
den: REAL ~ -2*val.x+12*val.y+3;
uPrime: REAL ~ 4*val.x/den;
vPrime: REAL ~ 9*val.y/den;
RETURN[[L: xxx, u: xxx, v: xxx]];
};