MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today
Asked by Farhad M on 19 Aug 2012

Let f(x,y,z) is a set of discrete data in 3D, namely a volume image. I want to map this image onto a deformed configuration, i.e.:

g(xx,yy,zz)=f(x-u,y-v,z-w) where u,v and w are functions of space coordinates (x,y,z).

I don't want to use interpolation for this purpose as it leads to phase errors attributed to interpolation methods. Therefore, I need to perform perfect resampling in Fourier space, for instance.

In 2D, I do this as follows:

DIM=size(f);

F=fftshift(fft2(f));

[x,y]=meshgrid(.5:DIM(2)-.5,.5:DIM(1)-.5);

g=0;

dim0=ceil((DIM+1)/2);

for i=(1:DIM(1))-dim0(1)

for j=(1:DIM(2))-dim0(2)

m=i+dim0(1);

n=j+dim0(2);

g=g+F(m,n)*exp(2*pi*1i*(j*(x-u)/DIM(2)+i*(y-v)/DIM(1)));

end end

g=real(f);

However, when it comes to 3D, this will take extremely large computational time, which is almost impossible to run even on strong computers.

I am no expert in non-uniform discrete fourier transform (NDFT), anyone knows if this problem can be solved using NDFT more efficiently, or any other better idea which leads to perfect resampling in 3D?

Your ideas are very much appreciated

## 1 Comment

Direct link to this comment:http://www.mathworks.de/matlabcentral/answers/46287#comment_95333

Sorry, in the last line: