%对比两张图片的Ycbcr分量的差值并形成直方图，改个路径就能用，图片分辨率要一致

close all;
clear all;
clc;
I1=imread('E:\test\resources\image\1.jpg');
I2=imread('E:\test\resources\image\2.jpg');

ycbcr1 = rgb2ycbcr(I1);
ycbcr2 = rgb2ycbcr(I2);

% 提取色度分量,Y（亮度）,Cb 和 Cr（色差）：色度或色差分量。Cb 表示蓝色和亮度之间的差异，而 Cr 表示红色和亮度之间的差异。
Y1 = ycbcr1(:, :, 1);
cb1 = ycbcr1(:, :, 2);
cr1 = ycbcr1(:, :, 3);
Y2 = ycbcr2(:, :, 1);
cb2 = ycbcr2(:, :, 2);
cr2 = ycbcr1(:, :, 3);

%显示彩色图像
figure;
subplot(141);%拿subplot(221)来说，就是一个2*2的矩阵画布，1代表图片处于第一个位置：
imshow(I1);
%R分量灰度图
subplot(142);imshow(Y1);
%G分量灰度图
subplot(143);imshow(cb1);
%B分量灰度图
subplot(144);imshow(cr1);

figure;
subplot(141);%拿subplot(221)来说，就是一个2*2的矩阵画布，1代表图片处于第一个位置：
imshow(I2);
%R分量灰度图
subplot(142);imshow(Y2);
%G分量灰度图
subplot(143);imshow(cb2);
%B分量灰度图
subplot(144);imshow(cr2);

Y_fault=abs(Y1-Y2);
cb_fault=abs(cb1-cb2);
cr_fault=abs(cr1-cr2);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(Y1);
%显示红色分辨率下的直方图
subplot(132);
imhist(cb1);
%显示红色分辨率下的直方图
subplot(133);
imhist(cr1);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(Y2);
%显示红色分辨率下的直方图
subplot(132);
imhist(cb2);
%显示红色分辨率下的直方图
subplot(133);
imhist(cr2);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(Y_fault);
%显示红色分辨率下的直方图
subplot(132);
imhist(cb_fault);
%显示红色分辨率下的直方图
subplot(133);
imhist(cr_fault);