rgb转换为lab源代码

在处理图像的过程中常常需要对图像的不同格式进行处理。
例如需要将rgb转换为lab，它们分属于不同的颜色空间。
一般来说rgb中，r、g、b的取值范围为（0，255）或者在（0，1）之间。
而lab中，l的范围为（0，100），a、b范围为（-120，120）。
而在转化过程中，通常需要将rgb先转化为xyz控件，然后再转换为lab空间。
因此需要用到两个函数，分别为，rgb2xyz,rgb2lab
下面是一段代码：
rgb2lab,来自：http://server.cs.ucf.edu/~vision/Code/ColorSpaces/rgb2lab.m

1. functiony=rgb2lab(im)
2. %rgb2lab(image)
3. %ConvertedeRGBparaLAB
5. %Defineasconstantes
6. Xn=95.13;
7. Yn=100.0;
8. Zn=108.86;
10. disp('IniciandoaconversaodeRGBparaLAB');
11. %ConvertedeRGBparaXYZ
12. xyz=rgb2xyz(im);
13. bx=xyz(:,:,1);
14. by=xyz(:,:,2);
15. bz=xyz(:,:,3);
17. %CalculaX/Xn,Y/YneZ/Zn
18. bx=bx/Xn;
19. by=by/Yn;
20. bz=bz/Zn;
22. %Determinaoselem.maioresouiguaisa8.856E-3
23. bin=(bx>=8.856E-3);
24. f_x=(bin.*bx).^(1/3);
25. bin=(by>=8.856E-3);
26. yyn=(bin.*by).^(1/3);
27. f_y=yyn;
28. bin=(bz>=8.856E-3);
29. f_z=(bin.*bz).^(1/3);
31. l=(116*yyn)-16;
33. %Determinaoselem.menoresque8.856E-3
34. bin=(bx<8.856E-3);
35. f_x=f_x+(7.787*(bin.*bx)+(16/116));
36. bin=(by<8.856E-3);
37. yyn=(bin.*by);
38. f_y=f_y+(7.787*yyn+(16/116));
39. bin=(bz<8.856E-3);
40. f_z=f_z+(7.787*(bin.*bz)+(16/116));
43. %CalculaL,AeB
44. l=l+(903.3*(yyn.^(1/3)));
45. a=(500*l).*(f_x-f_y);
46. b=(200*l).*(f_y-f_z);
48. y(:,:,1)=l;
49. y(:,:,2)=a;
50. y(:,:,3)=b;

另外一端代码，rgb2xyz，来自：http://scien.stanford.edu/class/psych221/projects/08/JaeMoPark/appendix/rgb2xyz.m

1. function[xyz]=rgb2xyz(rgb)
3. a=double(rgb);
4. mat=[.49.31.20;.17697.81240.01063;.00.01.99];
5. mat=(1/.17697)*mat;
7. R=a(:,:,1);
8. G=a(:,:,2);
9. B=a(:,:,3);
10. [m,n]=size(R);
11. RGB=[R(:)G(:)B(:)]';
12. XYZ=mat*RGB;
14. X=XYZ(1,:);x=reshape(X,m,n);
15. Y=XYZ(2,:);y=reshape(Y,m,n);
16. Z=XYZ(3,:);z=reshape(Z,m,n);
17. xyz=zeros(size(a));
18. xyz(:,:,1)=x;
19. xyz(:,:,2)=y;
20. xyz(:,:,3)=z;
22. end