[求助]从RGB空间到LAB空间的转换程序
 

求助，从RGB空间到LAB空间的转换程序
以下是一些简单的算法思想，希望谁能帮帮我
从RGB到XYZ空间，从XYZ空间到LMS空间，再取对数，最后从LMS到LAB空间

我正在做黑白图像的彩色化处理，其中有段这个程序，是我自己编的，但不清楚是否正确，请帮忙一起看看，谢谢!
[X,map]=imread('source.bmp');
r=X(:,:,1);
g=X(:,:,2);
b=X(:,:,3);
L=0.3811*r+0.5783*g+0.0402*b;
M=0.1967*r+0.7244*g+0.0782*b;
S=0.0241*r+0.1288*g+0.8444*b;
L=double(L)/255;
M=double(M)/255;
S=double(S)/255;
L1=log(L);
M1=log(M);
S1=log(S);
l=0.5774*L1+0.5774*M1+0.5774*S1;
a=0.4082*L1+0.4082*M1-0.8165*S1;
b=0.7071*L1-0.7071*M1;
lab=cat(3,l,a,b);