File Coverage

color_space.pd

Criterion	Covered	Total	%
statement	45	45	100.0
branch	1	2	50.0
condition			n/a
subroutine	13	13	100.0
pod	0	8	0.0
total	59	68	86.7

pp_add_exported(
  'rgb_to_xyz', 'xyz_to_rgb',
  'xyz_to_lab', 'lab_to_xyz',
  'rgb_to_lch', 'lch_to_rgb',
  'rgb_to_lab', 'lab_to_rgb',
  'add_rgb_space',
);

$PDL::Graphics::ColorSpace::VERSION = '0.206';

pp_setversion("'$PDL::Graphics::ColorSpace::VERSION'");

pp_addpm({At=>'Top'}, <<'EOD');

=encoding utf8

=head1 NAME

PDL::Graphics::ColorSpace - colour-space conversions for PDL

=head1 SYNOPSIS

    use PDL::LiteF;
    use PDL::IO::Pic;
    use PDL::Graphics::ColorSpace;

    my $image_rgb = PDL->rpic('photo.jpg') if PDL->rpiccan('JPEG');

    # convert RGB value from [0,255] to [0,1]
    $image_rgb = $image_rgb->double / 255;

    my $image_xyz = $image_rgb->rgb_to_xyz( 'sRGB' );

Or

    my $image_xyz = rgb_to_xyz( $image_rgb, 'sRGB' );

=head1 DESCRIPTION

Does image color space conversions such as RGB to XYZ and Lab to LCH. Derived from Graphics::ColorObject (Izvorski & Reibenschuh, 2005) but since it's implemented in C and PDL, it runs *much* faster.

Often the conversion will return out-of-gamut values. Retaining out-of-gamut values allows chained conversions to be lossless and reverse conversions to produce the original values. You can clip the values to be within-gamut if necessary. Please check the specific color space for the gamut range.

=head1 COLOR SPACES

=head2 RGB

An RGB color space is any additive color space based on the RGB color model. A particular RGB color space is defined by the three chromaticities of the red, green, and blue additive primaries, and can produce any chromaticity that is the triangle defined by those primary colors. The complete specification of an RGB color space also requires a white point chromaticity and a gamma correction curve.

For more info on the RGB color space, see L.

This module expects and produces RGB values normalized to be in the range of [0,1]. If you have / need integer value between [0,255], divide or multiply the values by 255.

=head2 CMYK

CMYK refers to the four inks used in some color printing: cyan, magenta, yellow, and key (black). The CMYK model works by partially or entirely masking colors on a lighter, usually white, background. The ink reduces the light that would otherwise be reflected. Such a model is called subtractive because inks "subtract" brightness from white.

In additive color models such as RGB, white is the "additive" combination of all primary colored lights, while black is the absence of light. In the CMYK model, it is the opposite: white is the natural color of the paper or other background, while black results from a full combination of colored inks. To save money on ink, and to produce deeper black tones, unsaturated and dark colors are produced by using black ink instead of the combination of cyan, magenta and yellow.

For more info, see L.

=head2 HSL

Hue, Saturation and Luminance (or brightness).

The HSL color space defines colors more naturally: Hue specifies the base color, the other two values then let you specify the saturation of that color and how bright the color should be.

Hue is specified here as degrees ranging from 0 to 360. There are 6 base colors:

    0    red
    60    yellow
    120    green
    180    cyan
    240    blue
    300    magenta
    360    red

Saturation specifies the distance from the middle of the color wheel. So a saturation value of 0 (0%) means "center of the wheel", i.e. a grey value, whereas a saturation value of 1 (100%) means "at the border of the wheel", where the color is fully saturated.

Luminance describes how "bright" the color is. 0 (0%) means 0 brightness and the color is black. 1 (100%) means maximum brightness and the color is white.

For more info, see L.

=head2 XYZ and xyY

The CIE XYZ color space was derived the CIE RGB color space. XYZ are three hypothetical primaries. Y means brightness, Z is quasi-equal to blue stimulation, and X is a mix which looks like red sensitivity curve of cones.  All visible colors can be represented by using only positive values of X, Y, and Z. The main advantage of the CIE XYZ space (and any color space based on it) is that this space is completely device-independent.

For more info, see L.

=head2 Lab

A Lab color space is a color-opponent space with dimension L for lightness and a and b for the color-opponent dimensions, based on nonlinearly compressed CIE XYZ color space coordinates. It's derived from the "master" space CIE 1931 XYZ color space but is more perceptually uniform than XYZ. The Lab space is relative to the white point of the XYZ data they were converted from. Lab values do not define absolute colors unless the white point is also specified.

