Aliasing and Antialiasing. What is Aliasing? What is Aliasing? What is Aliasing?

Size: px
Start display at page:

Download "Aliasing and Antialiasing. What is Aliasing? What is Aliasing? What is Aliasing?"

Transcription

1 What is Aliasing? Errors and Artifacts arising during rendering, due to the conversion from a continuously defined illumination field to a discrete raster grid of pixels 1 2 What is Aliasing? What is Aliasing? 3 4

2 What is Aliasing? Effects of Aliasing 5 6 Effects of Aliasing Effects of Aliasing 7 8

3 Effects of Aliasing Anti-aliasing 9 10 Area Sampling Techniques Anti-aliasing Techniques Prefiltering (unweighted/weighted area sampling) Postfiltering (supersampling, jittering) 11 12

4 Area Sampling Techniques 13 Area Sampling Techniques Area Sampling Techniques Area Sampling Techniques 16

5 Unweighted Area Sampling Pixel intensity is varied in proportion to the area of the pixel intercepted by the primitive. Unweighted equivalent to a box filter of unit height over pixel. Weighted Area Sampling Properties Intensity of pixel decreases as the distance between the pixel center and primitive increases. A primitive cannot influence a pixel s intensity if it does not intersect it. Equal areas (intersected) contribute equal intensity not a desirable property. 17 Equal areas can contribute unequally in terms of pixel intensity. Areas closer to the pixel center contribute more. Essentially results in filtering with a mask that is centered over the pixel with decreasing radial influence. Cone filters are a compromise between computational expense and optimality. 18 Postfiltering Techniques Supersampling (Regular Sampling) Very expensive. Not very satisfactory

6 Regular vs. Jittered Sampling 21 Filtering 23 Filtering Example Filtering 22 24

7 Filtering Example 25 Filtering Example Filtering Example 26 Aliasing from a Sampling Theory Viewpoint Sampling(Spatial Domain) 27 28

8 Frequency Domain Image is a spatial signal Sampling(Spatial Domain) X axis (position): frequency Y axis (height): strength of each frequency Examples: sine wave: impulse, square wave: infinite train of impulses How do we get to the Frequency Domain? What does the Fourier Transform Do to A Spatial Signal? Use the Fourier Transform Let φ(x) be a continuous function of a real variable x. Then I{φ(x)} = φ(ω) = φ(x)e j2πωx dx is the Fourier Transform of φ(x), with j = 1 and, I 1 {φ(ω)} = φ(x) = is the Inverse Fourier Transform. φ(x) is continuous and integrable φ(ω) is integrable x (spatial domain), ω (frequency domain) φ(ω)e j2πωx dω 31 ITCS 4120/5120 Signal in frequency domain is an 32 integration of individual Aliasingsinusoids. and Antialiasing

9 How does this related to Graphics? Sampling Theorem Continuous-time signal can be completely recovered from its samples iff the sampling rate is greater than twice the maximum frequency present in the signal. Claude Shannon Also known as the Nyquist rate Images are just a 2D signal and jagged edges are due to the pixel sampling rate not being high enough to capture the real signal Nyquist Rate Nyquist Rate:Undersampling The lower signal is undersampled and results in an aliased wave (dotted curve)

10 Comb Function Comb Function(contd) Application: Used to digitize continuous functions. Series of impulses (delta functions) Identity element of convolution: reproduces an indentical copy of the function f(x) FT of a comb function is another comb function 37 Multiplying f(x) with a comb in image space convolving their Fourier transforms, resulting in multiple identical copies of I{f(x)} Can result in aliasing if copies overlap Maximum allowable frequency is the Nyquist Frequency, which is half the sampling frequency. 38 Reconstruction Example(Adequate Sampling) Reconstruction Example(Inadequate Sampling) 39 40

11 Box Filter Pyramid Filter Reconstruction filter for nearest neighbor interpolation. Resampling images/volumes to a higher resolution using nearest neighbor values. FT of a box filter is the Sinc function ( sinπx πx ) Large side lobes continuing at regular intervals will cause aliasing. Aliasing in images manifests itself as jaggies 41 Reconstruction filter used in linear interpolation Computationally more expensive, but more accurate FT is much better behaved (side lobes much smaller) Less tendency to produce aliasing 42 Gaussian Filter The optimal filter in terms of avodiding side lobes FT of a Gaussian is another Gaussian Widely used to blur images and the basis for scale space 43

Filters. Materials from Prof. Klaus Mueller

Filters. Materials from Prof. Klaus Mueller Filters Materials from Prof. Klaus Mueller Think More about Pixels What exactly a pixel is in an image or on the screen? Solid square? This cannot be implemented A dot? Yes, but size matters Pixel Dots

More information

Antialiasing and Related Issues

Antialiasing and Related Issues Antialiasing and Related Issues OUTLINE: Antialiasing Prefiltering, Supersampling, Stochastic Sampling Rastering and Reconstruction Gamma Correction Antialiasing Methods To reduce aliasing, either: 1.

