Image processing in the spatial and frequency domain
|
|
- Peregrine Foster
- 5 years ago
- Views:
Transcription
1 Digital Media Lab Applied physics and electronics Ulrik Söderström Image processing in the spatial and frequency domain Fourier transform, Filtering and Image enhancement
2 Part 1. The Fourier Transform 1.1. Fourier analysis in two dimensions We are going to look at some simple signals in both the spatial- and the frequency domain. In this exercise we are using gray scale images with the size 64 x 64 pixels. Each pixel has a value between 0 and 1. Sometimes you have to rescale the resulting image to come back to values between one and zero. One way to get information about pixel values is the function "hist", which calculates the frequency for each pixel value. Let's start with a "one dimensional" sine wave. In the following code a 64 x 64 pixel image with two periods of a horizontal sine wave is created. The image is shown with "imshow", so we have to rescale the values to lie between 0 and 1 before displaying it. The signal is centred on the x-axis before the Fourier transform is calculated. (fft, fft2, fftshift) clear; %the image will be in format Q x Q Q =64; %number of periods period=2; %orginal image for m=1:1:q for n=1:1:q orig_image(m,n)= *sin(2.*pi.*(n./(q./period))); %the "0.5" constants give values between 0 and 1 for use with "imshow" subplot(1,2,1),imshow(orig_image);title('original image'); %delete the added scaling values before FFT, gives values symmetrically around the x- axis A=orig_image-.5; B=abs(fftshift(fft2(A))); subplot(1,2,2),imshow(b);title('magnitude of FFT'); Try to explain the result! What happens with the spectrum is the sine wave is vertical? Create a black image of size 64 x 64 with a 9 x 9 white square in the middle! In this case you might have to rescale the frequency image before you show the result. Try to use a log-transform (log(abs(fftshift(fft2(bild)))+1)), or rescale the image with the formula (bild2=(bildmin(min(bild)))./(max(max(bild))-min(min(bild)))). Is the spectrum looking like you expected? Describe the result Repeat the previous experiment but start with an image with zeros all over except one pixel with value 1.
3 1.2. The importance of object size Go back to problem but make the square 15 x 15 instead of 9 x How is spectrum changed? Explain the correspondence between size in spatial and frequency plane Fetch a "real" image (you can find smregit.bmp on the home page) and show the corresponding frequency domain image. Try to explain the result! (imread, ind2gray)
4 Part 2 - Filtering 2.1. Convolution (filtering in the spatial domain) If the filtering function is known and you want to calculate a specific out-signal from the in-signal, you can use two methods. In the spatial domain convolution is used. You apply convolution to the insignal and the impulse response of the filter. Usually the in-signal is a long sequence and the impulse response is a short sequence - usually called the convolution kernel. (conv, conv2) Create an image with a 15 x 15 square in the middle, like the one used in 1.2. Create the following convolution kernel with the code below: /9 filter1=ones(3, 3)./9; Calculate the output image through convolution How is the image changed in the filtering process? Compare the spectra for the images! Change the size of the filter to size 7 x 7. This is a mean filter, remember to divide each pixel in the filter with the sum of ones in the kernel. Repeat the previous exercise and compare the result! What importance does the size of the kernel have? The following filter is a Laplace filter: filter2=[0 1 0; 1-4 1; 0 1 0]; How is the image and its spectrum affected by this kind of filter? 2.2. Filtering in the frequency domain The other method of filtering is filtering in the frequency domain. The transform of the image is multiplied with a filter that attenuates certain frequencies. The filter can either be created directly in the frequency domain or be the transform of a filter created in the spatial domain. In this case the Fourier transform of the image is multiplied with the Fourier transform of the impulse response - the transfer function. Transform both the image and the 3x3 averaging filter to the frequency domain: [M1 N1]=size(bild); [M2 N2]=size(filter); M=M1+M2-1; %The images has to be zero-padded to this size! N=N1+N2-1;
5 fftbild=fftshift(fft2(bild, M, N)); % calculate the FFT, fftfilter=fftshift(fft2(filter, M, N)); %automatic zero-pad B1=log(abs(fftbild)+1); B1=(B1-min(min(B1)))./(max(max(B1))-min(min(B1))); C1=log(abs(fftfilter)+1); C1=(C1-min(min(C1)))./(max(max(C1))-min(min(C1))); figure subplot(1,2,1),imshow(b1);title('original'); subplot(1,2,2),imshow(c1);title('filter'); A new, filtered, image (in the frequency domain) is created through multiplication: nybild=fftbild.*fftfilter; %element wise multiplication D1=log(abs(nybild)+1); D1=(D1-min(min(D1)))./(max(max(D1))-min(min(D1))); figure subplot(1,2,1),imshow(d1);title('filtered image, magnitude as a function of frequency'); The last step is to do the inverse transform: nybild=fftshift(nybild); %shift back recreate=ifft2(nybild); subplot(1,2,2),imshow(recreate);title('filtered image'); Describe the function of the filter. How is the image affected? Compare with the result in What is the difference? See also what happens with a 7 x 7 filter. Compare with 2.1.2! Change the kernel to a Laplace filter and repeat the filtering. Compare with Create a filter directly in the frequency domain. For example, you can create an ideal low pass filter: [M, N]=size(fftbild); radius=15; filter=zeros(m, N); for i=1:m for j=1:n if (((i-(m+1)/2).^2+(j-(n+1)/2).^2).^0.5)<=radius filter(i,j)=1; Remember to do the zero-padding (double the size of the image in the transformation process). After the filtering you can go back to the original size with "imcrop". Try the filter above. What happens if the radius of the filter is changed? As you see - ideal filter will produce ripple.
6 2.2.6 Change the filter to get the "opposite" effect! If you have used the filter above, change it to a highpass filter. Describe the effect of this filter Discuss advantages and disadvantages with filtering in spatial and i frequency domain. You can get some hints on page 166 in the book
7 Part 3 - Image enhancement 3.1.Sharpening of images With a high pass filter it is possible to emphasize the contour of an object. This is usually a useful first step in further automatic manipulations of the image, but the result does not look "better" for the human eye. It is more common that you want to "sharpen" the image by enhancing its contours. This can be done by combining the original image with a high pass filtered version of it. On the home page you can find the image "debbie4.bmp" (imread, ind2gray). A high pass filter can be combined with the original image through the following filter: = Convolute the image with the filter. Show the result and describe the effect Histogram equalization An ATM machine in Umeå has been busted and a lot of money has been stolen. An image of the thief was captured by the surveillance system, but unfortunately it was damaged when the film was developed. The image seems to see totally gray ("robber.bmp"). A policeman interested in image processing suggests that histogram equalization should help in identifying the thief. Help the police (histeq) Show the result - is it possible to identify the thief? Report A report should be written and it shall contain: Task description Solution Answer to all the questions and the choices you have made based on these answers Figures Code Reflections on the task
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 informationImage Enhancement in the Spatial Domain Low and High Pass Filtering
Image Enhancement in the Spatial Domain Low and High Pass Filtering Topics Low Pass Filtering Averaging Median Filter High Pass Filtering Edge Detection Line Detection Low Pass Filtering Low pass filters
More informationSharpening Spatial Filters ( high pass)
Sharpening Spatial Filters ( high pass) Previously we have looked at smoothing filters which remove fine detail Sharpening spatial filters seek to highlight fine detail Remove blurring from images Highlight
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 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 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
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 informationMATLAB 6.5 Image Processing Toolbox Tutorial
MATLAB 6.5 Image Processing Toolbox Tutorial The purpose of this tutorial is to gain familiarity with MATLAB s Image Processing Toolbox. This tutorial does not contain all of the functions available in
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 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 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 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 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 informationIt is the speed and discrete nature of the FFT that allows us to analyze a signal's spectrum with MATLAB.
MATLAB Addendum on Fourier Stuff 1. Getting to know the FFT What is the FFT? FFT = Fast Fourier Transform. The FFT is a faster version of the Discrete Fourier Transform(DFT). The FFT utilizes some clever
More informationDigital Image Processing COSC 6380/4393
Digital Image Processing COSC 638/4393 Lecture 9 Sept 26 th, 217 Pranav Mantini Slides from Dr. Shishir K Shah and Frank (Qingzhong) Liu, S. Narasimhan HISTOGRAM SHAPING We now describe methods for histogram
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 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 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 informationAchim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University
Achim J. Lilienthal Mobile Robotics and Olfaction Lab, Room T29, Mo, -2 o'clock AASS, Örebro University (please drop me an email in advance) achim.lilienthal@oru.se 4.!!!!!!!!! Pre-Class Reading!!!!!!!!!
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 informationLog Booklet for EE2 Experiments
Log Booklet for EE2 Experiments Vasil Zlatanov DFT experiment Exercise 1 Code for sinegen.m function y = sinegen(fsamp, fsig, nsamp) tsamp = 1/fsamp; t = 0 : tsamp : (nsamp-1)*tsamp; y = sin(2*pi*fsig*t);
More informationUnderstanding Digital Signal Processing
Understanding Digital Signal Processing Richard G. Lyons PRENTICE HALL PTR PRENTICE HALL Professional Technical Reference Upper Saddle River, New Jersey 07458 www.photr,com Contents Preface xi 1 DISCRETE
More informationThe Discrete Fourier Transform. Claudia Feregrino-Uribe, Alicia Morales-Reyes Original material: Dr. René Cumplido
The Discrete Fourier Transform Claudia Feregrino-Uribe, Alicia Morales-Reyes Original material: Dr. René Cumplido CCC-INAOE Autumn 2015 The Discrete Fourier Transform Fourier analysis is a family of mathematical
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 informationLABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS
LABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS INTRODUCTION The objective of this lab is to explore many issues involved in sampling and reconstructing signals, including analysis of the frequency
More information02/02/10. Image Filtering. Computer Vision CS 543 / ECE 549 University of Illinois. Derek Hoiem
2/2/ Image Filtering Computer Vision CS 543 / ECE 549 University of Illinois Derek Hoiem Questions about HW? Questions about class? Room change starting thursday: Everitt 63, same time Key ideas from last
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 informationReading Instructions Chapters for this lecture. Computer Assisted Image Analysis Lecture 2 Point Processing. Image Processing
1/34 Reading Instructions Chapters for this lecture 2/34 Computer Assisted Image Analysis Lecture 2 Point Processing Anders Brun (anders@cb.uu.se) Centre for Image Analysis Swedish University of Agricultural
More informationDepartment of Electronic Engineering NED University of Engineering & Technology. LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202)
Department of Electronic Engineering NED University of Engineering & Technology LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202) Instructor Name: Student Name: Roll Number: Semester: Batch:
More informationLaboration Exercises in Digital Signal Processing
Laboration Exercises in Digital Signal Processing Mikael Swartling Department of Electrical and Information Technology Lund Institute of Technology revision 215 Introduction Introduction The traditional
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 informationElectrical & Computer Engineering Technology
Electrical & Computer Engineering Technology EET 419C Digital Signal Processing Laboratory Experiments by Masood Ejaz Experiment # 1 Quantization of Analog Signals and Calculation of Quantized noise Objective:
More informationLecture 3 Complex Exponential Signals
Lecture 3 Complex Exponential Signals Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/1 1 Review of Complex Numbers Using Euler s famous formula for the complex exponential The
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 informationThe Scientist and Engineer's Guide to Digital Signal Processing By Steven W. Smith, Ph.D.
The Scientist and Engineer's Guide to Digital Signal Processing By Steven W. Smith, Ph.D. Home The Book by Chapters About the Book Steven W. Smith Blog Contact Book Search Download this chapter in PDF
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 15 Image Processing 14/04/15 http://www.ee.unlv.edu/~b1morris/ee482/
More informationDigital Image Processing
Digital Image Processing Part : Image Enhancement in the Spatial Domain AASS Learning Systems Lab, Dep. Teknik Room T9 (Fr, - o'clock) achim.lilienthal@oru.se Course Book Chapter 3-4- Contents. Image Enhancement
More informationAdvanced Audiovisual Processing Expected Background
Advanced Audiovisual Processing Expected Background As an advanced module, we will not cover introductory topics in lecture. You are expected to already be proficient with all of the following topics,
More informationChapter 3 Image Enhancement in the Spatial Domain. Chapter 3 Image Enhancement in the Spatial Domain
It makes all the difference whether one sees darkness through the light or brightness through the shadows. - David Lindsay 3.1 Background 76 3.2 Some Basic Gray Level Transformations 78 3.3 Histogram Processing
More informationTDI2131 Digital Image Processing
TDI2131 Digital Image Processing Image Enhancement in Spatial Domain Lecture 3 John See Faculty of Information Technology Multimedia University Some portions of content adapted from Zhu Liu, AT&T Labs.
More informationColor Space 1: RGB Color Space. Color Space 2: HSV. RGB Cube Easy for devices But not perceptual Where do the grays live? Where is hue and saturation?
Color Space : RGB Color Space Color Space 2: HSV RGB Cube Easy for devices But not perceptual Where do the grays live? Where is hue and saturation? Hue, Saturation, Value (Intensity) RBG cube on its vertex
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 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 SIGNAL PROCESSING TOOLS VERSION 4.0
(Digital Signal Processing Tools) Indian Institute of Technology Roorkee, Roorkee DIGITAL SIGNAL PROCESSING TOOLS VERSION 4.0 A Guide that will help you to perform various DSP functions, for a course in
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 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 informationLAB 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 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 informationChapter 5 Window Functions. periodic with a period of N (number of samples). This is observed in table (3.1).
Chapter 5 Window Functions 5.1 Introduction As discussed in section (3.7.5), the DTFS assumes that the input waveform is periodic with a period of N (number of samples). This is observed in table (3.1).
More informationImage Enhancement. Image Enhancement
SPATIAL FILTERING g h * h g FREQUENCY DOMAIN FILTERING G H. F F H G Copright RMR / RDL - 999. PEE53 - Processamento Digital de Imagens LOW PASS FILTERING attenuate or eliminate high-requenc components
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 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 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 informationFFT Convolution. The Overlap-Add Method
CHAPTER 18 FFT Convolution This chapter presents two important DSP techniques, the overlap-add method, and FFT convolution. The overlap-add method is used to break long signals into smaller segments for
More informationIIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters
IIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters (ii) Ability to design lowpass IIR filters according to predefined specifications based on analog
More informationPart 2: Fourier transforms. Key to understanding NMR, X-ray crystallography, and all forms of microscopy
Part 2: Fourier transforms Key to understanding NMR, X-ray crystallography, and all forms of microscopy Sine waves y(t) = A sin(wt + p) y(x) = A sin(kx + p) To completely specify a sine wave, you need
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 informationLast Lecture. Lecture 2, Point Processing GW , & , Ida-Maria Which image is wich channel?
Last Lecture Lecture 2, Point Processing GW 2.6-2.6.4, & 3.1-3.4, Ida-Maria Ida.sintorn@it.uu.se Digitization -sampling in space (x,y) -sampling in amplitude (intensity) How often should you sample in
More informationECEn 487 Digital Signal Processing Laboratory. Lab 3 FFT-based Spectrum Analyzer
ECEn 487 Digital Signal Processing Laboratory Lab 3 FFT-based Spectrum Analyzer Due Dates This is a three week lab. All TA check off must be completed by Friday, March 14, at 3 PM or the lab will be marked
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 informationImage Processing Lecture 4
Image Enhancement Image enhancement aims to process an image so that the output image is more suitable than the original. It is used to solve some computer imaging problems, or to improve image quality.
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 informationDFT: Discrete Fourier Transform & Linear Signal Processing
DFT: Discrete Fourier Transform & Linear Signal Processing 2 nd Year Electronics Lab IMPERIAL COLLEGE LONDON Table of Contents Equipment... 2 Aims... 2 Objectives... 2 Recommended Textbooks... 3 Recommended
More informationMemorial University of Newfoundland Faculty of Engineering and Applied Science. Lab Manual
Memorial University of Newfoundland Faculty of Engineering and Applied Science Engineering 6871 Communication Principles Lab Manual Fall 2014 Lab 1 AMPLITUDE MODULATION Purpose: 1. Learn how to use Matlab
More informationLab 3 FFT based Spectrum Analyzer
ECEn 487 Digital Signal Processing Laboratory Lab 3 FFT based Spectrum Analyzer Due Dates This is a three week lab. All TA check off must be completed prior to the beginning of class on the lab book submission
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 Image Processing. Lecture # 3 Image Enhancement
Digital Image Processing Lecture # 3 Image Enhancement 1 Image Enhancement Image Enhancement 3 Image Enhancement 4 Image Enhancement Process an image so that the result is more suitable than the original
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 informationComputer Vision. Intensity transformations
Computer Vision Intensity transformations Filippo Bergamasco (filippo.bergamasco@unive.it) http://www.dais.unive.it/~bergamasco DAIS, Ca Foscari University of Venice Academic year 2016/2017 Introduction
More informationImage Processing. Chapter(3) Part 2:Intensity Transformation and spatial filters. Prepared by: Hanan Hardan. Hanan Hardan 1
Image Processing Chapter(3) Part 2:Intensity Transformation and spatial filters Prepared by: Hanan Hardan Hanan Hardan 1 Image Enhancement? Enhancement تحسين الصورة : is to process an image so that the
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 informationDIGITAL SIGNAL PROCESSING WITH VHDL
DIGITAL SIGNAL PROCESSING WITH VHDL GET HANDS-ON FROM THEORY TO PRACTICE IN 6 DAYS MODEL WITH SCILAB, BUILD WITH VHDL NUMEROUS MODELLING & SIMULATIONS DIRECTLY DESIGN DSP HARDWARE Brought to you by: Copyright(c)
More information30 lesions. 30 lesions. false positive fraction
Solutions to the exercises. 1.1 In a patient study for a new test for multiple sclerosis (MS), thirty-two of the one hundred patients studied actually have MS. For the data given below, complete the two-by-two
More informationOFDM Channel Modeling for WiMAX
OFDM Channel Modeling for WiMAX April 27, 2007 David Doria Goals: To develop a simplified model of a Rayleigh fading channel Apply this model to an OFDM system Implement the above in network simulation
More informationImplementation of Image Deblurring Techniques in Java
Implementation of Image Deblurring Techniques in Java Peter Chapman Computer Systems Lab 2007-2008 Thomas Jefferson High School for Science and Technology Alexandria, Virginia January 22, 2008 Abstract
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 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 informationSante FFT Imaging Copyright 2018 Santesoft, all rights reserved
Sante FFT Imaging Copyright 2018 Santesoft, all rights reserved Table of Contents About the program... 2 System Requirements... 2 The Fourier transform... 3 The user interface... 5 Customize the toolbar...
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 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 informationDIGITAL FILTERS. !! Finite Impulse Response (FIR) !! Infinite Impulse Response (IIR) !! Background. !! Matlab functions AGC DSP AGC DSP
DIGITAL FILTERS!! Finite Impulse Response (FIR)!! Infinite Impulse Response (IIR)!! Background!! Matlab functions 1!! Only the magnitude approximation problem!! Four basic types of ideal filters with magnitude
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 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 informationECE 421 Introduction to Signal Processing
ECE 421 Introduction to Signal Processing Dror Baron Assistant Professor Dept. of Electrical and Computer Engr. North Carolina State University, NC, USA Digital Filter Design [Reading material: Chapter
More informationInterference in stimuli employed to assess masking by substitution. Bernt Christian Skottun. Ullevaalsalleen 4C Oslo. Norway
Interference in stimuli employed to assess masking by substitution Bernt Christian Skottun Ullevaalsalleen 4C 0852 Oslo Norway Short heading: Interference ABSTRACT Enns and Di Lollo (1997, Psychological
More informationSignals A Preliminary Discussion EE442 Analog & Digital Communication Systems Lecture 2
Signals A Preliminary Discussion EE442 Analog & Digital Communication Systems Lecture 2 The Fourier transform of single pulse is the sinc function. EE 442 Signal Preliminaries 1 Communication Systems and
More informationTopic 6. The Digital Fourier Transform. (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith)
Topic 6 The Digital Fourier Transform (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith) 10 20 30 40 50 60 70 80 90 100 0-1 -0.8-0.6-0.4-0.2 0 0.2 0.4
More informationACM 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 informationChapter 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 informationLec 05 - Linear Filtering & Edge Detection
ECE 484 Digital Image Processing Lec 05 - Linear Filtering & Edge Detection Zhu Li Dept of CSEE, UMKC Office: FH560E, Email: lizhu@umkc.edu, Ph: x 2346. http://l.web.umkc.edu/lizhu Z. Li, ECE 484 Digital
More informationDigital Signal Processing ETI
2012 Digital Signal Processing ETI265 2012 Introduction In the course we have 2 laboratory works for 2012. Each laboratory work is a 3 hours lesson. We will use MATLAB for illustrate some features in digital
More informationFiltering. Image Enhancement Spatial and Frequency Based
Filtering Image Enhancement Spatial and Frequency Based Brent M. Dingle, Ph.D. 2015 Game Design and Development Program Mathematics, Statistics and Computer Science University of Wisconsin - Stout Lecture
More informationDigital Signal Processing ETI
2011 Digital Signal Processing ETI265 2011 Introduction In the course we have 2 laboratory works for 2011. Each laboratory work is a 3 hours lesson. We will use MATLAB for illustrate some features in digital
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 informationCSE 564: Scientific Visualization
CSE 564: Scientific Visualization Lecture 5: Image Processing Klaus Mueller Stony Brook University Computer Science Department Klaus Mueller, Stony Brook 2003 Image Processing Definitions Purpose: - enhance
More informationNotes on Fourier transforms
Fourier Transforms 1 Notes on Fourier transforms The Fourier transform is something we all toss around like we understand it, but it is often discussed in an offhand way that leads to confusion for those
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 informationSignal Processing for Speech Applications - Part 2-1. Signal Processing For Speech Applications - Part 2
Signal Processing for Speech Applications - Part 2-1 Signal Processing For Speech Applications - Part 2 May 14, 2013 Signal Processing for Speech Applications - Part 2-2 References Huang et al., Chapter
More informationDISCRETE 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 informationGraphs of sin x and cos x
Graphs of sin x and cos x One cycle of the graph of sin x, for values of x between 0 and 60, is given below. 1 0 90 180 270 60 1 It is this same shape that one gets between 60 and below). 720 and between
More information