C2: Medical Image Processing Linwei Wang
|
|
- Bethanie Dennis
- 6 years ago
- Views:
Transcription
1 C2: Medical Image Processing Linwei Wang
2 Content Enhancement Improve visual quality of the image When the image is too dark, too light, or has low contrast Highlight certain features of the image Restoration Reduce blur or noise in a degraded image Morphological image processing
3 Image Enhancement - Content Enhancement Spatial-domain image enhancement Transformations Single-point transformations Multi-frame transformations Filtering Smoothing filters Sharpening filters Nonlinear filtering Frequency-domain image enhancement Smoothing filters Sharpening filters
4 Restoration Image Restoration - Content Noise models Noise-only spatial filters Inverse Filter Wiener Filter (minimum mean square error filtering) Constrained least square filtering Geometric mean filter Total-variation image de-noising (*)
5 Morphological Processing- Content Morphological processing Dilation & erosion Opening & closing Basic morphological algorithms Boundary extraction Region filling Thinning Thickening Skeletons Pruning
6 Medical Images As Matrices
7 Medical Image Enhancement Spatial Domain
8 Single-Point Transformation One pixel in the input à one pixel in the output Only improve visual content, does not add information content No new information is obtained; might lose information if multiple pixels produces the same output Objective information is of secondary importance; subjective information is key to present to physicians & improve their understanding
9 Single-Point Transformation Map input gray levels to output gray levels through a look-up table (LUT) Window Scaling Selecting only necessary information for analyzing Visually or algorithmically simplify the task Fast speed The transformation function can be defined graphically by the user (by looking-up tables)
10 Basic Single-Point Transformations Linear Negative Logarithmic Log Inverse-log Nega%ve Nth root Log Iden%ty Power-law N-th power Nth power Inverse- log N-th root r: input gray level s: output gray level s =!(r)
11 Image Negatives Suited for enhancing white or gray detail embedded in dark regions of an image, especially when the black areas are dominant in size s = L!1! r s = L!1! r Mammogram with lesion nega%ve DIP Fig 3.4, Courtesy of GE medical system
12 Log Transformation s = clog(1+ r), r! 0 Expand the values of dark pixels & compress the higher-level values Maps a narrow range of low gray-level values into a wider range of output level Maps a wide range of high gray-level values into a narrower range of output level Compresses the dynamic range of images with large variations in pixel values Log Inverse- log
13 Power-Law Transformation Gamma-correction s = cr!, c,! > 0 Contrast manipulation Υ <1: ~ log: Expand the dark; images are brighter; lower contrast Υ >1: ~inverse log: Compress the dark; images are darker; higher contrast
14 Example original r=0.6 r=0.4 r=0.3 Reduced contrast washed- out
15 Example original r=0.6 r=0.4
16 Piecewise-Linear Transformation Contrast stretching To increase the dynamic range of the gray levels Can be arbitrarily complex; consider more user inputs
17 Image Statistics Global image statistics is useful for designing an image-specific LUT for the entire image Local image statistics is useful for adaptive contrast enhancement, with varying LUT across image Neighborhood statistics is useful for noise processing and filtering, not contrast enhancement Global Local Neighborhood
18 Single-Point Transformation Histogram Equalization Histogram equalization A standard technique for guiding single-point transformation Histograms provide a wealth of information about the image Dark Low contrast Bright Well readable
19 Histogram of an image f of size m by n: p(f k ) = n k /(mn) Let f be a random variable in [0, L-1]: p(f k ): probability density function (pdf) k Cumulative density function (cdf): monotonic function! j=0 Goal: to achieve uniform distribution of the pdf Let g be the output image after histogram equalization is applied on f: g = T(f) p(g) = p( f ) df Single-Point Transformation: Histogram Equalization dg = p( f ) dg df n j mn g = CDF( f ) = f # p(w)dw $ dg = p( f ) $ p(g) =1!" df
20 Single-Point Transformation: Histogram Equalization Discrete cases: images pdf: cdf: p( f k ) = p( f = k) = n k MN, CDF( f k ) =! p( f i ) = Output image g: Input image with arbitrary pdf -> Output image with uniform pdf Stretching the histogram for better contrast Straightforward; invertible Not computational expensive Indiscriminate k i=0 k! i=0 0! k < L n i MN g k = CDF( f k ) = k! j=0 n j MN
21 Example Single-Point Transformation Histogram Equalization 8 by 8 small image Frequency (histogram) CDF T(i) = round( CDF(i)!CDF min mn! CDF min " (L!1))
22 Single-Point Transformation Histogram Equalization Adaptive histogram equalization Compute several histograms for local sections of the image Regional contrast is increased to its maximum. Can produce significant noise, especially in homogeneous regions Contrast limited adaptive histogram equalization (CLAHE) Limit the regional contrast to a desired level Limit the maximum number of counts for each grayscale to a clip limit ; the extra counts beyond this limit are uniformly distributed among grayscale with counts less than the clip limit
23 Single-Point Transformation Histogram Equalization Early detection of ischemia in brain injury A: Directly after an injury B: After histogram equalization: Early diagnosis of ischemia C: Confirmed by CT scan 4 days later: treatment would have been too late
24 Multi-Image Transformation Single pixels from several input images à single pixel in one output image
25 Subtraction Multi-Image Transformation: Subtraction Angiography Why subtraction in addition to contrast media? Digital subtraction angiography (DSA) Both images show exactly the same area of the body Produces image with better contrast and quality; requires much less contrast medium. Can be applied to a sequence of images M j = I j!1 + I j!2 +!I j!n f j = I j! M j Information with temporal frequencies below 1/(n t) are removed, including stationary objects (0 Hz) or slow movements caused by breathing (low frequency)
26 Image Filtering Filtering: contextual/neighborhood operation (input several pixels in a neighborhood à output a single pixel) Linear filter: Additive Uniform!( f + g) =!( f )+!(g)!(" f ) = "!( f ), "! R Constancy in relation to a shift Usually operate according to convolution rule Nonlinear filter Filtering is useful for noise reduction!( f! h ) = [!( f )]! h Note: Results might need to be scaled (look-up table) to make sure the pixel values are integers within certain level
27 Image Filtering 1D Convolution: Continuous Discrete Example: $ # ( f! g)(x) = f (y)g(x " y)dy "# # $ ( f! g)[n] = f [m]g[n " m] m="#
28 Image Filtering 2D Convolution Continuous Discrete $ # $ # f (x, y)! g(x, y) = f (!,")g(x "!, y "")d! d" "# "# # f [m, n]! g[m, n] = f (i, j)g(m " i, n " j) # $ $ j="# i="# Example:
29 Low-Pass Filters: Mean Filters Low-Pass Filters The kernel (convolution mask) is filled with only positive values of coefficients Averaging filter / Mean filter! # # # " $ & & & %! # # # " $ & & & %
30 Low-Pass Filters: Gaussian Filters Gaussian kernel (smoothing) 1D: 1 2D: G(x) = G(x, y) = 2!" e! x 2 2" 2 x 1 2 +y 2 2!" e! 2" 2 Discrete approximation Smoothing / blurring images, similar to mean filter (weighted mean) Reduce the contrast of images Eliminate some material information about the edges or the texture of selected areas or structures
31 High-Pass Filters To isolate elements responsible for rapid brightness changes, e.g., contours and edges Often directional First-order derivative based Differentiation / finite difference Image Gradient Second-order derivative based Laplacian
32 High-Pass Filters: Differentiation Recall:!f!x = lim!x"0 f (x +!x, y)# f (x, y)!x " Possible kernels : directional discretize $ " $$$!f!x % f (x i#1, y)# f (x i+1, y) 2&x = ' x ( f " $ $ $ # ! % ' ' ' & " $ $ $ # 0 0 0! % ' ' ' & " $ $ $ # !1 0 % ' ' ' &
33 High-Pass Filters: Differentiation Responds strongly to noise Image noise usually results in pixels very different from neighbors Solutions Smoothing before differential filters Do we need to apply two filters (convolutions)? No, we can apply the derivative of Gaussian filter, because convolutions are associative! x "(G " f ) = (! x "G)" f! x "G! y "G
34 High-Pass Filters: Image Gradient Recall: let f(x,y) be a function of two variables Gradient at point (x,y)!f = [ The absolute value is large at boundaries "f "x "f "y % ] #! $ f ( x ' * &! y $ f )!f = ( "f "x )2 + ( "f "y )2 # (! x $ f ) 2 + (! y $ f ) 2 The direction corresponds to the steepest ascent Normally orthogonal to image boundaries small gradient
35 High-Pass Filters: Image Gradient Sobol derivative kernels # 1 0 "1 & % (! x : % 2 0 "2 ( % 1 0 "1 ( $ ' Directional variations à Edge detection # & % (!y :% ( % "1 "2 "1 ( $ '
36 High-Pass Filters: 2-nd Derivative 2 nd derivatives of images (divergence of gradient)!f = "2 f "x 2 + "2 f "y 2 = [ " "x " "y ] [ "f "x "f "y ] f = # #f Rota%onally Invariant
37 High-Pass Filters: 2-nd Derivative Using Laplacian for image enhancement #% g(x, y) = f (x, y)! "2 f (x, y) $ &% f (x, y)+ " 2 f (x, y) If center coefficient of Laplacian mask is nega%ve If center coefficient of Laplacian mask is posi%ve 2 nd derivative vs. gradient: 2 nd derivative: Superior in enhancing fine details Produce noisier results than the gradient Gradient: Has a stronger response in areas of significant gray-level transitions The response to fine details and noises are lower
38 Laplacian of Gaussian (LoG) Images should be smoothed a bit first! "G LoG(x, y) =! 1!" (1! x2 + y 2 )e! x 2 +y 2 2" 2 4 2" 2
39 Unsharp Masking Subtract a blurred version of an image from the image itself f s (x, y) = f (x, y)! f smoothed (x, y) High-boost filtering g(x, y) = wf (x, y)! f smoothed (x, y) = (w!1) f (x, y)+ f s (x, y) w: weight Laplacian is just one specific type of standard unsharp masking with w-1 =1 With increasing value of w, sharpening process becomes less and less important
40 Unsharp Masking Parameters Mask size: determine objects to be enhanced Small mask size: enhancement of small objects (vessels), high frequency components Large mask size: enhancement of large objects (vertebra), low frequency components Weight: determine strength of enhancement Unsharp mask kernal can be seen as a difference of two Gaussians (DoG) with 1 >> 2
41 Combining Spatial Enhancement Methods Laplacian of the image to obtain fine details Rather noisy: how to remove? Use a smoothed version of the gradient of the original image a. original b. Laplacian c. Laplacian sharpened (a+b)
42 d Sobel e smoothed sobel f laplacian sharpened smoothed (c) by e Combining Spatial Enhancement Methods To obtain a smoothed gradient To multiply it by the Laplacian image (works as a mask) Preserve details while reducing noise in the relatively flat area Roughly as combining the best features of the Laplacian and the gradient
43 Combining Spatial Enhancement Methods Add back to the original image for a sharpened image Power-law transformation to increase the dynamic range of the sharpened images g sharpened h Power- law transforma%on (r=0.5)
44 Combining Spatial Enhancement Methods a b c d e f g h
45 Combining Spatial Enhancement Summary Methods Utilize Laplacian to highlight fine details Utilize the gradient to enhance prominent features Smoothed version of the gradient is used to mask the Laplacian image to reduce noises Single-point transformation is used to increase the dynamic range of the gray levels
46 Filtering in MatLab Most of the filters / kernels mentioned earlier can be constructed in matlab by fspecial Average Gaussian Laplacian LoG Sobel Unsharp
47 Medical Image Enhancement Frequency Domain
48 Filtering in Frequency Domain Power spectrum of an image
49 Filtering in Frequency Domain Review: Discrete Fourier transform (DFT): F(u,v) F(u, v) = 1 MN N!1 M!1 "" y=0 x=0 xu! j2! ( M f (x, y)e + yv N ) Power spectrum: F(u,v) 2 High frequency Low frequency u v
50 Filtering in Frequency Domain Review 2: If F(u,v) = F [f(x,y)], H(u, v) = F [h(x,y)] Then F [f(x,y)*h(x,y)] = F(u,v) H(u,v) (element-by-element) FT of convolution of the image f(x,y) and the kernel h(x,y) is equal to the multiplication of each FT
51 Filtering in Frequency Domain Computational complexity: DFT: N-point data N 2 ultiplication, N(N-1) additionà O(N 2 ) F(u) = 1 N!1 xu! j2! ( N " f (x)e ) u = 0,1,..., N!1 N FFT: dividing into odd-number and even-number x=0 F(u) = 1 N = 1 N = 1 N n=0 N!1 " x=0 xu! j2! ( N f (x)e ) N /2!1 " f (2n) e! j2! ( 2nu N ) + 1 N N /2!1 " f (2n) e! j2! ( n=0 nu N /2 ) +e N /2!1 " f (2n +1) e! j2! ( n=0! j2! ( u N ) 1 N N /2!1 (2n+1)u ) N " f (2n +1) e! j2! ( n=0 nu N /2 ) = DFT N /2 ( f (0), f (2),... f (N! 2))+W N n DFT N /2 ( f (1), f (3),... f (N!1))
52 Filtering in Frequency Domain Computational complexity DFT: O(N 2 ) FFT: O(NlogN)
53 Filtering in Frequency Domain Flowchart of filtering in frequency domain Input image Real domain FFT Filtering (mul%plica%on) frequency domain Inverse FT Output image Real domain
54 Filtering in Frequency Domain 2D Filters in the frequency domain LPF HPF BPF f 0, f 1, f 2 are cut- off frequencies
55 Filtering in Frequency Domain X Original image A`er FFT
56 Smoothing Frequency-Domain Filters Attenuating a specified range of high-frequency components Ideal LPF! H(u, v) = " 1 # 0 If D(u,v) <= D 0 Otherwise Blurring: center component Ringing: concentric components
57 Smoothing Frequency-Domain Filters Butterworth LPF 1 H(u, v) = 1+ (D(u, v) / D 0 ) 2n n = 1, 2, 5, 20 n = 1: no ringing, LP effect is low n =2: imperceptible ringing but much less pronounced than ILPF n >2: more and more ringing effect BLPF of order 2 is a good compromise between effective lowpass filtering and acceptable ringing characteristics
58 Smoothing Frequency-Domain Filters Gassian lowpass filter H(u, v) = e!d2 (u,v)/2! 2 σ = D 0 is the cutoff frequency No ringing, suitable for applications where artifact is not acceptable (such as medical images!) Less smoothing compared to 2-order BLPF at the same cut-off frequency: BLPF is more suitable in applications where tight control of the transition around the cut-off frequency is needed
59 Sharpening Frequency-Domain Filters H hp (u, v) =1! H lp (u, v) Ideal HPF! H(u, v) = " 0 # 1 If D(u,v) <= D 0 Otherwise Ringing: concentric components Butterworth HPF Smoother than IHPF Gaussian HPF H(u, v) =1! e!d2 (u,v)/2! 2 1 H(u, v) = 1+ (D 0 / D(u, v)) 2n Smoother than IHPF & BHPF Can be constructed as DoG: allow more control over filter shape
60 Sharpening Frequency-Domain Filters Unsharp Masking & high-boost filtering f hb (x, y) = wf (x, y)! f lp (x, y) = (w!1) f (x, y)+ f hp (x, y) H hp (u, v) =1! H lp (u, v) H hb (u, v) = (w!1)+ H hp (u, v) w "1 high-frequency emphasis filtering H hfe (u, v) = a + bh hp (u, v), a! 0, b > a a: offset b: emphasizing high-frequency components
61 Nonlinear Filters Order-statistics filters Based on the ordering (ranking) the pixels in the neighborhood & replacing the value of the center pixel with the ranking results Median filter: for the processed point of the output image, the median filter selects a middle (median) value of the sorted values of the point and its neighbors in the input image Considered very radical Usually does not reduce the sharpness of object edges Does not introduce new values into the image, so no scaling is needed Unavoidable erosion of particularly small fragments Time-consuming Most frequently used to eliminate salt-and-pepper noise (impulse noise)
62 Geometric Mean Filters ˆf (x, y) = [ " (s,t)!s x,y g(s,t)] 1 MN " Achieve smoothing comparable to the arithmetic mean filter, but tends to lose less image details
Vivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT.
Vivekananda Collegee of Engineering & Technology Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT Dept. Prepared by Harivinod N Assistant Professor, of Computer Science and Engineering,
More informationCHAPTER 3 IMAGE ENHANCEMENT IN THE SPATIAL DOMAIN
CHAPTER 3 IMAGE ENHANCEMENT IN THE SPATIAL DOMAIN CHAPTER 3: IMAGE ENHANCEMENT IN THE SPATIAL DOMAIN Principal objective: to process an image so that the result is more suitable than the original image
More informationChapter 3: Intensity Transformations and Spatial Filtering
Chapter 3: Intensity Transformations and Spatial Filtering 3.1 Background 3.2 Some basic intensity transformation functions 3.3 Histogram processing 3.4 Fundamentals of spatial filtering 3.5 Smoothing
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More informationPoint and Spatial Processing
Filtering 1 Point and Spatial Processing Spatial Domain g(x,y) = T[ f(x,y) ] f(x,y) input image g(x,y) output image T is an operator on f Defined over some neighborhood of (x,y) can operate on a set of
More informationEEM 463 Introduction to Image Processing. Week 3: Intensity Transformations
EEM 463 Introduction to Image Processing Week 3: Intensity Transformations Fall 2013 Instructor: Hatice Çınar Akakın, Ph.D. haticecinarakakin@anadolu.edu.tr Anadolu University Enhancement Domains Spatial
More informationDigital Image Processing, 2nd ed. Digital Image Processing, 2nd ed. The principal objective of enhancement
Chapter 3 Image Enhancement in the Spatial Domain The principal objective of enhancement to process an image so that the result is more suitable than the original image for a specific application. Enhancement
More informationDigital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering
Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing
More informationImage Enhancement in Spatial Domain (Chapter 3)
Image Enhancement in Spatial Domain (Chapter 3) Yun Q. Shi shi@njit.edu Fall 11 Mask/Neighborhood Processing ECE643 2 1 Point Processing ECE643 3 Image Negatives S = (L 1) - r (3.2-1) Point processing
More informationImage restoration. Lecture 14. Milan Gavrilovic Centre for Image Analysis Uppsala University
Image restoration Lecture 14 Milan Gavrilovic milan@cb.uu.se Centre for Image Analysis Uppsala University Computer Assisted Image Analysis 2009-05-08 M. Gavrilovic (Uppsala University) L14 Image restoration
More informationEELE 5310: Digital Image Processing. Lecture 2 Ch. 3. Eng. Ruba A. Salamah. iugaza.edu
EELE 5310: Digital Image Processing Lecture 2 Ch. 3 Eng. Ruba A. Salamah Rsalamah @ iugaza.edu 1 Image Enhancement in the Spatial Domain 2 Lecture Reading 3.1 Background 3.2 Some Basic Gray Level Transformations
More informationEELE 5310: Digital Image Processing. Ch. 3. Eng. Ruba A. Salamah. iugaza.edu
EELE 531: Digital Image Processing Ch. 3 Eng. Ruba A. Salamah Rsalamah @ iugaza.edu 1 Image Enhancement in the Spatial Domain 2 Lecture Reading 3.1 Background 3.2 Some Basic Gray Level Transformations
More informationDigital Image Processing. Lecture 6
Digital Image Processing Lecture 6 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Fall 2016 Image Enhancement In The Frequency Domain Outline Jean Baptiste Joseph
More informationImage Enhancement: To improve the quality of images
Image Enhancement: To improve the quality of images Examples: Noise reduction (to improve SNR or subjective quality) Change contrast, brightness, color etc. Image smoothing Image sharpening Modify image
More informationOutlines. Medical Image Processing Using Transforms. 4. Transform in image space
Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca Outlines Image Quality Gray value transforms Histogram processing Transforms
More informationIntroduction to Digital Image Processing
Fall 2005 Image Enhancement in the Spatial Domain: Histograms, Arithmetic/Logic Operators, Basics of Spatial Filtering, Smoothing Spatial Filters Tuesday, February 7 2006, Overview (1): Before We Begin
More informationWhat will we learn? Neighborhood processing. Convolution and correlation. Neighborhood processing. Chapter 10 Neighborhood Processing
What will we learn? Lecture Slides ME 4060 Machine Vision and Vision-based Control Chapter 10 Neighborhood Processing By Dr. Debao Zhou 1 What is neighborhood processing and how does it differ from point
More informationFiltering and Enhancing Images
KECE471 Computer Vision Filtering and Enhancing Images Chang-Su Kim Chapter 5, Computer Vision by Shapiro and Stockman Note: Some figures and contents in the lecture notes of Dr. Stockman are used partly.
More informationClassification of image operations. Image enhancement (GW-Ch. 3) Point operations. Neighbourhood operation
Image enhancement (GW-Ch. 3) Classification of image operations Process of improving image quality so that the result is more suitable for a specific application. contrast stretching histogram processing
More informationLecture 4: Spatial Domain Transformations
# Lecture 4: Spatial Domain Transformations Saad J Bedros sbedros@umn.edu Reminder 2 nd Quiz on the manipulator Part is this Fri, April 7 205, :5 AM to :0 PM Open Book, Open Notes, Focus on the material
More informationSYDE 575: Introduction to Image Processing
SYDE 575: Introduction to Image Processing Image Enhancement and Restoration in Spatial Domain Chapter 3 Spatial Filtering Recall 2D discrete convolution g[m, n] = f [ m, n] h[ m, n] = f [i, j ] h[ m i,
More informationImage Restoration and Reconstruction
Image Restoration and Reconstruction Image restoration Objective process to improve an image, as opposed to the subjective process of image enhancement Enhancement uses heuristics to improve the image
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spatial Domain Filtering http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background Intensity
More informationEECS 556 Image Processing W 09. Image enhancement. Smoothing and noise removal Sharpening filters
EECS 556 Image Processing W 09 Image enhancement Smoothing and noise removal Sharpening filters What is image processing? Image processing is the application of 2D signal processing methods to images Image
More informationDigital Image Processing
Digital Image Processing Jen-Hui Chuang Department of Computer Science National Chiao Tung University 2 3 Image Enhancement in the Spatial Domain 3.1 Background 3.4 Enhancement Using Arithmetic/Logic Operations
More informationImage Restoration and Reconstruction
Image Restoration and Reconstruction Image restoration Objective process to improve an image Recover an image by using a priori knowledge of degradation phenomenon Exemplified by removal of blur by deblurring
More informationLinear Operations Using Masks
Linear Operations Using Masks Masks are patterns used to define the weights used in averaging the neighbors of a pixel to compute some result at that pixel Expressing linear operations on neighborhoods
More informationImage Processing. Traitement d images. Yuliya Tarabalka Tel.
Traitement d images Yuliya Tarabalka yuliya.tarabalka@hyperinet.eu yuliya.tarabalka@gipsa-lab.grenoble-inp.fr Tel. 04 76 82 62 68 Noise reduction Image restoration Restoration attempts to reconstruct an
More informationSampling and Reconstruction
Sampling and Reconstruction Sampling and Reconstruction Sampling and Spatial Resolution Spatial Aliasing Problem: Spatial aliasing is insufficient sampling of data along the space axis, which occurs because
More informationBiomedical Image Analysis. Spatial Filtering
Biomedical Image Analysis Contents: Spatial Filtering The mechanics of Spatial Filtering Smoothing and sharpening filters BMIA 15 V. Roth & P. Cattin 1 The Mechanics of Spatial Filtering Spatial filter:
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 03 Image Processing Basics 13/01/28 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationLecture 4 Image Enhancement in Spatial Domain
Digital Image Processing Lecture 4 Image Enhancement in Spatial Domain Fall 2010 2 domains Spatial Domain : (image plane) Techniques are based on direct manipulation of pixels in an image Frequency Domain
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
More informationImage Processing Lecture 10
Image Restoration Image restoration attempts to reconstruct or recover an image that has been degraded by a degradation phenomenon. Thus, restoration techniques are oriented toward modeling the degradation
More informationLecture 6: Edge Detection
#1 Lecture 6: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Options for Image Representation Introduced the concept of different representation or transformation Fourier Transform
More informationBabu Madhav Institute of Information Technology Years Integrated M.Sc.(IT)(Semester - 7)
5 Years Integrated M.Sc.(IT)(Semester - 7) 060010707 Digital Image Processing UNIT 1 Introduction to Image Processing Q: 1 Answer in short. 1. What is digital image? 1. Define pixel or picture element?
More informationLecture 4. Digital Image Enhancement. 1. Principle of image enhancement 2. Spatial domain transformation. Histogram processing
Lecture 4 Digital Image Enhancement 1. Principle of image enhancement 2. Spatial domain transformation Basic intensity it tranfomation ti Histogram processing Principle Objective of Enhancement Image enhancement
More informationBiometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong)
Biometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong) References: [1] http://homepages.inf.ed.ac.uk/rbf/hipr2/index.htm [2] http://www.cs.wisc.edu/~dyer/cs540/notes/vision.html
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationDigital Image Processing. Image Enhancement in the Spatial Domain (Chapter 4)
Digital Image Processing Image Enhancement in the Spatial Domain (Chapter 4) Objective The principal objective o enhancement is to process an images so that the result is more suitable than the original
More informationSECTION 5 IMAGE PROCESSING 2
SECTION 5 IMAGE PROCESSING 2 5.1 Resampling 3 5.1.1 Image Interpolation Comparison 3 5.2 Convolution 3 5.3 Smoothing Filters 3 5.3.1 Mean Filter 3 5.3.2 Median Filter 4 5.3.3 Pseudomedian Filter 6 5.3.4
More informationChapter4 Image Enhancement
Chapter4 Image Enhancement Preview 4.1 General introduction and Classification 4.2 Enhancement by Spatial Transforming(contrast enhancement) 4.3 Enhancement by Spatial Filtering (image smoothing) 4.4 Enhancement
More information1.Some Basic Gray Level Transformations
1.Some Basic Gray Level Transformations We begin the study of image enhancement techniques by discussing gray-level transformation functions.these are among the simplest of all image enhancement techniques.the
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
More informationChapter - 2 : IMAGE ENHANCEMENT
Chapter - : IMAGE ENHANCEMENT The principal objective of enhancement technique is to process a given image so that the result is more suitable than the original image for a specific application Image Enhancement
More informationDigital Image Processing. Image Enhancement - Filtering
Digital Image Processing Image Enhancement - Filtering Derivative Derivative is defined as a rate of change. Discrete Derivative Finite Distance Example Derivatives in 2-dimension Derivatives of Images
More informationFiltering Images. Contents
Image Processing and Data Visualization with MATLAB Filtering Images Hansrudi Noser June 8-9, 010 UZH, Multimedia and Robotics Summer School Noise Smoothing Filters Sigmoid Filters Gradient Filters Contents
More informationImage Enhancement in Spatial Domain. By Dr. Rajeev Srivastava
Image Enhancement in Spatial Domain By Dr. Rajeev Srivastava CONTENTS Image Enhancement in Spatial Domain Spatial Domain Methods 1. Point Processing Functions A. Gray Level Transformation functions for
More informationCS4442/9542b Artificial Intelligence II prof. Olga Veksler
CS4442/9542b Artificial Intelligence II prof. Olga Veksler Lecture 2 Computer Vision Introduction, Filtering Some slides from: D. Jacobs, D. Lowe, S. Seitz, A.Efros, X. Li, R. Fergus, J. Hayes, S. Lazebnik,
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationCS334: Digital Imaging and Multimedia Edges and Contours. Ahmed Elgammal Dept. of Computer Science Rutgers University
CS334: Digital Imaging and Multimedia Edges and Contours Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What makes an edge? Gradient-based edge detection Edge Operators From Edges
More informationPSD2B Digital Image Processing. Unit I -V
PSD2B Digital Image Processing Unit I -V Syllabus- Unit 1 Introduction Steps in Image Processing Image Acquisition Representation Sampling & Quantization Relationship between pixels Color Models Basics
More informationDigital Image Analysis and Processing
Digital Image Analysis and Processing CPE 0907544 Image Enhancement Part I Intensity Transformation Chapter 3 Sections: 3.1 3.3 Dr. Iyad Jafar Outline What is Image Enhancement? Background Intensity Transformation
More informationIMAGE ENHANCEMENT IN THE SPATIAL DOMAIN
1 Image Enhancement in the Spatial Domain 3 IMAGE ENHANCEMENT IN THE SPATIAL DOMAIN Unit structure : 3.0 Objectives 3.1 Introduction 3.2 Basic Grey Level Transform 3.3 Identity Transform Function 3.4 Image
More informationBasic relations between pixels (Chapter 2)
Basic relations between pixels (Chapter 2) Lecture 3 Basic Relationships Between Pixels Definitions: f(x,y): digital image Pixels: q, p (p,q f) A subset of pixels of f(x,y): S A typology of relations:
More informationChapter 3 Image Enhancement in the Spatial Domain
Chapter 3 Image Enhancement in the Spatial Domain Yinghua He School o Computer Science and Technology Tianjin University Image enhancement approaches Spatial domain image plane itsel Spatial domain methods
More informationPoint operation Spatial operation Transform operation Pseudocoloring
Image Enhancement Introduction Enhancement by point processing Simple intensity transformation Histogram processing Spatial filtering Smoothing filters Sharpening filters Enhancement in the frequency domain
More informationIntensity Transformations and Spatial Filtering
77 Chapter 3 Intensity Transformations and Spatial Filtering Spatial domain refers to the image plane itself, and image processing methods in this category are based on direct manipulation of pixels in
More informationImage processing. Reading. What is an image? Brian Curless CSE 457 Spring 2017
Reading Jain, Kasturi, Schunck, Machine Vision. McGraw-Hill, 1995. Sections 4.2-4.4, 4.5(intro), 4.5.5, 4.5.6, 5.1-5.4. [online handout] Image processing Brian Curless CSE 457 Spring 2017 1 2 What is an
More informationIntensity Transformation and Spatial Filtering
Intensity Transformation and Spatial Filtering Outline of the Lecture Introduction. Intensity Transformation Functions. Piecewise-Linear Transformation Functions. Introduction Definition: Image enhancement
More informationPROCESS > SPATIAL FILTERS
83 Spatial Filters There are 19 different spatial filters that can be applied to a data set. These are described in the table below. A filter can be applied to the entire volume or to selected objects
More informationComputer Vision and Graphics (ee2031) Digital Image Processing I
Computer Vision and Graphics (ee203) Digital Image Processing I Dr John Collomosse J.Collomosse@surrey.ac.uk Centre for Vision, Speech and Signal Processing University of Surrey Learning Outcomes After
More informationFundamentals of Digital Image Processing
\L\.6 Gw.i Fundamentals of Digital Image Processing A Practical Approach with Examples in Matlab Chris Solomon School of Physical Sciences, University of Kent, Canterbury, UK Toby Breckon School of Engineering,
More informationUnit - I Computer vision Fundamentals
Unit - I Computer vision Fundamentals It is an area which concentrates on mimicking human vision systems. As a scientific discipline, computer vision is concerned with the theory behind artificial systems
More informationImage restoration. Restoration: Enhancement:
Image restoration Most images obtained by optical, electronic, or electro-optic means is likely to be degraded. The degradation can be due to camera misfocus, relative motion between camera and object,
More informationCS4442/9542b Artificial Intelligence II prof. Olga Veksler
CS4442/9542b Artificial Intelligence II prof. Olga Veksler Lecture 8 Computer Vision Introduction, Filtering Some slides from: D. Jacobs, D. Lowe, S. Seitz, A.Efros, X. Li, R. Fergus, J. Hayes, S. Lazebnik,
More information2D Image Processing INFORMATIK. Kaiserlautern University. DFKI Deutsches Forschungszentrum für Künstliche Intelligenz
2D Image Processing - Filtering Prof. Didier Stricker Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de 1 What is image filtering?
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 informationImage gradients and edges April 11 th, 2017
4//27 Image gradients and edges April th, 27 Yong Jae Lee UC Davis PS due this Friday Announcements Questions? 2 Last time Image formation Linear filters and convolution useful for Image smoothing, removing
More informationDigital Image Fundamentals
Digital Image Fundamentals Image Quality Objective/ subjective Machine/human beings Mathematical and Probabilistic/ human intuition and perception 6 Structure of the Human Eye photoreceptor cells 75~50
More informationUNIT - 5 IMAGE ENHANCEMENT IN SPATIAL DOMAIN
UNIT - 5 IMAGE ENHANCEMENT IN SPATIAL DOMAIN Spatial domain methods Spatial domain refers to the image plane itself, and approaches in this category are based on direct manipulation of pixels in an image.
More informationImage gradients and edges April 10 th, 2018
Image gradients and edges April th, 28 Yong Jae Lee UC Davis PS due this Friday Announcements Questions? 2 Last time Image formation Linear filters and convolution useful for Image smoothing, removing
More informationDigital Image Processing
Digital Image Processing Image Restoration and Reconstruction (Noise Removal) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering 2 Image Restoration
More informationLecture 5: Frequency Domain Transformations
#1 Lecture 5: Frequency Domain Transformations Saad J Bedros sbedros@umn.edu From Last Lecture Spatial Domain Transformation Point Processing for Enhancement Area/Mask Processing Transformations Image
More informationCS534: Introduction to Computer Vision Edges and Contours. Ahmed Elgammal Dept. of Computer Science Rutgers University
CS534: Introduction to Computer Vision Edges and Contours Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What makes an edge? Gradient-based edge detection Edge Operators Laplacian
More informationDigital Image Processing
Digital Image Processing Part 2: Image Enhancement in the Spatial Domain AASS Learning Systems Lab, Dep. Teknik Room T1209 (Fr, 11-12 o'clock) achim.lilienthal@oru.se Course Book Chapter 3 2011-04-06 Contents
More informationLecture 7: Most Common Edge Detectors
#1 Lecture 7: Most Common Edge Detectors Saad Bedros sbedros@umn.edu Edge Detection Goal: Identify sudden changes (discontinuities) in an image Intuitively, most semantic and shape information from the
More informationNoise Model. Important Noise Probability Density Functions (Cont.) Important Noise Probability Density Functions
Others -- Noise Removal Techniques -- Edge Detection Techniques -- Geometric Operations -- Color Image Processing -- Color Spaces Xiaojun Qi Noise Model The principal sources of noise in digital images
More informationSegmentation and Grouping
Segmentation and Grouping How and what do we see? Fundamental Problems ' Focus of attention, or grouping ' What subsets of pixels do we consider as possible objects? ' All connected subsets? ' Representation
More informationIMAGE ENHANCEMENT in SPATIAL DOMAIN by Intensity Transformations
It makes all the difference whether one sees darkness through the light or brightness through the shadows David Lindsay IMAGE ENHANCEMENT in SPATIAL DOMAIN by Intensity Transformations Kalyan Kumar Barik
More informationIT Digital Image ProcessingVII Semester - Question Bank
UNIT I DIGITAL IMAGE FUNDAMENTALS PART A Elements of Digital Image processing (DIP) systems 1. What is a pixel? 2. Define Digital Image 3. What are the steps involved in DIP? 4. List the categories of
More informationVU Signal and Image Processing. Image Restoration. Torsten Möller + Hrvoje Bogunović + Raphael Sahann
052600 VU Signal and Image Processing Image Restoration Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at raphael.sahann@univie.ac.at vda.cs.univie.ac.at/teaching/sip/17s/
More informationLecture 4: Image Processing
Lecture 4: Image Processing Definitions Many graphics techniques that operate only on images Image processing: operations that take images as input, produce images as output In its most general form, an
More informationIslamic University of Gaza Faculty of Engineering Computer Engineering Department
Islamic University of Gaza Faculty of Engineering Computer Engineering Department EELE 5310: Digital Image Processing Spring 2011 Date: May 29, 2011 Time : 120 minutes Final Exam Student Name: Student
More informationSharpening through spatial filtering
Sharpening through spatial filtering Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 2017 2018 Sharpening The term sharpening is referred
More informationAn educational software package for digital image processing
An educational software package for digital image processing Juan Manuel ~amirez', Pilar Gomez-~il~, Luis G. Guerrero', David ~aez' & Ernesto ~orantes' '~e~artamento de Ingenieria Electrdnica, Universidad
More informationImage Processing Fundamentals. Nicolas Vazquez Principal Software Engineer National Instruments
Image Processing Fundamentals Nicolas Vazquez Principal Software Engineer National Instruments Agenda Objectives and Motivations Enhancing Images Checking for Presence Locating Parts Measuring Features
More informationBiomedical Image Analysis. Point, Edge and Line Detection
Biomedical Image Analysis Point, Edge and Line Detection Contents: Point and line detection Advanced edge detection: Canny Local/regional edge processing Global processing: Hough transform BMIA 15 V. Roth
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationComputer Vision I - Basics of Image Processing Part 1
Computer Vision I - Basics of Image Processing Part 1 Carsten Rother 28/10/2014 Computer Vision I: Basics of Image Processing Link to lectures Computer Vision I: Basics of Image Processing 28/10/2014 2
More informationImage gradients and edges
Image gradients and edges April 7 th, 2015 Yong Jae Lee UC Davis Announcements PS0 due this Friday Questions? 2 Last time Image formation Linear filters and convolution useful for Image smoothing, removing
More informationINTENSITY TRANSFORMATION AND SPATIAL FILTERING
1 INTENSITY TRANSFORMATION AND SPATIAL FILTERING Lecture 3 Image Domains 2 Spatial domain Refers to the image plane itself Image processing methods are based and directly applied to image pixels Transform
More informationImage Processing. Daniel Danilov July 13, 2015
Image Processing Daniel Danilov July 13, 2015 Overview 1. Principle of digital images and filters 2. Basic examples of filters 3. Edge detection and segmentation 1 / 25 Motivation For what image processing
More informationFourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)
Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) 1D: Common Transform Pairs Summary source FT Properties: Convolution See book DIP 4.2.5:
More informationReview for Exam I, EE552 2/2009
Gonale & Woods Review or Eam I, EE55 /009 Elements o Visual Perception Image Formation in the Ee and relation to a photographic camera). Brightness Adaption and Discrimination. Light and the Electromagnetic
More informationLecture Image Enhancement and Spatial Filtering
Lecture Image Enhancement and Spatial Filtering Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu September 29, 2005 Abstract Applications of
More informationIntroduction to Digital Image Processing
Introduction to Digital Image Processing Ranga Rodrigo June 9, 29 Outline Contents Introduction 2 Point Operations 2 Histogram Processing 5 Introduction We can process images either in spatial domain or
More informationCoE4TN3 Medical Image Processing
CoE4TN3 Medical Image Processing Image Restoration Noise Image sensor might produce noise because of environmental conditions or quality of sensing elements. Interference in the image transmission channel.
More informationJNTUWORLD. 4. Prove that the average value of laplacian of the equation 2 h = ((r2 σ 2 )/σ 4 ))exp( r 2 /2σ 2 ) is zero. [16]
Code No: 07A70401 R07 Set No. 2 1. (a) What are the basic properties of frequency domain with respect to the image processing. (b) Define the terms: i. Impulse function of strength a ii. Impulse function
More informationIn this lecture. Background. Background. Background. PAM3012 Digital Image Processing for Radiographers
PAM3012 Digital Image Processing for Radiographers Image Enhancement in the Spatial Domain (Part I) In this lecture Image Enhancement Introduction to spatial domain Information Greyscale transformations
More information