ELEC Dr Reji Mathew Electrical Engineering UNSW
|
|
- Roxanne Phillips
- 5 years ago
- Views:
Transcription
1 ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW
2 Filter Design Circularly symmetric 2-D low-pass filter Pass-band radial frequency: ω p Stop-band radial frequency: ω s 1 δ p Pass-band tolerances: δ p Stop-band tolerances: δ s How to design this low pass filter? Windowing? ω p ω s ω δ s
3 Windowing Design a 2-D FIR filter Start with a desired frequency response h d (ω) Where h d (ω) satisfies the requirements (stop & pass band) Perhaps an ideal filter h d (ω) Then take the inverse DSFT to obtain h d [n] h d [n] may have infinite support Apply a spatial window w[n] to limit support to a finite region n n
4 Windowing Windowing in the spatial domain: affects filter characteristics in the frequency domain Does h(ω) which is the DSFT of h n still satisfy the original constrains? Multiplication of spatial signals equates to convolution of the corresponding DSFT of the signals Ideally - window will have as narrow a bandwidth as possible That is w(ω) should decay rapidly with increasing ω However, narrow band windows have large spatial support Constrained by the desired support size
5 Windowing Remember: we want a filter with zero phase h ω = h(ω) w(ω) and h d (ω) will both need to be real-valued So h ω will also be real-valued (zero phase) Follows on from the convolution equation Convenient: express window as spatially continuous function w x And w n = w x x=n We assume that the window is sufficiently narrow band Essentially bandlimited to the Nyquist region, ω [ π, π] 2 So FT of w x is identical to the DSFT of w n Therefore, we are only concerned with the properties of the spatially continuous function w x
6 Windowing Simplest case: window is a separable product of 1D functions Determines the extent of the window This trivial window is rarely used in practice because its frequency response decays only slowly.
7 Windowing Smooth function: means high frequencies components are avoided minimising the effective bandwidth of w(ω) one period of a cosine raised with its negative peaks just touch zero ( raised cosine ) t Multiply a window: To Restrict an impulse response that is infinite t
8 Windowing
9 Windowing For some applications, superior results can be obtained by using non-separable windows. When designing a circularly symmetric filter, a circularly symmetric window function is generally preferable w(s) = w t t= s Where w t can be the 1D window functions we discussed earlier.
10 Filter Design Frequency Sampling Find a finite set of filter tap values Which minimizes a weighted combination of squared differences between the desired and actual frequency responses Can uniformly distribute ω k in the region ω [ π, π] 2 Alternatively, can place more samples in the neighbourhood of critical frequency regions such as (pass and stop) band edges
11 Frequency Sampling Want to ensure that the designed filter h n is zero phase Remember that zero phase implies a symmetric filter. h n = h n The sin terms are cancelled out. Note summation over
12 Frequency Sampling We wish to minimize the weighted squared error expression weight which reflects the relative importance of matching the desired response at frequency ω k Summation for k = 1 K
13 Frequency Sampling
14 Frequency Sampling Minimization of the weighted squared error expression is a classic linear algebra problem whose solution is well known: Can adjust the weights and number of samples (K) Reduce the density of samples in critical frequency regions by increasing weights associated with a smaller number of samples
15 Filter Design Transformation Method (McClellan Transformation) 1D filter design problem is much easier Methods exist for designing optimum 1D zero phase FIR filters subject to constraints on tolerances in the pass- and stop-bands Map the 2-D design problem into a 1-D design problem Without losing the finite impulse response property. The method is relevant only for zero phase designs. An approximate circularly symmetric function specifies the 2D to 1D mapping
16 Transformation Method Frequency response of a 1-D zero phase filter with 2N + 1 taps Taking the DTFT of zero phase filter h n, defined for -N to N Remember that zero phase implies a symmetric filter. h n = h n The sin terms are cancelled out. Note summation from 1 to N
17 Transformation Method Frequency response of a 1-D zero phase filter with 2N + 1 taps Taking the DTFT of zero phase filter h n, defined for -N to N By application of Chebychev polynomials
18 Transformation Method
19 Transformation Method x = cos ω C n = cos nω C 0 = cos 0ω = 1 C 1 = cos ω C 2 = cos 2ω = 2 cos ω cos ω 1 = 2cos 2 ω 1 C 3 = cos 3ω = 2 cos ω cos 2ω cos ω = 4cos 3 ω 3 cos ω cos nω can be expressed as an nth-degree polynomial in cos ω
20 Transformation Method Frequency response of a 1-D zero phase filter with 2N + 1 taps Taking the DTFT of zero phase filter h n, defined for -N to N By application of Chebychev polynomials Coefficients b n can be obtained from a n Next: Frequency mapping from ω 1, ω 2 cos (ω)
21 Transformation Method μ ω is the DSFT of a finite support sequence μ[n] Essentially μ[n] controls the mapping (kernel) h[n] has finite support Summation of n-fold convolution of finite support sequence μ[n]
22 Transformation Method Given μ[n], what can we say about the Region of Support for the final 2D filter h[n]? And what is DSFT of μ[n], which governs the mapping ω 1, ω 2 cos (ω)?
23 Transformation Method Example: for a 7 tap 1D filter N=3 (-3 0 3) So μ n convolved three times (due to μ ω 3 ) Since μ n has a region of support [ 1, 1] 2 The support of h[n] will be [ 3, 3] 2
24 Transformation Method DSFT of μ n The mapping represents quite a good approximation to a circularly symmetric function. (plot next slide)
25
26 Transformation Method Steps to follow for the Frequency Transform Method 1. Select a transformation kernel, μ n. Possible to design your own. The original McClellan kernel is a common choice. 2. Map the 2-D frequency domain specifications on h ω back to specification on the 1-D kernel, h ω. Stop band requirement for the 2D filter: ω ω s Find the largest ω s such that μ ω cos (ω s ) for all ω ω s ω 1 tolerances map directly from 2-D to the 1-D ω ω 2
27 Transformation Method Steps to follow for the Frequency Transform Method 1. Select a transformation kernel, μ n. Possible to design your own. The original McClellan kernel is a common choice. 2. Map the 2-D frequency domain specifications on h ω back to specification on the 1-D kernel, h ω. 3. Design a suitable 1-D kernel using well established methods. Thereby determine the b n coefficients. ω 1 ω ω 2
28 Transformation Method Input Image x n μ n μ n μ n b[0] b[1] b[2] b[3] Example: N = 3 Weighted summation of convolutions of the mapping kernel. Weights b n as determined from the 1D filter design
29 Filtering Examples Gaussian Filters Normalisation chosen so that Spatially continuous impulse response has unit DC gain
30 Gaussian Filters In most cases the off-diagonal term, σ 1,2, is zero So that h(s) is a separable Gaussian function whose contours have an elliptical cross-section (circular if σ 1 = σ 2 ). s 2 s 1
31 Gaussian Filters The σ 1,2 term may be used to design arbitrarily rotated elliptical cross-sections. s 2 s 1
32 Gaussian Filters Gaussian filters have the property that their Fourier transform is also a Gaussian function
33 Gaussian Filters Then: DSFT of the sampled impulse response, h n, will be almost identical to the spatially continuous Fourier transform of h(s). DSFT of the sampled impulse response is Gaussian and also that h n = 1 n Since the spatial impulse response decays rapidly as n becomes large, can simply truncate the impulse response (i.e. apply a rectangular window) to some suitably large region of support.
34 Gaussian Filters Important properties Gaussian filters - their Fourier transform is also a Gaussian function Both spatial & frequency domain representations are real & positive. For simplicity - Ignore the possibility that σ 1,2 might not be zero, Horizontal and vertical spread of the Gaussian function is given by σ 1 and σ 2 Horizontal and vertical spread in the frequency domain is given by 1 σ 1 and 1 σ 2 So the product of the spatial & frequency spreads is always exactly 1.
35 Moving Average Filters Low complexity Low pass, high pass and band pass implementations A2 B2 A1 B1 (B1 A1 + 1)(B2 A2 + 1) multiplications & additions for each pixel (B1 A1 + 1) + (B2 A2 + 1) multiplications & additions per output pixel
36 Moving Average Filters Low-pass filter, whose frequency response is a sinc function Filter is separable into the row and column filters Can be implemented independent of filter s region of support?
37 Moving Average Filters x(0,0) x(0,1) x(0,2) x(0,3) x(0,4) x(0,5) x(0,6) x(0,7) x(1,0) x(1,1) x(1,2) x(1,3) x(1,4) x(1,5) x(1,6) x(1,7) x(2,0) x(2,1) x(2,2) x(2,3) x(2,4) x(2,5) x(2,6) x(2,7) x(3,0) x(3,1) x(3,2) x(3,3) x(3,4) x(3,5) x(3,6) x(3,7) x(4,0) x(4,1) x(4,2) x(4,3) x(4,4) x(4,5) x(4,6) x(4,7) Each subsequent scaled output value requires only two additions
38 Moving Average Filters
39 Unsharp Masking One of the oldest image processing operations Ad-hoc technique: artificially increasing sharpness of an image x n
40 Image Acquisition Models Scene illuminant reflects off scene surfaces; the resulting radiant flux passes through an optical system, is integrated over sensor surfaces, converted to an electrical signal and quantized
41 Image Acquisition Models Scene surface reflectivities as a function of wavelength λ, and position s = (s 1, s 2 ) Radiant power spectral density ( scene radiance ) at the surface of the sensor Spectral power density of the illuminant as a function of position and wavelength Note: Variation in s and λ Example of variations in λ : Outdoor and indoor images
42 Image Acquisition Models We ignore surface orientation and its impact upon scene radiance. Also, we ignore fluorescence.
43 Optics Add the impact of optics to our simple model Resulting scene radiance at the surface of the sensor Point spread function (PSF) of the optical system In general, the optical PSF depends on both wavelength λ and spatial position z (ie varies across the sensor). We ignore the dependence of h o z,λ(s) on z. optics are not ideal, but are at least linear
44 Optics This is the convolution integral. Linear systems which are modeled by convolution are known as Linear Shift Invariant (LSI) systems. LSI property is valid only when we ignore the dependence of h o z,λ(s) on z. LSI model is sufficiently accurate in the so-called centre-field
45 Optics Away from centre field a number of aberrations appear Example: Chromatic aberrations h o z,λ(s) introduces shifts which depend upon both position z and wavelength λ. Different wavelengths (colours) experience different geometric distortions do not properly line up blurring
46 Optics Simplifications ignore spatial variations of the optical system PSF restrict to centre field Even with these simplifications the PSF is not simply an impulse non-ideal causes of these non-idealities are refractive and diffractive limitations Refractive: lens imperfections Diffractive: finite lens size wavelength dependent Narrow aperture: comparable in size to incident light wavelength
47 Sensor Integration Integration across all wavelength Spectral sensitivity 2D sequence of pixel values Integration over the aperture of the sensor Ideal sensor
48 Sensor Integration
49 Sensing Devices CCD photons converted into electrons: accumulated in charge wells collection of 1-D analog shift registers entire array of pixels must be shifted difficult to retrieve a subset of the pixels clocking the out pixels : Consumes a large amount of power
50 Sensing Devices CMOS Most popular even at high quality end Sensing element CMOS transistor with a transparent gate Or a photodiode connected to a CMOS amplifier At each photoreceptor is a CMOS circuit, which amplify the signal & provide addressing capability Pixels in CMOS sensors addressed in a manner similar random access memory Low power, inexpensive manufacturing System-on-a-chip: A-D conversion, image processing on the same chip
51 Sensing Devices CCD and CMOS typically respond to infra-red light Infra-red is outside the visible spectrum Therefore require infra-red absorbtion filters in the optical path
52 Combining Optical and Aperture Effects Combine the effects of optical and sensor integration into a single PSF h v λ (s) Integral over wavelength Convolution with optical PSF Multiplication with the aperture
53 Combining Optical and Aperture Effects
54 Combining Optical and Aperture Effects
55 Combining Optical and Aperture Effects System PSF is obtained by convolving the optical PSF with a flipped version of the sensor integration aperture. Sensor samples are then obtained by convolving radiant scene intensity at the sensor surface by the system PSF impulsively sampling the resulting spatially and integrating over the spectral response characteristic of the sensor
56 Combining Optical and Aperture Effects What s wrong if the system PSF is an impulse? What s wrong if the system PSF is excessively smooth with extent larger than the separation between pixels (excessively smooth disc)?
Chapter 2 Fourier Integral Representation of an Optical Image
Chapter 2 Fourier Integral Representation of an Optical This chapter describes optical transfer functions. The concepts of linearity and shift invariance were introduced in Chapter 1. This chapter continues
More informationDesign of FIR Filters
Design of FIR Filters Elena Punskaya www-sigproc.eng.cam.ac.uk/~op205 Some material adapted from courses by Prof. Simon Godsill, Dr. Arnaud Doucet, Dr. Malcolm Macleod and Prof. Peter Rayner 1 FIR as a
More informationDigital Filters IIR (& Their Corresponding Analog Filters) Week Date Lecture Title
http://elec3004.com Digital Filters IIR (& Their Corresponding Analog Filters) 2017 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date
More informationLecture Notes 10 Image Sensor Optics. Imaging optics. Pixel optics. Microlens
Lecture Notes 10 Image Sensor Optics Imaging optics Space-invariant model Space-varying model Pixel optics Transmission Vignetting Microlens EE 392B: Image Sensor Optics 10-1 Image Sensor Optics Microlens
More informationToday. 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 informationDeconvolution , , Computational Photography Fall 2018, Lecture 12
Deconvolution http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2018, Lecture 12 Course announcements Homework 3 is out. - Due October 12 th. - Any questions?
More informationImage 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 informationSUPER RESOLUTION INTRODUCTION
SUPER RESOLUTION Jnanavardhini - Online MultiDisciplinary Research Journal Ms. Amalorpavam.G Assistant Professor, Department of Computer Sciences, Sambhram Academy of Management. Studies, Bangalore Abstract:-
More informationTSBB09 Image Sensors 2018-HT2. Image Formation Part 1
TSBB09 Image Sensors 2018-HT2 Image Formation Part 1 Basic physics Electromagnetic radiation consists of electromagnetic waves With energy That propagate through space The waves consist of transversal
More informationFourier 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 informationF I R Filter (Finite Impulse Response)
F I R Filter (Finite Impulse Response) Ir. Dadang Gunawan, Ph.D Electrical Engineering University of Indonesia The Outline 7.1 State-of-the-art 7.2 Type of Linear Phase Filter 7.3 Summary of 4 Types FIR
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Multi-Resolution Processing Gaussian Pyramid Starting with an image x[n], which we will also label x 0 [n], Construct a sequence of progressively lower
More informationBiomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar
Biomedical Signals Signals and Images in Medicine Dr Nabeel Anwar Noise Removal: Time Domain Techniques 1. Synchronized Averaging (covered in lecture 1) 2. Moving Average Filters (today s topic) 3. Derivative
More informationΕισαγωγική στην Οπτική Απεικόνιση
Εισαγωγική στην Οπτική Απεικόνιση Δημήτριος Τζεράνης, Ph.D. Εμβιομηχανική και Βιοϊατρική Τεχνολογία Τμήμα Μηχανολόγων Μηχανικών Ε.Μ.Π. Χειμερινό Εξάμηνο 2015 Light: A type of EM Radiation EM radiation:
More informationDigital Filters FIR and IIR Systems
Digital Filters FIR and IIR Systems ELEC 3004: Systems: Signals & Controls Dr. Surya Singh (Some material adapted from courses by Russ Tedrake and Elena Punskaya) Lecture 16 elec3004@itee.uq.edu.au http://robotics.itee.uq.edu.au/~elec3004/
More informationECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction
More informationImage acquisition. In both cases, the digital sensing element is one of the following: Line array Area array. Single sensor
Image acquisition Digital images are acquired by direct digital acquisition (digital still/video cameras), or scanning material acquired as analog signals (slides, photographs, etc.). In both cases, the
More informationIMAGE ENHANCEMENT IN SPATIAL DOMAIN
A First Course in Machine Vision IMAGE ENHANCEMENT IN SPATIAL DOMAIN By: Ehsan Khoramshahi Definitions The principal objective of enhancement is to process an image so that the result is more suitable
More informationIMAGE FORMATION. Light source properties. Sensor characteristics Surface. Surface reflectance properties. Optics
IMAGE FORMATION Light source properties Sensor characteristics Surface Exposure shape Optics Surface reflectance properties ANALOG IMAGES An image can be understood as a 2D light intensity function f(x,y)
More informationSensors and Sensing Cameras and Camera Calibration
Sensors and Sensing Cameras and Camera Calibration Todor Stoyanov Mobile Robotics and Olfaction Lab Center for Applied Autonomous Sensor Systems Örebro University, Sweden todor.stoyanov@oru.se 20.11.2014
More informationModule 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 informationObjectives. Presentation Outline. Digital Modulation Lecture 03
Digital Modulation Lecture 03 Inter-Symbol Interference Power Spectral Density Richard Harris Objectives To be able to discuss Inter-Symbol Interference (ISI), its causes and possible remedies. To be able
More informationCorso di DATI e SEGNALI BIOMEDICI 1. Carmelina Ruggiero Laboratorio MedInfo
Corso di DATI e SEGNALI BIOMEDICI 1 Carmelina Ruggiero Laboratorio MedInfo Digital Filters Function of a Filter In signal processing, the functions of a filter are: to remove unwanted parts of the signal,
More informationBe aware that there is no universal notation for the various quantities.
Fourier Optics v2.4 Ray tracing is limited in its ability to describe optics because it ignores the wave properties of light. Diffraction is needed to explain image spatial resolution and contrast and
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More informationDeconvolution , , Computational Photography Fall 2017, Lecture 17
Deconvolution http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2017, Lecture 17 Course announcements Homework 4 is out. - Due October 26 th. - There was another
More informationLane Detection in Automotive
Lane Detection in Automotive Contents Introduction... 2 Image Processing... 2 Reading an image... 3 RGB to Gray... 3 Mean and Gaussian filtering... 5 Defining our Region of Interest... 6 BirdsEyeView Transformation...
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More informationPROBLEM SET 6. Note: This version is preliminary in that it does not yet have instructions for uploading the MATLAB problems.
PROBLEM SET 6 Issued: 2/32/19 Due: 3/1/19 Reading: During the past week we discussed change of discrete-time sampling rate, introducing the techniques of decimation and interpolation, which is covered
More informationPrinciples of Baseband Digital Data Transmission
Principles of Baseband Digital Data Transmission Prof. Wangrok Oh Dept. of Information Communications Eng. Chungnam National University Prof. Wangrok Oh(CNU) / 3 Overview Baseband Digital Data Transmission
More informationME scope Application Note 01 The FFT, Leakage, and Windowing
INTRODUCTION ME scope Application Note 01 The FFT, Leakage, and Windowing NOTE: The steps in this Application Note can be duplicated using any Package that includes the VES-3600 Advanced Signal Processing
More informationFinal Exam Practice Questions for Music 421, with Solutions
Final Exam Practice Questions for Music 4, with Solutions Elementary Fourier Relationships. For the window w = [/,,/ ], what is (a) the dc magnitude of the window transform? + (b) the magnitude at half
More informationUSE OF FT IN IMAGE PROCESSING IMAGE PROCESSING (RRY025)
IMAGE PROCESSIG (RRY25) USE OF FT I IMAGE PROCESSIG Optics- originofimperfectionsinimagingsystems(limited resolution/blurring related to 2D FTs)- need to understand using Continuous FT. Sampling -Capturecontinuousimageontoasetofdiscrete
More informationDigital Image Processing
Digital Image Processing Part 2: Image Enhancement Digital Image Processing Course Introduction in the Spatial Domain Lecture AASS Learning Systems Lab, Teknik Room T26 achim.lilienthal@tech.oru.se Course
More informationDESIGN NOTE: DIFFRACTION EFFECTS
NASA IRTF / UNIVERSITY OF HAWAII Document #: TMP-1.3.4.2-00-X.doc Template created on: 15 March 2009 Last Modified on: 5 April 2010 DESIGN NOTE: DIFFRACTION EFFECTS Original Author: John Rayner NASA Infrared
More informationImage Filtering. Median Filtering
Image Filtering Image filtering is used to: Remove noise Sharpen contrast Highlight contours Detect edges Other uses? Image filters can be classified as linear or nonlinear. Linear filters are also know
More informationFast MTF measurement of CMOS imagers using ISO slantededge methodology
Fast MTF measurement of CMOS imagers using ISO 2233 slantededge methodology M.Estribeau*, P.Magnan** SUPAERO Integrated Image Sensors Laboratory, avenue Edouard Belin, 34 Toulouse, France ABSTRACT The
More informationFourier transforms, SIM
Fourier transforms, SIM Last class More STED Minflux Fourier transforms This class More FTs 2D FTs SIM 1 Intensity.5 -.5 FT -1.5 1 1.5 2 2.5 3 3.5 4 4.5 5 6 Time (s) IFT 4 2 5 1 15 Frequency (Hz) ff tt
More informationSolution Set #2
05-78-0 Solution Set #. For the sampling function shown, analyze to determine its characteristics, e.g., the associated Nyquist sampling frequency (if any), whether a function sampled with s [x; x] may
More informationConvolution Pyramids. Zeev Farbman, Raanan Fattal and Dani Lischinski SIGGRAPH Asia Conference (2011) Julian Steil. Prof. Dr.
Zeev Farbman, Raanan Fattal and Dani Lischinski SIGGRAPH Asia Conference (2011) presented by: Julian Steil supervisor: Prof. Dr. Joachim Weickert Fig. 1.1: Gradient integration example Seminar - Milestones
More informationDetermining MTF with a Slant Edge Target ABSTRACT AND INTRODUCTION
Determining MTF with a Slant Edge Target Douglas A. Kerr Issue 2 October 13, 2010 ABSTRACT AND INTRODUCTION The modulation transfer function (MTF) of a photographic lens tells us how effectively the lens
More informationRadiometry I: Illumination. cs348b Matt Pharr
Radiometry I: Illumination cs348b Matt Pharr Administrivia Extra copies of lrt book Bug fix for assignment 1 polynomial.h file Onward To The Physical Description of Light Four key quantities Power Radiant
More informationThe Formation of an Aerial Image, part 2
T h e L i t h o g r a p h y T u t o r (April 1993) The Formation of an Aerial Image, part 2 Chris A. Mack, FINLE Technologies, Austin, Texas In the last issue, we began to described how a projection system
More informationPolarization Experiments Using Jones Calculus
Polarization Experiments Using Jones Calculus Reference http://chaos.swarthmore.edu/courses/physics50_2008/p50_optics/04_polariz_matrices.pdf Theory In Jones calculus, the polarization state of light is
More information2) How fast can we implement these in a system
Filtration Now that we have looked at the concept of interpolation we have seen practically that a "digital filter" (hold, or interpolate) can affect the frequency response of the overall system. We need
More informationCameras. CSE 455, Winter 2010 January 25, 2010
Cameras CSE 455, Winter 2010 January 25, 2010 Announcements New Lecturer! Neel Joshi, Ph.D. Post-Doctoral Researcher Microsoft Research neel@cs Project 1b (seam carving) was due on Friday the 22 nd Project
More informationMidterm 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 informationLecture 3 Review of Signals and Systems: Part 2. EE4900/EE6720 Digital Communications
EE4900/EE6720: Digital Communications 1 Lecture 3 Review of Signals and Systems: Part 2 Block Diagrams of Communication System Digital Communication System 2 Informatio n (sound, video, text, data, ) Transducer
More information>>> from numpy import random as r >>> I = r.rand(256,256);
WHAT IS AN IMAGE? >>> from numpy import random as r >>> I = r.rand(256,256); Think-Pair-Share: - What is this? What does it look like? - Which values does it take? - How many values can it take? - Is it
More informationdigital film technology Scanity multi application film scanner white paper
digital film technology Scanity multi application film scanner white paper standing the test of time multi application film scanner Scanity >>> In the last few years, both digital intermediate (DI) postproduction
More informationMotivation: Image denoising. How can we reduce noise in a photograph?
Linear filtering Motivation: Image denoising How can we reduce noise in a photograph? Moving average Let s replace each pixel with a weighted average of its neighborhood The weights are called the filter
More informationLecture 3, Multirate Signal Processing
Lecture 3, Multirate Signal Processing Frequency Response If we have coefficients of an Finite Impulse Response (FIR) filter h, or in general the impulse response, its frequency response becomes (using
More informationSubband coring for image noise reduction. Edward H. Adelson Internal Report, RCA David Sarnoff Research Center, Nov
Subband coring for image noise reduction. dward H. Adelson Internal Report, RCA David Sarnoff Research Center, Nov. 26 1986. Let an image consisting of the array of pixels, (x,y), be denoted (the boldface
More informationLENSES. INEL 6088 Computer Vision
LENSES INEL 6088 Computer Vision Digital camera A digital camera replaces film with a sensor array Each cell in the array is a Charge Coupled Device light-sensitive diode that converts photons to electrons
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain convolution, convolution theorem, cross-correlation Revision:.3, dated: December 7, 5 Tomáš Svoboda Czech Technical University, Faculty of Electrical Engineering Center
More informationObservational Astronomy
Observational Astronomy Instruments The telescope- instruments combination forms a tightly coupled system: Telescope = collecting photons and forming an image Instruments = registering and analyzing the
More informationOptics of Wavefront. Austin Roorda, Ph.D. University of Houston College of Optometry
Optics of Wavefront Austin Roorda, Ph.D. University of Houston College of Optometry Geometrical Optics Relationships between pupil size, refractive error and blur Optics of the eye: Depth of Focus 2 mm
More informationSIGNALS AND SYSTEMS LABORATORY 13: Digital Communication
SIGNALS AND SYSTEMS LABORATORY 13: Digital Communication INTRODUCTION Digital Communication refers to the transmission of binary, or digital, information over analog channels. In this laboratory you will
More informationDigital Communication System
Digital Communication System Purpose: communicate information at required rate between geographically separated locations reliably (quality) Important point: rate, quality spectral bandwidth, power requirements
More informationPractical Image and Video Processing Using MATLAB
Practical Image and Video Processing Using MATLAB Chapter 10 Neighborhood processing What will we learn? What is neighborhood processing and how does it differ from point processing? What is convolution
More informationMotivation: Image denoising. How can we reduce noise in a photograph?
Linear filtering Motivation: Image denoising How can we reduce noise in a photograph? Moving average Let s replace each pixel with a weighted average of its neighborhood The weights are called the filter
More informationImage Deblurring. This chapter describes how to deblur an image using the toolbox deblurring functions.
12 Image Deblurring This chapter describes how to deblur an image using the toolbox deblurring functions. Understanding Deblurring (p. 12-2) Using the Deblurring Functions (p. 12-5) Avoiding Ringing in
More informationELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet
ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet Lecture 10: Summary Taneli Riihonen 16.05.2016 Lecture 10 in Course Book Sanjit K. Mitra, Digital Signal Processing: A Computer-Based Approach, 4th
More informationECE 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 informationImaging Optics Fundamentals
Imaging Optics Fundamentals Gregory Hollows Director, Machine Vision Solutions Edmund Optics Why Are We Here? Topics for Discussion Fundamental Parameters of your system Field of View Working Distance
More informationCS3291: 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 informationSingle Slit Diffraction
PC1142 Physics II Single Slit Diffraction 1 Objectives Investigate the single-slit diffraction pattern produced by monochromatic laser light. Determine the wavelength of the laser light from measurements
More informationDigital 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 informationPart I Feature Extraction (1) Image Enhancement. CSc I6716 Spring Local, meaningful, detectable parts of the image.
CSc I6716 Spring 211 Introduction Part I Feature Extraction (1) Zhigang Zhu, City College of New York zhu@cs.ccny.cuny.edu Image Enhancement What are Image Features? Local, meaningful, detectable parts
More information1.Discuss the frequency domain techniques of image enhancement in detail.
1.Discuss the frequency domain techniques of image enhancement in detail. Enhancement In Frequency Domain: The frequency domain methods of image enhancement are based on convolution theorem. This is represented
More informationSignals and Filtering
FILTERING OBJECTIVES The objectives of this lecture are to: Introduce signal filtering concepts Introduce filter performance criteria Introduce Finite Impulse Response (FIR) filters Introduce Infinite
More informationImage Enhancement. DD2423 Image Analysis and Computer Vision. Computational Vision and Active Perception School of Computer Science and Communication
Image Enhancement DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 15, 2013 Mårten Björkman (CVAP)
More informationDSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters
Islamic University of Gaza OBJECTIVES: Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters To demonstrate the concept
More informationEnhanced Sample Rate Mode Measurement Precision
Enhanced Sample Rate Mode Measurement Precision Summary Enhanced Sample Rate, combined with the low-noise system architecture and the tailored brick-wall frequency response in the HDO4000A, HDO6000A, HDO8000A
More informationFilip Malmberg 1TD396 fall 2018 Today s lecture
Today s lecture Local neighbourhood processing Convolution smoothing an image sharpening an image And more What is it? What is it useful for? How can I compute it? Removing uncorrelated noise from an image
More informationProjection. Announcements. Müller-Lyer Illusion. Image formation. Readings Nalwa 2.1
Announcements Mailing list (you should have received messages) Project 1 additional test sequences online Projection Readings Nalwa 2.1 Müller-Lyer Illusion Image formation object film by Pravin Bhat http://www.michaelbach.de/ot/sze_muelue/index.html
More informationPerformance Factors. Technical Assistance. Fundamental Optics
Performance Factors After paraxial formulas have been used to select values for component focal length(s) and diameter(s), the final step is to select actual lenses. As in any engineering problem, this
More informationDigital Processing of
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 informationAnalysis 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 informationPHYS 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 informationMultispectral. imaging device. ADVANCED LIGHT ANALYSIS by. Most accurate homogeneity MeasureMent of spectral radiance. UMasterMS1 & UMasterMS2
Multispectral imaging device Most accurate homogeneity MeasureMent of spectral radiance UMasterMS1 & UMasterMS2 ADVANCED LIGHT ANALYSIS by UMaster Ms Multispectral Imaging Device UMaster MS Description
More informationDigital Signal Processing
Digital Signal Processing Lecture 9 Discrete-Time Processing of Continuous-Time Signals Alp Ertürk alp.erturk@kocaeli.edu.tr Analog to Digital Conversion Most real life signals are analog signals These
More information8.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 informationSAMPLING 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 informationSECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS
RADT 3463 - COMPUTERIZED IMAGING Section I: Chapter 2 RADT 3463 Computerized Imaging 1 SECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS RADT 3463 COMPUTERIZED IMAGING Section I: Chapter 2 RADT
More informationVisible Light Communication-based Indoor Positioning with Mobile Devices
Visible Light Communication-based Indoor Positioning with Mobile Devices Author: Zsolczai Viktor Introduction With the spreading of high power LED lighting fixtures, there is a growing interest in communication
More informationSupplementary Figure 1. Effect of the spacer thickness on the resonance properties of the gold and silver metasurface layers.
Supplementary Figure 1. Effect of the spacer thickness on the resonance properties of the gold and silver metasurface layers. Finite-difference time-domain calculations of the optical transmittance through
More informationUnit 1: Image Formation
Unit 1: Image Formation 1. Geometry 2. Optics 3. Photometry 4. Sensor Readings Szeliski 2.1-2.3 & 6.3.5 1 Physical parameters of image formation Geometric Type of projection Camera pose Optical Sensor
More informationCSE 473/573 Computer Vision and Image Processing (CVIP)
CSE 473/573 Computer Vision and Image Processing (CVIP) Ifeoma Nwogu inwogu@buffalo.edu Lecture 4 Image formation(part I) Schedule Last class linear algebra overview Today Image formation and camera properties
More informationComputer 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 informationImplementation of Digital Signal Processing: Some Background on GFSK Modulation
Implementation of Digital Signal Processing: Some Background on GFSK Modulation Sabih H. Gerez University of Twente, Department of Electrical Engineering s.h.gerez@utwente.nl Version 5 (March 9, 2016)
More informationLab 3.0. Pulse Shaping and Rayleigh Channel. Faculty of Information Engineering & Technology. The Communications Department
Faculty of Information Engineering & Technology The Communications Department Course: Advanced Communication Lab [COMM 1005] Lab 3.0 Pulse Shaping and Rayleigh Channel 1 TABLE OF CONTENTS 2 Summary...
More informationSignals 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 informationThe Noise about Noise
The Noise about Noise I have found that few topics in astrophotography cause as much confusion as noise and proper exposure. In this column I will attempt to present some of the theory that goes into determining
More informationCS 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 informationFrequency 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 informationFIR FILTER DESIGN USING A NEW WINDOW FUNCTION
FIR FILTER DESIGN USING A NEW WINDOW FUNCTION Mahroh G. Shayesteh and Mahdi Mottaghi-Kashtiban, Department of Electrical Engineering, Urmia University, Urmia, Iran Sonar Seraj System Cor., Urmia, Iran
More informationdigital film technology Resolution Matters what's in a pattern white paper standing the test of time
digital film technology Resolution Matters what's in a pattern white paper standing the test of time standing the test of time An introduction >>> Film archives are of great historical importance as they
More informationOptical transfer function shaping and depth of focus by using a phase only filter
Optical transfer function shaping and depth of focus by using a phase only filter Dina Elkind, Zeev Zalevsky, Uriel Levy, and David Mendlovic The design of a desired optical transfer function OTF is a
More informationReflectors vs. Refractors
1 Telescope Types - Telescopes collect and concentrate light (which can then be magnified, dispersed as a spectrum, etc). - In the end it is the collecting area that counts. - There are two primary telescope
More information