图像YCbCr与RGB空间相互转换与MATLAB实现

Halcom

编写RGB到YCbCr颜色空间的函数程序如下：

1. function im1 = rgb2ycbcr(im)
   
2. % 转化矩阵
   
3. T = [65.481,128.553,24.966;
   
4.     -37.797,-74.203,112.0;
   
5.     112.0,-93.786,-18.214];  % 矩阵系数
   
6. Ta = [16;128;128];
   
7. 
   
8. R = im(:,:,1);  % R
   
9. G = im(:,:,2);  % G
   
10. B = im(:,:,3);  % B
    
11. R = im2double(R); % 转化为double类型
    
12. G = im2double(G); % 转化为double类型
    
13. B = im2double(B); % 转化为double类型
    
14. 
    
15. Y = Ta(1,1) + T(1,1).* R + T(1,2).*G + T(1,3).*B;
    
16. Cb = Ta(2,1) + T(2,1).* R + T(2,2).*G + T(2,3).*B;
    
17. Cr = Ta(3,1) + T(3,1).* R + T(3,2).*G + T(3,3).*B;
    
18. im1(:,:,1) = Y;
    
19. im1(:,:,2) = Cb;
    
20. im1(:,:,3) = Cr;
    
21. im1 = uint8(im1);  % 类型转换
    
22. 调用该函数，程序如下：
    
23. clc,clear,close all  % 清理命令区、清理工作区、关闭显示图形
    
24. warning off       % 消除警告
    
25. feature jit off      % 加速代码运行
    
26. im = imread('coloredChips.png');
    
27. figure('color',[1,1,1])
    
28. subplot(121),imshow(im,[]);title('RGB')
    
29. im1 = rgb2ycbcr(im);    % RGB转化为YCbCr
    
30. subplot(122),imshow(im1,[]);title('YCbCr')
    

复制代码

编写YCbCr到RGB颜色空间的函数程序如下：

1. function im1 = ycbcr2rgb(im)
   
2. % 转化矩阵
   
3. T = [65.481,128.553,24.966;
   
4.     -37.797,-74.203,112.0;
   
5.     112.0,-93.786,-18.214];  % 矩阵系数
   
6. Ta = [16;128;128];
   
7. 
   
8. TT = inv(T);   % 求逆矩阵
   
9. Tb = TT*Ta;
   
10. 
    
11. Y = im(:,:,1);   % Y
    
12. Cb = im(:,:,2);  % Cb
    
13. Cr = im(:,:,3);  % Cr
    
14. Y = im2double(Y);   % 转化为double类型
    
15. Cb = im2double(Cb); % 转化为double类型
    
16. Cr = im2double(Cr);  % 转化为double类型
    
17. 
    
18. R = -Tb(1,1) + TT(1,1).* Y + TT(1,2).*Cb + TT(1,3).*Cr;
    
19. G = -Tb(2,1) + TT(2,1).* Y + TT(2,2).*Cb + TT(2,3).*Cr;
    
20. B = -Tb(3,1) + TT(3,1).* Y + TT(3,2).*Cb + TT(3,3).*Cr;
    
21. R = mat2gray(R); % 灰度值 转化为 0 - 1之间
    
22. G = mat2gray(G); % 灰度值 转化为 0 - 1之间
    
23. B = mat2gray(B); % 灰度值 转化为 0 - 1之间
    
24. im1(:,:,1) = R;  
    
25. im1(:,:,2) = G;
    
26. im1(:,:,3) = B;
    
27. im1 = im2uint8(im1);  % 类型转换
    
28. 调用该函数，程序如下：
    
29. %% YCbCr -->  RGB
    
30. clc,clear,close all  % 清理命令区、清理工作区、关闭显示图形
    
31. warning off       % 消除警告
    
32. feature jit off      % 加速代码运行
    
33. im = imread('coloredChips.png');
    
34. im1 = rgb2ycbcr(im);    % RGB转化为YCbCr
    
35. figure('color',[1,1,1])
    
36. subplot(121),imshow(im1,[]);title('YCbCr')
    
37. im2 = ycbcr2rgb(im1);
    
38. subplot(122),imshow(im2,[]);title('RGB')

复制代码