From 3de1f2611036dc0faebe7c45f6b68332f4dd3148 Mon Sep 17 00:00:00 2001 From: avinal <185067@nith.ac.in> Date: Thu, 11 Feb 2021 23:52:15 +0530 Subject: [PATCH] high pass filter added --- digital_image_processing/lab/lab2/highpass.m | 33 ++++++++++++++++++++ 1 file changed, 33 insertions(+) create mode 100644 digital_image_processing/lab/lab2/highpass.m diff --git a/digital_image_processing/lab/lab2/highpass.m b/digital_image_processing/lab/lab2/highpass.m new file mode 100644 index 0000000..6dcbb3c --- /dev/null +++ b/digital_image_processing/lab/lab2/highpass.m @@ -0,0 +1,33 @@ +input_image = imread('download.jfif'); +image = rgb2gray(input_image); +[M, N] = size(image); + +FT_img = fft2(double(image)); + +D0 = 5; + + +u = 0:(M-1); +idx = find(u>M/2); +u(idx) = u(idx)-M; +v = 0:(N-1); +idy = find(v>N/2); +v(idy) = v(idy)-N; + + +[V, U] = meshgrid(v, u); + + +D = sqrt(U.^2+V.^2); + + +H = double(D > D0); + + +G = H.*FT_img; + +output_image = real(ifft2(double(G))); + + +subplot(2, 1, 1), imshow(image), +subplot(2, 1, 2), imshow(output_image, [ ]);