Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Sense of zeropadding in a time domain

I have the task related to Radon transform which contains a subtask which uses resampling by means of DFT.

Let's consider the non-periodical discretized signal (Fig.1) (for example the string of pixels) having 515 pixels length. In my implementation for resampling contains following steps:

  1. Cyclic left shift (Fig.2).
  2. Add zeros to the center in order to the length of signal became 2^n (in our case 1024-515 = 509 zeros we must have add) (Fig. 3).
  3. Get DFT from this signal (Fig. 4).
  4. Cyclic right shift. (for shifting 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 we must perform cyclic shift of the signal and add zeros exactly in the center? (I assumed what this made the signal periodic) Zeropadding makes interpolation DFT spectrum, is it correct? (I asked and someone says what it is not quite so) Maybe someone can explain in simple way what happens with signal after zeropadding.

I have made some experiments in a Matlab and found out that any other sequence of actions can not give required result.

Now let's consider two cases:

a) (THIS CORRECT VARIANT) We has the non-periodical discretized signal (for example the string of pixels) which will be cyclic shifted to left and filled zeros in the center after that will be obtained DFT from this and to shift it back. enter image description here

b) We has the non-periodical discretized signal (for example the set strings of pixels) which will be filled zeros from left and right after that will be obtained DFT from this.

enter image description here

What the difference these DFT spectrums?

I have read some books but not found the answer of this zeropadding's case. It seems this can be found only by own experience.

Answer in book:

A. C. Kak and Malcolm Slaney, Principles of Computerized Tomographic Imaging, Society of Industrial and Applied Mathematics, 2001 on page 25

like image 877
Roman Shkarin Avatar asked Jun 10 '12 08:06

Roman Shkarin


1 Answers

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

Circular shifting in the time domain corresponds to a "phase twist" in the frequency domain; each bin has a complex rotation applied to it. I have no idea why you've been asked to do that in your application!

like image 192
Oliver Charlesworth Avatar answered Nov 10 '22 14:11

Oliver Charlesworth