parallel computing for an image stack

Dear all,

I've got an image stack with the size of 1024 pixels *1024 pixels *1600 frames. I call it M(x,y,z). I'd like to do some computations like smooth and FFT for all the pixels simultaneously. I don't want to use parfor on CPU cores. How can I realize it using GPU?

Answer 1

[Edric Ellis]

I'm not sure exactly what you're after here, but note that the gpuArray implementation of fft2 automatically applies to "pages" of a 3-dimensional array, like this:

x = ones(2, 2, 3, "gpuArray") .* reshape(1:3, 1, 1, 3)

x = 
x(:,:,1) =

     1     1
     1     1


x(:,:,2) =

     2     2
     2     2


x(:,:,3) =

     3     3
     3     3

fft2(x)

ans = 
ans(:,:,1) =

     4     0
     0     0


ans(:,:,2) =

     8     0
     0     0


ans(:,:,3) =

    12     0
     0     0