More information

Sampling and Reconstruction

Sampling and Reconstruction Sampling and Reconstruction Peter Rautek, Eduard Gröller, Thomas Theußl Institute of Computer Graphics and Algorithms Vienna University of Technology Motivation Theory and practice of sampling and reconstruction

More information

Image Processing. Image Processing. What is an Image? Image Resolution. Overview. Sources of Error. Filtering Blur Detect edges

Image Processing. Image Processing. What is an Image? Image Resolution. Overview. Sources of Error. Filtering Blur Detect edges Thomas Funkhouser Princeton University COS 46, Spring 004 Quantization Random dither Ordered dither Floyd-Steinberg dither Pixel operations Add random noise Add luminance Add contrast Add saturation ing

More information

Sampling Theory. CS5625 Lecture Steve Marschner. Cornell CS5625 Spring 2016 Lecture 7

Sampling Theory. CS5625 Lecture Steve Marschner. Cornell CS5625 Spring 2016 Lecture 7 Sampling Theory CS5625 Lecture 7 Sampling example (reminder) When we sample a high-frequency signal we don t get what we expect result looks like a lower frequency not possible to distinguish between this

More information

Computer Graphics (Fall 2011) Outline. CS 184 Guest Lecture: Sampling and Reconstruction Ravi Ramamoorthi

Computer Graphics (Fall 2011) Outline. CS 184 Guest Lecture: Sampling and Reconstruction Ravi Ramamoorthi Computer Graphics (Fall 2011) CS 184 Guest Lecture: Sampling and Reconstruction Ravi Ramamoorthi Some slides courtesy Thomas Funkhouser and Pat Hanrahan Adapted version of CS 283 lecture http://inst.eecs.berkeley.edu/~cs283/fa10

More information

CS 775: Advanced Computer Graphics. Lecture 12 : Antialiasing

CS 775: Advanced Computer Graphics. Lecture 12 : Antialiasing CS 775: Advanced Computer Graphics Lecture 12 : Antialiasing Antialiasing How to prevent aliasing? Prefiltering Analytic Approximate Postfiltering Supersampling Stochastic Supersampling Antialiasing Textures

More information

06: Thinking in Frequencies. CS 5840: Computer Vision Instructor: Jonathan Ventura

