# Kernel Convolution in Frequency Domain - Cyclic Padding

I don't know whether this is the right place to post this, but I suppose it is.

I know that frequency multiplication = circular convolution in time space for discrete signals (vectors).

I also know that "the convolution theorem yields the desired linear convolution result only if $$x(n)$$ and $$h(n)$$ are padded with zeros prior to the DFT such that their respective lengths are $$N_x+N_h-1$$, essentially zeroing out all circular artifacts."

and everything works with vectors.. but my goal is circular convolution with matrices as in this paper:

Victor Podlozhnyuk (nVidia) - FFT Based 2D Convolution

If you watch the first two figures (figure 1 and 2) you'll see that the kernel is padded in a weird way I've never seen before, what's this?

Figures 1 and 2 are not showing any padding whatsoever. The larger matrix is the data (probably image) matrix, not a padded kernel matrix. The figures are simply showing how the circular aspect of the convolution works in 2 dimensions.

• How can I find more information on how the circular convolution works with matrices? Jul 4 '12 at 7:08
• @JohnPell I would concentrate on really understanding one-dimensional convolution, and then the generalization to two dimensions is pretty simple. Starting with two-dimensions, though, I think makes understanding convolution more difficult. Pretty much any introduction to signal processing will talk about convolution. Jul 5 '12 at 2:59
• Regarding circular convolution there aren't many resources talking about it in 2D, you can find a lot of stuff about 1D though Jul 5 '12 at 6:10

In my StackExchange Signal Processing Q38542 GitHub Repository (See Applying Image Filtering (Circular Convolution) in Frequency Domain) you will be able to see a code which implements 2D Circular Convolution both in Spatial and Frequency Domain.

Pay attention to the function CircularExtension2D().
This function align the axis origin between the image and the kernel before working in the Frequency Domain.
Remember that for Discrete Signals the implicit assumption on signals, In frequency Domain analysis, is being periodic (Circular).
This is basically what's in FFT Based 2D Cyclic Convolution.

In the discrete case one could indeed apply Circular Convolution by element wise multiplication in the Frequency Domain.

With proper padding one could apply linear convolution using circular convolution hence Linear Convolution can also be achieved using multiplication in the Frequency Domain.

See: