FFT on series of images

Hi there!

I wish to perform an FFT on 10,000 images. Each pixel in the image, together with the same pixel from all the other images forms the series of data to be FFT. Is this possible? How do you import data from 10,000 files recursively?

I had code to do this is an older version of Matlab but now it's giving me the error "Undefined operator '==' for input arguments of type 'table'."

L = length (R);
Fs = 458.72
phase = 0
Frequency = 7
ZeroRef = R-ave;
Ref = sin(2*pi*Frequency*Time+phase);
Ref2 = cos(2*pi*Frequency*Time+phase);
PSD = Ref.*ZeroRef;
PSD2 = Ref2.*ZeroRef;
Y = fft(PSD,L,1);
Y2 = abs(Y/L);
Y1 = Y2(1:L/2+1);             % L/2+1 when even, L/2+1.5 when odd %
Y1(2:end-1) = 2*Y1(2:end-1);

Z = fft(PSD2,L,1);
Z2 = abs(Z/L);
Z1 = Z2(1:L/2+1);             % L/2+1 when even, L/2+1.5 when odd %
Z1(2:end-1) = 2*Z1(2:end-1);

A = fft(ZeroRef,L,1);
A2 = abs(A/L);
A1 = A2(1:L/2+1);             % L/2+1 when even, L/2+1.5 when odd %
A1(2:end-1) = 2*A1(2:end-1);

f = Fs*(0:(L/2))/L;

Amplitude
Amplitude*2
PhaseDiff

It would be good if the amplitude from each pixel can be arranged back into its pixel position, with another file for phase difference at each pixel.

Thank you in advance for any and all pointers.

[Image Analyst]

Which variable is your image?

[Jaslyn Tan]

I'm not sure I understand your question but image is a 128x128 array of tab-separated integer values. There are 10,000 of these arrays. If each image was a array in x, y direction, each image would be in the z axis.