29/8/25 11:21 AM D:\B1 BIP\P4.
m 1 of 2
clc;
clear all;
close all;
im1=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\9.
tif');
im2=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\10.
tif');
im3=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\11.
tif');
im4=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\12.
tif');
im5=imread ('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\13.
tif');
% image histogram
figure,subplot(4,2,1),imshow(im1),title('dark image');
subplot(4,2,2),imhist(im1),title('dark image histogram');
subplot(4,2,3),imshow(im2),title('white image');
subplot(4,2,4),imhist(im2),title('white image histogram');
subplot(4,2,5),imshow(im3),title('low contrast image');
subplot(4,2,6),imhist(im3),title('low contrast image histogram');
subplot(4,2,7),imshow(im4),title('high contrast image');
subplot(4,2,8),imhist(im4),title('high contrast image histogram');
% % histogram equalisation
g1 = histeq(im1);
figure,subplot(4,3,1),imshow(im1);title('dark image histogram');
subplot(4,3,2),imshow(g1),title('histogram equalised image');
subplot(4,3,3),imhist(g1);
g2 = histeq(im2);
subplot(4,3,4),imshow(im2),title('white image');
subplot(4,3,5),imshow(g2),title('histogram equalised image');
subplot(4,3,6),imhist(g2);
g3 = histeq(im3);
subplot(4,3,7),imshow(im3),title('low contrast image');
subplot(4,3,8),imshow(g3),title('histogram equalised image');
subplot(4,3,9),imhist(g3);
g4 = histeq(im4);
subplot(4,3,10),imshow(im4),title('high contrast image');
subplot(4,3,11),imshow(g4),title('histogram equalised image');
subplot(4,3,12),imhist(g4);
% original image and histogram
figure,subplot(4,2,1),imshow(im5);title('original image')
subplot(4,2,2),imhist(im5);
hnorm = imhist(im5)./numel(im5);
cdf1 = cumsum(hnorm);
x=[Link];
subplot(4,2,3),plot(x,cdf1);
g5 = histeq(im5);
subplot(4,2,4),imshow(g5);title('histogram equalized image')
subplot(4,2,5),imhist(g5);
%
% %Specify New Histogram
29/8/25 11:21 AM D:\B1 BIP\P4.m 2 of 2
x = 0:255; y = 1:256;
y(1:20) = 0:100000/(20-1):100000;
y(21:30) = 100000:(5000-100000)/(20-11):5000;
y(31:100) = 5000:(-5000)/(100-31):0;
y(101:200) = 0:5000/(200-101):5000;
y(201:256) = 5000:(-5000)/(256-201):0;
hsum = sum(y); y = y./hsum; y = y*length(im5(:));
x1=[Link];
subplot(4,2,6),plot(x1,y);
% % histogram matching
g21 = histeq(im5,y);
subplot(4,2,7),imshow(g21);title('histogram matching image')
subplot(4,2,8),imhist(g21);
