1 / 31

Sampling and Aliasing

Sampling and Aliasing. Gilad Lerman Math 5467 (stealing slides from Gonzalez & Woods, and Efros). The Sampling Theorem. Theorem: If f is in L 1 (  ) & supported on [- B 0 , B 0 ], then Recall Proof: We view as (2 B 0 )-periodic function with coefficients:

arnoldo
Download Presentation

Sampling and Aliasing

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Sampling and Aliasing Gilad Lerman Math 5467 (stealing slides from Gonzalez & Woods, and Efros)

  2. The Sampling Theorem Theorem: If f is in L1() & supported on [-B0, B0], then Recall Proof: We view as (2B0)-periodic function with coefficients: At last, find f using IFT and using FS of

  3. More on the Sampling Theorem Frequency band: Time: Note: Theorem holds for B>B0. Indeed, then If B<B0, the above equation is not true for all 

  4. Sampling Theorem (meaning) • Interpretation: If a function f(t) contains no frequencies higher than W cps, it is completely determined by giving its ordinates at a series of points spaced 1/(2W) seconds apart • Remark: For L1 function a freq. = W is fine but for more general functions we need > W…

  5. Simple Example (not L1) Assume a cosine (it is not L1() but will be instrumental) Freq: a (“& -a”), Freq Band: =[-a,a], Time: 1/(2a) Here one needs B>B0 (B=B0 doesn’t work) Example: for all 3 functions freq: 0.5, time: 1 The sampled function has different aliases…

  6. Aliasing • If the sampling condition is not satisfied, frequencies will overlap (high freq → low freq) • The reconstructed signal is said to be an alias of the original signal

  7. Input signal: Matlab output: x = 0:.05:5; imagesc(sin((2.^x).*x)) Example: Increased Frequency Related Image: Picket fence receding Into the distance will produce aliasing…

  8. One more example at the Fourier domain

  9. Aliasing in Images (Fourier domain)

  10. Good and Bad Sampling • Good sampling: • Sample often or, • Sample wisely • Bad sampling: • see aliasing in action!

  11. Texture makes its worse(high frequencies)

  12. Even worse for synthetic images Slide by Steve Seitz

  13. Really bad in video Slide by Paul Heckbert

  14. Wheels of Wagons in Westerns

  15. Moiré pattern • Definition: Interference pattern created, e.g., when two grids are overlaid at an angle, or when they have slightly different mesh sizes. • In images produced e.g., when scanning a halftone picture or due to undersampling a fine regular pattern.

  16. Moiré pattern due to undersampling Original image downsampled image

  17. Antialiasing • What can be done? Sampling rate ≥ 2 * max frequency in the image • Raise sampling rate by oversampling • Sample at k times the resolution • continuous signal: easy • discrete signal: need to interpolate • 2. Lower the max frequency by prefiltering • Smooth the signal enough • Works on discrete signals • 3. Improve sampling quality with better sampling • Nyquist is best case! • Stratified sampling • Importance sampling • Relies on domain knowledge

  18. Gaussian pre-filtering G 1/8 G 1/4 Gaussian 1/2 • Solution: filter the image, then subsample • Filter size should double for each ½ size reduction.

  19. Subsampling with Gaussian pre-filtering Gaussian 1/2 G 1/4 G 1/8

  20. Compare with... 1/2 1/4 (2x zoom) 1/8 (4x zoom)

  21. Correcting some Moiré patterns

  22. Rethinking of the Cooley-Tukey FFT • Step 1(top to bottom): Create two subsampled signals (even and odd coordinates) • Note that the two subsampled signals are associated with half bands in the frequency domain (Shannon) • Step 2 (bottom-up): Combine the two signals by the formulas: • Interpretation: combining the two half bands in the right way (in frequency domain) to exactly recover the signal

More Related