man Graphics::ColorObject () - convert between color spaces

NAME

Graphics::ColorObject - convert between color spaces

SYNOPSIS

  use Graphics::ColorObject;

  # rgb to hsv
  $color = Graphics::ColorObject->new_RGB([$r, $g, $b]);
  ($h, $s, $v) = @{ $color->as_HSV() };

  # one rgb space to another (NTSC to PAL)
  $color = Graphics::ColorObject->new_RGB([$r, $g, $b], space=>'NTSC');
  ($r, $g, $b) = @{ $color->as_RGB(space=>'PAL') };

ABSTRACT

Use this module to convert between all the common color spaces. As a pure Perl module, it is not very fast, and so it you want to convert entire images quickly, this is probably not what you want. The emphasis is on completeness and accurate conversion.

Supported color spaces are: RGB (including sRGB, Apple, Adobe, CIE Rec 601, CIE Rec 709, CIE ITU, and about a dozen other RGB spaces), CMY, CMYK, HSL, HSV, XYZ, xyY, Lab, LCHab, Luv, LCHuv, YPbPr, YCbCr, YUV, YIQ, PhotoYCC.

Conversion between different RGB working spaces, and between different white-points, is fully supported.

DESCRIPTION

For any supported color space XXX, there is one constructor new_XXX that creates a color using data in that color space, and one method as_XXX that returns the current color as expressed in that color space. For example, for RGB there is new_RGB and as_RGB. The color data is always passed as an array reference to a three-element array (four-element in the case of CMYK). Thus, to convert from RGB to HSL, you can use:

  $color = Graphics::ColorObject->new_RGB([$r, $g, $b]);
  ($h, $s, $l) = @{ $color->as_HSL() };

The constructor can always take a hash of optional arguments in addition to the color value, namely the working RGB space and the white point. For example:

  $color = Graphics::ColorObject->new_RGB([$r, $g, $b], space=>'Adobe', white_point=>'D65');

For a list of all supported color spaces, call Graphics::ColorObject->list_colorspaces(). For a list of all RGB working spaces and of all white points that this module supports, call Graphics::ColorObject->list_rgb_spaces() and Graphics::ColorObject->list_white_points().

If not specified, the working RGB space will be sRGB. Many non-RGB conversions also rely on an implicit RGB space, and passing an RGB space as an option (either to the constructor or later) will have an effect on the values.

VARIOUS NOTES AND GOTCHAS

Most conversions will return out-of-gamut values if necessary, because that way they are lossless and can be chained in calculations, or reversed to produce the original values. Many conversion methods will take an optional boolean clip parameter to restrict the returned values to be within gamut:

  ($r, $g, $b) = @{ $color->as_RGB(space=>'sRGB', clip=>1) };

Currently clipping is supported in RGB, RGB-derived (HSL, CMY) and chroma-luma separated (YUV, etc) spaces, but not in XYZ-derived spaces. The only way to check whether a value is within gamut is to convert it with and without the clip option and compare the two results. An RGB value is within gamut simply if R, G and B are between 0 and 1, but other spaces can be much harder to check.

RGB values are non-linear (gamma-adjusted) floating-point values scaled in the range from 0 to 1. If you want integer values in the range 0..255, use the new_RGB255/as_RGB255 functions instead. If you want linear RGB (not gamma-adjusted) use RGB_to_linear_RGB([$r, CW$g, CW$b]).

Functions that use an angle value always express it in degrees from 0 to 360. That includes the hue H in HSL, HSV, LCHab and LCHuv. Use rad2deg and deg2rad from Math::Trig to convert to/from degrees if necessary.

There is some confusion in the naming of YUV and related (Y-something-something) colorspaces. Most of the time when YUV or YCC is used in software, for example in JPEG and MPEG2, that is actually YCbCr, a chroma-luma separated space with integer values of Y in the range [16..235], Cb and Cr in [16..240]. JPEG uses a modified YCbCr with values in [0..255] (which is not implemented in this module). As used here, YUV is a floating-point representation of the analog signal in PAL TV, YIQ is the same for NTSC TV, YPbPr is component analog video, and PhotoYCC or YCC is the Kodak PhotoCD standard.

The set_white_point() function can take arbitrary temperatures as well as the predefined standard illuminants. The valid range of temperatures is from 4000K to 25000K.

RECOMMENDATIONS

Aside from converting from one space to another, what colorspace is the best one to use for a particular task? This section attempts to answer that question.

For generic RGB values, use sRGB (which is the default).

For 2D effects filters, use Lab (or LCHab).

For adjustment of brightness, saturation and hue, use LCHab or LSHab.

For compression, use YCbCr, or use YPbPr and convert to integer values in a way that makes sense in your application.

For representing data as colors, use Lab (straight lines between points in Lab are more-or-less uniform gradients, unlike straight lines in RGB, for example).

UPGRADING FROM 0.3a2 AND OLDER VERSIONS

Version 0.4 and later are a complete rewrite from the previous major version, 0.3a2. The API is completely changed. The old API should be emulated exactly in all cases. Please test any code that uses this module when upgrading. If you encounter any strange behavior, please downgrade to 0.3a2 and email me a bug report. Additionally, the exact values returned by some functions may be slightly different, this is not a bug - the new values are (more) correct.

METHODS

