matlab练习程序,RGB2YCbCr、YCbCr2RGB

  这个YCbCr是从YUV派生出来的,比较适合处理数字图像,JPEG压缩就是在这个彩色空间中处理的。变换公式

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))