Can someone break down this code and tell me what it does so I can better understand it?

So far I know the following:

This code first asked the user to choose an image file of choice. The rgb2gray command converts the selected image to grey scale. In the main for loop, the Discrete Fourier Transform of the image is calculated. Once the transform has been done calculating. It is changed to a colormap to reflect the magnitudes of the transform. The absolute value is then plotted on a graph.

Here is the code:

clc;
clear all;
close all;

[filename, pathname, filterindex] = uigetfile( ...
{  '*.jpg','JPEG (*.jpg)'; ...
   '*.bmp','Windows Bitmap (*.bmp)'; ...
   '*.fig','Figures (*.fig)'; ...
   '*.*',  'All Files (*.*)'}, ...
   'Choose image(s) to be processed', ...
   'MultiSelect', 'off');

if filterindex==0, break;end
filename=cellstr(filename);
im2= imread(horzcat(pathname,char(filename)));

im1=rgb2gray(im2);
im3=im2double(im1);
[n,m]=size(im3);

c1=0;
h = waitbar(0,'Calculating XXX please wait......'); 
k=1;l=1;
for l=0:1:m-1
    for k=0:1:n-1
        for x=0:1:n-1
            for y=0:1:m-1
                a=x+1;b=y+1;
                c= im3(a,b) * exp(-1i*2*pi*(k*x/n + l*y/m));
                c1=c1+c;
            end
        end
        aa=l+1;bb=k+1;
        im(bb,aa)=c1;
        c1=0;
    end
    waitbar(l / m);
end
ims = im*255;
close(h)

imshow(ims);title('XXX plot');
figure
imshow(log(abs(ims)),[-1 5]); colormap(jet); colorbar;title('absolute value of XXX plot');

[Rik]

Why did you delete the old question? I'm not going to invest time in writing code for it to be ignored and the question being deleted.

[Westin Messer]

I'm sorry I shouldn't have deleted it. I didn't want people to think I was posting a duplicate question. This time I just wanted to understand what was going on in the code instead of trying to make it run.

[Stephen23]

@Westin Messer: it is often best to keep a discussion in one thread: it makes it easier for everyone to keep track of what information and clarifications have been given.

[Westin Messer]

That makes sense, I'll defiantly do that in the future. Sorry guys!

[Rena Berman]

(Answers dev) I undeleted the question.

[John BG]

good, can we now focus on the question?

[Westin Messer]

Yes please!

Answer 1

[John BG]

Hi Westin Messer

this is John BG <mailto:[e-mail] [e-mail]>

1. instead of all those for loops, for a 2D fft, would you consider using the command fft2 instead?

clc;clear all;close all;

im1=imread('atrium.jpg');
figure;imshow(im1);im2=rgb2gray(im1);
figure;imshow(im2);im3=im2double(im2);
figure;imshow(im3);

ims=fft2(im3);

figure;imshow(abs(ims));title('this is FFT2D(im) plot');

figure;imshow(10*log10(abs(ims)));
colormap(jet); colorbar;title('absolute value of FFT(im) plot');

2.-

there's something wrong with the way you are attempting to use the waitbars.

If I were you I would leave this cosmetics feature for the end of the coding process. Just a hint, if you produce a handle to the waitbar named 'h' then you should be using that handle in the gradual update of waitbar, but you are currently not using such handle, in the update I mean.

Contrary to many other MATLAB commands, the handle goes at the end of input fields of waitbar.

3.-

Instead of log have you tried 10*log10?

figure;surf(abs(ims),'EdgeColor','none');
figure;surf(10*log10(abs(ims)),'EdgeColor','none');

4.-

At Mathworks they allow us to read the invoking part of the function fft2, key in the following

type fft2

Basically fft2 uses fft recursively.

fft2 explained here

https://uk.mathworks.com/help/matlab/ref/fft2.html?s_tid=srchtitle

the command fft2 performs a spectral sorting of the image components, in X dimension and Y dimension.

A vertical curtain has all horizontal variations and no vertical variations.

A horizontal curtain has all vertical variations and no horizontal variations.

In the equivalent SIMULINK block for FFT 2D there is a short explanation

https://uk.mathworks.com/help/vision/ref/2dfft.html?s_tid=srchtitle

.

Herr Messer

If you find this answer useful would you please be so kind to consider marking my answer as Accepted Answer?

To any other reader, if you find this answer useful please consider clicking on the thumbs-up vote link

thanks in advance

John BG

<mailto:[e-mail] [e-mail]>

[John BG]

Danke Schön