FourierTools.jl | Working with the Frequency Space
2021-07-30 , Purple

FourierTools.jl aims at simplifying work in Fourier/Frequency space without loosing efficiency.
We provide several convenient wrappers to speed-up the common fft(fftshift(x)) pattern.
This package also brings functionality to up and downsample signals through sinc interpolation.
Furthermore, based on FFTs it provides shearing, rotation, convolution and (sub) pixel shift functions which can be applied to N-dimensional data efficiently.


Fourier space is commonly used for convolution operations, as the Fast Fourier Transformation (FFT) is, as its name may suggest, O(N log N) fast. The FFT algorithm typically produces data at a mangled form that makes it difficult to directly apply functions to. fftshift is a way to deal with this but involves data copies.
Based on the packages ShiftedArrays.jl and PaddedViews.jl, the FourierTools.jl package implements views to the results of the FFTW routines fft and rfft and their inverse ifft and irfft including the respective fftshift operations but implemented as views rather than copying data. The indexing is, in notable difference to FFTViews.jl kept as ordinary arrays are indexed. This helps with the seamless integration across packages.
To implement an FFT-based resample operation of real-valued data, a new view, derived from AbstractArray is introduced, handling potential copy and addition operations for even-sized arrays to enforce the real-valuedness of the corresponding real space data (select_region_ft). In the community it has been discussed, whether such an operation is necessary. Referring to this discussion, we argue that the Fourier-space operations cannot be replaced by casting to real, since the latter violates Parseval's theorem.
In addition to the resample operation FourierTools.jl also provides a tool for sub-pixel shifting based FFTs. Further algorithms like shearing, sub-pixel shifting and rotation can be also implemented via the use of the Fourier shift theorem and due to the generality of the FFT these can be applied to N-dimensional datasets efficiently.

Felix Wechsler studied Physics and Informatics at the Technical University of Munich in Germany. For his master studies in Photonics he moved to Jena (city of light). Currently he finishes his master thesis in the field of computational microscopy at the Biomedical Imaging Group of Leibniz Institute of Photonic Technology under the joint supervision of Rainer Heintzmann and Ivo Ihrke.

This speaker also appears in:

My research focuses on developing tools to image cellular function at high resolution. We develop techniques to measure multidimensional information in small biological objects such as cells, cellular organelles or other small structures of interest. Computer-based reconstruction methods, especially inverse-modeling based in automatic differentiation are of primary interest.

This speaker also appears in: