算法：图像彩色空间的转换(RGB2YCbCr、YCbCr2RGB)

clear all;

close all;

clc;

img=imread('lena_color.jpg');

%img=mat2gray(img);  %任意区间映射到[0,1];

[m n dim]=size(img);

imshow(img);

%%图像的RGB

img=double(img);

R=img(:,:,1);

G=img(:,:,2);

B=img(:,:,3);

%%RGB2YCbCr

Y=zeros(m,n);    %亮度

Cb=zeros(m,n);   %彩度

Cr=zeros(m,n);   %浓度

matrix=[0.299 0.587 0.114;

       -0.1687 -0.3313 0.5;

       0.5 -0.4187 -0.0813];

for i=1:m

  for j=1:n

       tmp=matrix*[R(i,j) G(i,j) B(i,j)]';

       Y(i,j)=tmp(1);

       Cb(i,j)=tmp(2)+128;

       Cr(i,j)=tmp(3)+128;

  end

end

%%YCbCr2RGB

matrix=inv(matrix);

for i=1:m

  for j=1:n

       tmp=matrix*[Y(i,j) Cb(i,j)-128 Cr(i,j)-128]';

       R(i,j)=tmp(1);

       G(i,j)=tmp(2);

       B(i,j)=tmp(3);

  end

end

%%如果正反变换都没错的话，那么图像是不变的

img(:,:,1)=R;

img(:,:,2)=G;

img(:,:,3)=B;

figure;

imshow(uint8(img))