On Wed, May 09, 2007 at 04:40:36PM +0200, fred wrote: > I use an FFT 2D on a matrix to compute a convolution like this: > > a = fft2(input_data) > b = fft2(output_data) > c = real(ifft2(a*b)) Also take a look at scipy.signal.fftconvolve. Cheers Stéfan