File Exchange

image thumbnail

Low Freq 2D Fourier Transform

version 1.3 (2.06 KB) by

Produces high resolution 2D Fourier Transform for low spatial frequencies.

1 Download

Updated

View License

Useage case:
You are only interested in the low spatial frequencies of your matrix, but you want a high resolution transform. Using fft2 with lots of zero-padding you get a decent resolution but it takes too long. What can you do? Try sft2_low!

Algorithm:
This function uses the basic Fourier transform rather than the FFT, but because it is only evaluating it at a small number of points the total execution time is considerably less (it can be 20x faster or more).

The code is fully vectorised, so it's pretty quick even though it's doing a lot of arithmetic. At some point I may generalise the code to allow you to specify your exact frequencies of interest.

Warning:
The output should be correct, but I recommend testing it and comparing the values with fft2. Please report problems and suggestions in the comments.

Comments and Ratings (2)

Ali

Ali (view profile)

Wow. This is even better than the first version and even faster.

Ali

Ali (view profile)

This is awesome. Thank you.

Updates

1.3

improved comments

1.1

Code is now optimized to work specifically on real-valued inputs. It also now optionally gives the frequencies as outputs.

MATLAB Release
MATLAB 8.0 (R2012b)
Acknowledgements

Inspired by: Numerical Methods for Physics

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video