For more info, see L.

=head2 LCH

This is possibly a little easier to comprehend than the Lab colour space, with which it shares several features. It is more correctly known as  L*C*H*.  Essentially it is in the form of a sphere. There are three axes; L* and C* and H°.

The L* axis represents Lightness. This is vertical; from 0, which has no lightness (i.e. absolute black), at the bottom; through 50 in the middle, to 100 which is maximum lightness (i.e. absolute white) at the top.

The C* axis represents Chroma or "saturation". This ranges from 0 at the centre of the circle, which is completely unsaturated (i.e. a neutral grey, black or white) to 100 or more at the edge of the circle for very high Chroma (saturation) or "colour purity".

If we take a horizontal slice through the centre, we see a coloured circle. Around the edge of the circle we see every possible saturated colour, or Hue. This circular axis is known as H° for Hue. The units are in the form of degrees° (or angles), ranging from 0 (red) through 90 (yellow), 180 (green), 270 (blue) and back to 0.

For more info, see L.

=head1 OPTIONS

Some conversions require specifying the RGB space which includes gamma curve and white point definitions. Supported RGB spaces include (aliases in square brackets):

    Adobe RGB (1998) [Adobe]
    Apple RGB [Apple]
    BestRGB
    Beta RGB
    BruceRGB
    CIE
    ColorMatch
    DonRGB4
    ECI
    Ekta Space PS5
    NTSC [601] [CIE Rec 601]
    PAL/SECAM [PAL] [CIE ITU]
    ProPhoto
    SMPTE-C [SMPTE]
    WideGamut
    sRGB [709] [CIE Rec 709]
    lsRGB

You can also add custom RGB space definitions via the function add_rgb_space.
Alternatively, as of 0.202 you can directly supply the function with a
hash-ref in the same format.

=head1 CONVERSIONS

Some conversions, if not already included as functions, can be achieved
by chaining existing functions. For example, LCH to HSV conversion can
be achieved by chaining lch_to_rgb and rgb_to_hsv:

    my $hsv = rgb_to_hsv( lch_to_rgb( $lch, 'sRGB' ), 'sRGB' );

To generate a local diagram of what conversions are available between
formats, using GraphViz, you can use this script:

  use blib;
  use PDL::Graphics::ColorSpace;
  print "digraph {\n";
  print join(' -> ', split /_to_/), "\n"
    for grep !/^_/ && /_to_/, @PDL::Graphics::ColorSpace::EXPORT_OK;
  print "}\n";
  # then:
  perl scriptname >d.dot; dot -Tsvg d.dot >d.svg; display d.svg

(As of 0.202, this is everything to and from C, plus C <->
C <-> C)

=cut

use strict;
use warnings;

use Carp;
use PDL::LiteF;
use PDL::Graphics::ColorSpace::RGBSpace;

my $RGB_SPACE = $PDL::Graphics::ColorSpace::RGBSpace::RGB_SPACE;

EOD

pp_addhdr('
#include
#include "color_space.h"  /* Local decs */
'
);

pp_def('rgb_to_cmyk',
    Pars => 'rgb(c=3); [o]cmyk(d=4)',
    HandleBad => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (d) %{
                $SETBAD(cmyk());
            %}
            /* skip to the next cmyk triple */
        }
        else,) {
            rgb2cmyk($P(rgb), $P(cmyk));
        }
    ',

    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an RGB color triple to an CMYK color quadruple.

The first dimension of the ndarrays holding the rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...). The first dimension of the ndarrays holding the cmyk values must be size 4.

=for usage

Usage:

    my $cmyk = rgb_to_cmyk( $rgb );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated C, M, Y, and K values will all be marked as bad.

=cut

BADDOC
);


pp_def('cmyk_to_rgb',
    Pars => 'cmyk(d=4); [o]rgb(c=3)',
    HandleBad => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(cmyk(d=>0)) || $ISBAD(cmyk(d=>1)) || $ISBAD(cmyk(d=>2)) || $ISBAD(cmyk(d=>3))) {
            loop (c) %{
                $SETBAD(rgb());
            %}
            /* skip to the next cmyk triple */
        }
        else,) {
            cmyk2rgb($P(cmyk), $P(rgb));
        }
    ',
    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an CMYK color quadruple to an RGB color triple

