Hi Ahereza,
From what I gather you are trying to compress images of same dimensions for PSNR analysis but are getting images of different dimensions.
To achieve the desired results, you can go through the following steps along with the attached code snippets:
1. Implement a ‘try-catch’ block to handle potential errors when loading the image file, ensuring that the script does not crash if the file is missing, or the path is incorrect.
try
a = imread('img.png');
catch
error('File "img.png" does not exist.');
end
2. Convert the image to grayscale using ‘rgb2gray’ function and store the dimensions in ‘row’ and ‘col’ variables.
grayIm = rgb2gray(a);
[row, col] = size(grayIm);
3. Apply Fourier Transformation using ‘fft2’ function converting it into frequency domain and sort the Fourier coefficients by magnitude to facilitate thresholding based on size.
% Compute the 2D Fourier Transform
A = fft2(grayIm);
% Sort the Fourier coefficients by magnitude
Btsort = sort(abs(A(:)));
4. Iterate over pre-defined compression levels, using a threshold to retain significant Fourier coefficients. Then apply ‘ifft2’ to reconstruct the image from these retained components.
% Loop over different levels of compression
for keep = [0.99, 0.05, 0.01, 0.002]
thresh = Btsort(floor((1 - keep) * length(Btsort)));
ind = abs(A) > thresh;
count = row * col - sum(ind(:));
Alow = A .* ind;
per = 100 * sum(ind(:)) / (row * col);
Blow = ifft2(Alow);
5. Extract the real part of the reconstructed image, convert it into ‘uint8’ format and display the compressed image.
Blow = uint8(real(Blow));
figure;
imshow(Blow);
title([num2str(per) '% of FFT basis']);
For better understanding, you can refer to the output attached below:
For more information, kindly go through the documentations linked below:
- ‘fft2’ function: https://www.mathworks.com/help/matlab/ref/fft2.html
- ‘ifft2’ function: https://www.mathworks.com/help/matlab/ref/ifft2.html
- ‘rgb2gray’ function: https://www.mathworks.com/help/matlab/ref/rgb2gray.html
Happy Coding!
