这里的变换公式是模拟领域的公式,这里虽然是数字图像,不过也可以这样写,YCbCr采样另外的公式。
1 clear all; 2 close all; 3 clc; 4 5 img=imread('lena_color.jpg'); 6 img=mat2gray(img); %任意区间映射到[0,1]; 7 [m n dim]=size(img); 8 imshow(img); 9 %%图像的RGB 10 R=img(:,:,1); 11 G=img(:,:,2); 12 B=img(:,:,3); 13 14 %%RGB2YUV 15 Y=zeros(m,n); %亮度 16 U=zeros(m,n); %彩度 17 V=zeros(m,n); %浓度 18 matrix=[0.299 0.587 0.114; 19 -0.14713 -0.28886 0.436; 20 0.615 -0.51498 -0.10001]; 21 for i=1:m 22 for j=1:n 23 tmp=matrix*[R(i,j) G(i,j) B(i,j)]'; 24 Y(i,j)=tmp(1); 25 U(i,j)=tmp(2); 26 V(i,j)=tmp(3); 27 end 28 end 29 30 %%YUV2RGB 31 matrix=inv(matrix); 32 for i=1:m 33 for j=1:n 34 tmp=matrix*[Y(i,j) U(i,j) V(i,j)]'; 35 R(i,j)=tmp(1); 36 G(i,j)=tmp(2); 37 B(i,j)=tmp(3); 38 end 39 end 40 41 %%如果正反变换都没错的话,那么图像是不变的 42 img(:,:,1)=R; 43 img(:,:,2)=G; 44 img(:,:,3)=B; 45 figure; 46 imshow(img)