The first dimension of the ndarrays holding the cmyk values must be size 4, i.e. the dimensions must look like (4, m, n, ...). The first dimension of the ndarray holding the rgb values must be 3.

=for usage

Usage:

    my $rgb = cmyk_to_rgb( $cmyk );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the C, M, Y, or K quantities, the output ndarray will be marked as bad and the associated R, G, and B color values will all be marked as bad.

=cut

BADDOC
);


pp_def('rgb_to_hsl',
    Pars => 'rgb(c=3); [o]hsl(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (c) %{
                $SETBAD(hsl());
            %}
            /* skip to the next hsl triple */
        }
        else,) {
            rgb2hsl($P(rgb), $P(hsl));
        }
    ',

    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an RGB color triple to an HSL color triple.

The first dimension of the ndarrays holding the hsl and rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $hsl = rgb_to_hsl( $rgb );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated H, S, and L values will all be marked as bad.

=cut

BADDOC
);


pp_def('hsl_to_rgb',
    Pars => 'hsl(c=3); [o]rgb(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(hsl(c=>0)) || $ISBAD(hsl(c=>1)) || $ISBAD(hsl(c=>2))) {
            loop (c) %{
                $SETBAD(rgb());
            %}
            /* skip to the next hsl triple */
        }
        else,) {
            hsl2rgb($P(hsl), $P(rgb));
        }
    ',
    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an HSL color triple to an RGB color triple

The first dimension of the ndarrays holding the hsl and rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $rgb = hsl_to_rgb( $hsl );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the H, S, or V quantities, the output ndarray will be marked as bad and the associated R, G, and B color values will all be marked as bad.

=cut

BADDOC
);


pp_def('rgb_to_hsv',
    Pars => 'rgb(c=3); [o]hsv(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (c) %{
                $SETBAD(hsv());
            %}
            /* skip to the next hsv triple */
        }
        else,) {
            rgb2hsv($P(rgb), $P(hsv));
        }
    ',

    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an RGB color triple to an HSV color triple.

The first dimension of the ndarrays holding the hsv and rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $hsv = rgb_to_hsv( $rgb );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated H, S, and V values will all be marked as bad.

=cut

BADDOC
);


pp_def('hsv_to_rgb',
    Pars => 'hsv(c=3); [o]rgb(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(hsv(c=>0)) || $ISBAD(hsv(c=>1)) || $ISBAD(hsv(c=>2))) {
            loop (c) %{
                $SETBAD(rgb());
            %}
            /* skip to the next rgb triple */
        }
        else,) {
            hsv2rgb($P(hsv), $P(rgb));
        }
    ',

    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an HSV color triple to an RGB color triple

The first dimension of the ndarrays holding the hsv and rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $rgb = hsv_to_rgb( $hsv );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the H, S, or V quantities, the output ndarray will be marked as bad and the associated R, G, and B color values will all be marked as bad.

=cut

BADDOC
);


pp_def('xyY_to_xyz',
    Pars => 'xyY(c=3); [o]xyz(c=3)',
    Doc => 'Internal function for white point calculation. Use it if you must.',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(xyY(c=>0)) || $ISBAD(xyY(c=>1)) || $ISBAD(xyY(c=>2))) {
            loop (c) %{
                $SETBAD(xyz());
            %}
            /* skip to the next hsl triple */
        }
        else,) {
            xyY2xyz($P(xyY), $P(xyz));
        }
    ',
);


pp_def('_rgb_to_xyz',
    Pars => 'rgb(c=3); gamma(); m(i=3,rows=3); [o]xyz(c=3)',
    HandleBad => 1,
    Inplace => ['rgb'],
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (c) %{
                $SETBAD(xyz());
            %}
            /* skip to the next xyz triple */
        }
        else,) {
            rgb2xyz($P(rgb), $gamma(), &$m(i=>0,rows=>0), &$m(i=>0,rows=>1), &$m(i=>0,rows=>2), $P(xyz));
        }
    ',
    Doc => undef,
);


pp_def('_xyz_to_rgb',
    Pars => 'xyz(c=3); gamma(); mstar(i=3,rows=3); [o]rgb(c=3)',
    HandleBad => 1,
    Inplace => ['xyz'],
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(xyz(c=>0)) || $ISBAD(xyz(c=>1)) || $ISBAD(xyz(c=>2))) {
            loop (c) %{
                $SETBAD(rgb());
            %}
            /* skip to the next rgb triple */
        }
        else,) {
            xyz2rgb($P(xyz), $gamma(), &$mstar(i=>0,rows=>0), &$mstar(i=>0,rows=>1), &$mstar(i=>0,rows=>2), $P(rgb));
        }
    ',
    Doc => undef,
);


