matlab-对比两张图片的CIElab分量的差值并形成直方图

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

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

lab1 = rgb2lab(I1);
lab2 = rgb2lab(I2);

% 提取色度分量,L（亮度）、a（从绿色到红色的颜色分量）和b（从蓝色到黄色的颜色分量）

L1=lab1(:, :, 1);
A1=lab1(:, :, 2);
B1=lab1(:, :, 3);

L2=lab2(:,:,1);
A2=lab2(:,:,2);
B2=lab2(:,:,3);

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

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

L_fault=abs(L1-L2);
A_fault=abs(A1-A2);
B_fault=abs(B1-B2);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(L1);
%显示红色分辨率下的直方图
subplot(132);
imhist(A1);
%显示红色分辨率下的直方图
subplot(133);
imhist(B1);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(L2);
%显示红色分辨率下的直方图
subplot(132);
imhist(A2);
%显示红色分辨率下的直方图
subplot(133);
imhist(B2);

figure;
%显示红色分辨率下的直方图
subplot(131);
imhist(L_fault);
%显示红色分辨率下的直方图
subplot(132);
imhist(A_fault);
%显示红色分辨率下的直方图
subplot(133);
imhist(B_fault);