The meaning of numbering in the time domain

I have a task related to the Radon transform, which contains a subtask that uses re-sampling using DFT.

Consider a non-periodic sampled signal (Fig. 1) (for example, a row of pixels) having a length of 515 pixels. My implementation for re-fetching contains the following steps:

  • A cyclic shift to the left (Fig. 2).
  • Add zeros to the center so that the signal length becomes 2 ^ n (in our case, 1024-515 = 509 zeros, which we must add) (Fig. 3).
  • Get the DFT from this signal (Figure 4).
  • Cyclical shift to the right. (to shift the low frequencies to the center) (Fig. 5)

Fig. 1 Original image

Fig. 2 Cyclic left shift

Fig. 3 Zeropadded

Fig. 4 DFT Spectrum

fig. 5 DFT back shifted

The main question:

Why do we need to cycle the signal and add zeros exactly to the center? (I assumed that this made the signal periodic) Zeropadding does the interpolation DFT spectrum, is it correct? (I asked, and someone says that this is not entirely true) Maybe someone can just explain what happens to the signal after the zero update.

I conducted several experiments in Matlab and found that any other sequence of actions cannot give the desired result.

Now consider two cases:

a) (THIS IS THE RIGHT OPTION) We have a non-periodic discretized signal (for example, a row of pixels) that will be cyclically shifted to the left and filled with zeros in the center, after which we will receive the DFT from this and shift it back. enter image description here

b) We have a non-periodic sampled signal (for example, given pixel strings) that will be filled with zeros on the left and right after the DFT is received after that.

enter image description here

What is the difference between these DFT spectra?

I read several books, but did not find an answer to this case for numbering. It seems that this can only be found from personal experience.

The answer in the book:

a. C. Kak and Malcolm Slaney, Principles of Computed Tomographic Imaging, Society for Industrial and Applied Mathematics, 2001 on p. 25

+6
source share
2 answers

Zero-padding in the time domain corresponds to interpolation in the frequency domain.

The circular displacement in the time domain corresponds to a “phase rotation” in the frequency domain; each bit has a complex rotation applied to it. I have no idea why you were asked to do this in your application!

+5
source

The offset of the data points and zero fill the exact center of the FFT aperture has the property that all even (symmetrical) signals in the original data window end in the real component of the complex FFT result and all odd signals end in the imaginary component. for example, the parity to odd ratio is maintained, which allows the phase to be interpolated. The ability to interpolate the phase is important in the case of an FFT with zero fill, since zero fill also interpolates the spectral value.

If you do not center the zero padding, the phase must be “spun” in the frequency domain before any additional interpolation can produce reasonable phase results.

+1
source

Source: https://habr.com/ru/post/917712/


All Articles