06: Thinking in Frequencies. CS 5840: Computer Vision Instructor: Jonathan Ventura 06: Thinking in Frequencies CS 5840: Computer Vision Instructor: Jonathan Ventura Decomposition of Functions Taylor series: Sum of polynomials f(x) =f(a)+f 0 (a)(x a)+ f 00 (a) 2! (x a) 2 + f 000 (a) (x

More information

CS4495/6495 Introduction to Computer Vision. 2C-L3 Aliasing

CS4495/6495 Introduction to Computer Vision. 2C-L3 Aliasing CS4495/6495 Introduction to Computer Vision 2C-L3 Aliasing Recall: Fourier Pairs (from Szeliski) Fourier Transform Sampling Pairs FT of an impulse train is an impulse train Sampling and Aliasing Sampling

More information

Image Scaling. This image is too big to fit on the screen. How can we reduce it? How to generate a halfsized

Image Scaling. This image is too big to fit on the screen. How can we reduce it? How to generate a halfsized Resampling Image Scaling This image is too big to fit on the screen. How can we reduce it? How to generate a halfsized version? Image sub-sampling 1/8 1/4 Throw away every other row and column to create

More information

Image Sampling. Moire patterns. - Source: F. Durand

Image Sampling. Moire patterns. -  Source: F. Durand Image Sampling Moire patterns Source: F. Durand - http://www.sandlotscience.com/moire/circular_3_moire.htm Any questions on project 1? For extra credits, attach before/after images how your extra feature

More information

Sampling and reconstruction. CS 4620 Lecture 13

Sampling and reconstruction. CS 4620 Lecture 13 Sampling and reconstruction CS 4620 Lecture 13 Lecture 13 1 Outline Review signal processing Sampling Reconstruction Filtering Convolution Closely related to computer graphics topics such as Image processing

More information

Sampling and Pyramids

Sampling and Pyramids Sampling and Pyramids 15-463: Rendering and Image Processing Alexei Efros with lots of slides from Steve Seitz Today Sampling Nyquist Rate Antialiasing Gaussian and Laplacian Pyramids 1 Fourier transform

More information

Module 3 : Sampling and Reconstruction Problem Set 3

Module 3 : Sampling and Reconstruction Problem Set 3 Module 3 : Sampling and Reconstruction Problem Set 3 Problem 1 Shown in figure below is a system in which the sampling signal is an impulse train with alternating sign. The sampling signal p(t), the Fourier

More information

Sampling and Reconstruction

Sampling and Reconstruction Sampling and reconstruction COMP 575/COMP 770 Fall 2010 Stephen J. Guy 1 Review What is Computer Graphics? Computer graphics: The study of creating, manipulating, and using visual images in the computer.

More information

ECE 484 Digital Image Processing Lec 09 - Image Resampling

ECE 484 Digital Image Processing Lec 09 - Image Resampling ECE 484 Digital Image Processing Lec 09 - Image Resampling Zhu Li Dept of CSEE, UMKC Office: FH560E, Email: lizhu@umkc.edu, Ph: x 2346. http://l.web.umkc.edu/lizhu slides created with WPS Office Linux

More information

קורס גרפיקה ממוחשבת 2008 סמסטר ב' Image Processing 1 חלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור

קורס גרפיקה ממוחשבת 2008 סמסטר ב' Image Processing 1 חלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור קורס גרפיקה ממוחשבת 2008 סמסטר ב' Image Processing 1 חלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור What is an image? An image is a discrete array of samples representing a continuous

More information

Sampling and reconstruction

Sampling and reconstruction Sampling and reconstruction Week 10 Acknowledgement: The course slides are adapted from the slides prepared by Steve Marschner of Cornell University 1 Sampled representations How to store and compute with

More information

Image Processing. What is an image? קורס גרפיקה ממוחשבת 2008 סמסטר ב' Converting to digital form. Sampling and Reconstruction.

Image Processing. What is an image? קורס גרפיקה ממוחשבת 2008 סמסטר ב' Converting to digital form. Sampling and Reconstruction. Amplitude 5/1/008 What is an image? An image is a discrete array of samples representing a continuous D function קורס גרפיקה ממוחשבת 008 סמסטר ב' Continuous function Discrete samples 1 חלק מהשקפים מעובדים

More information

Announcements. Image Processing. What s an image? Images as functions. Image processing. What s a digital image?

Announcements. Image Processing. What s an image? Images as functions. Image processing. What s a digital image? Image Processing Images by Pawan Sinha Today s readings Forsyth & Ponce, chapters 8.-8. http://www.cs.washington.edu/education/courses/49cv/wi/readings/book-7-revised-a-indx.pdf For Monday Watt,.3-.4 (handout)

More information

Sampling and Signal Processing

Sampling and Signal Processing Sampling and Signal Processing Sampling Methods Sampling is most commonly done with two devices, the sample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquires a continuous-time signal

More information

Sampling and reconstruction

Sampling and reconstruction Sampling and reconstruction CS 5625 Lecture 6 Lecture 6 1 Sampled representations How to store and compute with continuous functions? Common scheme for representation: samples write down the function s

More information

2D Discrete Fourier Transform

2D Discrete Fourier Transform 2D Discrete Fourier Transform In these lecture notes the figures have been removed for copyright reasons. References to figures are given instead, please check the figures yourself as given in the course

More information

Antialiasing & Compositing

Antialiasing & Compositing Antialiasing & Compositing CS4620 Lecture 14 Cornell CS4620/5620 Fall 2013 Lecture 14 (with previous instructors James/Bala, and some slides courtesy Leonard McMillan) 1 Pixel coverage Antialiasing and

More information

Sampling and Reconstruction of Analog Signals

Sampling and Reconstruction of Analog Signals Sampling and Reconstruction of Analog Signals Chapter Intended Learning Outcomes: (i) Ability to convert an analog signal to a discrete-time sequence via sampling (ii) Ability to construct an analog signal

More information

SAMPLING THEORY. Representing continuous signals with discrete numbers

SAMPLING THEORY. Representing continuous signals with discrete numbers SAMPLING THEORY Representing continuous signals with discrete numbers Roger B. Dannenberg Professor of Computer Science, Art, and Music Carnegie Mellon University ICM Week 3 Copyright 2002-2013 by Roger

More information

Prof. Feng Liu. Fall /04/2018

Prof. Feng Liu. Fall /04/2018 Prof. Feng Liu Fall 2018 http://www.cs.pdx.edu/~fliu/courses/cs447/ 10/04/2018 1 Last Time Image file formats Color quantization 2 Today Dithering Signal Processing Homework 1 due today in class Homework

More information

Sampling and pixels. CS 178, Spring Marc Levoy Computer Science Department Stanford University. Begun 4/23, finished 4/25.

Sampling and pixels. CS 178, Spring Marc Levoy Computer Science Department Stanford University. Begun 4/23, finished 4/25. Sampling and pixels CS 178, Spring 2013 Begun 4/23, finished 4/25. Marc Levoy Computer Science Department Stanford University Why study sampling theory? Why do I sometimes get moiré artifacts in my images?

More information

Last Lecture. photomatix.com

Last Lecture. photomatix.com Last Lecture photomatix.com Today Image Processing: from basic concepts to latest techniques Filtering Edge detection Re-sampling and aliasing Image Pyramids (Gaussian and Laplacian) Removing handshake

More information

Exercise Problems: Information Theory and Coding

Exercise Problems: Information Theory and Coding Exercise Problems: Information Theory and Coding Exercise 9 1. An error-correcting Hamming code uses a 7 bit block size in order to guarantee the detection, and hence the correction, of any single bit

More information

ECE 2111 Signals and Systems Spring 2012, UMD Experiment 9: Sampling

ECE 2111 Signals and Systems Spring 2012, UMD Experiment 9: Sampling ECE 2111 Signals and Systems Spring 2012, UMD Experiment 9: Sampling Objective: In this experiment the properties and limitations of the sampling theorem are investigated. A specific sampling circuit will

More information

!"!#"#$% Lecture 2: Media Creation. Some materials taken from Prof. Yao Wang s slides RECAP

!!##$% Lecture 2: Media Creation. Some materials taken from Prof. Yao Wang s slides RECAP Lecture 2: Media Creation Some materials taken from Prof. Yao Wang s slides RECAP #% A Big Umbrella Content Creation: produce the media, compress it to a format that is portable/ deliverable Distribution:

More information

Lecture Schedule: Week Date Lecture Title

Lecture Schedule: Week Date Lecture Title http://elec3004.org Sampling & More 2014 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 1 2-Mar Introduction 3-Mar

More information

Fourier Transform. Any signal can be expressed as a linear combination of a bunch of sine gratings of different frequency Amplitude Phase

Fourier Transform. Any signal can be expressed as a linear combination of a bunch of sine gratings of different frequency Amplitude Phase Fourier Transform Fourier Transform Any signal can be expressed as a linear combination of a bunch of sine gratings of different frequency Amplitude Phase 2 1 3 3 3 1 sin 3 3 1 3 sin 3 1 sin 5 5 1 3 sin

More information

Image Pyramids. Sanja Fidler CSC420: Intro to Image Understanding 1 / 35

Image Pyramids. Sanja Fidler CSC420: Intro to Image Understanding 1 / 35 Image Pyramids Sanja Fidler CSC420: Intro to Image Understanding 1 / 35 Finding Waldo Let s revisit the problem of finding Waldo This time he is on the road template (filter) image Sanja Fidler CSC420:

More information

Image Filtering and Gaussian Pyramids

Image Filtering and Gaussian Pyramids Image Filtering and Gaussian Pyramids CS94: Image Manipulation & Computational Photography Alexei Efros, UC Berkeley, Fall 27 Limitations of Point Processing Q: What happens if I reshuffle all pixels within

More information

Signals and Systems. Lecture 13 Wednesday 6 th December 2017 DR TANIA STATHAKI

Signals and Systems. Lecture 13 Wednesday 6 th December 2017 DR TANIA STATHAKI Signals and Systems Lecture 13 Wednesday 6 th December 2017 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Continuous time versus discrete time Continuous time

More information

image Scanner, digital camera, media, brushes,

image Scanner, digital camera, media, brushes, 118 Also known as rasterr graphics Record a value for every pixel in the image Often created from an external source Scanner, digital camera, Painting P i programs allow direct creation of images with

More information

Sampling of Continuous-Time Signals. Reference chapter 4 in Oppenheim and Schafer.

Sampling of Continuous-Time Signals. Reference chapter 4 in Oppenheim and Schafer. Sampling of Continuous-Time Signals Reference chapter 4 in Oppenheim and Schafer. Periodic Sampling of Continuous Signals T = sampling period fs = sampling frequency when expressing frequencies in radians

More information

Last Lecture. photomatix.com

Last Lecture. photomatix.com Last Lecture photomatix.com HDR Video Assorted pixel (Single Exposure HDR) Assorted pixel Assorted pixel Pixel with Adaptive Exposure Control light attenuator element detector element T t+1 I t controller

More information

CS 465 Prelim 1. Tuesday 4 October hours. Problem 1: Image formats (18 pts)

CS 465 Prelim 1. Tuesday 4 October hours. Problem 1: Image formats (18 pts) CS 465 Prelim 1 Tuesday 4 October 2005 1.5 hours Problem 1: Image formats (18 pts) 1. Give a common pixel data format that uses up the following numbers of bits per pixel: 8, 16, 32, 36. For instance,

More information

Image Resizing. Reminder no class next Tuesday, but your problem set is due. 9/19/08 Comp 665 Image Resizing 1

Image Resizing. Reminder no class next Tuesday, but your problem set is due. 9/19/08 Comp 665 Image Resizing 1 Image Resizing Narbonic Shaenon Garrity http://www.narbonic.com Reminder no class next Tuesday, but your problem set is due. 9/19/08 Comp 665 Image Resizing 1 Magnifica/on = Reconstruc/on Conceptually

More information

Fourier Transforms in Radio Astronomy

Fourier Transforms in Radio Astronomy Fourier Transforms in Radio Astronomy Kavilan Moodley, UKZN Slides taken from N Gupta s lectures: SKA School 2013 van-cittert Zernike theorem Extended, quasi-monochromatic, incoherent source X (l,m) Y

More information

+ a(t) exp( 2πif t)dt (1.1) In order to go back to the independent variable t, we define the inverse transform as: + A(f) exp(2πif t)df (1.

+ a(t) exp( 2πif t)dt (1.1) In order to go back to the independent variable t, we define the inverse transform as: + A(f) exp(2πif t)df (1. Chapter Fourier analysis In this chapter we review some basic results from signal analysis and processing. We shall not go into detail and assume the reader has some basic background in signal analysis

More information

CoE4TN4 Image Processing. Chapter 4 Filtering in the Frequency Domain

CoE4TN4 Image Processing. Chapter 4 Filtering in the Frequency Domain CoE4TN4 Image Processing Chapter 4 Filtering in the Frequency Domain Fourier Transform Sections 4.1 to 4.5 will be done on the board 2 2D Fourier Transform 3 2D Sampling and Aliasing 4 2D Sampling and

More information

Frequency Domain Enhancement

Frequency Domain Enhancement Tutorial Report Frequency Domain Enhancement Page 1 of 21 Frequency Domain Enhancement ESE 558 - DIGITAL IMAGE PROCESSING Tutorial Report Instructor: Murali Subbarao Written by: Tutorial Report Frequency

More information

8.2 IMAGE PROCESSING VERSUS IMAGE ANALYSIS Image processing: The collection of routines and

8.2 IMAGE PROCESSING VERSUS IMAGE ANALYSIS Image processing: The collection of routines and 8.1 INTRODUCTION In this chapter, we will study and discuss some fundamental techniques for image processing and image analysis, with a few examples of routines developed for certain purposes. 8.2 IMAGE

More information

EE521 Analog and Digital Communications

EE521 Analog and Digital Communications EE51 Analog and Digital Communications January 5, 006 Instructor: James K Beard, PhD Office: Ft. Washington TBA Email:. jkbeard@temple.edu, jkbeard@comcast.net Office Hours: Wednesdays 3:00 PM to 4:30

More information

Image Processing for Mechatronics Engineering For senior undergraduate students Academic Year 2017/2018, Winter Semester

Image Processing for Mechatronics Engineering For senior undergraduate students Academic Year 2017/2018, Winter Semester Image Processing for Mechatronics Engineering For senior undergraduate students Academic Year 2017/2018, Winter Semester Lecture 6: Image Acquisition and Digitization 14.10.2017 Dr. Mohammed Abdel-Megeed

More information

In The Name of Almighty. Lec. 2: Sampling

In The Name of Almighty. Lec. 2: Sampling In The Name of Almighty Lec. 2: Sampling Lecturer: Hooman Farkhani Department of Electrical Engineering Islamic Azad University of Najafabad Feb. 2016. Email: H_farkhani@yahoo.com A/D and D/A Conversion

More information

Image Interpolation. Image Processing

Image Interpolation. Image Processing Image Interpolation Image Processing Brent M. Dingle, Ph.D. 2015 Game Design and Development Program Mathematics, Statistics and Computer Science University of Wisconsin - Stout public domain image from

More information

Thinking in Frequency

Thinking in Frequency Thinking in Frequency Computer Vision Brown James Hays Slides: Hoiem, Efros, and others Recap of Wednesday linear filtering convolution differential filters filter types boundary conditions. Review: questions

More information

Texture mapping from 0 to infinity

Texture mapping from 0 to infinity Announcements CS4620/5620: Lecture 24 HW 3 out Barycentric coordinates for Problem 1 Texture Mapping 1 2 Texture mapping from 0 to infinity When you go close... When viewed from a distance Aliasing! 3

More information

Topic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music)

Topic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music) Topic 2 Signal Processing Review (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music) Recording Sound Mechanical Vibration Pressure Waves Motion->Voltage Transducer

More information

Image features: Histograms, Aliasing, Filters, Orientation and HOG. D.A. Forsyth

Image features: Histograms, Aliasing, Filters, Orientation and HOG. D.A. Forsyth Image features: Histograms, Aliasing, Filters, Orientation and HOG D.A. Forsyth Simple color features Histogram of image colors in a window Opponent color representations R-G B-Y=B-(R+G)/2 Intensity=(R+G+B)/3

More information

Fourier analysis of images

Fourier analysis of images Fourier analysis of images Intensity Image Fourier Image Slides: James Hays, Hoiem, Efros, and others http://sharp.bu.edu/~slehar/fourier/fourier.html#filtering Signals can be composed + = http://sharp.bu.edu/~slehar/fourier/fourier.html#filtering

More information

Recall. Sampling. Why discrete time? Why discrete time? Many signals are continuous-time signals Light Object wave CCD

Recall. Sampling. Why discrete time? Why discrete time? Many signals are continuous-time signals Light Object wave CCD Recall Many signals are continuous-time signals Light Object wave CCD Sampling mic Lens change of voltage change of voltage 2 Why discrete time? With the advance of computer technology, we want to process

More information

ACM Fast Image Convolutions. by: Wojciech Jarosz

ACM Fast Image Convolutions. by: Wojciech Jarosz ACM SIGGRAPH@UIUC Fast Image Convolutions by: Wojciech Jarosz Image Convolution Traditionally, image convolution is performed by what is called the sliding window approach. For each pixel in the image,

More information

Filters. Motivating Example. Tracking a fly, oh my! Moving Weighted Average Filter. General Picture

Filters. Motivating Example. Tracking a fly, oh my! Moving Weighted Average Filter. General Picture Motivating Example Filters Consider we are tracking a fly Sensor reports the fly s position several times a second Some noise in the sensor Goal: reconstruct the fly s actual path Problem: can t rely on

More information

Sistemas de Aquisição de Dados. Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 3-29 de Setembro

Sistemas de Aquisição de Dados. Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 3-29 de Setembro Sistemas de Aquisição de Dados Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 3-29 de Setembro Aliasing Example fsig=101khz fsig=899 khz All sampled signals are equal! fsig=1101 khz 2 How to

More information

Vision Review: Image Processing. Course web page:

Vision Review: Image Processing. Course web page: Vision Review: Image Processing Course web page: www.cis.udel.edu/~cer/arv September 7, Announcements Homework and paper presentation guidelines are up on web page Readings for next Tuesday: Chapters 6,.,

More information

Midterm is on Thursday!

Midterm is on Thursday! Midterm is on Thursday! Project presentations are May 17th, 22nd and 24th Next week there is a strike on campus. Class is therefore cancelled on Tuesday. Please work on your presentations instead! REVIEW

More information

To Do. Advanced Computer Graphics. Image Compositing. Digital Image Compositing. Outline. Blue Screen Matting

To Do. Advanced Computer Graphics. Image Compositing. Digital Image Compositing. Outline. Blue Screen Matting Advanced Computer Graphics CSE 163 [Spring 2018], Lecture 5 Ravi Ramamoorthi http://www.cs.ucsd.edu/~ravir To Do Assignment 1, Due Apr 27. This lecture only extra credit and clear up difficulties Questions/difficulties

More information

Image Processing for feature extraction

Image Processing for feature extraction Image Processing for feature extraction 1 Outline Rationale for image pre-processing Gray-scale transformations Geometric transformations Local preprocessing Reading: Sonka et al 5.1, 5.2, 5.3 2 Image

More information

CS6670: Computer Vision Noah Snavely. Administrivia. Administrivia. Reading. Last time: Convolution. Last time: Cross correlation 9/8/2009

CS6670: Computer Vision Noah Snavely. Administrivia. Administrivia. Reading. Last time: Convolution. Last time: Cross correlation 9/8/2009 CS667: Computer Vision Noah Snavely Administrivia New room starting Thursday: HLS B Lecture 2: Edge detection and resampling From Sandlot Science Administrivia Assignment (feature detection and matching)

More information

Islamic University of Gaza. Faculty of Engineering Electrical Engineering Department Spring-2011

Islamic University of Gaza. Faculty of Engineering Electrical Engineering Department Spring-2011 Islamic University of Gaza Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#4 Sampling and Quantization OBJECTIVES: When you have completed this assignment,

More information

Chapter 2: Digitization of Sound

Chapter 2: Digitization of Sound Chapter 2: Digitization of Sound Acoustics pressure waves are converted to electrical signals by use of a microphone. The output signal from the microphone is an analog signal, i.e., a continuous-valued

More information

System Identification & Parameter Estimation

System Identification & Parameter Estimation System Identification & Parameter Estimation Wb2301: SIPE lecture 4 Perturbation signal design Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE 3/9/2010 Delft University of Technology

More information

Biomedical Instrumentation B2. Dealing with noise

Biomedical Instrumentation B2. Dealing with noise Biomedical Instrumentation B2. Dealing with noise B18/BME2 Dr Gari Clifford Noise & artifact in biomedical signals Ambient / power line interference: 50 ±0.2 Hz mains noise (or 60 Hz in many data sets)

More information

Analysis and design of filters for differentiation

Analysis and design of filters for differentiation Differential filters Analysis and design of filters for differentiation John C. Bancroft and Hugh D. Geiger SUMMARY Differential equations are an integral part of seismic processing. In the discrete computer

More information

Module 3: Video Sampling Lecture 18: Filtering operations in Camera and display devices. The Lecture Contains: Effect of Temporal Aperture:

Module 3: Video Sampling Lecture 18: Filtering operations in Camera and display devices. The Lecture Contains: Effect of Temporal Aperture: The Lecture Contains: Effect of Temporal Aperture: Spatial Aperture: Effect of Display Aperture: file:///d /...e%20(ganesh%20rana)/my%20course_ganesh%20rana/prof.%20sumana%20gupta/final%20dvsp/lecture18/18_1.htm[12/30/2015

More information

Next Classes. Spatial frequency Fourier transform and frequency domain. Reminder: Textbook. Frequency view of filtering Hybrid images Sampling

Next Classes. Spatial frequency Fourier transform and frequency domain. Reminder: Textbook. Frequency view of filtering Hybrid images Sampling Salvador Dali, 1976 Next Classes Spatial frequency Fourier transform and frequency domain Frequency view of filtering Hybrid images Sampling Reminder: Textbook Today s lecture covers material in 3.4 Slide:

More information

MULTIMEDIA SYSTEMS

MULTIMEDIA SYSTEMS 1 Department of Computer Engineering, g, Faculty of Engineering King Mongkut s Institute of Technology Ladkrabang 01076531 MULTIMEDIA SYSTEMS Pakorn Watanachaturaporn, Ph.D. pakorn@live.kmitl.ac.th, pwatanac@gmail.com

More information

On the rendering of synthetic images with specific point spread functions F. van den Bergh

On the rendering of synthetic images with specific point spread functions F. van den Bergh On the rendering of synthetic images with specific point spread functions F. van den Bergh Remote Sensing Research Unit, Meraka Institute CSIR, PO Box 395, Pretoria South Africa, 0001 fvdbergh@csir.co.za

More information

The Intuitions of Signal Processing (for Motion Editing)

The Intuitions of Signal Processing (for Motion Editing) The Intuitions of Signal Processing (for Motion Editing) This chapter will be an appendix of the book Motion Capture and Motion Editing: Bridging Principle and Practice, by Jung, Fischer, Gleicher, and

More information

Computer Vision, Lecture 3

Computer Vision, Lecture 3 Computer Vision, Lecture 3 Professor Hager http://www.cs.jhu.edu/~hager /4/200 CS 46, Copyright G.D. Hager Outline for Today Image noise Filtering by Convolution Properties of Convolution /4/200 CS 46,

More information

Midterm Examination CS 534: Computational Photography

Midterm Examination CS 534: Computational Photography Midterm Examination CS 534: Computational Photography November 3, 2015 NAME: SOLUTIONS Problem Score Max Score 1 8 2 8 3 9 4 4 5 3 6 4 7 6 8 13 9 7 10 4 11 7 12 10 13 9 14 8 Total 100 1 1. [8] What are

More information

Human Reconstruction of Digitized Graphical Signals

Human Reconstruction of Digitized Graphical Signals Proceedings of the International MultiConference of Engineers and Computer Scientists 8 Vol II IMECS 8, March -, 8, Hong Kong Human Reconstruction of Digitized Graphical s Coskun DIZMEN,, and Errol R.

More information

Experiment 8: Sampling

Experiment 8: Sampling Prepared By: 1 Experiment 8: Sampling Objective The objective of this Lab is to understand concepts and observe the effects of periodically sampling a continuous signal at different sampling rates, changing

More information

Table of contents. Vision industrielle 2002/2003. Local and semi-local smoothing. Linear noise filtering: example. Convolution: introduction

Table of contents. Vision industrielle 2002/2003. Local and semi-local smoothing. Linear noise filtering: example. Convolution: introduction Table of contents Vision industrielle 2002/2003 Session - Image Processing Département Génie Productique INSA de Lyon Christian Wolf wolf@rfv.insa-lyon.fr Introduction Motivation, human vision, history,

More information

DISCRETE FOURIER TRANSFORM AND FILTER DESIGN

DISCRETE FOURIER TRANSFORM AND FILTER DESIGN DISCRETE FOURIER TRANSFORM AND FILTER DESIGN N. C. State University CSC557 Multimedia Computing and Networking Fall 2001 Lecture # 03 Spectrum of a Square Wave 2 Results of Some Filters 3 Notation 4 x[n]

More information

Photo-Consistent Motion Blur Modeling for Realistic Image Synthesis

Photo-Consistent Motion Blur Modeling for Realistic Image Synthesis Photo-Consistent Motion Blur Modeling for Realistic Image Synthesis Huei-Yung Lin and Chia-Hong Chang Department of Electrical Engineering, National Chung Cheng University, 168 University Rd., Min-Hsiung

More information

PHYS 352. FFT Convolution. More Advanced Digital Signal Processing Techniques

PHYS 352. FFT Convolution. More Advanced Digital Signal Processing Techniques PHYS 352 More Advanced Digital Signal Processing Techniques FFT Convolution take a chunk of your signal (say N=128 samples) apply FFT to it multiply the frequency domain signal by your desired transfer

More information

Chapter 6 CONTINUOUS-TIME, IMPULSE-MODULATED, AND DISCRETE-TIME SIGNALS. 6.6 Sampling Theorem 6.7 Aliasing 6.8 Interrelations

Chapter 6 CONTINUOUS-TIME, IMPULSE-MODULATED, AND DISCRETE-TIME SIGNALS. 6.6 Sampling Theorem 6.7 Aliasing 6.8 Interrelations Chapter 6 CONTINUOUS-TIME, IMPULSE-MODULATED, AND DISCRETE-TIME SIGNALS 6.6 Sampling Theorem 6.7 Aliasing 6.8 Interrelations Copyright c 2005- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org

More information

Sampling and Reconstruction

Sampling and Reconstruction Sampling and Reconstruction Many slides from Steve Marschner 15-463: Computational Photography Alexei Efros, CMU, Fall 211 Sampling and Reconstruction Sampled representations How to store and compute with

More information

CS3291: Digital Signal Processing

CS3291: Digital Signal Processing CS39 Exam Jan 005 //08 /BMGC University of Manchester Department of Computer Science First Semester Year 3 Examination Paper CS39: Digital Signal Processing Date of Examination: January 005 Answer THREE

More information

CSCI 1290: Comp Photo

CSCI 1290: Comp Photo CSCI 29: Comp Photo Fall 28 @ Brown University James Tompkin Many slides thanks to James Hays old CS 29 course, along with all of its acknowledgements. Things I forgot on Thursday Grads are not required

More information

SIGMA-DELTA CONVERTER

SIGMA-DELTA CONVERTER SIGMA-DELTA CONVERTER (1995: Pacífico R. Concetti Western A. Geophysical-Argentina) The Sigma-Delta A/D Converter is not new in electronic engineering since it has been previously used as part of many

More information

LAB 2: Sampling & aliasing; quantization & false contouring

LAB 2: Sampling & aliasing; quantization & false contouring CEE 615: Digital Image Processing Spring 2016 1 LAB 2: Sampling & aliasing; quantization & false contouring A. SAMPLING: Observe the effects of the sampling interval near the resolution limit. The goal

More information

Transforms and Frequency Filtering

Transforms and Frequency Filtering Transforms and Frequency Filtering Khalid Niazi Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University 2 Reading Instructions Chapter 4: Image Enhancement in the Frequency

More information

ELECTRONIC HOLOGRAPHY

ELECTRONIC HOLOGRAPHY ELECTRONIC HOLOGRAPHY CCD-camera replaces film as the recording medium. Electronic holography is better suited than film-based holography to quantitative applications including: - phase microscopy - metrology

More information

Design of practical color filter array interpolation algorithms for digital cameras

Design of practical color filter array interpolation algorithms for digital cameras Design of practical color filter array interpolation algorithms for digital cameras James E. Adams, Jr. Eastman Kodak Company, Imaging Research and Advanced Development Rochester, New York 14653-5408 ABSTRACT

More information

XXXX - ANTI-ALIASING AND RESAMPLING 1 N/08/08

XXXX - ANTI-ALIASING AND RESAMPLING 1 N/08/08 INTRODUCTION TO GRAPHICS Anti-Aliasing and Resampling Information Sheet No. XXXX The fundamental fundamentals of bitmap images and anti-aliasing are a fair enough topic for beginners and it s not a bad

More information

Sampling Rate = Resolution Quantization Level = Color Depth = Bit Depth = Number of Colors

Sampling Rate = Resolution Quantization Level = Color Depth = Bit Depth = Number of Colors ITEC2110 FALL 2011 TEST 2 REVIEW Chapters 2-3: Images I. Concepts Graphics A. Bitmaps and Vector Representations Logical vs. Physical Pixels - Images are modeled internally as an array of pixel values

More information

Sampling and Reconstruction

Sampling and Reconstruction Sampling and Reconstruction Salvador Dali, Dali from the Back Painting Gala from the Back Eternalized by Six Virtual Corneas Provisionally Reflected by Six Real Mirrors Many slides from Steve Marschner,

More information

Dynamically Reparameterized Light Fields & Fourier Slice Photography. Oliver Barth, 2009 Max Planck Institute Saarbrücken

Dynamically Reparameterized Light Fields & Fourier Slice Photography. Oliver Barth, 2009 Max Planck Institute Saarbrücken Dynamically Reparameterized Light Fields & Fourier Slice Photography Oliver Barth, 2009 Max Planck Institute Saarbrücken Background What we are talking about? 2 / 83 Background What we are talking about?

More information

Digital Processing of Continuous-Time Signals

Digital Processing of Continuous-Time Signals Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital

More information

Anti aliasing and Graphics Formats

Anti aliasing and Graphics Formats Anti aliasing and Graphics Formats Eric C. McCreath School of Computer Science The Australian National University ACT 0200 Australia ericm@cs.anu.edu.au Overview 2 Nyquist sampling frequency supersampling

More information

Image and Video Processing

Image and Video Processing Image and Video Processing () Image Representation Dr. Miles Hansard miles.hansard@qmul.ac.uk Segmentation 2 Today s agenda Digital image representation Sampling Quantization Sub-sampling Pixel interpolation

More information

Today. Defocus. Deconvolution / inverse filters. MIT 2.71/2.710 Optics 12/12/05 wk15-a-1

Today. Defocus. Deconvolution / inverse filters. MIT 2.71/2.710 Optics 12/12/05 wk15-a-1 Today Defocus Deconvolution / inverse filters MIT.7/.70 Optics //05 wk5-a- MIT.7/.70 Optics //05 wk5-a- Defocus MIT.7/.70 Optics //05 wk5-a-3 0 th Century Fox Focus in classical imaging in-focus defocus

More information