色空間の変換（sRGBからLab, XYZ, Yxy）

function [XYZ, xyY, Lab]= mySRGBtoColData(SRGBdata, illumiType)

% 引数SRGBdataは、3つの行ベクトルからなり、それぞれの行がRGBの値（0から255）を取ります。

% SRGBのデータは、D65の光源を想定したデータです。

% 引数illumiTypeは、C光源（'C'）、D50光源（'D50'）、D65光源（'D65'）のいずれかです。

% これら以外の光源を使う場合には、http://www.brucelindbloom.com/index.html?Eqn_ChromAdapt.html

% を参考にソースを書き換えてください。

origRGB = SRGBGammaUncorrect(SRGBdata); % ガンマ補正前のデータに戻す（戻り値は、0から1の範囲のデータであることに注意）

XYZdata = SRGBPrimaryToXYZ(origRGB); % 補正前のRGBの値をXYZに変換。

% 光源が変化する場合には、Bradford変換などを用いて、XYZの色データを変換する必要がある。

% もともとのSRGBのデータは光源がD65であることに注意。

switch illumiType

    case 'C'

        bradfordData = [ % D65 -> C

            1.0097785 0.0070419 0.0127971

            0.0123113 0.9847094 0.0032962

            0.0038284 -0.0072331 1.0891639];

        illuminant = transpose([0.98074 1.00000 1.18232]); % illuminant C

    case 'D50'

        bradfordData = [ % D65 -> D50

            1.0478112 0.0228866 -0.0501270

            0.0295424 0.9904844 -0.0170491

            -0.0092345 0.0150436 0.7521316];

    case 'D65'

        bradfordData = eye(3); % 変換の必要なし。

    otherwise

        error('Invalid illumiType!');

        %return;

end;

XYZdata = bradfordData * XYZdata; % Bradford変換

Lab = XYZToLab(XYZdata, illuminant); % XYZ -> Lab空間への変換

XYZ = XYZdata .* 100;

xyY = XYZToxyY(XYZ); % XYZ -> xyY空間への変換

XYZToLab

clear all;

dataFileName = 'rgbData.csv'; % RGBのカンマ区切りのデータ

illumiType = 'D50'; % Labの光源の情報。一般的にはD50だが、C または D65 にも対応

illumiType = upper(illumiType); % 小文字から大文字へ

colData = csvread(dataFileName); % CSVファイルから読み取り

[pathstr, fileName, exqt] = fileparts(dataFileName); % 拡張子を取り除いたファイル名を取得。

[XYZ, xyY, Lab]= mySRGBtoColData(transpose(colData), illumiType);

csvwrite([fileName '-Lab-' illumiType '.csv'], transpose(Lab));

csvwrite([fileName '-XYZ-' illumiType '.csv'], transpose(XYZ)); % bradford変換後であることに注意