pp_def('_xyz_to_lab',
    Pars => 'xyz(c=3); w(d=2);  [o]lab(c=3)',
    Doc => undef,
    HandleBad => 1,
    Inplace => ['xyz'],
    GenericTypes=>['D'],
    Code => '
        /* construct white point */
        double xyY[3] = { $w(d=>0), $w(d=>1), 1.0 };
        double xyz_white[3];
        xyY2xyz( &xyY[0], &xyz_white[0] );

        threadloop %{
            /* First check for bad values */
            PDL_IF_BAD(if ($ISBAD(xyz(c=>0)) || $ISBAD(xyz(c=>1)) || $ISBAD(xyz(c=>2))) {
                loop (c) %{
                    $SETBAD(lab());
                %}
                /* skip to the next xyz triple */
            }
            else,) {
                xyz2lab( $P(xyz), &xyz_white[0], $P(lab) );
            }
        %}
    ',
);


pp_def('_lab_to_xyz',
    Pars => 'lab(c=3); w(d=2);  [o]xyz(c=3)',
    Doc => undef,

    HandleBad => 1,
    Inplace => ['lab'],
    GenericTypes=>['D'],
    Code => '
        /* construct white point */
        double xyY[3] = { $w(d=>0), $w(d=>1), 1.0 };
        double xyz_white[3];
        xyY2xyz( &xyY[0], &xyz_white[0] );

        threadloop %{
            /* First check for bad values */
            PDL_IF_BAD(if ($ISBAD(lab(c=>0)) || $ISBAD(lab(c=>1)) || $ISBAD(lab(c=>2))) {
                loop (c) %{
                    $SETBAD(xyz());
                %}
                /* skip to the next lab triple */
            }
            else,) {
                lab2xyz( $P(lab), &xyz_white[0], $P(xyz) );
            }
        %}
    ',
);


pp_def('lab_to_lch',
    Pars => 'lab(c=3); [o]lch(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(lab(c=>0)) || $ISBAD(lab(c=>1)) || $ISBAD(lab(c=>2))) {
            loop (c) %{
                $SETBAD(lch());
            %}
            /* skip to the next lch triple */
        }
        else,) {
            lab2lch( $P(lab), $P(lch) );
        }
    ',
    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an Lab color triple to an LCH color triple.

The first dimension of the ndarrays holding the lab values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $lch = lab_to_lch( $lab );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the L, a, or b values the output ndarray will be marked as bad and the associated L, C, and H values will all be marked as bad.

=cut

BADDOC
);


pp_def('lch_to_lab',
    Pars => 'lch(c=3); [o]lab(c=3)',
    HandleBad => 1,
    Inplace => 1,
    GenericTypes=>['D'],
    Code => '
        /* First check for bad values */
        PDL_IF_BAD(if ($ISBAD(lch(c=>0)) || $ISBAD(lch(c=>1)) || $ISBAD(lch(c=>2))) {
            loop (c) %{
                $SETBAD(lab());
            %}
            /* skip to the next lab triple */
        }
        else,) {
            lch2lab( $P(lch), $P(lab) );
        }
    ',
    Doc => <<'DOCUMENTATION',

=pod

=for ref

Converts an LCH color triple to an Lab color triple.

The first dimension of the ndarrays holding the lch values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for usage

Usage:

    my $lab = lch_to_lab( $lch );

=cut

DOCUMENTATION
    BadDoc => <

=for bad

If C encounters a bad value in any of the L, C, or H values the output ndarray will be marked as bad and the associated L, a, and b values will all be marked as bad.

=cut

BADDOC
);


pp_def('rgb_to_linear',
    Pars => 'rgb(c=3); gamma(); [o]out(c=3)',
    HandleBad => 1,
    Inplace => ['rgb'],
    GenericTypes=>['D'],
    Code => '
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (c) %{ $SETBAD(out()); %}
        }
        else,) {
            rgb2linear( $P(rgb), $gamma(), $P(out) );
        }
    ',
    Doc => <<'DOCUMENTATION',
=for ref

Converts an RGB color triple (presumably with gamma) to an RGB color triple
with linear values.

=for usage

