29/8/25 11:22 AM D:\B1 BIP\P5.
m 1 of 2
clc;
clear all;
close all;
% im=imread('C:\Documents and Settings\All Users\Documents\My Pictures\Sample
Pictures\Water lilies.jpg');
im1=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\14.
tif');
im2=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\15.
tif');
im3=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\16.
tif');
im4=imread('D:\Biomedical Image Processing\sample images\DIP3E_Original_Images_CH03\14.
tif');
% x=rgb2gray(im);
x=im1;
figure(1),subplot(2,3,1),imshow(x),title('original image');
% spatial low pass filtering using different sizes masks
% ns=input('input low pass filter mask size (odd number): ');
B=ones(3,3)/(3*3);
y=filter2(B,x);
subplot(2,3,2),imagesc(y),colormap('gray'),title('low passed image filter mask=[3 3]');
c=ones(5,5)/(5*5);
y1=filter2(c,x);
subplot(2,3,3),imagesc(y1),colormap('gray'),title('low passed image filter mask=[5 5]');
d1=ones(9,9)/(9*9);
y2=filter2(d1,x);
subplot(2,3,4),imagesc(y2),colormap('gray'),title('low passed image filter mask=[9 9]');
d2=ones(15,15)/(15*15);
y2=filter2(d2,x);
subplot(2,3,5),imagesc(y2),colormap('gray'),title('low passed image filter mask=[15
15]');
d3=ones(35,35)/(35*35);
y2=filter2(d3,x);
subplot(2,3,6),imagesc(y2),colormap('gray'),title('low passed image filter mask=[35
35]');
% low pass filter application
d4=ones(15,15)/(15*15);
y21=filter2(d4,im3);
y21=uint8(y21);
figure(2),subplot(1,3,1),imshow(im3),colormap('gray'),title('original image]');
subplot(1,3,2),imshow(y21),colormap('gray'),title('low passed image filter mask=[15
15]');
% binary image
y23=im2bw(y21,0.3);
subplot(1,3,3),imshow(y23),title('binary image]');
% median filtering
% y4=imnoise(im4,'salt & pepper');
y5=medfilt2(im4);
figure(3),subplot(2,1,1),imshow(im4),colormap('gray'),title('image with salt & pepper
noise ')
subplot(2,1,2),imshow(y5),colormap('gray'),title('median filtered image');
29/8/25 11:22 AM D:\B1 BIP\P5.m 2 of 2
% min filter
y6 = ordfilt2(x,1,ones(5,5));
figure(4),subplot(2,1,1),imshow(x),title('original image ')
subplot(2,1,2),imshow(y6),title('min filtered image');
% max filter
y7 = ordfilt2(x,25,ones(5,5));
figure(5),subplot(2,1,1),imshow(x),title('original image ')
subplot(2,1,2),imshow(y7),title('max filtered image');
