Frequency Domain Interpolation: Convolution with Sinc Function
The "discrete sinc kernel" is also known as a Dirichlet Kernel.  It looks like a you got a 2D (u and v) going there.  I'm not that familiar with that, but we just did a pretty good number on the 1D case.

This post explains the inverse DFT can be interpreted as "epicycle drawing":

How to get Fourier coefficients to draw any shape using DFT?

This post then follows up with the relationship between the sampling, epicycle drawing, upsampling, and the normalized sinc function:

Absolute convergence of periodic sinc interpolation

The conclusions are:

Suppose a continuous band limited function (with band limit $L$) is sampled, then reconstructed from the samples.

1) If L &lt; N/2 the function can be unambiguously reconstructed.  The results from an infinite series sinc reconstruction (Whittaker-Shannon) and taking a DFT, zero padding at the Nyquist and then the inverse DFT to upsample, will produce the same path.

2) If L = N/2, the Nyquist bin must be split between N/2 and -N/2 for the sinc reconstruction to match the DFT reconstruction.  It may or not may not match the original function, but it will match at all the sample points.

3) If L > N/2, and N is even, the Nyquist bin must still be split for the sinc reconstruction to match the DFT reconstruction, but neither will match the original function, yet both will still match at all the sample points.

The relationship between the sinc function and the Direchlet Kernel is this:

1) The sinc function is the limit of the Dirichlet kernel as the sample count goes to infinity

2) For odd N, the Dirichlet kernel is an infinite sum of sinc functions.  For even N, it is an adjusted one.

See the posts for details and discussion.



What I forgot to overtly mention is that the Dirichlet Kernel formula is directly derived from, and equivalent to, doing the FFT, zero pad at Nyquist (split Nyquist on even N) then inverse DFT for resampling or coefficients for Fourier Series.  Your link in the comment is about the reverse case, but since the forward DFT and inverse DFT are mirror image twins (especially when $1/\sqrt{N}$ normalized), so what gets said in one direction has to be true in the other.  This means if you zero pad a signal, and have an even number of points, one of the outside points needs to be halved and the other half appended at the other end before zero padding outward, leaving your data centered.

I already explained the relationship between the Sinc and the Dirichlet, just above.  The math can be found in the second link in the section of my answer that has $y_m = Y(w)$ in it.  Towards the bottom.