sub new_XYZ { my ($pkgname, CW$xyz, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = CW$xyz; return CW$this; } sub new_xyY { my ($pkgname, CW$xyy, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &xyY_to_XYZ($xyy); return CW$this; } sub new_RGB { my ($pkgname, CW$rgb, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &RGB_to_XYZ($rgb, CW$this->{space}); return CW$this; } sub new_RGB255 { my ($pkgname, CW$rgb255, CW%opts) = CW@_; return &new_RGB($pkgname, &RGB255_to_RGB($rgb255), CW%opts); }

$color = Graphics::ColorObject->new_RGBhex($rgbhex) =cut

sub new_RGBhex { my ($pkgname, CW$rgbhex, CW%opts) = CW@_; return &new_RGB($pkgname, &RGBhex_to_RGB($rgbhex), CW%opts); } sub new_Lab { my ($pkgname, CW$lab, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &Lab_to_XYZ($lab, CW$this->get_XYZ_white()); return CW$this; } sub new_LCHab { my ($pkgname, CW$lch, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &Lab_to_XYZ(&LCHab_to_Lab($lch), CW$this->get_XYZ_white()); return CW$this; } sub new_Luv { my ($pkgname, CW$luv, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &Luv_to_XYZ($luv, CW$this->get_XYZ_white()); return CW$this; } sub new_LCHuv { my ($pkgname, CW$lch, CW%opts) = CW@_; my CW$this = &new($pkgname, CW%opts); CW$this->{xyz} = &Luv_to_XYZ(&LCHuv_to_Luv($lch), CW$this->get_XYZ_white()); return CW$this; } sub new_HSL { my ($pkgname, CW$hsl, CW%opts) = CW@_; return &new_RGB($pkgname, &HSL_to_RGB($hsl), CW%opts); } sub new_HSV { my ($pkgname, CW$hsv, CW%opts) = CW@_; return &new_RGB($pkgname, &HSV_to_RGB($hsv), CW%opts); } sub new_CMY { my ($pkgname, CW$cmy, CW%opts) = CW@_; return &new_RGB($pkgname, &CMY_to_RGB($cmy), CW%opts); } sub new_CMYK { my ($pkgname, CW$cmyk, CW%opts) = CW@_; return &new_RGB($pkgname, &CMY_to_RGB(&CMYK_to_CMY($cmyk)), CW%opts); } sub new_YPbPr { my ($pkgname, CW$ypbpr, CW%opts) = CW@_; return &new_RGB($pkgname, &YPbPr_to_RGB($ypbpr), space => 'NTSC'); # force NTSC } sub new_YCbCr { my ($pkgname, CW$ycbcr, CW%opts) = CW@_; return &new_RGB($pkgname, &YCbCr_to_RGB($ycbcr), space => 'NTSC'); # force NTSC } sub new_YUV { my ($pkgname, CW$yuv, CW%opts) = CW@_; return &new_RGB($pkgname, &YUV_to_RGB($yuv), space => 'NTSC'); # force NTSC } sub new_YIQ { my ($pkgname, CW$yiq, CW%opts) = CW@_; return &new_RGB($pkgname, &YIQ_to_RGB($yiq), space => 'NTSC'); # force NTSC } sub new_PhotoYCC { my ($pkgname, CW$ycc, CW%opts) = CW@_; return &new_RGB($pkgname, &PhotoYCC_to_RGB($ycc), space => 'sRGB'); # force sRGB } sub as_XYZ { my ($this, CW%opts) = CW@_; my CW$xyz = CW$this->{xyz}; if ($opts{clip}) { # TODO check this is correct my ($Xw, CW$Yw, CW$Zw) = @{ CW$this->get_XYZ_white() }; CW$xyz = &_generic_clip($xyz, [[0,$Xw], [0,$Yw], [0,$Zw]]); } return CW$xyz; } sub as_RGB { my ($this, CW%opts) = CW@_; my CW$space = CW$opts{space} || CW$this->{space}; my CW$rgb = &XYZ_to_RGB($this->{xyz}, CW$space); if ($opts{clip}) { CW$rgb = &_generic_clip($rgb, [[0,1], [0,1], [0,1]]); }; return CW$rgb; } sub as_RGB255 { my ($this) = CW@_; # always clipped return &RGB_to_RGB255($this->as_RGB()); } sub as_RGBhex { my ($this) = CW@_; # always clipped return &RGB_to_RGBhex($this->as_RGB()); } sub as_xyY { my ($this, CW%opts) = CW@_; my CW$xyy = &XYZ_to_xyY($this->{xyz}, CW$this->get_XYZ_white()); return CW$xyy; } sub as_Lab { my ($this) = CW@_; my CW$lab = &XYZ_to_Lab($this->{xyz}, CW$this->get_XYZ_white()); return CW$lab; } sub as_LCHab { my ($this) = CW@_; my CW$lchab = &Lab_to_LCHab( &XYZ_to_Lab($this->{xyz}, CW$this->get_XYZ_white()) ); return CW$lchab; } sub as_Luv { my ($this) = CW@_; my CW$luv = &XYZ_to_Luv($this->{xyz}, CW$this->get_XYZ_white()); return CW$luv; } sub as_LCHuv { my ($this) = CW@_; my CW$lchuv = &Luv_to_LCHuv( &XYZ_to_Luv($this->{xyz}, CW$this->get_XYZ_white()) ); return CW$lchuv; } sub as_HSL { my ($this, CW%opts) = CW@_; my CW$hsl = &RGB_to_HSL( CW$this->as_RGB() ); if ($opts{clip}) { CW$hsl = &_generic_clip($hsl, [[0,360], [0,1], [0,1]]); }; return CW$hsl; } sub as_HSV { my ($this, CW%opts) = CW@_; my CW$hsv = &RGB_to_HSV( CW$this->as_RGB() ); if ($opts{clip}) { CW$hsv = &_generic_clip($hsv, [[0,360], [0,1], [0,1]]); }; return CW$hsv; } sub as_CMY { my ($this, CW%opts) = CW@_; my CW$cmy = &RGB_to_CMY( CW$this->as_RGB() ); if ($opts{clip}) { CW$cmy = &_generic_clip($cmy, [[0,1], [0,1], [0,1]]); }; return CW$cmy; } sub as_CMYK { my ($this) = CW@_; my CW$cmyk = &CMY_to_CMYK( &RGB_to_CMY( CW$this->as_RGB() ) ); # TODO clip return CW$cmyk; } sub as_YPbPr { my ($this, CW%opts) = CW@_; my CW$ypbpr = &RGB_to_YPbPr( CW$this->as_RGB( space => 'NTSC' ) ); if ($opts{clip}) { CW$ypbpr = &_generic_clip($ypbpr, [[0,1], [-0.5,0.5], [-0.5,0.5]]); }; return CW$ypbpr; } sub as_YCbCr { my ($this, CW%opts) = CW@_; my CW$ycbcr = &RGB_to_YCbCr( CW$this->as_RGB( space => 'NTSC' ) ); if ($opts{clip}) { CW$ycbcr = &_generic_clip($ycbcr, [[16,235], [16,239], [16,239]]); }; # TODO round to integers return CW$ycbcr; } sub as_YUV { my ($this) = CW@_; my CW$yuv = &RGB_to_YUV( CW$this->as_RGB( space => 'NTSC' ) ); return CW$yuv; } sub as_YIQ { my ($this) = CW@_; my CW$yiq = &RGB_to_YIQ( CW$this->as_RGB( space => 'NTSC' ) ); return CW$yiq; } sub as_PhotoYCC { my ($this) = CW@_; my CW$ycc = &RGB_to_PhotoYCC( CW$this->as_RGB( space => 'sRGB' ) ); return CW$ycc; }

# returns the XYZ value of the white point actually used (always defined, default is D65) sub get_XYZ_white { my ($this, CW%opts) = CW@_; my CW$white_point = CW$opts{white_point} || CW$this->{white_point} || &_get_RGB_space_by_name( CW$opts{space} || CW$this->{space} )->{white_point};

        $white_point = &_check_white_point($white_point);

        my $xy = $WHITE_POINTS{ $white_point };

        my ($x, $y) = @{ $xy };
        return &xyY_to_XYZ([$x, $y, 1.0]);
        #return &RGB_to_XYZ([1, 1, 1], $this->{space});
}
# returns the name of the white point actually used # FIXME should be always defined sub get_white_point { my ($this) = CW@_; return CW$this->{white_point}; } sub set_white_point { my ($this, CW$white_point) = CW@_;

        $white_point = &_check_white_point($white_point);

        if (&_check_white_point($this->{white_point}) ne $white_point)
        {
                $this->{xyz} = &XYZ_change_white_point($this->{xyz}, $this->get_XYZ_white(), $this->get_XYZ_white($white_point));
                $this->{white_point} = $white_point;
        }

        return $this;
}
# FIXME should be always defined sub get_rgb_space { my ($this) = CW@_; return CW$this->{space}; } sub set_rgb_space { my ($this, CW$space) = CW@_; my CW$s = &_get_RGB_space_by_name($space); if ($this->get_white_point() ne CW$s->{white_point}) { CW$this->set_white_point($s->{white_point}); } CW$this->{space} = CW$space; return CW$this; } sub copy { my ($this) = CW@_; my CW$copy = +{ xyz => CW$this->{xyz}, space => CW$this->{space}, white_point => CW$this->{white_point} }; bless CW$copy, ref CW$this; return CW$copy; }

if ($color->equals($color2)) { ... } Checks if another color is the same as this one. Optionally takes an accuracy parameter which is the distance between the two colors as measured by the city-block metric in XYZ space (default accuracy is 0.01%). =cut

sub equals { my ($this, CW$other, CW%opts) = CW@_; CW$other = CW$other->copy(); CW$other->set_white_point($this->{white_point}); CW$other->set_rgb_space($this->{space}); my CW$accuracy = CW$opts{accuracy} || 0.0001; if (&_delta_v3($this->{xyz}, CW$other->{xyz}) < CW$accuracy) { return 1; } else { return 0; } } sub difference { my ($this, CW$other) = CW@_; return CW$this->difference_CIE1976($other); }

# reference: http://www.brucelindbloom.com/index.html?Eqn_DeltaE_CIE76.html sub difference_CIE1976 { my ($this, CW$other) = CW@_;

        my ($L1, $a1, $b1) = @{ $this->as_Lab() };
        my ($L2, $a2, $b2) = @{ $other->as_Lab() };

        my $deltaE = sqrt(&_sqr($L1-$L2) + &_sqr($a1-$a2) + &_sqr($b1-$b2));

        return $deltaE;
}

sub difference_CIE1994 { # TODO }

# reference: http://www.brucelindbloom.com/index.html?Eqn_DeltaE_CMC.html sub difference_CMC { my ($this, CW$other, CW%opts) = CW@_;

        my $l = $opts{l} || 1;
        my $c = $opts{c} || 1;

        my ($L1, $a1, $b1) = @{ $this->as_Lab() };
        my ($L2, $a2, $b2) = @{ $other->as_Lab() };

        my $C1 = sqrt($a1*$a1 + $b1*$b1);
        my $C2 = sqrt($a2*$a2 + $b2*$b2);

        my $dH = sqrt(&_sqr($a1-$a2) + &_sqr($b1-$b2) - &_sqr($C1-$C2));

        my $SL = ($L1 < 16 ? 
                          0.511 : 
                          0.040975 * $L1 / ( 1 + 0.01765 * $L1)
                          );

        my $SC = 0.638 + 0.0638 * $C1 / ( 1 + 0.0131 * $C1 );

        my $F = sqrt(pow($C1, 4) / ( pow($C1, 4) + 1900 ));

        my $H1 = atan2($b1, $a1);

        my $T = ((deg2rad(1) <= $H1 && $H1 <= deg2rad(3)) ? 
                         0.56 + abs(0.2 * cos($H1 + deg2rad(1))) :
                         0.36 + abs(0.4 * cos($H1 + deg2rad(3)))
                         );

        my $SH = $SC * ($F*$T - $F + 1);

        my $deltaE = sqrt(&_sqr(($L1 - $L2)/($l * $SL)) + 
                                          &_sqr(($C1 - $C2)/($c * $SC)) + 
                                          &_sqr($dH/$SH)
                                          );

        return $deltaE;
}

@colorspaces = &Graphics::ColorObject->list_colorspaces() Returns a list of all supported colorspaces. =cut

sub list_colorspaces { return qw(RGB XYZ xyY Lab LCHab Luv LCHuv HSL HSV CMY CMYK YCbCr YPbPr YUV YIQ); # PhotoYCC }

@rgb_spaces = &Graphics::ColorObject->list_rgb_spaces() Returns a list of all supported RGB spaces. Some items are aliases, so the same space may be listed more than once under different names. =cut

sub list_rgb_spaces { return sort keys CW%RGB_SPACES; }

@white_points = &Graphics::ColorObject->list_white_points() Returns a list of all supported white points. =cut

sub list_white_points { return sort keys CW%WHITE_POINTS; }

############# non-OO interface ###########

sub RGB_to_RGB255 { my ($rgb) = CW@_; my ($r, CW$g, CW$b) = @{$rgb}; if ($r < 0) { CW$r = 0; } elsif ($r > 1) { CW$r = 1; } if ($g < 0) { CW$g = 0; } elsif ($g > 1) { CW$g = 1; } if ($b < 0) { CW$b = 0; } elsif ($b > 1) { CW$b = 1; } # FIXME use round, not sprintf return [ sprintf('%.0f', 255*$r), sprintf('%.0f', 255*$g), sprintf('%.0f', 255*$b) ]; }

sub RGB255_to_RGB { my ($rgb255) = CW@_; my ($r, CW$g, CW$b) = @{$rgb255}; return [ CW$r/255, CW$g/255, CW$b/255 ]; }

sub RGBhex_to_RGB { my ($rgbhex) = CW@_; my ($r, CW$g, CW$b); if ($rgbhex =~ m!^\#([0-9a-fA-F]{6})!) { CW$rgbhex = CW$1; } if ($rgbhex =~ m!^[0-9a-fA-F]{6}$!) { CW$r=hex(substr($rgbhex,0,2)); CW$g=hex(substr($rgbhex,2,2)); CW$b=hex(substr($rgbhex,4,2)); } return &RGB255_to_RGB([$r, CW$g, CW$b]); # return &RGB255_to_RGB([ unpack(C*,pack(N,hex($rgbhex)<<8)) ]); }

sub RGB_to_RGBhex { my ($rgb) = CW@_; my CW$rgb255 = &RGB_to_RGB255($rgb); return sprintf('%02X%02X%02X', @{$rgb255}); }

sub RGB_to_XYZ { my ($rgb, CW$space) = CW@_; my CW$s = &_get_RGB_space_by_name($space); my CW$rgb_lin = &RGB_to_linear_RGB($rgb, CW$space); my CW$xyz = &_mult_v3_m33($rgb_lin, CW$s->{m}); return ($xyz); }

sub XYZ_to_RGB { my ($xyz, CW$space) = CW@_; my CW$s = &_get_RGB_space_by_name($space); my CW$rgb_lin = &_mult_v3_m33($xyz, CW$s->{mstar}); my CW$rgb = &linear_RGB_to_RGB($rgb_lin, CW$space); return ($rgb); }

sub XYZ_to_Lab { my ($xyz, CW$xyz_white) = CW@_; my ($X, CW$Y, CW$Z) = @{$xyz}; my ($Xw, CW$Yw, CW$Zw) = @{$xyz_white}; my ($L, CW$a, CW$b);

        my $epsilon =  0.008856;
        my $kappa = 903.3;

        my ($fx, $fy, $fz);
        my ($xr, $yr, $zr) = ( $X /  $Xw, 
                                                   $Y /  $Yw, 
                                                   $Z /  $Zw );

        if ($xr > $epsilon) { $fx = pow($xr, 1/3); } else { $fx = ($kappa*$xr + 16)/116; }
        if ($yr > $epsilon) { $fy = pow($yr, 1/3); } else { $fy = ($kappa*$yr + 16)/116; }
        if ($zr > $epsilon) { $fz = pow($zr, 1/3); } else { $fz = ($kappa*$zr + 16)/116; }

        $L = 116 * $fy - 16;
        $a = 500 * ($fx - $fy);
        $b = 200 * ($fy - $fz);

        return [ $L, $a, $b ];
}

sub Lab_to_XYZ { my ($lab, CW$xyz_white) = CW@_; my ($L, CW$a, CW$b) = @{$lab}; my ($Xw, CW$Yw, CW$Zw) = @{$xyz_white}; my ($X, CW$Y, CW$Z);

        my $epsilon =  0.008856;
        my $kappa = 903.3;

        my ($fx, $fy, $fz);
        my ($xr, $yr, $zr);

        if ($L > $kappa*$epsilon) { $yr = pow( ($L + 16)/116, 3 ); } else { $yr = $L / $kappa; }
        if ( $yr > $epsilon ) { $fy =  ($L + 16)/116; } else { $fy  =  ($kappa*$yr + 16)/116; }

        $fx = ($a / 500) + $fy;
        $fz = $fy - ($b / 200);

        if (pow($fx, 3) > $epsilon) { $xr = pow($fx, 3); } else { $xr = (116 * $fx - 16)/$kappa; }
        if (pow($fz, 3) > $epsilon) { $zr = pow($fz, 3); } else { $zr = (116 * $fz - 16)/$kappa; }
        if ($L > $kappa*$epsilon) { $yr = pow(($L + 16)/116, 3);  } else { $yr = $L/$kappa; }

        $X = $xr * $Xw;
        $Y = $yr * $Yw;
        $Z = $zr * $Zw;

        return [ $X, $Y, $Z ];
}

sub RGB_to_Lab { my ($rgb, CW$space) = CW@_; my CW$xyz_white = &RGB_to_XYZ([ 1.0, 1.0, 1.0 ], CW$space); my CW$xyz = &RGB_to_XYZ($rgb, CW$space);

        return &XYZ_to_Lab($xyz, $xyz_white);
}

sub Lab_to_RGB { my ($lab, CW$space) = CW@_; my CW$xyz_white = &RGB_to_XYZ([ 1.0, 1.0, 1.0 ], CW$space); my CW$xyz = &Lab_to_XYZ($lab, CW$xyz_white);

        return &XYZ_to_RGB($xyz, $space);
}

sub XYZ_to_Luv { my ($xyz, CW$xyz_white) = CW@_; my ($X, CW$Y, CW$Z) = @{$xyz}; my ($Xw, CW$Yw, CW$Zw) = @{$xyz_white}; my ($L, CW$u, CW$v);

        my $epsilon =  0.008856;
        my $kappa = 903.3;

        my ($yr) = ( $Y /  $Yw );

        if ($yr > $epsilon) { $L = 116 * pow($yr, 1/3) - 16; }
        else { $L = $kappa*$yr; }

        my ($up, $vp);
        my ($upw, $vpw);

        ($upw, $vpw) = ( 4 * $Xw / ( $Xw + 15 * $Yw + 3 * $Zw ),
                                                9 * $Yw / ( $Xw + 15 * $Yw + 3 * $Zw ) );

        if (! ($X == 0 && $Y == 0 && $Z == 0))
        {
                ($up, $vp) = ( 4 * $X / ( $X + 15 * $Y + 3 * $Z ),
                                           9 * $Y / ( $X + 15 * $Y + 3 * $Z ) );
        }
        else
        {
                ($up, $vp) = ($upw, $vpw);
        }

        ($u, $v) = ( 13 * $L * ($up - $upw),
                                 13 * $L * ($vp - $vpw) );

        return [ $L, $u, $v ];
}

sub Luv_to_XYZ { my ($luv, CW$xyz_white) = CW@_; my ($L, CW$u, CW$v) = @{$luv}; my ($Xw, CW$Yw, CW$Zw) = @{$xyz_white}; my ($X, CW$Y, CW$Z);

        my $epsilon =  0.008856;
        my $kappa = 903.3;

        if ($L > $kappa*$epsilon) { $Y = pow( ($L + 16)/116, 3 ); } else { $Y = $L / $kappa; }

        my ($upw, $vpw) = ( 4 * $Xw / ( $Xw + 15 * $Yw + 3 * $Zw ),
                                                9 * $Yw / ( $Xw + 15 * $Yw + 3 * $Zw ) );

        if (! ($L == 0 && $u == 0 && $v == 0))
        {
                my $a = (1/3)*( ((52 * $L) / ($u + 13 * $L * $upw)) - 1 );
                my $b = -5 * $Y;
                my $c = -1/3;
                my $d = $Y * ( ((39 * $L) / ($v + 13 * $L * $vpw)) - 5 );

                $X = ($d - $b)/($a - $c);
                $Z = $X * $a + $b;
        }
        else
        {
                ($X, $Z) = (0.0, 0.0);
        }

        return [ $X, $Y, $Z ];
}

sub Luv_to_LCHuv { my ($luv) = CW@_; my ($L, CW$u, CW$v) = @{$luv}; my ($C, CW$H);

        $C = sqrt( $u*$u + $v*$v );
        $H = atan2( $v, $u );
        $H = rad2deg($H);

        return [ $L, $C, $H ];
}

sub LCHuv_to_Luv { my ($lch) = CW@_; my ($L, CW$C, CW$H) = @{$lch}; my ($u, CW$v);

        $H = deg2rad($H);
        my $th = tan($H);
        $u = $C / sqrt( $th * $th + 1 );
        $v = sqrt($C*$C - $u*$u);

        #$H = $H - 2*pi*int($H / 2*pi); # convert H to 0..2*pi - this seems to be wrong
        if ($H < 0) { $H = $H + 2*pi; }
        if ($H > pi/2 && $H < 3*pi/2) { $u = - $u; }
        if ($H > pi) { $v = - $v; }

        return [ $L, $u, $v ];
}

sub XYZ_to_xyY { my ($xyz, CW$xyz_white) = CW@_; my ($X, CW$Y, CW$Z) = @{$xyz}; my ($Xw, CW$Yw, CW$Zw) = @{$xyz_white}; my ($x, CW$y);

        if (! ($X == 0 && $Y == 0 && $Z == 0))
        {
                $x = $X / ($X + $Y + $Z);
                $y = $Y / ($X + $Y + $Z);
        }
        else
        {
                $x = $Xw / ( $Xw + $Yw + $Zw );
                $y = $Yw / ( $Xw + $Yw + $Zw );
        }

        return [ $x, $y, $Y ];
}

sub xyY_to_XYZ { my ($xyy) = CW@_; my ($x, CW$y, CW$Y) = @{$xyy}; my ($X, CW$Z);

        if (! ($y == 0))
        {
                $X = $x * $Y / $y;
                $Z = (1 - $x - $y) * $Y / $y;
        }
        else
        {
                $X = 0; $Y = 0; $Z = 0;
        }

        return [ $X, $Y, $Z ];
}

sub Lab_to_LCHab { my ($lab) = CW@_; my ($L, CW$a, CW$b) = @{$lab}; my ($C, CW$H);

        $C = sqrt( $a*$a + $b*$b );
        $H = atan2( $b, $a );
        $H = rad2deg($H);

        return [ $L, $C, $H ];
}

sub LCHab_to_Lab { my ($lch) = CW@_; my ($L, CW$C, CW$H) = @{$lch}; my ($a, CW$b);

        $H = deg2rad($H);
        my $th = tan($H);
        $a = $C / sqrt( $th * $th + 1 );
        $b = sqrt($C*$C - $a*$a);

        #$H = $H - 2*pi*int($H / 2*pi); # convert H to 0..2*pi - this seems to be wrong
        if ($H < 0) { $H = $H + 2*pi; }
        if ($H > pi/2 && $H < 3*pi/2) { $a = - $a; }
        if ($H > pi) { $b = - $b; }

        return [ $L, $a, $b ];
}

sub RGB_to_linear_RGB { my ($rgb, CW$space) = CW@_; my ($R, CW$G, CW$B) = @{$rgb};

        my $s = &_get_RGB_space_by_name($space);
        if ($s->{gamma} eq 'sRGB') # handle special sRGB gamma curve
        {
                if ( abs($R) <= 0.04045 ) { $R = $R / 12.92; }
                else { $R = &_apow( ( $R + 0.055 ) / 1.055 , 2.4 ); }

                if ( abs($G) <= 0.04045 ) { $G = $G / 12.92; }
                else { $G = &_apow( ( $G + 0.055 ) / 1.055 , 2.4 ); }

                if ( abs($B) <= 0.04045 ) { $B = $B / 12.92; }
                else { $B = &_apow( ( $B + 0.055 ) / 1.055 , 2.4 ); }
        }
        else 
        {
                $R = &_apow($R, $s->{gamma});
                $G = &_apow($G, $s->{gamma});
                $B = &_apow($B, $s->{gamma});
        }

        return [ $R, $G, $B ];
}

sub linear_RGB_to_RGB { my ($rgb, CW$space) = CW@_; my ($R, CW$G, CW$B) = @{$rgb};

        my $s = &_get_RGB_space_by_name($space);
        if ($s->{gamma} eq 'sRGB') # handle special sRGB gamma curve
        {
                if ( abs($R) <= 0.0031308 ) { $R = 12.92 * $R; }
                else { $R = 1.055 * &_apow($R, 1/2.4) - 0.055; };

                if ( abs($G) <= 0.0031308 ) { $G = 12.92 * $G; }
                else { $G = 1.055 * &_apow($G, 1/2.4) - 0.055; }

                if ( abs($B) <= 0.0031308 ) { $B = 12.92 * $B; }
                else { $B = 1.055 * &_apow($B, 1/2.4) - 0.055; }
        }
        else 
        {
                $R = &_apow($R, 1/$s->{gamma});
                $G = &_apow($G, 1/$s->{gamma});
                $B = &_apow($B, 1/$s->{gamma});
        }

        return [ $R, $G, $B ];
}

# reference: http://en.wikipedia.org/wiki/YIQ sub RGB_to_YIQ { my ($rgb) = CW@_; # input should be CIE Rec 601/NTSC non-linear rgb my CW$m = [[0.299 , 0.587 , 0.114 ], [0.59590059, -0.27455667, -0.32134392], [0.21153661, -0.52273617, 0.31119955]];

        my $yiq = &_mult_m33_v3($m, $rgb);
        return $yiq;
}

sub YIQ_to_RGB { my ($yiq) = CW@_; my CW$mstar = [[ 1.0 , 0.95598634, 0.6208248 ], [ 1.0 , -0.27201283, -0.64720424], [ 1.0 , -1.1067402 , 1.7042305 ]];

        my $rgb = &_mult_m33_v3($mstar, $yiq);
        return $rgb; # result is NTSC non-linear rgb
}

sub RGB_to_YUV { my ($rgb) = CW@_; # input should be CIE Rec 601/NTSC non-linear rgb my CW$m = [[ 0.299 , 0.587 , 0.114 ], [-0.14714119, -0.28886916, 0.43601035 ], [ 0.61497538, -0.51496512, -0.10001026 ]];

        my $yuv = &_mult_m33_v3($m, $rgb);
        return $yuv;
}

sub YUV_to_RGB { my ($yuv) = CW@_; my CW$mstar = [[ 1.0, 0.0 , 1.139883 ], [ 1.0, -0.39464233, -0.58062185], [ 1.0, 2.0320619 , 0.0 ]];

        my $rgb = &_mult_m33_v3($mstar, $yuv);
        return $rgb; # result is NTSC non-linear rgb
}

# reference: http://www.poynton.com/notes/colour_and_gamma/ColorFAQ.txt # Y is [0..1], Pb and Pr are [-0.5..0.5] sub RGB_to_YPbPr { my ($rgb) = CW@_; # input should be CIE Rec 601/NTSC non-linear rgb my CW$m = [[ 0.299 , 0.587 , 0.114 ], [-0.168736,-0.331264, 0.5 ], [ 0.5 ,-0.418688,-0.081312]];

        my $ypbpr = &_mult_m33_v3($m, $rgb);
        return $ypbpr;
}

sub YPbPr_to_RGB { my ($ypbpr) = CW@_; my CW$mstar = [[ 1.0 , 0.0 , 1.402 ], [ 1.0 ,-0.344136,-0.714136], [ 1.0 , 1.772 , 0.0 ]];

        my $rgb = &_mult_m33_v3($mstar, $ypbpr);
        return $rgb; # result is NTSC non-linear rgb
}

# Y is [16..235], Cb and Cr are [16..239] sub RGB_to_YCbCr { my ($rgb) = CW@_; # input should be NTSC non-linear rgb my CW$m = [[ 65.481, 128.553, 24.966], [ -37.797, -74.203, 112.0 ], [ 112.0 , -93.786, -18.214]];

        my $ycbcr = &_add_v3( &_mult_m33_v3($m, $rgb), [ 16, 128, 128 ] );
        return $ycbcr;
}

sub YCbCr_to_RGB { my ($ycbcr) = CW@_; my CW$mstar = [[ 0.00456621, 0.0 , 0.00625893], [ 0.00456621,-0.00153632,-0.00318811], [ 0.00456621, 0.00791071, 0.0 ]];

        my $rgb = &_mult_m33_v3($mstar, &_add_v3($ycbcr, [-16, -128, -128]));
        return $rgb;
}

# reference: http://wwwde.kodak.com/global/en/professional/products/storage/pcd/techInfo/pcd-045.jhtml sub RGB_to_PhotoYCC { my ($rgb) = CW@_; # input should be CIE Rec 709 non-linear rgb my CW$m = [[ 0.299 , 0.587 , 0.114 ], [-0.299 , -0.587 , 0.866 ], [ 0.701 , -0.587 , -0.114 ]]; my CW$ycc = &_add_v3([0, 156, 137], &_mult_m33_v3([[255/1.402, 0, 0], [0, 111.40, 0], [0, 0, 135.64]], &_mult_m33_v3($m, CW$rgb))); return CW$ycc; }

sub PhotoYCC_to_RGB { my ($ycc) = CW@_; my CW$mstar = [[ 1.0 , 0.0 , 1.0 ], [ 0.99603657, -0.19817126, -0.50936968], [ 1.0204082 , 1.0204082 , 0.0 ]];

        my $rgb = &_mult_m33_v3($mstar, 
                                                        &_mult_m33_v3([[1/(255/1.402), 0, 0], [0, 1/111.40, 0], [0, 0, 1/135.64]], 
                                                                                  &_add_v3([0, -156, -137], $ycc)));
        return $rgb; # result is CIE 709 non-linear rgb
}

sub RGB_to_HSV { my ($rgb) = CW@_; my ($r, CW$g, CW$b)=@{$rgb}; my ($h, CW$s, CW$v);

        my $min= &_min($r, $g, $b);
        my $max= &_max($r, $g, $b);

        $v = $max;                              
        my $delta = $max - $min;

        if( $delta != 0 )
        {
                $s = $delta / $max;
        }
        else
        {
                $s = 0;
                $h = 0;
                return [ $h, $s, $v];
        }

        if( $r == $max )
        {
                $h = ( $g - $b ) / $delta; 
        }
        elsif ( $g == $max )
        {
                $h = 2 + ( $b - $r ) / $delta; 
        }
        else # if $b == $max
        {
                $h = 4 + ( $r - $g ) / $delta;
        }

        $h *= 60;
        if( $h < 0 ) { $h += 360; }
        return [ $h, $s, $v ];
}

sub HSV_to_RGB { my ($hsv) = CW@_; my ($h, CW$s, CW$v)=@{$hsv}; my ($r, CW$g, CW$b);

        # force $h to 0 <= $h < 360
        # FIXME should not loop, looks infinite
        while ($h < 0) { $h += 360; }
        while ($h >= 360) { $h -= 360; }

        $h /= 60;                       ## sector 0 to 5
        my $i = POSIX::floor( $h );
        my $f = $h - $i;                   ## fractional part of h
        my $p = $v * ( 1 - $s );
        my $q = $v * ( 1 - $s * $f );
        my $t = $v * ( 1 - $s * ( 1 - $f ) );

        if($i == 0)
        {
                $r = $v;
                $g = $t;
                $b = $p;
        }
        elsif($i == 1)
        {
                $r = $q;
                $g = $v;
                $b = $p;
        }
        elsif($i == 2)
        {
                $r = $p;
                $g = $v;
                $b = $t;
        }
        elsif($i == 3)
        {
                $r = $p;
                $g = $q;
                $b = $v;
        }
        elsif($i == 4)
        {
                $r = $t;
                $g = $p;
                $b = $v;
        }
        else # if $i == 5
        {
                $r = $v;
                $g = $p;
                $b = $q;
        }

        return [ $r, $g, $b ];
}

sub RGB_to_HSL { my ($rgb) = CW@_; my ($r,$g,$b)=@{$rgb};

        my ($h, $s, $v) = @{ &RGB_to_HSV($rgb) };

        my $min= &_min($r, $g, $b);
        my $max= &_max($r, $g, $b);
        my $delta = $max - $min;

        my $l = ($max+$min)/2.0;

        if( $delta == 0 )
        {
                return [0, 0, $l];
        }
        else
        {
                if($l <= 0.5)
                {
                        $s = $delta/($max+$min); # FIXME possible divide-by-zero
                }
                else
                {
                        $s = $delta/(2-$max-$min); # FIXME possible divide-by-zero
                }
        }
        return [$h, $s, $l];
}

sub HSL_to_RGB { my ($hsl) = CW@_; my ($h, CW$s, CW$l) = @{$hsl}; my ($r, CW$g, CW$b); my ($p1, CW$p2);

        if( $l <= 0.5 )
        {
                $p1 = $l * (1-$s);
                $p2 = 2*$l - $p1;
        }
        else
        {
                $p2 = $l + $s - ($l*$s);
                $p1 = 2*$l - $p2;
        }

        $r = &_rgbquant($p1, $p2, $h+120);
        $g = &_rgbquant($p1, $p2, $h);
        $b = &_rgbquant($p1, $p2, $h-120);

        return [ $r, $g, $b ];
}

sub _rgbquant { my ($q1, CW$q2, CW$h) = CW@_;

        # force $h to 0 <= $h < 360
        # FIXME should not loop, looks infinite
        while ($h < 0) { $h += 360; }
        while ($h >= 360) { $h -= 360; }

        if ($h < 60)
        {
                return ($q1 + (($q2-$q1)*$h/60) );
        }
        elsif ($h < 180)
        {
                return ($q2);
        }
        elsif ($h < 240)
        {
                return ($q1 + (($q2-$q1)*(240-$h)/60) );
        }
        else
        {
                return ($q1);
        }
}

sub RGB_to_CMY { my ($rgb) = CW@_; return [ map { 1 - CW$_ } @{$rgb} ]; }

sub CMY_to_RGB { my ($cmy) = CW@_; return [ map { 1 - CW$_ } @{$cmy} ]; }

sub CMY_to_CMYK { my ($cmy) = CW@_; my CW$k = &_min(@{$cmy}); return [ (map { CW$_-$k } @{$cmy}),$k ]; }

sub CMYK_to_CMY { my ($cmyk) = CW@_; my ($c, CW$m, CW$y, CW$k) = @{$cmyk}; return [ CW$c+$k, CW$m+$k, CW$y+$k ]; }

sub XYZ_change_white_point { my ($xyz, CW$xyz_old_white_point, CW$xyz_new_white_point) = CW@_;

        # matrices for Bradford color-adaptation
        my $ma = [[ 0.8951,   -0.7502,    0.0389 ],  
                          [ 0.2664,    1.7135,   -0.0685 ], 
                          [ -0.1614,   0.0367,    1.0296 ]];

        my $ma_star =  [[ 0.986993,  0.432305, -0.008529 ],
                                        [-0.147054,  0.518360,  0.040043 ],
                                        [ 0.159963,  0.049291,  0.968487 ]];

        # cone = cone response domain value (rho, ypsilon, beta)
        my $cone_old = &_mult_v3_m33($xyz_old_white_point, $ma);
        my $cone_new = &_mult_v3_m33($xyz_new_white_point, $ma);

        my $q = [[ $cone_new->[0]/$cone_old->[0], 0, 0 ],
                         [ 0, $cone_new->[1]/$cone_old->[1], 0 ],
                         [ 0, 0, $cone_new->[2]/$cone_old->[2] ]];

        my $m = &_mult_m33_m33($ma, &_mult_m33_m33($q, $ma_star));

        my $xyz_new = &_mult_v3_m33($xyz, $m);

        return $xyz_new;
}

# reference: http://www.brucelindbloom.com/index.html?Eqn_T_to_xy.html sub white_point_from_temperature { my ($temp) = CW@_; my ($x, CW$y);

        if ($temp < 4000 || $temp > 25000)
        {
                carp "color temperature out of range: $temp, should be between 4000 and 25000 Kelvin";
        }

        if ($temp <= 7000)
        {
                $x = -4.6070e+9 / ($temp*$temp*$temp) + 
                        2.9678e+6   / ($temp*$temp) +
                        0.09911e+3  / $temp +
                        0.244063;
        }
        else # $temp > 7000
        {
                $x = -2.0064e+9 / ($temp*$temp*$temp) + 
                        1.9018e+6   / ($temp*$temp) +
                        0.24748e+3  / $temp +
                        0.237040;
        }

        $y = -3.0 * $x * $x + 2.87 * $x - 0.275;

        return [ $x, $y ];
}

######### private utility functions ########

sub _get_RGB_space_by_name { my ($space) = CW@_; # FIXME the logic here is a bit convoluted, this could be cleaned up a lot

        if (! defined $space)
        {
                # carp("no rgb space specified in operation that requires it, defaulting to sRGB");
                $space = 'sRGB';
        }
        elsif (! $RGB_SPACES{ $space })
        {
                carp("rgb space not found: ".$space.", defaulting to sRGB");
                $space = 'sRGB';
        }

        my $s = $RGB_SPACES{$space};
        if ($s && ! ref $s)
        {
                $s = $RGB_SPACES{$s}; # follow aliases
        }

        return $s;
}

sub _check_white_point { my ($white_point) = CW@_;

        if (! defined $white_point)
        {
                # carp("no white point specified in operation that requires it, defaulting to D65");
                $white_point = 'D65';
        }
        elsif ($white_point =~ m!^(\d+)K$!)
        {
                my $temperature = $1;
                #$white_point = $temperature.'K'; # already in that form
                $WHITE_POINTS{ $white_point } = &white_point_from_temperature($temperature);
        }
        elsif (! $WHITE_POINTS{ $white_point })
        {
                carp("white point not found: ". $white_point.", defaulting to D65");
                $white_point = 'D65';
        }

        return $white_point;
}

sub _mult_v3_m33 { my ($v, CW$m) = CW@_; my CW$vout = [ ( CW$v->[0] * CW$m->[0]->[0] + CW$v->[1] * CW$m->[1]->[0] + CW$v->[2] * CW$m->[2]->[0] ), ( CW$v->[0] * CW$m->[0]->[1] + CW$v->[1] * CW$m->[1]->[1] + CW$v->[2] * CW$m->[2]->[1] ), ( CW$v->[0] * CW$m->[0]->[2] + CW$v->[1] * CW$m->[1]->[2] + CW$v->[2] * CW$m->[2]->[2] ) ]; return CW$vout; }

sub _mult_m33_v3 { my ($m, CW$v) = CW@_; my CW$vout = [ ( CW$v->[0] * CW$m->[0]->[0] + CW$v->[1] * CW$m->[0]->[1] + CW$v->[2] * CW$m->[0]->[2] ), ( CW$v->[0] * CW$m->[1]->[0] + CW$v->[1] * CW$m->[1]->[1] + CW$v->[2] * CW$m->[1]->[2] ), ( CW$v->[0] * CW$m->[2]->[0] + CW$v->[1] * CW$m->[2]->[1] + CW$v->[2] * CW$m->[2]->[2] ) ]; return CW$vout; }

sub _mult_m33_m33 { my ($m, CW$n) = CW@_; my CW$q = []; foreach my CW$i (0..2) { foreach my CW$j (0..2) { foreach my CW$k (0..2) { CW$q->[$i]->[$j] += CW$m->[$i]->[$k] * CW$n->[$k]->[$j]; } } } return CW$q; }

sub _add_v3 { my ($a, CW$b) = CW@_; my CW$c = [ CW$a->[0] + CW$b->[0], CW$a->[1] + CW$b->[1], CW$a->[2] + CW$b->[2] ]; return CW$c; }

sub _pow_v3 { my ($v3, CW$c) = CW@_; my CW$v3out = [ pow($v3->[0], CW$c), pow($v3->[1], CW$c), pow($v3->[2], CW$c) ]; return CW$v3out; }

sub _delta_v3 { my ($a3, CW$b3) = CW@_; return ( abs($a3->[0] - CW$b3->[0]) + abs($a3->[1] - CW$b3->[1]) + abs($a3->[2] - CW$b3->[2]) ); }

sub _generic_clip { my ($v3, CW$c32) = CW@_; if ($v3->[0] < CW$c32->[0]->[0]) { CW$v3->[0] = CW$c32->[0]->[0]; } if ($v3->[0] > CW$c32->[0]->[1]) { CW$v3->[0] = CW$c32->[0]->[1]; } if ($v3->[1] < CW$c32->[1]->[0]) { CW$v3->[1] = CW$c32->[1]->[0]; } if ($v3->[1] > CW$c32->[1]->[1]) { CW$v3->[1] = CW$c32->[1]->[1]; } if ($v3->[2] < CW$c32->[2]->[0]) { CW$v3->[2] = CW$c32->[2]->[0]; } if ($v3->[2] > CW$c32->[2]->[1]) { CW$v3->[2] = CW$c32->[2]->[1]; } return CW$v3; }

sub _apow { my ($v, CW$p) = CW@_; return ($v >= 0 ? pow($v, CW$p) : -pow(-$v, CW$p)); }

sub _sqr { my ($v) = CW@_; return CW$v*$v; }

sub _is_zero { my ($v) = CW@_; return (abs($v) < 0.000001); }

sub _min { my CW$min = shift(@_); foreach my CW$v (@_) { if ($v <= CW$min) { CW$min = CW$v; } }; return CW$min; }

sub _max { my CW$max = shift(@_); foreach my CW$v (@_) { if ($v >= CW$max) { CW$max = CW$v; } }; return CW$max; }

######### colorspace tables ########

# reference: http://www.brucelindbloom.com/Eqn_RGB_XYZ_Matrix.html # All the rgb spaces that this module knows about. # Key is the name, value is either another name (i.e. this is an alias), or a hashref containg a white point, gamma, a conversion matrix m for rgb-to-xyz and a reverse matrix mstar for xyz-to-rgb transformations our CW%RGB_SPACES = ( 'Adobe' => 'Adobe RGB (1998)', 'Adobe RGB (1998)' => { white_point => 'D65', gamma => 2.2, m => [ [ 0.5767001212121210, 0.2973609999999999, 0.0270328181818181 ], [ 0.1855557042253521, 0.6273550000000000, 0.0706878873239437 ], [ 0.1882125000000000, 0.0752850000000000, 0.9912525000000000 ] ], mstar => [ [ 2.0414778828777158, -0.9692568708746859, 0.0134454339800522 ], [ -0.5649765261191881, 1.8759931170154693, -0.1183725462165374 ], [ -0.3447127732462102, 0.0415556248231326, 1.0152620834741313 ] ], }, 'Apple' => 'Apple RGB', 'Apple RGB' => { white_point => 'D65', gamma => 1.8, m => [ [ 0.4496948529411764, 0.2446340000000000, 0.0251829117647059 ], [ 0.3162512941176471, 0.6720340000000000, 0.1411836134453782 ], [ 0.1845208571428572, 0.0833320000000000, 0.9226042857142855 ] ], mstar => [ [ 2.9517603398020569, -1.0851001264872848, 0.0854802409232915 ], [ -1.2895090072470441, 1.9908397072633022, -0.2694550155056003 ], [ -0.4738802866606785, 0.0372022452865781, 1.0911301341384845 ] ], }, 'BestRGB' => { white_point => 'D50', gamma => 2.2, m => [ [ 0.6326700260082926, 0.2284570000000000, 0.0000000000000000 ], [ 0.2045557161290322, 0.7373519999999999, 0.0095142193548387 ], [ 0.1269951428571429, 0.0341910000000000, 0.8156995714285713 ] ], mstar => [ [ 1.7552588897490133, -0.5441338472581142, 0.0063467101890703 ], [ -0.4836782739368681, 1.5068795234848715, -0.0175760572028268 ], [ -0.2529998994965047, 0.0215528345168675, 1.2256901641540674 ] ], }, 'Beta RGB' => { white_point => 'D50', gamma => 2.2, m => [ [ 0.6712546349614399, 0.3032730000000001, 0.0000000000000001 ], [ 0.1745833659117997, 0.6637859999999999, 0.0407009558998808 ], [ 0.1183817187500000, 0.0329410000000000, 0.7845011448863635 ] ], mstar => [ [ 1.6832246105012654, -0.7710229999344457, 0.0400016919321019 ], [ -0.4282356869228009, 1.7065573340451357, -0.0885384492378917 ], [ -0.2360181522709381, 0.0446899574535591, 1.2723768250932299 ] ], }, 'BruceRGB' => { white_point => 'D65', gamma => 2.2, m => [ [ 0.4673842424242424, 0.2409950000000000, 0.0219086363636363 ], [ 0.2944540307692308, 0.6835539999999999, 0.0736135076923076 ], [ 0.1886300000000000, 0.0754520000000000, 0.9934513333333335 ] ], mstar => [ [ 2.7456543761403882, -0.9692568108426551, 0.0112706581772173 ], [ -1.1358911781912031, 1.8759930008236942, -0.1139588771251973 ], [ -0.4350565642146659, 0.0415556222493375, 1.0131069405965349 ] ], }, 'CIE' => { white_point => 'E', gamma => 2.2, m => [ [ 0.4887167547169811, 0.1762040000000000, 0.0000000000000000 ], [ 0.3106804602510461, 0.8129850000000002, 0.0102048326359833 ], [ 0.2006041111111111, 0.0108110000000000, 0.9898071111111111 ] ], mstar => [ [ 2.3706802022946527, -0.5138847730830187, 0.0052981111618865 ], [ -0.9000427625776859, 1.4253030498717687, -0.0146947611471193 ], [ -0.4706349622815629, 0.0885813466699250, 1.0093845871252884 ] ], }, 'ColorMatch' => { white_point => 'D50', gamma => 1.8, m => [ [ 0.5093438823529410, 0.2748840000000000, 0.0242544705882353 ], [ 0.3209073388429752, 0.6581320000000002, 0.1087821487603307 ], [ 0.1339700000000000, 0.0669850000000000, 0.6921783333333333 ] ], mstar => [ [ 2.6422872594587332, -1.1119754096457255, 0.0821692807629542 ], [ -1.2234269646206919, 2.0590166676215107, -0.2807234418494614 ], [ -0.3930142794480749, 0.0159613695164458, 1.4559774449385248 ] ], }, 'DonRGB4' => { white_point => 'D50', gamma => 2.2, m => [ [ 0.6457719999999998, 0.2783499999999999, 0.0037113333333334 ], [ 0.1933510457516340, 0.6879700000000001, 0.0179861437908497 ], [ 0.1250971428571429, 0.0336800000000000, 0.8035085714285716 ] ], mstar => [ [ 1.7603878846606116, -0.7126289975811030, 0.0078207770365325 ], [ -0.4881191497764036, 1.6527436537605511, -0.0347412748629646 ], [ -0.2536122811541382, 0.0416715470705678, 1.2447804103656714 ] ], }, 'ECI' => { white_point => 'D50', gamma => 1.8, m => [ [ 0.6502045454545454, 0.3202500000000000, -0.0000000000000001 ], [ 0.1780773380281691, 0.6020710000000000, 0.0678389859154930 ], [ 0.1359382500000000, 0.0776790000000000, 0.7573702500000002 ] ], mstar => [ [ 1.7827609790470664, -0.9593624312689213, 0.0859317810050046 ], [ -0.4969845184555761, 1.9477964513641737, -0.1744675553737970 ], [ -0.2690099687053119, -0.0275807381172883, 1.3228286288043098 ] ], }, 'Ekta Space PS5' => { white_point => 'D50', gamma => 2.2, m => [ [ 0.5938923114754098, 0.2606289999999999, 0.0000000000000000 ], [ 0.2729799428571429, 0.7349460000000001, 0.0419969142857143 ], [ 0.0973500000000000, 0.0044250000000000, 0.7832250000000001 ] ], mstar => [ [ 2.0043787360968186, -0.7110290170493107, 0.0381257297502959 ], [ -0.7304832564783660, 1.6202136618008882, -0.0868766628736253 ], [ -0.2450047962579189, 0.0792227384931296, 1.2725243569115190 ] ], }, '601' => 'NTSC', 'CIE Rec 601' => 'NTSC', 'NTSC' => { white_point => 'C', gamma => 2.2, m => [ [ 0.6067337272727271, 0.2988389999999999, -0.0000000000000001 ], [ 0.1735638169014085, 0.5868110000000000, 0.0661195492957747 ], [ 0.2001125000000000, 0.1143500000000000, 1.1149125000000002 ] ], mstar => [ [ 1.9104909450902432, -0.9843106185066585, 0.0583742441336926 ], [ -0.5325921048972800, 1.9984488315135187, -0.1185174047562849 ], [ -0.2882837998985277, -0.0282979742694222, 0.8986095763610844 ] ], }, 'CIE ITU' => 'PAL/SECAM', 'PAL' => 'PAL/SECAM', 'PAL/SECAM' => { white_point => 'D65', gamma => 2.2, m => [ [ 0.4305861818181819, 0.2220210000000001, 0.0201837272727273 ], [ 0.3415450833333333, 0.7066450000000000, 0.1295515833333333 ], [ 0.1783350000000000, 0.0713340000000000, 0.9392309999999999 ] ], mstar => [ [ 3.0631308078036081, -0.9692570313532748, 0.0678676345258901 ], [ -1.3932854294802033, 1.8759934276211896, -0.2288214781555966 ], [ -0.4757879688629482, 0.0415556317034429, 1.0691933898259074 ] ], }, 'ProPhoto' => { white_point => 'D50', gamma => 1.8, m => [ [ 0.7976742857142858, 0.2880400000000000, 0.0000000000000000 ], [ 0.1351916830080914, 0.7118740000000000, 0.0000000000000000 ], [ 0.0314760000000000, 0.0000860000000000, 0.8284380000000000 ] ], mstar => [ [ 1.3459444124134017, -0.5445989438461810, -0.0000000000000000 ], [ -0.2556077203964527, 1.5081675237232912, -0.0000000000000000 ], [ -0.0511118080787822, 0.0205351443915685, 1.2070909349884964 ] ], }, 'SMPTE' => 'SMPTE-C', 'SMPTE-C' => { white_point => 'D65', gamma => 2.2, m => [ [ 0.3935554411764707, 0.2123950000000001, 0.0187407352941176 ], [ 0.3652524201680672, 0.7010489999999999, 0.1119321932773109 ], [ 0.1916597142857142, 0.0865560000000000, 0.9582985714285710 ] ], mstar => [ [ 3.5056956039694129, -1.0690641158576772, 0.0563116543373650 ], [ -1.7396380462846184, 1.9778095119692913, -0.1969933651732733 ], [ -0.5440105230649496, 0.0351719640259221, 1.0500467308790999 ] ], }, '709' => 'sRGB', 'CIE Rec 709' => 'sRGB', 'sRGB' => { white_point => 'D65', gamma => 'sRGB', # 2.4, m => [ [ 0.4124237575757575, 0.2126560000000000, 0.0193323636363636 ], [ 0.3575789999999999, 0.7151579999999998, 0.1191930000000000 ], [ 0.1804650000000000, 0.0721860000000000, 0.9504490000000001 ] ], mstar => [ [ 3.2407109439941704, -0.9692581090654827, 0.0556349466243886 ], [ -1.5372603195869781, 1.8759955135292130, -0.2039948042894247 ], [ -0.4985709144606416, 0.0415556779089489, 1.0570639858633826 ] ], }, 'WideGamut' => { white_point => 'D50', gamma => 2.2, m => [ [ 0.7161035660377360, 0.2581870000000001, 0.0000000000000000 ], [ 0.1009296246973366, 0.7249380000000000, 0.0517812857142858 ], [ 0.1471875000000000, 0.0168750000000000, 0.7734375000000001 ] ], mstar => [ [ 1.4628087611158722, -0.5217931929785991, 0.0349338148323482 ], [ -0.1840625990709008, 1.4472377239217746, -0.0968919015161355 ], [ -0.2743610287417160, 0.0677227300206644, 1.2883952872306403 ] ], } );

# reference: http://www.aim-dtp.net/aim/technology/cie_xyz/cie_xyz.htm # reference: Wyszecki, G. and Stiles, W. S. Color Science Concepts and Methods, Wiley (2000). ISBN 0471399183 # based on CIE1931 (2 degree FOV) our CW%WHITE_POINTS = ( 'A' => [ 0.44757, 0.40745 ], # Tungsten lamp 2856K 'D50' => [ 0.34567, 0.35850 ], # Bright tungsten 'B' => [ 0.34842, 0.35161 ], # CIE Std illuminant B 'D55' => [ 0.33242, 0.34743 ], # Cloudy daylight 'E' => [ 0.333333, 0.333333 ], # Normalized reference source 'D65' => [ 0.312713, 0.329016 ], # Daylight 6504K 'C' => [ 0.310063, 0.316158 ], # North daylight 6774K 'D75' => [ 0.29902, 0.31485 ], # 7500K 'D93' => [ 0.28480, 0.29320 ], # old CRT monitors 'F2' => [ 0.37207, 0.37512 ], # Cool white fluorescent 4200K 'F7' => [ 0.31285, 0.32918 ], # Narrow daylight fluorescent 6500K 'F11' => [ 0.38054, 0.37691 ], # Narrow white fluorescent );

EXPORT

None by default. The 'all' tag causes the non-object-oriented interface to be exported, and you get all the XXX_to_YYY functions, for example RGB_to_XYZ. Please note that some of these functions need extra arguments in addition to the color value to be converted.

BUGS

Backwards compatibility with versions before 0.4 is not very well tested.

This module will produce results that are, in some cases, different from other software. Most of the time that is not a bug in this module, but rather a case where the other software uses an approximate (trading accuracy for speed) algorithm. That is particularly true for YUV and related conversions which are often implemented using integer-math approximations. As far as possible, this module produces results which are exact according to the definitions in the relevant CIE/ITU or other standards.

Some color transformations are not exactly reversible. In particular, conversions between different white points are almost but not exactly reversible. This is not a bug.

There is no way to choose any other color-adaptation algorithm than the Bradford algorithm. That is probably ok since the Bradford algorithm is better than other algorithms (such as Von Kries or simple scaling).

There is no way to choose a RGB space other than the built-in ones.

Support for CMYK is very basic, it relies on assumptions that completely do not work in the physical world of subtractive pigment mixtures. If you tried to convert an image to CMYK directly for printing using these functions, the results will not be very good, to say the least.

TODO

Add clipping to gamut for every color space.

Choose between several clipping algorithms (nearest, luminance-preserving, hue-preserving).

Add a simpler way to check whether something is within gamut.

Add user-defined RGB spaces.

Calculate RGB matrices from chromaticity coordinates.

Only load non-RGB matrices once at startup.

Add colorspaces: uvw, YOZ, RYB, others?

Add RGB spaces: ROMM, others?

Convert arrays of colors efficiently (maybe someday in C).

SEE ALSO

The Color FAQ by Charles Poynton is one of the definitive references on the subject: http://www.poynton.com/notes/colour_and_gamma/ColorFAQ.txt

Bruce Lindbloom's web site contains a tremendous amount of information on color: http://www.brucelindbloom.com/index.html?Math.html

AUTHOR

Alex Izvorski, <izv@dslextreme.com>

Alfred Reibenschuh <alfredreibenschuh@yahoo.com> was the original author for versions up to 0.3a2.

Many thanks to:

Alfred Reibenschuh <alfredreibenschuh@yahoo.com> for the previous versions of Graphics::ColorObject, and for the HSL/HSV/CMYK code.

Bruce Lindbloom <info@brucelindbloom.com> for providing a wealth of information on color space conversion and color adaptation algorithms, and for the precalculated RGB conversion matrices.

Charles Poynton <colorfaq@poynton.com> for the Color FAQ.

Timo Autiokari <timo.autiokari@aim-dtp.net> for information on white points.

COPYRIGHT AND LICENSE

Copyright 2003-2005 by Alex Izvorski

Portions Copyright 2001-2003 by Alfred Reibenschuh

This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself.