Digital Image Processing. Frequency Domain Filtering
|
|
- Simon Day
- 5 years ago
- Views:
Transcription
1 Digital Image Processing Frequency Domain Filtering
2 DFT Matlab demo clear all; close all; a=imread('testpat1.png');b=imdouble(a); figure;imshow(b); Fb = fft(b);fbshift=fftshift(fb); figure;imshow(log(abs(fbshift) ),[]); FMask=zeros(56,56);FMask([96:160],[96:160])=1.0; Fbband=Fbshift.*FMask; figure;imshow(log(abs(fbband) ),[]); Fbband=ifftshift(Fbband);bandrebuilt=ifft(Fbband); figure;imshow(bandrebuilt); FMask=zeros(56,56);FMask([64:19],[64:19])=1.0; Fbband=Fbshift.*FMask; figure;imshow(log(abs(fbband) ),[]); Fbband=ifftshift(Fbband);bandrebuilt=ifft(Fbband); figure;imshow(bandrebuilt); FMask=zeros(56,56);FMask([3:4],[3:4])=1.0; Fbband=Fbshift.*FMask; figure;imshow(log(abs(fbband) ),[]); Fbband=ifftshift(Fbband);bandrebuilt=ifft(Fbband); figure;imshow(bandrebuilt);
3 FT Signal decomposition(coefficient computation) IFT Signal reconstruction
4 IFT IFT
5 The central part of FT, i.e. the low frequency components are responsible for the general gray-level appearance of an image. The high frequency components of FT are responsible for the detail information of an image.
6 Image Frequency Domain (log magnitude) v Detail u General appearance
7 5 % 10 % 0 % 50 %
8 Image Filtering Image filtering techniques: Spatial domain methods Frequency domain methods Spatial (time) domain techniques are techniques that operate directly on pixels. Frequency domain techniques are based on modifying the Fourier transform of an image.
9 Frequency Domain Filtering
10 Frequency Domain Filtering Edges and sharp transitions (e.g., noise) in an image contribute significantly to highfrequency content of FT. Low frequency contents in the FT are responsible to the general appearance of the image over smooth areas. Blurring (smoothing) is achieved by attenuating range of high frequency components of FT.
11 Convolution in Time Domain g(x,y)=h(x,y) f(x,y) M 1 M 1 g( x, y) = h( x', y' ) f( x x', y y' ) x' = 0 y' = 0 f( x, y)* h( x, y) f(x,y) is the input image g(x,y) is the filtered image h(x,y): impulse response
12 Convolution Theorem G(u,v)=F(u,v) H(u,v) g(x,y)= f(x,y) h(x,y) Multiplication in Frequency Domain Convolution in Time Domain Filtering in Frequency Domain with H(u,v) is equivalent to filtering in Spatial Domain with h(x,y).
13 blue line = sum of 3 sinusoids (0, 50, and 80 Hz) + random noise red line = sum of 3 sinusoids without noise blue line = sum of 3 sinusoids after filtering in time domain 1x average [ ] / 5 blue line = sum of 3 sinusoids after filtering in frequency domain cut-off 90 Hz
14 Convolution Property of the Fourier Transform Let functions f ( r, c) and g( r, c) have Fourier Transforms F( u, v) and G( u, v). Then, F { f Moreover, F { f g} = g} = F G. F G. * = convolution = multiplication The The Fourier Transform of of a convolution equals the the product of of the the Fourier Transforms. Similarly, the the Fourier Transform of of a multiplication is is the the convolution of of the the Fourier Transforms
15 Convolution via Fourier Transform Image & Mask Transforms Pixel-wise Product Inverse Transform
16 How to Convolve via FT in Matlab 1. Read the image from a file into a variable, say I.. Read in or create the convolution mask, h. 3. Compute the sum of the mask: s = sum(sum(h)); 4. If s == 0, set s = 1; 5. Create: H = zeros(size(i)); 6. Copy h into the middle of H. 7. Shift H into position: H = ifftshift(h); 8. Take the D FT of I and H: FI=fft(I); FH=fft(H); 9. Pointwise multiply the FTs: FJ=FI.*FH; 10. Compute the inverse FT: J = real(ifft(fj)); 11. Normalize the result: J = J/s; The The mask mask is is usually usually 1-band 1-band For For color color images images you you may may need need to to do do each each step step for for each each band band separately.
17 Coordinate Origin of the FFT Center Center = (floor(r/)+1, (floor(r/)+1, floor(c/)+1) floor(c/)+1) Even Odd Even Odd Image Origin Image Origin Weight Matrix Origin Weight Matrix Origin After FFT shift After FFT shift After IFFT shift After IFFT shift
18 Matlab s fftshift and ifftshift J = fftshift(i): I (1,1) J ( R/ +1, C/ +1) I = ifftshift(j): J ( R/ +1, C/ +1) I (1,1) where x = floor(x) = the largest integer smaller than x.
19 Algorithm Complexity We can compute the DFT directly using the formula An N point DFT would require N floating point multiplications per output point Since there are N output points, the computational complexity of the DFT is N 4 N 4 =4x10 9 for N=56 Bad news! Many hours on a workstation
20 Algorithm Complexity The FFT algorithm was developed in the 60 s for seismic exploration Reduced the DFT complexity to N log N N log N~10 6 for N=56 A few seconds on a workstation F( uv, ) = F( u, v) * original image Fourier log magnitude Fourier phase
21 Examples of Filters Frequency domain Gaussian lowpass filter Gaussian highpass filter Spatial domain
22 Blurring: Averaging / Lowpass Filtering Blurring results from: Pixel averaging in the spatial domain: Each pixel in the output is a weighted average of its neighbors. Is a convolution whose weight matrix sums to 1. Lowpass filtering in the frequency domain: High frequencies are diminished or eliminated Individual frequency components are multiplied by a nonincreasing function of ω such as 1/ω = 1/ (u +v ).
23 Ideal Lowpass Filter 理想低通滤波器 Image Image size: size: 51x51 51x51 FD FD filter filter radius: radius: Multiply by this, or convolve by this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
24 Ideal Lowpass Filter Image Image size: size: 51x51 51x51 FD FD filter filter radius: radius: 8 Multiply by this, or convolve by this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
25 Power Spectrum and Phase of an Image Consider the the image below: Original Image Image Power Power Spectrum Phase Phase
26 Ideal Lowpass Filter Image Image size: size: 51x51 51x51 FD FD filter filter radius: radius: Original Image Image Power Power Spectrum Ideal Ideal LPF LPF in in FD FD
27 Ideal Lowpass Filter Image Image size: size: 51x51 51x51 FD FD filter filter radius: radius: Filtered Image Image Filtered Power Power Spectrum Original Image Image
28 Ideal low-pass filter (ILPF) H( u, v) 1 = 0 D( u, v) D D ( u, v) > D 0 0 1/ Duv (, ) = [( u M/ ) + ( v N/ ) ] (M/,N/): center in frequency domain. D 0 is called the cutoff frequency( 截止频率 ).
29 FT Ideal in frequency domain means non-ideal in spatial domain, vice versa. ringing and blurring
30 Shape of ILPF Frequency domain h(x,y) Spatial domain
31 Fourier transform basis functions Approximating a square wave as the sum of sine waves.
32 F(u,v) H(u,v) f(x,y) h(x,y) The limiting case Frequency domain Spatial Domain H(u,v)=1 h(x,y)= And vice versa
33 Treating the ringing 譬如人染沉疴, 当先用糜粥以饮之, 和药以服之 ; 待其腑脏调和, 形体渐安, 然后用肉食以补之, 猛药以治之 : 则病根尽去, 人得全生也 若不待气脉和缓, 便投以猛药厚味, 欲求安保, 诚为难矣 三国演义 第四十三回
34 Butterworth Lowpass Filters (BLPF) Smooth transfer function, no sharp discontinuity, no clear cutoff frequency. H ( u, v) = n 1+ 1 D( u, v) D0 1
35 Butterworth Lowpass Filters (BLPF) h(x) n=1 n= n=5 n=0 Order=1 Order= Order=5 Order=0 H ( u, v) = n 1+ 1 D( u, v) D0 Spatial representations of BLPF with different orders Notice the ringing of different orders
36 No serious ringing artifacts
37 Gaussian Lowpass Filters (GLPF) Smooth transfer function, smooth impulse response, no ringing H ( u, v) = e D ( u, v) D 0
38 GLPF Frequency domain Gaussian lowpass filter Spatial domain
39 No ringing artifacts
40 Examples of Lowpass Filtering
41 Examples of Lowpass Filtering Low-pass filter H(u,v) Original image and its FT Filtered image and its FT
42 Sharpening: Differencing / Highpass Filtering Sharpening results from adding to the image, a copy of itself that has been: Pixel-differenced in the spatial domain: Each pixel in the output is a difference between itself and a weighted average of its neighbors. Is a convolution whose weight matrix sums to 0. Highpass filtered in the frequency domain: High frequencies are enhanced or amplified. Individual frequency components are multiplied by an increasing function of ω such as αω = α (u +v ), where α is a constant.
43 Ideal Highpass Filter 理想高通滤波器 Image Image size: size: 51x51 51x51 FD FD notch notch radius: radius: Multiply Multiply by by this, this, or or convolve convolve by by this this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
44 Ideal Highpass Filter Image Image size: size: 51x51 51x51 FD FD notch notch radius: radius: Original Image Image Power Power Spectrum Ideal Ideal HPF HPF in in FD FD
45 Ideal Highpass Filter Image Image size: size: 51x51 51x51 FD FD notch notch radius: radius: Filtered Image Image * * Filtered Power Power Spectrum Original Image Image
46 High-pass Filters 高通滤波器 H hp (u,v)=1-h lp (u,v) Ideal: Butterworth: Gaussian: n v D u D v u H 0 ), ( 1 1 ), ( + = ) = 0 1, ( v H u 0 ), ( D v u D > 0 ), ( D v u D 0 )/, ( 1 ), ( D v u D e v u H =
47
48 h(x,y) High-pass Filters
49 Ideal High-pass Filtering ringing artifacts
50 Butterworth High-pass Filtering
51 Gaussian High-pass Filtering
52 Filtered image and its FT Gaussian High-pass Filtering Original image Gaussian filter H(u,v)
53 Laplacian in Frequency Domain ) ( ), ( 1 v u v u H + = ), ( ) ( ] ), ( ), ( [ v u F v u y y x f x y x f + = + I y f x f f + = Laplacian operator Spatial domain Frequency domain
54 The Uncertainty Relation 不确定性关系 space frequency If Δx Δ y is the extent of FT FT the object in space and if Δu Δv is its extent in frequency then, space frequency Δx Δ y Δu Δv 1 16π FT FT A small object in in space has has a large frequency extent and and vice-versa.
55 Ideal Filters Do Not Produce Ideal Results IFT IFT A sharp cutoff in in the the frequency domain causes ringing in in the the spatial domain.
56 Ideal Filters Do Not Produce Ideal Results Ideal LPF LPF Blurring the the image above with with an an ideal lowpass filter distorts the the results with with ringing or or ghosting.
57 Optimal Filter: The Gaussian IFT IFT The The Gaussian filter optimizes the the uncertainty relation. It It provides the the sharpest cutoff with with the the least ringing.
58 One-Dimensional Gaussian g( x) = 1 σ π e ( x μ ) σ
59 Two-Dimensional Gaussian ) ( ) ( 1 ) ( ) ( 1 ) ( ) ( ), ( c r c r r c c r c c r r c r y x y r e e c g r g c r g σ σ μ σ μ σ π σ σ σ μ σ μ π σ σ + = = = ) ( ) ( 1 ), ( σ μ μ π σ + = c r e c r g If μ and σ are different for r & c If μ and σ are different for r & c or if μ and σ are the same for r & c. or if μ and σ are the same for r & c. r c R = 51, C = 51 μ = 57, σ = 64
60 Optimal Filter: The Gaussian Gaussian LPF LPF With a gaussian lowpass filter, the the image above is is blurred without ringing or or ghosting.
61 Gaussian Lowpass Filter 高斯低通滤波器 Image Image size: size: 51x51 51x51 SD SD filter filter sigma sigma = = 8 Multiply Multiply by by this, this, or or convolve convolve by by this this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
62 Gaussian Lowpass Filter Image Image size: size: 51x51 51x51 SD SD filter filter sigma sigma = = Multiply Multiply by by this, this, or or convolve convolve by by this this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
63 Gaussian Lowpass Filter Image Image size: size: 51x51 51x51 SD SD filter filter sigma sigma = = 8 Original Image Image Power Power Spectrum Gaussian LPF LPF in in FD FD
64 Gaussian Lowpass Filter Image Image size: size: 51x51 51x51 SD SD filter filter sigma sigma = = 8 Filtered Image Image Filtered Power Power Spectrum Original Image Image
65 Comparison of Ideal and Gaussian Lowpass Filters Ideal Ideal LPF LPF Original Image Image Gaussian LPF LPF
66 Gaussian Highpass Filter 高斯高通滤波器 Image Image size: size: 51x51 51x51 FD FD notch notch sigma sigma = = 8 Multiply Multiply by by this, this, or or convolve convolve by by this this Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
67 Gaussian Highpass Filter Image Image size: size: 51x51 51x51 FD FD notch notch sigma sigma = = 8 Original Image Image Power Power Spectrum Gaussian HPF HPF in in FD FD
68 Gaussian Highpass Filter Image Image size: size: 51x51 51x51 FD FD notch notch sigma sigma = = 8 Filtered Image Image * * Filtered Power Power Spectrum Original Image Image
69 Comparison of Ideal and Gaussian Highpass Filters * * signed signed image; image; 0 mapped mapped to to Ideal Ideal HPF HPF * * Original Image Image Gaussian HPF HPF * *
70 Another Highpass Filter * * signed signed image; image; 0 mapped mapped to to original image image filter filter power power spectrum filtered image image * *
71 Ideal Bandpass Filter 理想带通滤波器 * * signed signed image; image; 0 mapped mapped to to original image image filter filter power power spectrum filtered image image * *
72 Gaussian Bandpass Filter 高斯带通滤波器 Image Image size: size: 51x51 51x51 sigma sigma = = --sigma sigma = = 8 Fourier Domain Rep. Rep. Spatial Spatial Representation Central Profile Profile
73 Gaussian Bandpass Filter Image Image size: size: 51x51 51x51 sigma sigma = = --sigma sigma = = 8 Original Image Image Power Power Spectrum Gaussian BPF BPF in in FD FD
74 Gaussian Bandpass Filter Image Image size: size: 51x51 51x51 sigma sigma = = --sigma sigma = = 8 Filtered Image Image * * Filtered Power Power Spectrum Original Image Image * * signed signed image; image; 0 mapped mapped to to 18 18
75 Comparison of Ideal and Gaussian Bandpass Filters * * signed signed image; image; 0 mapped mapped to to Ideal Ideal BPF BPF * * Original Image Image Gaussian BPF BPF * *
76 DFT Matlab demo clear all;close all; iptsetpref('imshowborder','tight'); bw = imread('text.png'); imshow(bw); a = bw(3:45,88:98); figure, imshow(a); C = real(ifft(fft(bw).* fft(rot90(a,),56,56))); figure, imshow(c,[]) max(c(:)); %cys this call output 68 thresh = 60; % Set a threshold a little less than max. figure, imshow(c > thresh); % highlight pixels over threshold.
77 Homework VII Design your own DFT exploration experimence. It can be either of the following DFT and reconstruction DFT Real/Imagenary part, magnitude/phase Frequency domain filtering DFT based image analysis (Detection/Editing ) Submit your test images, codes, experiment results and documents.
78 End of Lecture
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 informationFourier Transforms and the Frequency Domain
Fourier Transforms and the Frequency Domain Lecture 11 Magnus Gedda magnus.gedda@cb.uu.se Centre for Image Analysis Uppsala University Computer Assisted Image Analysis 04/27/2006 Gedda (Uppsala University)
More informationDigital Image Processing
Digital Image Processing Filtering in the Frequency Domain (Application) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering 2 Periodicity of the
More informationCoE4TN4 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 informationDigital Image Processing. Filtering in the Frequency Domain (Application)
Digital Image Processing Filtering in the Frequency Domain (Application) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Periodicity of the DFT The range
More informationTransforms 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 informationTDI2131 Digital Image Processing
TDI131 Digital Image Processing Frequency Domain Filtering Lecture 6 John See Faculty of Information Technology Multimedia University Some portions of content adapted from Zhu Liu, AT&T Labs. Most figures
More informationLecture #10. EECS490: Digital Image Processing
Lecture #10 Wraparound and padding Image Correlation Image Processing in the frequency domain A simple frequency domain filter Frequency domain filters High-pass, low-pass Apodization Zero-phase filtering
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 informationMidterm Review. Image Processing CSE 166 Lecture 10
Midterm Review Image Processing CSE 166 Lecture 10 Topics covered Image acquisition, geometric transformations, and image interpolation Intensity transformations Spatial filtering Fourier transform and
More informationDIGITAL IMAGE PROCESSING UNIT III
DIGITAL IMAGE PROCESSING UNIT III 3.1 Image Enhancement in Frequency Domain: Frequency refers to the rate of repetition of some periodic events. In image processing, spatial frequency refers to the variation
More informationDigital Image Processing. Image Enhancement: Filtering in the Frequency Domain
Digital Image Processing Image Enhancement: Filtering in the Frequency Domain 2 Contents In this lecture we will look at image enhancement in the frequency domain Jean Baptiste Joseph Fourier The Fourier
More informationSmoothing frequency domain filters
Smoothing frequency domain filters Ideal Lowpass Filter (ILPF) ILPF is the simplest lowpass filter that cuts off all high frequency components of the DFT that are at a distance greater than a specified
More informationSmoothing frequency domain filters
Smoothing frequency domain filters Ideal Lowpass Filter (ILPF) ILPF is the simplest lowpass filter that cuts off all high frequency components of the DFT that are at a distance greater than a specified
More informationThe Fourier Transform
The Fourier Transform Introduction to Digital Signal Processing (886457) 6 1 / 56 Contents Introduction Fourier Transforms One-dimensional DFT Two-dimensional DFT Fourier Transforms Function in Octave
More informationImage acquisition. Midterm Review. Digitization, line of image. Digitization, whole image. Geometric transformations. Interpolation 10/26/2016
Image acquisition Midterm Review Image Processing CSE 166 Lecture 10 2 Digitization, line of image Digitization, whole image 3 4 Geometric transformations Interpolation CSE 166 Transpose these matrices
More informationLAB MANUAL SUBJECT: IMAGE PROCESSING BE (COMPUTER) SEM VII
LAB MANUAL SUBJECT: IMAGE PROCESSING BE (COMPUTER) SEM VII IMAGE PROCESSING INDEX CLASS: B.E(COMPUTER) SR. NO SEMESTER:VII TITLE OF THE EXPERIMENT. 1 Point processing in spatial domain a. Negation of an
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 informationImage Smoothening and Sharpening using Frequency Domain Filtering Technique
Volume 5, Issue 4, April (17) Image Smoothening and Sharpening using Frequency Domain Filtering Technique Swati Dewangan M.Tech. Scholar, Computer Networks, Bhilai Institute of Technology, Durg, India.
More informationFiltering in the spatial domain (Spatial Filtering)
Filtering in the spatial domain (Spatial Filtering) refers to image operators that change the gray value at any pixel (x,y) depending on the pixel values in a square neighborhood centered at (x,y) using
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 information8. Lecture. Image restoration: Fourier domain
8. Lecture Image restoration: Fourier domain 1 Structured noise 2 Motion blur 3 Filtering in the Fourier domain ² Spatial ltering (average, Gaussian,..) can be done in the Fourier domain (convolution theorem)
More informationFourier 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 informationIMAGE PROCESSING: AREA OPERATIONS (FILTERING)
IMAGE PROCESSING: AREA OPERATIONS (FILTERING) N. C. State University CSC557 Multimedia Computing and Networking Fall 2001 Lecture # 13 IMAGE PROCESSING: AREA OPERATIONS (FILTERING) N. C. State University
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 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 informationDiscrete Fourier Transform (DFT)
Amplitude Amplitude Discrete Fourier Transform (DFT) DFT transforms the time domain signal samples to the frequency domain components. DFT Signal Spectrum Time Frequency DFT is often used to do frequency
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 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 informationComputer 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 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 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 informationLaboratory Assignment 4. Fourier Sound Synthesis
Laboratory Assignment 4 Fourier Sound Synthesis PURPOSE This lab investigates how to use a computer to evaluate the Fourier series for periodic signals and to synthesize audio signals from Fourier series
More informationProf. 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 informationImage Enhancement using Histogram Equalization and Spatial Filtering
Image Enhancement using Histogram Equalization and Spatial Filtering Fari Muhammad Abubakar 1 1 Department of Electronics Engineering Tianjin University of Technology and Education (TUTE) Tianjin, P.R.
More informationInvestigation of Optimal Denoising Filter for MRI Images
International Journal of Applied Engineering Research ISSN 0973-456 Volume 13, Number 15 (018) pp. 164-171 Investigation of Optimal Denoising Filter for MRI Images Ch. Rajasekhara Rao, M N V S S Kumar,
More informationAutomatic processing to restore data of MODIS band 6
Automatic processing to restore data of MODIS band 6 --Final Project for ECE 533 Abstract An automatic processing to restore data of MODIS band 6 is introduced. For each granule of MODIS data, 6% of the
More informationImage Filtering. Reading Today s Lecture. Reading for Next Time. What would be the result? Some Questions from Last Lecture
Image Filtering HCI/ComS 575X: Computational Perception Instructor: Alexander Stoytchev http://www.cs.iastate.edu/~alex/classes/2007_spring_575x/ January 24, 2007 HCI/ComS 575X: Computational Perception
More informationCoE4TN4 Image Processing. Chapter 3: Intensity Transformation and Spatial Filtering
CoE4TN4 Image Processing Chapter 3: Intensity Transformation and Spatial Filtering Image Enhancement Enhancement techniques: to process an image so that the result is more suitable than the original image
More informationIntroduce cascaded first-order op-amp filters. Faculty of Electrical and Electronic Engineering
Yıldız Technical University Cascaded FirstOrder Filters Introduce cascaded first-order op-amp filters Faculty of Electrical and Electronic Engineering Lesson Objectives Introduce cascaded filters Introduce
More informationIMAGE PROCESSING (RRY025) THE CONTINUOUS 2D FOURIER TRANSFORM
IMAGE PROCESSING (RRY5) THE CONTINUOUS D FOURIER TRANSFORM INTRODUCTION A vital tool in image processing. Also a prototype of other image transforms, cosine, Wavelet etc. Applications Image Filtering -
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 informationSYLLABUS CHAPTER - 2 : INTENSITY TRANSFORMATIONS. Some Basic Intensity Transformation Functions, Histogram Processing.
Contents i SYLLABUS UNIT - I CHAPTER - 1 : INTRODUCTION TO DIGITAL IMAGE PROCESSING Introduction, Origins of Digital Image Processing, Applications of Digital Image Processing, Fundamental Steps, Components,
More informationSystem 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 informationImage Restoration. Lecture 7, March 23 rd, Lexing Xie. EE4830 Digital Image Processing
Image Restoration Lecture 7, March 23 rd, 2009 Lexing Xie EE4830 Digital Image Processing http://www.ee.columbia.edu/~xlx/ee4830/ thanks to G&W website, Min Wu and others for slide materials 1 Announcements
More informationExamples of image processing
Examples of image processing Example 1: We would like to automatically detect and count rings in the image 3 Detection by correlation Correlation = degree of similarity Correlation between f(x, y) and
More informationExercise 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 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 informationFinal Exam Solutions June 14, 2006
Name or 6-Digit Code: PSU Student ID Number: Final Exam Solutions June 14, 2006 ECE 223: Signals & Systems II Dr. McNames Keep your exam flat during the entire exam. If you have to leave the exam temporarily,
More informationHead, IICT, Indus University, India
International Journal of Emerging Research in Management &Technology Research Article December 2015 Comparison Between Spatial and Frequency Domain Methods 1 Anuradha Naik, 2 Nikhil Barot, 3 Rutvi Brahmbhatt,
More informationLecture - 10 Image Enhancement in the Frequency Domain
Lectre - Image Enhancement in the Freqenc Domain Cosimo Distante Backgrond An fnction that periodicall repeats itself can be epressed as the sm of sines and/or cosines of different freqencies each mltiplied
More informationImage filtering, image operations. Jana Kosecka
Image filtering, image operations Jana Kosecka - photometric aspects of image formation - gray level images - point-wise operations - linear filtering Image Brightness values I(x,y) Images Images contain
More informationMatlab (see Homework 1: Intro to Matlab) Linear Filters (Reading: 7.1, ) Correlation. Convolution. Linear Filtering (warm-up slide) R ij
Matlab (see Homework : Intro to Matlab) Starting Matlab from Unix: matlab & OR matlab nodisplay Image representations in Matlab: Unsigned 8bit values (when first read) Values in range [, 255], = black,
More informationEEL 6562 Image Processing and Computer Vision Image Restoration
DEPARTMENT OF ELECTRICAL & COMPUTER ENGINEERING EEL 6562 Image Processing and Computer Vision Image Restoration Rajesh Pydipati Introduction Image Processing is defined as the analysis, manipulation, storage,
More informationDigital Image Processing 3/e
Laboratory Projects for Digital Image Processing 3/e by Gonzalez and Woods 2008 Prentice Hall Upper Saddle River, NJ 07458 USA www.imageprocessingplace.com The following sample laboratory projects are
More informationChapter 2 Image Enhancement in the Spatial Domain
Chapter 2 Image Enhancement in the Spatial Domain Abstract Although the transform domain processing is essential, as the images naturally occur in the spatial domain, image enhancement in the spatial domain
More informationLecture 12: Image Processing and 2D Transforms
Lecture 12: Image Processing and 2D Transforms Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu October 18, 2005 Abstract The Fourier transform
More informationChapter 4 SPEECH ENHANCEMENT
44 Chapter 4 SPEECH ENHANCEMENT 4.1 INTRODUCTION: Enhancement is defined as improvement in the value or Quality of something. Speech enhancement is defined as the improvement in intelligibility and/or
More informationFFT analysis in practice
FFT analysis in practice Perception & Multimedia Computing Lecture 13 Rebecca Fiebrink Lecturer, Department of Computing Goldsmiths, University of London 1 Last Week Review of complex numbers: rectangular
More informationDigital Signal Processing
Digital Signal Processing System Analysis and Design Paulo S. R. Diniz Eduardo A. B. da Silva and Sergio L. Netto Federal University of Rio de Janeiro CAMBRIDGE UNIVERSITY PRESS Preface page xv Introduction
More informationFinal Exam Solutions June 7, 2004
Name: Final Exam Solutions June 7, 24 ECE 223: Signals & Systems II Dr. McNames Write your name above. Keep your exam flat during the entire exam period. If you have to leave the exam temporarily, close
More informationA.V.C. COLLEGE OF ENGINEERING DEPARTEMENT OF CSE CP7004- IMAGE PROCESSING AND ANALYSIS UNIT 1- QUESTION BANK
A.V.C. COLLEGE OF ENGINEERING DEPARTEMENT OF CSE CP7004- IMAGE PROCESSING AND ANALYSIS UNIT 1- QUESTION BANK STAFF NAME: TAMILSELVAN K UNIT I SPATIAL DOMAIN PROCESSING Introduction to image processing
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 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 informationLow Spatial Frequency Noise Reduction with Applications to Light Field Moment Imaging
Low Spatial Frequency Noise Reduction with Applications to Light Field Moment Imaging Christopher Madsen Stanford University cmadsen@stanford.edu Abstract This project involves the implementation of multiple
More informationFrequency Domain Representation of Signals
Frequency Domain Representation of Signals The Discrete Fourier Transform (DFT) of a sampled time domain waveform x n x 0, x 1,..., x 1 is a set of Fourier Coefficients whose samples are 1 n0 X k X0, X
More informationDigital Image Processing
Digital Image Processing Lecture # 5 Image Enhancement in Spatial Domain- I ALI JAVED Lecturer SOFTWARE ENGINEERING DEPARTMENT U.E.T TAXILA Email:: ali.javed@uettaxila.edu.pk Office Room #:: 7 Presentation
More informationImage Restoration. Lecture 7, March 23 rd, Lexing Xie. EE4830 Digital Image Processing
Image Restoration Lecture 7, March 23 rd, 2008 Lexing Xie EE4830 Digital Image Processing http://www.ee.columbia.edu/~xlx/ee4830/ thanks to G&W website, Min Wu and others for slide materials 1 Announcements
More informationSampling 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 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 information2D 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 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 informationProf. Vidya Manian Dept. of Electrical and Comptuer Engineering
Image Processing Intensity Transformations Chapter 3 Prof. Vidya Manian Dept. of Electrical and Comptuer Engineering INEL 5327 ECE, UPRM Intensity Transformations 1 Overview Background Basic intensity
More informationChapter 2: Signal Representation
Chapter 2: Signal Representation Aveek Dutta Assistant Professor Department of Electrical and Computer Engineering University at Albany Spring 2018 Images and equations adopted from: Digital Communications
More informationTable 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 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 informationVision 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 informationIntroduction to Wavelet Transform. Chapter 7 Instructor: Hossein Pourghassem
Introduction to Wavelet Transform Chapter 7 Instructor: Hossein Pourghassem Introduction Most of the signals in practice, are TIME-DOMAIN signals in their raw format. It means that measured signal is a
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 informationLecture 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 informationCSE 564: Visualization. Image Operations. Motivation. Provide the user (scientist, t doctor, ) with some means to: Global operations:
Motivation CSE 564: Visualization mage Operations Klaus Mueller Computer Science Department Stony Brook University Provide the user (scientist, t doctor, ) with some means to: enhance contrast of local
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad
INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 500 043 ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK Course Title Course Code Class Branch DIGITAL IMAGE PROCESSING A70436 IV B. Tech.
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Filter Design Circularly symmetric 2-D low-pass filter Pass-band radial frequency: ω p Stop-band radial frequency: ω s 1 δ p Pass-band tolerances: δ
More information06: 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 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 informationDefocusing and Deblurring by Using with Fourier Transfer
Defocusing and Deblurring by Using with Fourier Transfer AKIRA YANAGAWA and TATSUYA KATO 1. Introduction Image data may be obtained through an image system, such as a video camera or a digital still camera.
More informationAPPENDIX A to VOLUME A1 TIMS FILTER RESPONSES
APPENDIX A to VOLUME A1 TIMS FILTER RESPONSES A2 TABLE OF CONTENTS... 5 Filter Specifications... 7 3 khz LPF (within the HEADPHONE AMPLIFIER)... 8 TUNEABLE LPF... 9 BASEBAND CHANNEL FILTERS - #2 Butterworth
More informationLaboratory Assignment 5 Amplitude Modulation
Laboratory Assignment 5 Amplitude Modulation PURPOSE In this assignment, you will explore the use of digital computers for the analysis, design, synthesis, and simulation of an amplitude modulation (AM)
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 informationAnalysis of Image Enhancement Techniques Used in Remote Sensing Satellite Imagery
Analysis of Image Enhancement Techniques Used in Remote Sensing Satellite Imagery Kriti Bajpai MTech Student, Department of Computer Science Engineering, Gyan Ganga Institute of Technology & Sciences,
More informationCS534 Introduction to Computer Vision. Linear Filters. Ahmed Elgammal Dept. of Computer Science Rutgers University
CS534 Introduction to Computer Vision Linear Filters Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What are Filters Linear Filters Convolution operation Properties of Linear Filters
More informationCG401 Advanced Signal Processing. Dr Stuart Lawson Room A330 Tel: January 2003
CG40 Advanced Dr Stuart Lawson Room A330 Tel: 23780 e-mail: ssl@eng.warwick.ac.uk 03 January 2003 Lecture : Overview INTRODUCTION What is a signal? An information-bearing quantity. Examples of -D and 2-D
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 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 informationSignal processing preliminaries
Signal processing preliminaries ISMIR Graduate School, October 4th-9th, 2004 Contents: Digital audio signals Fourier transform Spectrum estimation Filters Signal Proc. 2 1 Digital signals Advantages of
More informationImage analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror
Image analysis CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror 1 Outline Images in molecular and cellular biology Reducing image noise Mean and Gaussian filters Frequency domain interpretation
More informationImage analysis. CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror
Image analysis CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror A two- dimensional image can be described as a function of two variables f(x,y). For a grayscale image, the value of f(x,y) specifies the brightness
More informationDigital Imaging Systems for Historical Documents
Digital Imaging Systems for Historical Documents Improvement Legibility by Frequency Filters Kimiyoshi Miyata* and Hiroshi Kurushima** * Department Museum Science, ** Department History National Museum
More informationNeurophysiology. The action potential. Why should we care? AP is the elemental until of nervous system communication
Neurophysiology Why should we care? AP is the elemental until of nervous system communication The action potential Time course, propagation velocity, and patterns all constrain hypotheses on how the brain
More informationECE 440L. Experiment 1: Signals and Noise (1 week)
ECE 440L Experiment 1: Signals and Noise (1 week) I. OBJECTIVES Upon completion of this experiment, you should be able to: 1. Use the signal generators and filters in the lab to generate and filter noise
More information