Gaussian_Low_Pass_Filter.m

function output_image = Gaussian_Low_Pass_Filter(input_image)
    if size(input_image, 3) == 3
        input_image = rgb2gray(input_image);
    end
    cutoff_frequency = 50;  % You can set this to any value or get it from GUI
    
    % Step 1: Obtain the size of the image
    [M, N] = size(input_image);
    
    % Step 2: Transform the image from spatial domain to frequency domain
    F = fft2(double(input_image));
    F_shifted = fftshift(F);
    
    % Create the Gaussian Low Pass Filter (GLPF)
    [X, Y] = meshgrid(1:N, 1:M);
    D = sqrt((X - N/2).^2 + (Y - M/2).^2);
    
    % Gaussian Low Pass Filter Equation
    sigma = cutoff_frequency / 2;  % You can adjust sigma based on cutoff frequency
    H = exp(- (D.^2) / (2 * sigma^2));  % Gaussian filter
    
    % Step 3: Apply the filter
    F_filtered = F_shifted .* H;
    
    % Step 4: Inverse Fourier Transform to get the filtered image
    F_inverse_shifted = ifftshift(F_filtered);
    output_image = real(ifft2(F_inverse_shifted));
    output_image=mat2gray(output_image);
end