Hello,
You can refer to the below code to perform FFT on series of images and to compute an average FFT per image.
% Read in the four images
img1 = imread('img.jpg');
img2 = imread('img.jpg');
img3 = imread('img.jpg');
img4 = imread('img.jpg');
% Store the images in a cell array
imgs = {img1, img2, img3, img4};
% Compute the FFT for each image and store the magnitudes in a cell array
ffts = cell(1, 4);
for i = 1:4
ffts{i} = abs(fft2(imgs{i})).^2; % Square the magnitudes to get power
end
% Compute the average FFT power across all images
avg_fft_power = mean(cat(3, ffts{:}),'all');
disp(['Average FFT power: ', num2str(avg_fft_power)]);
To compute average FFT power over all the frequencies of an image, You can follow the code shown below
% Read in the image
img = imread('img.jpg');
% Compute the FFT and power spectrum
fft_power = abs(fft2(img)).^2;
% Compute the average FFT power across all frequencies in the image
avg_fft_power = mean(fft_power, 'all');
% Display the average FFT power
disp(['Average FFT power: ', num2str(avg_fft_power)]);
