Signals and Sampling 7.1 Sampling Theory 7.2 Image Sampling Interface

Total Page:16

File Type:pdf, Size:1020Kb

Load more

Chapter 7 Signals and Sampling 7.1 Sampling Theory 7.2 Image Sampling Interface Chapter 7 of “Physically Based 7.3 Stratified Sampling Rendering” by Pharr&Humphreys 7.4 Low-Discrepancy Sampling 7.5 Best-Candidate Sampling Patterns 7.6 Image Reconstruction 2 Additional Reading Motivation Chapter 14.10 of “CG: Principles & Practice” by • Real World - continuous Foley, van Dam et al. • Digital (Computer) world - discrete Chapter 4, 5, 8, 9, 10 in “Principles of Digital Image Synthesis,” by A. Glassner • Typically we have to either: Chapter 4, 5, 6 of “Digital Image Warping” by – create discrete data from continuous or (e.g. Wolberg rendering/ray-tracing, illumination models, Chapter 2, 4 of “Discrete-Time Signal Processing” by morphing) Oppenheim, Shafer – manipulate discrete data (textures, surface description, image processing,tone mapping) Motivation Motivation • Artifacts occurring in sampling - aliasing: Engineering approach: – Jaggies nearest neighbor: – Moire Convolution with – Flickering small objects box filter – Sparkling highlights – Temporal strobing linear filter: • Preventing these artifacts - Antialiasing Convolution with tent filter Motivation- Graphics Motivation Engineering approach: Original “manipulated” • black-box (continuous) signal “Graphics” (continuous) signal “System” or Algorithm • discretization: sampling Reconstruction sampled Multiplication with filter “shah” function signal Convolution Convolution (2) • How can we characterize our “black box”? • Impulse: • We assume to have a “nice” box/algorithm: – linear – time-invariant • discrete impulse: • then it can be characterized through the response to an “impulse”: • Finite Impulse Response (FIR) vs. “System” or • Infinite Impulse Response (IIR) Algorithm Convolution (3) Convolution (4) • for a time-invariant system h[k-n] would • An arbitrary signal x[k] can be written as: be the impulse response to a delayed impulse d[k-n] • Let the impulse response be h[k]: • hence, if y[k] is the response of our system to the input x[k] (and we assume a linear system): “System” or ![k] h[k] IIR - N=inf. Algorithm FIR - N<inf. “System” or x[k] y[k] Algorithm Fourier Transforms Fourier Transforms (2) • Let’s look at a special input sequence: • Hence is an eigen-function and H(") its eigenvalue • then: • H(") is the Fourier-Transform of the h[n] and hence characterizes the underlying system in terms of frequencies • H(") is periodic with period 2# • H(") is decomposed into – phase (angle) response – magnitude response Properties Properties (2) • Linear • Parseval’s Theorem • scaling • convolution • Multiplication • preserves “Energy” - overall signal content • Differentiation • delay/shift Transforms Pairs Transform Pairs - Shah • Sampling = Multiplication with a Shah Fourier function: Transform sampling Average T 1/T Filter • multiplication in spatial domain = convolution in the frequency domain Box/Sinc • frequency replica of primary spectrum Filter (also called aliased spectra) Transforms Pairs (2) General Process Original function Linear Sampled function Acquisition Filter Gaussian Reconstructed Filter Reconstruction Function Re-sampled function Resampling derivative Filter How? - Reconstruction Sampling Theorem • A signal can be reconstructed from its Spatial Domain: Frequency Domain: samples without loss of information if the original signal has no frequencies above 1/2 of the sampling frequency • For a given bandlimited function, the rate at Mathematically: which it must be sampled is called the • Convolution: f(x)*h(x) • Multiplication: Nyquist frequency Evaluated at discrete online demo points (sum) Example Example 2D 1D Given Given Needed Needed Nearest neighbor Linear Interpolation General Process - Pre-Filtering Frequency Domain Original function Sampled function Original function Band-limited function Acquisition Pre-Filtering Reconstructed Sampled Acquisition Reconstruction Function Function Re-sampled function Reconstructed function Resampling Reconstruction Once Again ... Pipeline - Example Spatial domain Frequency domain Pre-filter Post-aliasing x * Pre-aliasingsampling Reconstruction filter * x Pipeline - Example (2) Pipeline - Example (3) Spatial domain Frequency domain Spatial domain Frequency domain * x * x x * Sources of Aliasing Aliasing • Non-bandlimited signal sampling • Low sampling rate (below Nyquist) sampling • Non perfect reconstruction Bandlimited Interpolation Spatial Domain: Frequency Domain: • convolution is exact • cut off freq. replica Derivatives Reconstruction Kernels Spatial Domain: Frequency Domain: • Nearest Neighbor (Box) • convolution is exact • cut off freq. replica • Linear • Sinc • Gaussian • Many others Spatial d. Frequency d. Ideal Reconstruction Ideal Reconstruction • Box filter in frequency domain = • Use the sinc function – to bandlimit the • Sinc Filter in spatial domain sampled signal and remove all copies of the • impossible to realize (really?) spectra introduced by sampling Pass-band stop- • But: band – The sinc has infinite extent and we must use Ideal filter simpler filters with finite extents. Smoothing – The windowed versions of sinc may introduce Practical ringing artifacts which are perceptually filter Post-aliasing objectionable. Reconstructing with Sinc Ideal Reconstruction – Realizable filters do not have sharp transitions; also have ringing in pass/stop bands # # # Low-pass band-pass high-pass filter filter filter Higher Dimensions? Possible Errors • Design typically in 1D • Post-aliasing – reconstruction filter passes frequencies beyond the • extensions to higher dimensions (typically): Nyquist frequency (of duplicated frequency spectrum) – separable filters => frequency components of the original signal appear in the reconstructed signal at different frequencies – radially symmetric filters • Smoothing – limited results – frequencies below the Nyquist frequency are attenuated • research topic • Ringing (overshoot) T – occurs when trying to sample/reconstruct discontinuity • Anisotropy ? – caused by not spherically symmetric filters Aliasing vs. Noise Antialiasing • Antialiasing = Preventing aliasing • 1. Analytically pre-filter the signal – Solvable for points, lines and polygons – Not solvable in general (e.g. procedurally defined images) • 2. Uniform supersampling and resample • 3. Nonuniform or stochastic sampling Distribution of Extrafoveal Uniform Supersampling Cones • Increasing the sampling rate moves each • Yellot theory (1983) copy of the spectra further apart, potentially – Aliases replaced by noise reducing the overlap and thus aliasing – Visual system less sensitive to high freq noise • Resulting samples must be resampled Monkey eye cone distribution Fourier Transform (filtered) to image sampling rate Non-Uniform Sampling - Non-Uniform Sampling - Intuition Patterns • Uniform sampling • Poisson – The spectrum of uniformly spaced samples is also a set – Pick n random points in sample space of uniformly spaced spikes – Multiplying the signal by the sampling pattern • Uniform Jitter corresponds to placing a copy of the spectrum at each – Subdivide sample space into n regions spike (in freq. space) – Aliases are coherent, and very noticeable • Poisson Disk • Non-uniform sampling – Pick n random points, but not too close – Samples at non-uniform locations have a different spectrum; a single spike plus noise – Sampling a signal in this way converts aliases into broadband noise – Noise is incoherent, and much less objectionable Poisson Disk Sampling Uniform Jittered Sampling Spatial Domain Fourier Domain Spatial Domain Fourier Domain Non-Uniform Sampling - Stratified Sampling Patterns • Spectral characteristics of these • Put at least one sample in each strata distributions: • Multiple samples in strata do no good – Poisson: completely uniform (white noise). • Also have samples far away from each High and low frequencies equally present other – Poisson disc: Pulse at origin (DC component of image), surrounded by empty ring (no low frequencies), surrounded by white noise • Graphics: jittering – Jitter: Approximates Poisson disc spectrum, but with a smaller empty disc. Stratification Stratification • OR – Split up the integration domain in N disjoint sub-domains or strata – Evaluate the integral in each of the sub- domains separately with one or more samples. • More precisely: More Jittered Sequences Jitter • Place samples in the grid • Perturb the samples up to 1/2 width or height Texture Example Multiple Dimensions • Too many samples Exact – 256 samples/pixel Jitter with 1 sample/pixel • 1D • 2D 3D 1 sample/pixel Jitter with 4 samples/pixel Jitter Problems Multiple Dimensions • How to deal with higher dimensions? • Make (separate) strata for each dimension – Curse of dimensionality • Randomly associate strata among each – D dimensions means ND “cells” (if we use a other separable extension) • Ensure good sample “distribution” • Solutions: – Example: 2D screen position; 2D lens position; – We can look at each dimension independently 1D time – We can either look in non-separable geometries – Latin Hypercube (or N-Rook) sampling Optimal sampling lattices Optimal sampling lattices • Dividing space up into equal cells doesn’t • We have to deal with different geometry have to be on a Cartesian lattices • 2D - hexagon • In fact - Cartesian is NOT the optimal way • 3D - truncated octahedron how to divide up space uniformly Cartesian Hexagonal Latin Hypercubes - N-Rooks Stratification - problems • Distributing n samples in D dimensions, • Clamping (LHS helps) } even if n is not a power of D • Could still have large
Recommended publications
  • Convolution! (CDT-14) Luciano Da Fontoura Costa

    Convolution! (CDT-14) Luciano Da Fontoura Costa

    Convolution! (CDT-14) Luciano da Fontoura Costa To cite this version: Luciano da Fontoura Costa. Convolution! (CDT-14). 2019. hal-02334910 HAL Id: hal-02334910 https://hal.archives-ouvertes.fr/hal-02334910 Preprint submitted on 27 Oct 2019 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Convolution! (CDT-14) Luciano da Fontoura Costa [email protected] S~aoCarlos Institute of Physics { DFCM/USP October 22, 2019 Abstract The convolution between two functions, yielding a third function, is a particularly important concept in several areas including physics, engineering, statistics, and mathematics, to name but a few. Yet, it is not often so easy to be conceptually understood, as a consequence of its seemingly intricate definition. In this text, we develop a conceptual framework aimed at hopefully providing a more complete and integrated conceptual understanding of this important operation. In particular, we adopt an alternative graphical interpretation in the time domain, allowing the shift implied in the convolution to proceed over free variable instead of the additive inverse of this variable. In addition, we discuss two possible conceptual interpretations of the convolution of two functions as: (i) the `blending' of these functions, and (ii) as a quantification of `matching' between those functions.
  • Eliminating Aliasing Caused by Discontinuities Using Integrals of the Sinc Function

    Eliminating Aliasing Caused by Discontinuities Using Integrals of the Sinc Function

    Sound production – Sound synthesis: Paper ISMR2016-48 Eliminating aliasing caused by discontinuities using integrals of the sinc function Fabián Esqueda(a), Stefan Bilbao(b), Vesa Välimäki(c) (a)Aalto University, Dept. of Signal Processing and Acoustics, Espoo, Finland, fabian.esqueda@aalto.fi (b)Acoustics and Audio Group, University of Edinburgh, United Kingdom, [email protected] (c)Aalto University, Dept. of Signal Processing and Acoustics, Espoo, Finland, vesa.valimaki@aalto.fi Abstract A study on the limits of bandlimited correction functions used to eliminate aliasing in audio signals with discontinuities is presented. Trivial sampling of signals with discontinuities in their waveform or their derivatives causes high levels of aliasing distortion due to the infinite bandwidth of these discontinuities. Geometrical oscillator waveforms used in subtractive synthesis are a common example of audio signals with these characteristics. However, discontinuities may also be in- troduced in arbitrary signals during operations such as signal clipping and rectification. Several existing techniques aim to increase the perceived quality of oscillators by attenuating aliasing suf- ficiently to be inaudible. One family of these techniques consists on using the bandlimited step (BLEP) and ramp (BLAMP) functions to quasi-bandlimit discontinuities. Recent work on antialias- ing clipped audio signals has demonstrated the suitability of the BLAMP method in this context. This work evaluates the performance of the BLEP, BLAMP, and integrated BLAMP functions by testing whether they can be used to fully bandlimit aliased signals. Of particular interest are cases where discontinuities appear past the first derivative of a signal, like in hard clipping. These cases require more than one correction function to be applied at every discontinuity.
  • Automatic Detection of Perceived Ringing Regions in Compressed Images

    Automatic Detection of Perceived Ringing Regions in Compressed Images

    International Journal of Electronics and Communication Engineering. ISSN 0974-2166 Volume 4, Number 5 (2011), pp. 491-516 © International Research Publication House http://www.irphouse.com Automatic Detection of Perceived Ringing Regions in Compressed Images D. Minola Davids Research Scholar, Singhania University, Rajasthan, India Abstract An efficient approach toward a no-reference ringing metric intrinsically exists of two steps: first detecting regions in an image where ringing might occur, and second quantifying the ringing annoyance in these regions. This paper presents a novel approach toward the first step: the automatic detection of regions visually impaired by ringing artifacts in compressed images. It is a no- reference approach, taking into account the specific physical structure of ringing artifacts combined with properties of the human visual system (HVS). To maintain low complexity for real-time applications, the proposed approach adopts a perceptually relevant edge detector to capture regions in the image susceptible to ringing, and a simple yet efficient model of visual masking to determine ringing visibility. Index Terms: Luminance masking, per-ceptual edge, ringing artifact, texture masking. Introduction In current visual communication systems, the most essential task is to fit a large amount of visual information into the narrow bandwidth of transmission channels or into a limited storage space, while maintaining the best possible perceived quality for the viewer. A variety of compression algorithms, for example, such as JPEG[1] and MPEG/H.26xhave been widely adopted in image and video coding trying to achieve high compression efficiency at high quality. These lossy compression techniques, however, inevitably result in various coding artifacts, which by now are known and classified as blockiness, ringing, blur, etc.
  • Generalizing Sampling Theory for Time-Varying Nyquist Rates Using Self-Adjoint Extensions of Symmetric Operators with Deficiency Indices (1,1) in Hilbert Spaces

    Generalizing Sampling Theory for Time-Varying Nyquist Rates Using Self-Adjoint Extensions of Symmetric Operators with Deficiency Indices (1,1) in Hilbert Spaces

    Generalizing Sampling Theory for Time-Varying Nyquist Rates using Self-Adjoint Extensions of Symmetric Operators with Deficiency Indices (1,1) in Hilbert Spaces by Yufang Hao A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in Applied Mathematics Waterloo, Ontario, Canada, 2011 c Yufang Hao 2011 I hereby declare that I am the sole author of this thesis. This is a true copy of the thesis, including any required final revisions, as accepted by my examiners. I understand that my thesis may be made electronically available to the public. ii Abstract Sampling theory studies the equivalence between continuous and discrete representa- tions of information. This equivalence is ubiquitously used in communication engineering and signal processing. For example, it allows engineers to store continuous signals as discrete data on digital media. The classical sampling theorem, also known as the theorem of Whittaker-Shannon- Kotel'nikov, enables one to perfectly and stably reconstruct continuous signals with a con- stant bandwidth from their discrete samples at a constant Nyquist rate. The Nyquist rate depends on the bandwidth of the signals, namely, the frequency upper bound. Intuitively, a signal's `information density' and ‘effective bandwidth' should vary in time. Adjusting the sampling rate accordingly should improve the sampling efficiency and information storage. While this old idea has been pursued in numerous publications, fundamental problems have remained: How can a reliable concept of time-varying bandwidth been defined? How can samples taken at a time-varying Nyquist rate lead to perfect and stable reconstruction of the continuous signals? This thesis develops a new non-Fourier generalized sampling theory which takes samples only as often as necessary at a time-varying Nyquist rate and maintains the ability to perfectly reconstruct the signals.
  • Discrete - Time Signals and Systems

    Discrete - Time Signals and Systems

    Discrete - Time Signals and Systems Sampling – II Sampling theorem & Reconstruction Yogananda Isukapalli Sampling at diffe- -rent rates From these figures, it can be concluded that it is very important to sample the signal adequately to avoid problems in reconstruction, which leads us to Shannon’s sampling theorem 2 Fig:7.1 Claude Shannon: The man who started the digital revolution Shannon arrived at the revolutionary idea of digital representation by sampling the information source at an appropriate rate, and converting the samples to a bit stream Before Shannon, it was commonly believed that the only way of achieving arbitrarily small probability of error in a communication channel was to 1916-2001 reduce the transmission rate to zero. All this changed in 1948 with the publication of “A Mathematical Theory of Communication”—Shannon’s landmark work Shannon’s Sampling theorem A continuous signal xt( ) with frequencies no higher than fmax can be reconstructed exactly from its samples xn[ ]= xn [Ts ], if the samples are taken at a rate ffs ³ 2,max where fTss= 1 This simple theorem is one of the theoretical Pillars of digital communications, control and signal processing Shannon’s Sampling theorem, • States that reconstruction from the samples is possible, but it doesn’t specify any algorithm for reconstruction • It gives a minimum sampling rate that is dependent only on the frequency content of the continuous signal x(t) • The minimum sampling rate of 2fmax is called the “Nyquist rate” Example1: Sampling theorem-Nyquist rate x( t )= 2cos(20p t ), find the Nyquist frequency ? xt( )= 2cos(2p (10) t ) The only frequency in the continuous- time signal is 10 Hz \ fHzmax =10 Nyquist sampling rate Sampling rate, ffsnyq ==2max 20 Hz Continuous-time sinusoid of frequency 10Hz Fig:7.2 Sampled at Nyquist rate, so, the theorem states that 2 samples are enough per period.
  • The Best Approximation of the Sinc Function by a Polynomial of Degree N with the Square Norm

    The Best Approximation of the Sinc Function by a Polynomial of Degree N with the Square Norm

    Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 2010, Article ID 307892, 12 pages doi:10.1155/2010/307892 Research Article The Best Approximation of the Sinc Function by a Polynomial of Degree n with the Square Norm Yuyang Qiu and Ling Zhu College of Statistics and Mathematics, Zhejiang Gongshang University, Hangzhou 310018, China Correspondence should be addressed to Yuyang Qiu, [email protected] Received 9 April 2010; Accepted 31 August 2010 Academic Editor: Wing-Sum Cheung Copyright q 2010 Y. Qiu and L. Zhu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The polynomial of degree n which is the best approximation of the sinc function on the interval 0, π/2 with the square norm is considered. By using Lagrange’s method of multipliers, we construct the polynomial explicitly. This method is also generalized to the continuous function on the closed interval a, b. Numerical examples are given to show the effectiveness. 1. Introduction Let sin cxsin x/x be the sinc function; the following result is known as Jordan inequality 1: 2 π ≤ sin cx < 1, 0 <x≤ , 1.1 π 2 where the left-handed equality holds if and only if x π/2. This inequality has been further refined by many scholars in the past few years 2–30. Ozban¨ 12 presented a new lower bound for the sinc function and obtained the following inequality: 2 1 4π − 3 π 2 π2 − 4x2 x − ≤ sin cx.
  • Oversampling of Stochastic Processes

    Oversampling of Stochastic Processes

    OVERSAMPLING OF STOCHASTIC PROCESSES by D.S.G. POLLOCK University of Leicester In the theory of stochastic differential equations, it is commonly assumed that the forcing function is a Wiener process. Such a process has an infinite bandwidth in the frequency domain. In practice, however, all stochastic processes have a limited bandwidth. A theory of band-limited linear stochastic processes is described that reflects this reality, and it is shown how the corresponding ARMA models can be estimated. By ignoring the limitation on the frequencies of the forcing function, in the process of fitting a conventional ARMA model, one is liable to derive estimates that are severely biased. If the data are sampled too rapidly, then maximum frequency in the sampled data will be less than the Nyquist value. However, the underlying continuous function can be reconstituted by sinc function or Fourier interpolation; and it can be resampled at a lesser rate corresponding to the maximum frequency of the forcing function. Then, there will be a direct correspondence between the parameters of the band-limited ARMA model and those of an equivalent continuous-time process; and the estimation biases can be avoided. 1 POLLOCK: Band-Limited Processes 1. Time-Limited versus Band-Limited Processes Stochastic processes in continuous time are usually modelled by filtered versions of Wiener processes which have infinite bandwidth. This seems inappropriate for modelling the slowly evolving trajectories of macroeconomic data. Therefore, we shall model these as processes that are limited in frequency. A function cannot be simultaneously limited in frequency and limited in time. One must choose either a band-limited function, which extends infinitely in time, or a time-limited function, which extents over an infinite range of frequencies.
  • CHAPTER 3 ADC and DAC

    CHAPTER 3 ADC and DAC

    CHAPTER 3 ADC and DAC Most of the signals directly encountered in science and engineering are continuous: light intensity that changes with distance; voltage that varies over time; a chemical reaction rate that depends on temperature, etc. Analog-to-Digital Conversion (ADC) and Digital-to-Analog Conversion (DAC) are the processes that allow digital computers to interact with these everyday signals. Digital information is different from its continuous counterpart in two important respects: it is sampled, and it is quantized. Both of these restrict how much information a digital signal can contain. This chapter is about information management: understanding what information you need to retain, and what information you can afford to lose. In turn, this dictates the selection of the sampling frequency, number of bits, and type of analog filtering needed for converting between the analog and digital realms. Quantization First, a bit of trivia. As you know, it is a digital computer, not a digit computer. The information processed is called digital data, not digit data. Why then, is analog-to-digital conversion generally called: digitize and digitization, rather than digitalize and digitalization? The answer is nothing you would expect. When electronics got around to inventing digital techniques, the preferred names had already been snatched up by the medical community nearly a century before. Digitalize and digitalization mean to administer the heart stimulant digitalis. Figure 3-1 shows the electronic waveforms of a typical analog-to-digital conversion. Figure (a) is the analog signal to be digitized. As shown by the labels on the graph, this signal is a voltage that varies over time.
  • Analysis of Ringing Artifact in Image Fusion Using Directional Wavelet

    Analysis of Ringing Artifact in Image Fusion Using Directional Wavelet

    Special Issue - 2021 International Journal of Engineering Research & Technology (IJERT) ISSN: 2278-0181 NTASU - 2020 Conference Proceedings Analysis of Ringing Artifact in Image Fusion Using Directional Wavelet Transforms y z x Ashish V. Vanmali Tushar Kataria , Samrudhha G. Kelkar , Vikram M. Gadre Dept. of Information Technology Dept. of Electrical Engineering Vidyavardhini’s C.O.E. & Tech. Indian Institute of Technology, Bombay Vasai, Mumbai, India – 401202 Powai, Mumbai, India – 400076 Abstract—In the field of multi-data analysis and fusion, image • Images taken from multiple sensors. Examples include fusion plays a vital role for many applications. With inventions near infrared (NIR) images, IR images, CT, MRI, PET, of new sensors, the demand of high quality image fusion fMRI etc. algorithms has seen tremendous growth. Wavelet based fusion is a popular choice for many image fusion algorithms, because of its We can broadly classify the image fusion techniques into four ability to decouple different features of information. However, categories: it suffers from ringing artifacts generated in the output. This 1) Component substitution based fusion algorithms [1]–[5] paper presents an analysis of ringing artifacts in application of 2) Optimization based fusion algorithms [6]–[10] image fusion using directional wavelets (curvelets, contourlets, non-subsampled contourlets etc.). We compare the performance 3) Multi-resolution (wavelets and others) based fusion algo- of various fusion rules for directional wavelets available in rithms [11]–[15] and literature. The experimental results suggest that the ringing 4) Neural network based fusion algorithms [16]–[19]. artifacts are present in all types of wavelets with the extent of Wavelet based multi-resolution analysis decouples data into artifact varying with type of the wavelet, fusion rule used and levels of decomposition.
  • Evaluating Oscilloscope Sample Rates Vs. Sampling Fidelity: How to Make the Most Accurate Digital Measurements

    Evaluating Oscilloscope Sample Rates Vs. Sampling Fidelity: How to Make the Most Accurate Digital Measurements

    AC 2011-2914: EVALUATING OSCILLOSCOPE SAMPLE RATES VS. SAM- PLING FIDELITY Johnnie Lynn Hancock, Agilent Technologies About the Author Johnnie Hancock is a Product Manager at Agilent Technologies Digital Test Division. He began his career with Hewlett-Packard in 1979 as an embedded hardware designer, and holds a patent for digital oscillo- scope amplifier calibration. Johnnie is currently responsible for worldwide application support activities that promote Agilent’s digitizing oscilloscopes and he regularly speaks at technical conferences world- wide. Johnnie graduated from the University of South Florida with a degree in electrical engineering. In his spare time, he enjoys spending time with his four grandchildren and restoring his century-old Victorian home located in Colorado Springs. Contact Information: Johnnie Hancock Agilent Technologies 1900 Garden of the Gods Rd Colorado Springs, CO 80907 USA +1 (719) 590-3183 johnnie [email protected] c American Society for Engineering Education, 2011 Evaluating Oscilloscope Sample Rates vs. Sampling Fidelity: How to Make the Most Accurate Digital Measurements Introduction Digital storage oscilloscopes (DSO) are the primary tools used today by digital designers to perform signal integrity measurements such as setup/hold times, rise/fall times, and eye margin tests. High performance oscilloscopes are also widely used in university research labs to accurately characterize high-speed digital devices and systems, as well as to perform high energy physics experiments such as pulsed laser testing. In addition, general-purpose oscilloscopes are used extensively by Electrical Engineering students in their various EE analog and digital circuits lab courses. The two key banner specifications that affect an oscilloscope’s signal integrity measurement accuracy are bandwidth and sample rate.
  • Enhancing ADC Resolution by Oversampling

    Enhancing ADC Resolution by Oversampling

    AVR121: Enhancing ADC resolution by oversampling 8-bit Features Microcontrollers • Increasing the resolution by oversampling • Averaging and decimation • Noise reduction by averaging samples Application Note 1 Introduction Atmel’s AVR controller offers an Analog to Digital Converter with 10-bit resolution. In most cases 10-bit resolution is sufficient, but in some cases higher accuracy is desired. Special signal processing techniques can be used to improve the resolution of the measurement. By using a method called ‘Oversampling and Decimation’ higher resolution might be achieved, without using an external ADC. This Application Note explains the method, and which conditions need to be fulfilled to make this method work properly. Figure 1-1. Enhancing the resolution. A/D A/D A/D 10-bit 11-bit 12-bit t t t Rev. 8003A-AVR-09/05 2 Theory of operation Before reading the rest of this Application Note, the reader is encouraged to read Application Note AVR120 - ‘Calibration of the ADC’, and the ADC section in the AVR datasheet. The following examples and numbers are calculated for Single Ended Input in a Free Running Mode. ADC Noise Reduction Mode is not used. This method is also valid in the other modes, though the numbers in the following examples will be different. The ADCs reference voltage and the ADCs resolution define the ADC step size. The ADC’s reference voltage, VREF, may be selected to AVCC, an internal 2.56V / 1.1V reference, or a reference voltage at the AREF pin. A lower VREF provides a higher voltage precision but minimizes the dynamic range of the input signal.
  • Fourier Analysis and Sampling Theory

    Fourier Analysis and Sampling Theory

    Reading Required: Shirley, Ch. 9 Recommended: Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Fourier analysis and Don P. Mitchell and Arun N. Netravali, “Reconstruction Filters in Computer Computer sampling theory Graphics ,” Computer Graphics, (Proceedings of SIGGRAPH 88). 22 (4), pp. 221-228, 1988. Brian Curless CSE 557 Fall 2009 1 2 What is an image? Images as functions We can think of an image as a function, f, from R2 to R: f(x,y) gives the intensity of a channel at position (x,y) Realistically, we expect the image only to be defined over a rectangle, with a finite range: • f: [a,b]x[c,d] Æ [0,1] A color image is just three functions pasted together. We can write this as a “vector-valued” function: ⎡⎤rxy(, ) f (,xy )= ⎢⎥ gxy (, ) ⎢⎥ ⎣⎦⎢⎥bxy(, ) We’ll focus in grayscale (scalar-valued) images for now. 3 4 Digital images Motivation: filtering and resizing In computer graphics, we usually create or operate What if we now want to: on digital (discrete)images: smooth an image? Sample the space on a regular grid sharpen an image? Quantize each sample (round to nearest enlarge an image? integer) shrink an image? If our samples are Δ apart, we can write this as: In this lecture, we will explore the mathematical underpinnings of these operations. f [n ,m] = Quantize{ f (n Δ, m Δ) } 5 6 Convolution Convolution in 2D One of the most common methods for filtering a In two dimensions, convolution becomes: function, e.g., for smoothing or sharpening, is called convolution.