Usage:

    my $rgb_linear = rgb_to_linear( $gammaed, 2.2 );
DOCUMENTATION
    BadDoc => <
=for bad

If C encounters a bad value in any of the R, G, or B
values the output ndarray will be marked as bad and the associated R,
G, and B values will all be marked as bad.
BADDOC
);


pp_def('rgb_from_linear',
    Pars => 'rgb(c=3); gamma(); [o]out(c=3)',
    HandleBad => 1,
    Inplace => ['rgb'],
    GenericTypes=>['D'],
    Code => '
        PDL_IF_BAD(if ($ISBAD(rgb(c=>0)) || $ISBAD(rgb(c=>1)) || $ISBAD(rgb(c=>2))) {
            loop (c) %{ $SETBAD(out()); %}
        }
        else,) {
            rgb2gamma( $P(rgb), $gamma(), $P(out) );
        }
    ',
    Doc => <<'DOCUMENTATION',
=for ref

Converts an RGB color triple (presumably linear) to an RGB color triple
with the specified gamma.

=for usage

Usage:

    my $gammaed = rgb_from_linear( $rgb_linear, 2.2 );
DOCUMENTATION
    BadDoc => <
=for bad

If C encounters a bad value in any of the R, G, or B
values the output ndarray will be marked as bad and the associated R,
G, and B values will all be marked as bad.
BADDOC
);


pp_addpm(<<'EOD');


=head2 rgb_to_xyz

=for ref

Converts an RGB color triple to an XYZ color triple.

The first dimension of the ndarrays holding the rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated X, Y, and Z values will all be marked as bad.

=for usage

Usage:

    my $xyz = rgb_to_xyz( $rgb, 'sRGB' );
    my $xyz = rgb_to_xyz( $rgb, \%rgb_spec );

=cut

*rgb_to_xyz = \&PDL::rgb_to_xyz;
sub PDL::rgb_to_xyz {
    my ($rgb, $space) = @_;
    my $spec = get_space($space);
    my $m = PDL->topdl( $spec->{m} );
    return _rgb_to_xyz( $rgb, $spec->{gamma}, $m );
}


=head2 xyz_to_rgb

=for ref

Converts an XYZ color triple to an RGB color triple.

The first dimension of the ndarrays holding the xyz and rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the X, Y, or Z values the output ndarray will be marked as bad and the associated R, G, and B values will all be marked as bad.

=for usage

Usage:

    my $rgb = xyz_to_rgb( $xyz, 'sRGB' );
    my $rgb = xyz_to_rgb( $xyz, \%rgb_spec );

=cut

*xyz_to_rgb = \&PDL::xyz_to_rgb;
sub PDL::xyz_to_rgb {
    my ($xyz, $space) = @_;
    my $spec = get_space($space);
    my $mstar = exists $spec->{mstar} ? PDL->topdl( $spec->{mstar} ) : PDL->topdl( $spec->{m} )->inv;
    return _xyz_to_rgb( $xyz, $spec->{gamma}, $mstar );
}


=head2 xyz_to_lab

=for ref

Converts an XYZ color triple to an Lab color triple.

The first dimension of the ndarrays holding the xyz values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the X, Y, or Z values the output ndarray will be marked as bad and the associated L, a, and b values will all be marked as bad.

=for usage

Usage:

    my $lab = xyz_to_lab( $xyz, 'sRGB' );
    my $lab = xyz_to_lab( $xyz, \%rgb_spec );

=cut

*xyz_to_lab = \&PDL::xyz_to_lab;
sub PDL::xyz_to_lab {
    my ($xyz, $space) = @_;
    my $spec = get_space($space);
    my $w = PDL->topdl($spec->{white_point});
    return _xyz_to_lab( $xyz, $w );
}


=head2 lab_to_xyz

=for ref

Converts an Lab color triple to an XYZ color triple.

The first dimension of the ndarrays holding the lab values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the L, a, or b values the output ndarray will be marked as bad and the associated X, Y, and Z values will all be marked as bad.

=for usage

Usage:

    my $xyz = lab_to_xyz( $lab, 'sRGB' );
    my $xyz = lab_to_xyz( $lab, \%rgb_spec );

=cut

*lab_to_xyz = \&PDL::lab_to_xyz;
sub PDL::lab_to_xyz {
    my ($lab, $space) = @_;
    my $spec = get_space($space);
    my $w = PDL->topdl($spec->{white_point});
    return _lab_to_xyz( $lab, $w );
}


=head2 rgb_to_lch

=for ref

Converts an RGB color triple to an LCH color triple.

The first dimension of the ndarrays holding the rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated L, C, and H values will all be marked as bad.

=for usage

Usage:

    my $lch = rgb_to_lch( $rgb, 'sRGB' );
    my $lch = rgb_to_lch( $rgb, \%rgb_spec );

=cut

*rgb_to_lch = \&PDL::rgb_to_lch;
sub PDL::rgb_to_lch {
    my ($rgb, $space) = @_;
    my $spec = get_space($space);
    my $lab = xyz_to_lab( rgb_to_xyz( $rgb, $spec ), $spec );
    return lab_to_lch( $lab );
}


=head2 lch_to_rgb

=for ref

Converts an LCH color triple to an RGB color triple.

The first dimension of the ndarrays holding the lch values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the L, C, or H values the output ndarray will be marked as bad and the associated R, G, and B values will all be marked as bad.

=for usage

Usage:

    my $rgb = lch_to_rgb( $lch, 'sRGB' );
    my $rgb = lch_to_rgb( $lch, \%rgb_spec );

=cut

*lch_to_rgb = \&PDL::lch_to_rgb;
sub PDL::lch_to_rgb {
    my ($lch, $space) = @_;
    my $spec = get_space($space);
    my $xyz = lab_to_xyz( lch_to_lab( $lch ), $spec );
    return xyz_to_rgb( $xyz, $spec );
}


=head2 rgb_to_lab

=for ref

Converts an RGB color triple to an LAB color triple.

The first dimension of the ndarrays holding the rgb values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the R, G, or B values the output ndarray will be marked as bad and the associated L, A, and B values will all be marked as bad.

=for usage

Usage:

    my $lab = rgb_to_lab( $rgb, 'sRGB' );
    my $lab = rgb_to_lab( $rgb, \%rgb_spec );

=cut

*rgb_to_lab = \&PDL::rgb_to_lab;
sub PDL::rgb_to_lab {
    my ($rgb, $space) = @_;
    my $spec = get_space($space);
    return xyz_to_lab( rgb_to_xyz( $rgb, $spec ), $spec );
}


=head2 lab_to_rgb

=for ref

Converts an LAB color triple to an RGB color triple.

The first dimension of the ndarrays holding the lab values must be size 3, i.e. the dimensions must look like (3, m, n, ...).

=for bad

If C encounters a bad value in any of the L, A, or B values the output ndarray will be marked as bad and the associated R, G, and B values will all be marked as bad.

=for usage

Usage:

    my $rgb = lab_to_rgb( $lab, 'sRGB' );
    my $rgb = lab_to_rgb( $lab, \%rgb_spec );

=cut

*lab_to_rgb = \&PDL::lab_to_rgb;
sub PDL::lab_to_rgb {
    my ($lab, $space) = @_;
    my $spec = get_space($space);
    return xyz_to_rgb( lab_to_xyz( $lab, $spec ), $spec );
}


=head2 add_rgb_space

Supports adding custom RGB space definitions. The C and C
can be supplied as PDL ndarrays if desired. As of 0.202, you don't need
to provide an C since the inverse of the C will be calculated
(once) as a default.

Usage:

    my %custom_space = (
        custom_1 => {
          'gamma' => '2.2',
          'm' => [
                   [
                     '0.467384242424242',
                     '0.240995',
                     '0.0219086363636363'
                   ],
                   [
                     '0.294454030769231',
                     '0.683554',
                     '0.0736135076923076'
                   ],
                   [
                     '0.18863',
                     '0.075452',
                     '0.993451333333334'
                   ]
                 ],
          'white_point' => [
                             '0.312713',
                             '0.329016'
                           ],
        },
        custom_2 => { ... },
    );

    add_rgb_space( \%custom_space );

    my $rgb = lch_to_rgb( $lch, 'custom_1' );

=cut

*add_rgb_space = \&PDL::Graphics::ColorSpace::RGBSpace::add_rgb_space;
*get_space = \&PDL::Graphics::ColorSpace::RGBSpace::get_space;


=head1 SEE ALSO

Graphics::ColorObject

=head1 AUTHOR

~~~~~~~~~~~~ ~~~~~ ~~~~~~~~ ~~~~~ ~~~ `` ><(((">

Copyright (C) 2012 Maggie J. Xiong

Original work sponsored by Shutterstock, LLC L

All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.

=cut

EOD

pp_done();