Medical Image Analysis
|
|
- Victor Joseph
- 5 years ago
- Views:
Transcription
1 Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 23rd January 2006
2 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of the paper... 5 Anisotropic Diffusion
3 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of the paper... 5 Anisotropic Diffusion
4 Edge Enhancement Inhomogeneous diffusion may actually enhance edges, for a certain choice of g(). 1D example: Let s(x) = I x, and φ(s) = g(s)s = g(i x)i x. The 1D heat equation becomes I t = x (g(i x)i x ) = x φ(s) = φ (s)i xx t (I x) is the rate of change of edge slope with respect to time. t (I x) = x (I t) = φ (s)i 2 xx + φ (s)i xxx
5 Edge Enhancement t (I x) = φ (s)i 2 xx + φ (s)i xxx For a step edge with I x > 0 look at the inflection point, p. Observe that I xx (p) = 0, and I xxx (p) < 0. I, I x, I xx, I xxx t (I x)(p) = φ (s)i xxx (p) The sign of this quantity depends only on φ (s).
6 Edge Enhancement t (I x) = φ (s)i 2 xx + φ (s)i xxx For a step edge with I x > 0 look at the inflection point, p. Observe that I xx (p) = 0, and I xxx (p) < 0. I, I x, I xx, I xxx t (I x)(p) = φ (s)i xxx (p) The sign of this quantity depends only on φ (s).
7 Edge Enhancement At the inflection point: t (I x)(p) = φ (s)i xxx (p) If φ (s) > 0, then t (I x)(p) < 0 (slope is decreasing). If φ (s) < 0, then t (I x)(p) > 0 (slope is increasing). Since φ(s) = g(s)s, selecting the function g(s) determines which edges are smoothed and which are sharpened.
8 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of the paper... 5 Anisotropic Diffusion
9 The function φ(s) = g(s)s φ(0) = 0 φ (s) > 0 for s < K φ (s) < 0 for s > K lim s φ(s) 0
10 The function g(s) Perona and Malik suggest two possible functions (s = I ): I ( ) g( I ) = e K 2 g( I ) = ( I (α > 0) K )1+α
11 Effect of varying K on g() 1 g( I ) = 1 + ( I (α > 0) K )1+α Figure: K = 2, 4, 6
12 Effect of varying α on g() 1 g( I ) = 1 + ( I (α > 0) K )1+α Figure: α = 1, 3, 5, 7, 9
13 Effect of varying K and α on c() Figure: I and I.
14 Effect of varying K on c() Figure: K = 3, 5, 10, 100.
15 Effect of varying α on c() Figure: α = 1, 2, 3, 5.
16 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results Discretized Inhomogeneous Heat Equation Perona-Malik Implementation Discrete Maximum Principle Adaptive Parameter Setting 4 The rest of the paper... 5 Anisotropic Diffusion
17 Discretized Inhomogeneous Heat Equation Explicit Formulation I t = c(x, y, t) 2 I + c I Using centered differences for the Laplacian and gradients: I t+1 x,y I t x,y λ = c x,y (I t x 1,y + It x+1,y + It x,y 1 + It x,y+1 4It x,y) + ( c x+1,y c x 1,y 2 + ( c x,y+1 c x,y 1 2 )( I x+1,y I x 1,y ) 2 )( I x,y+1 I x,y 1 ) 2
18 Discretized Inhomogeneous Heat Equation Implicit Formulation Same diagonal structure as homogeneous heat equation? Yes. Symmetric? No. Diagonal dominance? Data dependent.
19 Discretized Inhomogeneous Heat Equation Implicit Formulation Same diagonal structure as homogeneous heat equation? Yes. Symmetric? No. Diagonal dominance? Data dependent.
20 Discretized Inhomogeneous Heat Equation Implicit Formulation Same diagonal structure as homogeneous heat equation? Yes. Symmetric? No. Diagonal dominance? Data dependent.
21 Discretized Inhomogeneous Heat Equation Implicit Formulation Same diagonal structure as homogeneous heat equation? Yes. Symmetric? No. Diagonal dominance? Data dependent.
22 Perona-Malik Implementation Explicit Formulation I t = c(x, y, t) 2 I + c I By splitting the Laplacian and averaging the forward and backward differences in the gradient: I t+1 x,y I t x,y λ = c x,y [(I t x 1,y It x,y) + (I t x+1,y It x,y) + (I t x,y 1 It x,y) + (I t x,y+1 It x,y)] + c x (I x+1,y I x,y 2 + c y (I x,y+1 I x,y 2 + I x,y I x 1,y ) 2 + I x,y I x,y 1 ) 2
23 Perona-Malik Implementation Explicit Formulation I t = c(x, y, t) 2 I + c I I t+1 x,y I t x,y λ = (c x,y 1 c 2 x )(It x 1,y It x,y) + (c x,y + 1 c 2 x )(It x+1,y It x,y) + (c x,y 1 c 2 x )(It x,y 1 It x,y) + (c x,y + 1 c 2 x )(It x,y+1 It x,y)
24 Perona-Malik Implementation Explicit Formulation c x,y + 1 c 2 x c x,y 1 c 2 x c x+ 1 2,y c x 1 2,y c x+ 1 c x 1 Where s x,y = I(x, y). 2,y g( s x,y + s x+1,y ) 2 2,y g( s x,y + s x 1,y ) 2
25 Perona-Malik Implementation Explicit Formulation I t+1 x,y I t x,y λ = g( s x,y + s x 1,y )(Ix 1,y t 2 It x,y) + g( s x,y + s x+1,y )(Ix+1,y t 2 It x,y) + g( s x,y + s x,y 1 )(Ix,y 1 t 2 It x,y) + g( s x,y + s x,y+1 )(Ix,y+1 t 2 It x,y)
26 Perona-Malik Implementation Anisotropic Implementation Compute g() using the projection of the gradient along one direction. For example, in g( sx,y +s x+1,y 2 ), let s x,y = I (x, y) x s x+1,y = I (x + 1, y) x Computing s x,y using forward differences, and s x+1,y using backward differences s x,y = I x+1,y I x,y s x+1,y = I x+1,y I x,y, so g( sx,y +s x+1,y 2 ) = g( I(x + 1, y) I(x, y) ).
27 Perona-Malik Implementation Explicit Formulation I t+1 x,y I t x,y λ = g( I x 1,y I x,y )(I t x 1,y It x,y) + g( I x+1,y I x,y )(I t x+1,y It x,y) + g( I x,y 1 I x,y )(I t x,y 1 It x,y) + g( I x,y+1 I x,y )(I t x,y+1 It x,y) Notation: The authors use to denote finite differences. This is not the gradient operator ( ).
28 Perona-Malik Implementation N I i,j I i 1,j I i,j S I i,j I i+1,j I i,j E I i,j I i,j+1 I i,j W I i,j I i,j 1 I i,j Image neighborhood system c Ni,j = g( N I i,j ) c Si,j = g( S I i,j ) c Ei,j = g( E I i,j ) c Wi,j = g( W I i,j )
29 Perona-Malik Implementation Explicit Formulation The previous explicit formulation I t+1 x,y I t x,y λ = g( I x 1,y I x,y )(I t x 1,y It x,y) + g( I x+1,y I x,y )(I t x+1,y It x,y) + g( I x,y 1 I x,y )(I t x,y 1 It x,y) + g( I x,y+1 I x,y )(I t x,y+1 It x,y) can be rewritten as I t+1 x,y = I t x,y + λ(c Ni,j N I i,j + c Si,j S I i,j + c Ei,j E I i,j + c Wi,j W I i,j ) t
30 Discrete Maximum Principle Let I Mij = max(i, I N, I S, I E, I W ) Ix,y t+1 = Ix,y t + λ(c Ni,j N I i,j + c Si,j S I i,j + c Ei,j E I i,j + c Wi,j W I i,j ) Ix,y t+1 = Ix,y(1 t λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) + λ(c Ni,j I Ni,j + c Si,j I Si,j + c Ei,j I Ei,j + c Wi,j I Wi,j ) Since all c s and λ are positive and between 0,1: I t+1 x,y I t M ij (1 λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) + I t M ij (λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) = I t M ij
31 Discrete Maximum Principle Let I Mij = max(i, I N, I S, I E, I W ) Ix,y t+1 = Ix,y t + λ(c Ni,j N I i,j + c Si,j S I i,j + c Ei,j E I i,j + c Wi,j W I i,j ) Ix,y t+1 = Ix,y(1 t λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) + λ(c Ni,j I Ni,j + c Si,j I Si,j + c Ei,j I Ei,j + c Wi,j I Wi,j ) Since all c s and λ are positive and between 0,1: I t+1 x,y I t M ij (1 λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) + I t M ij (λ(c Ni,j + c Si,j + c Ei,j + c Wi,j )) = I t M ij
32 Adaptive Parameter Setting Set K every iteration Compute a histogram, f i, of I Find K such that 90% of the pixels have gradient magnitude < K. (If b i=1 f i 0.9n 2 then bin b corresponds to gradient magnitude K ).
33 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of the paper... 5 Anisotropic Diffusion
34 Edge Detection Due to Edge enhancement Edge localization edge detection algorithms will benefit from using this nonlinear diffusion process rather than using linear diffusion (Gaussian convolution).
35 Energy Minimization and Markov Random Field Models The anisotropic diffusion algorithm minimizes some energy function. Smoothness conditional dependence on nearest neighbors (Markovian property)
36 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of the paper... 5 Anisotropic Diffusion Introduction
37 Introduction The problem with inhomogeneous diffusion Noise on edges will not be reduced. Even with the anisotropic implementation of Perona-Malik, noise on diagonal edges will not be handled properly.
38 Introduction The problem with inhomogeneous diffusion The inhomogeneous heat equation I t = div(c(x, y, t) I) only controls the magnitude of intensity (or heat) flow. The anisotropic heat equation I t = div(d(x, y, t) I), where D(x, y, t) is a matrix-valued function, can also control the direction of intensity (or heat) flow.
39 Introduction Next: Wednesday: Susan Lemieux will continue the MRI acquisition lecture. Friday: Weickert A review of nonlinear diffusion filtering. The physical laws of heat flow and diffusion. Anisotropic diffusion filtering.
Medical Image Analysis
Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 20th January 2006 Outline 1 Discretizing the heat equation 2 Outline 1 Discretizing the heat
More informationLinear Diffusion. E9 242 STIP- R. Venkatesh Babu IISc
Linear Diffusion Derivation of Heat equation Consider a 2D hot plate with Initial temperature profile I 0 (x, y) Uniform (isotropic) conduction coefficient c Unit thickness (along z) Problem: What is temperature
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More informationLaplacian Filters. Sobel Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters
Sobel Filters Note that smoothing the image before applying a Sobel filter typically gives better results. Even thresholding the Sobel filtered image cannot usually create precise, i.e., -pixel wide, edges.
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationLecture 7: Edge Detection
#1 Lecture 7: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Definition of an Edge First Order Derivative Approximation as Edge Detector #2 This Lecture Examples of Edge Detection
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationNonlinear Diffusion. 1 Introduction: Motivation for non-standard diffusion
Nonlinear Diffusion These notes summarize the way I present this material, for my benefit. But everything in here is said in more detail, and better, in Weickert s paper. 1 Introduction: Motivation for
More informationIntroduction to Computer Vision
Introduction to Computer Vision Michael J. Black Sept 2009 Lecture 8: Pyramids and image derivatives Goals Images as functions Derivatives of images Edges and gradients Laplacian pyramids Code for lecture
More informationReview Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.
Spatial Filtering Dr. Praveen Sankaran Department of ECE NIT Calicut January 7, 203 Outline 2 Linear Nonlinear 3 Spatial Domain Refers to the image plane itself. Direct manipulation of image pixels. Figure:
More informationEdge Detection in Computer Vision Systems
1 CS332 Visual Processing in Computer and Biological Vision Systems Edge Detection in Computer Vision Systems This handout summarizes much of the material on the detection and description of intensity
More informationT H E S I S. Computer Engineering May Smoothing of Matrix-Valued Data. Thomas Brox
D I P L O M A Computer Engineering May 2002 T H E S I S Smoothing of Matrix-Valued Data Thomas Brox Computer Vision, Graphics, and Pattern Recognition Group Department of Mathematics and Computer Science
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
2/9/7 ECE592-064 Digital Image Processing and Introduction to Computer Vision Depart. of ECE, NC State University Instructor: Tianfu (Matt) Wu Spring 207. Recap Outline 2. Sharpening Filtering Illustration
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationIntroduction to Nonlinear Image Processing
Introduction to Nonlinear Image Processing 1 IPAM Summer School on Computer Vision July 22, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Mean and median 2 Observations
More informationPDE Based Image Diffusion and AOS
PDE Based Image Diffusion and AOS by Jarno Ralli April 27, 204 Abstract This technical paper shows how AOS Additive Operator Splitting scheme, together with TDMA Tri-Diagonal Matrix Algorithm, can be used
More informationCN780 Final Lecture. Low-Level Vision, Scale-Space, and Polyakov Action. Neil I. Weisenfeld
CN780 Final Lecture Low-Level Vision, Scale-Space, and Polyakov Action Neil I. Weisenfeld Department of Cognitive and Neural Systems Boston University chapter 14.2-14.3 May 9, 2005 p.1/25 Introduction
More informationTRACKING and DETECTION in COMPUTER VISION Filtering and edge detection
Technischen Universität München Winter Semester 0/0 TRACKING and DETECTION in COMPUTER VISION Filtering and edge detection Slobodan Ilić Overview Image formation Convolution Non-liner filtering: Median
More informationNON-LINEAR DIFFUSION FILTERING
NON-LINEAR DIFFUSION FILTERING Chalmers University of Technology Page 1 Summary Introduction Linear vs Nonlinear Diffusion Non-Linear Diffusion Theory Applications Implementation References Page 2 Introduction
More informationReading. 3. Image processing. Pixel movement. Image processing Y R I G Q
Reading Jain, Kasturi, Schunck, Machine Vision. McGraw-Hill, 1995. Sections 4.-4.4, 4.5(intro), 4.5.5, 4.5.6, 5.1-5.4. 3. Image processing 1 Image processing An image processing operation typically defines
More informationNonlinear diffusion filtering on extended neighborhood
Applied Numerical Mathematics 5 005) 1 11 www.elsevier.com/locate/apnum Nonlinear diffusion filtering on extended neighborhood Danny Barash Genome Diversity Center, Institute of Evolution, University of
More informationGraph Signal Processing for Image Compression & Restoration (Part II)
ene Cheung, Xianming Liu National Institute of Informatics 11 th July, 016 raph Signal Processing for Image Compression & Restoration (Part II). ICME'16 utorial 07/11/016 1 Outline (Part II) Image Restoration
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationEfficient Beltrami Filtering of Color Images via Vector Extrapolation
Efficient Beltrami Filtering of Color Images via Vector Extrapolation Lorina Dascal, Guy Rosman, and Ron Kimmel Computer Science Department, Technion, Institute of Technology, Haifa 32000, Israel Abstract.
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationLecture 8: Interest Point Detection. Saad J Bedros
#1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Review of Edge Detectors #2 Today s Lecture Interest Points Detection What do we mean with Interest Point Detection in an Image Goal:
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
More informationA Four-Pixel Scheme for Singular Differential Equations
A Four-Pixel Scheme for Singular Differential Equations Martin Welk 1, Joachim Weickert 1, and Gabriele Steidl 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer Science, Bldg. 7 Saarland
More informationComplex Diffusion Processes for Image Filtering
Complex Diffusion Processes for Image Filtering Guy Gilboa, Yehoshua Y. Zeevi, and Nir A. Sochen Department of Electrical Engineering, Technion - Israel Institute of Technology Technion City, Haifa 3,
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationConvolutional networks. Sebastian Seung
Convolutional networks Sebastian Seung Convolutional network Neural network with spatial organization every neuron has a location usually on a grid Translation invariance synaptic strength depends on locations
More informationRIGOROUS MATHEMATICAL INVESTIGATION OF A NONLINEAR ANISOTROPIC DIFFUSION-BASED IMAGE RESTORATION MODEL
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 129, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu RIGOROUS MATHEMATICAL
More informationCSE 473/573 Computer Vision and Image Processing (CVIP)
CSE 473/573 Computer Vision and Image Processing (CVIP) Ifeoma Nwogu inwogu@buffalo.edu Lecture 11 Local Features 1 Schedule Last class We started local features Today More on local features Readings for
More informationConsistent Positive Directional Splitting of Anisotropic Diffusion
In: Boštjan Likar (ed.): Proc. of Computer Vision Winter Workshop, Bled, Slovenia, Feb. 7-9, 2001, pp. 37-48. Consistent Positive Directional Splitting of Anisotropic Diffusion Pavel Mrázek and Mirko Navara
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sharpening, image derivatives, Laplacian, edges Revision: 1.2, dated: May 25, 2007 Tomáš Svoboda Czech Technical University, Faculty of Electrical Engineering Center
More informationI Chen Lin, Assistant Professor Dept. of CS, National Chiao Tung University. Computer Vision: 4. Filtering
I Chen Lin, Assistant Professor Dept. of CS, National Chiao Tung University Computer Vision: 4. Filtering Outline Impulse response and convolution. Linear filter and image pyramid. Textbook: David A. Forsyth
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationA BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT
A BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT Aleksandra Pižurica, Iris Vanhamel, Hichem Sahli, Wilfried Philips and Antonis Katartzis Ghent University,
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationLesson 04. KAZE, Non-linear diffusion filtering, ORB, MSER. Ing. Marek Hrúz, Ph.D.
Lesson 04 KAZE, Non-linear diffusion filtering, ORB, MSER Ing. Marek Hrúz, Ph.D. Katedra Kybernetiky Fakulta aplikovaných věd Západočeská univerzita v Plzni Lesson 04 KAZE ORB: an efficient alternative
More informationCS 3710: Visual Recognition Describing Images with Features. Adriana Kovashka Department of Computer Science January 8, 2015
CS 3710: Visual Recognition Describing Images with Features Adriana Kovashka Department of Computer Science January 8, 2015 Plan for Today Presentation assignments + schedule changes Image filtering Feature
More informationErgodicity for Infinite Dimensional Systems
London Mathematical Society Lecture Note Series. 229 Ergodicity for Infinite Dimensional Systems G. DaPrato Scuola Normale Superiore, Pisa J. Zabczyk Polish Academy of Sciences, Warsaw If CAMBRIDGE UNIVERSITY
More informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationImage Enhancement: Methods. Digital Image Processing. No Explicit definition. Spatial Domain: Frequency Domain:
Image Enhancement: No Explicit definition Methods Spatial Domain: Linear Nonlinear Frequency Domain: Linear Nonlinear 1 Spatial Domain Process,, g x y T f x y 2 For 1 1 neighborhood: Contrast Enhancement/Stretching/Point
More informationLecture 8: Interest Point Detection. Saad J Bedros
#1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Last Lecture : Edge Detection Preprocessing of image is desired to eliminate or at least minimize noise effects There is always tradeoff
More informationImage Processing. Waleed A. Yousef Faculty of Computers and Information, Helwan University. April 3, 2010
Image Processing Waleed A. Yousef Faculty of Computers and Information, Helwan University. April 3, 2010 Ch3. Image Enhancement in the Spatial Domain Note that T (m) = 0.5 E. The general law of contrast
More informationMultiple Random Variables
Multiple Random Variables This Version: July 30, 2015 Multiple Random Variables 2 Now we consider models with more than one r.v. These are called multivariate models For instance: height and weight An
More informationDetectors part II Descriptors
EECS 442 Computer vision Detectors part II Descriptors Blob detectors Invariance Descriptors Some slides of this lectures are courtesy of prof F. Li, prof S. Lazebnik, and various other lecturers Goal:
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationOld painting digital color restoration
Old painting digital color restoration Michail Pappas Ioannis Pitas Dept. of Informatics, Aristotle University of Thessaloniki GR-54643 Thessaloniki, Greece Abstract Many old paintings suffer from the
More informationGaussian derivatives
Gaussian derivatives UCU Winter School 2017 James Pritts Czech Tecnical University January 16, 2017 1 Images taken from Noah Snavely s and Robert Collins s course notes Definition An image (grayscale)
More informationImage Gradients and Gradient Filtering Computer Vision
Image Gradients and Gradient Filtering 16-385 Computer Vision What is an image edge? Recall that an image is a 2D function f(x) edge edge How would you detect an edge? What kinds of filter would you use?
More informationJoint Optimization of Segmentation and Appearance Models
Joint Optimization of Segmentation and Appearance Models David Mandle, Sameep Tandon April 29, 2013 David Mandle, Sameep Tandon (Stanford) April 29, 2013 1 / 19 Overview 1 Recap: Image Segmentation 2 Optimization
More informationCorners, Blobs & Descriptors. With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros
Corners, Blobs & Descriptors With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros Motivation: Build a Panorama M. Brown and D. G. Lowe. Recognising Panoramas. ICCV 2003 How do we build panorama?
More informationSolving linear systems (6 lectures)
Chapter 2 Solving linear systems (6 lectures) 2.1 Solving linear systems: LU factorization (1 lectures) Reference: [Trefethen, Bau III] Lecture 20, 21 How do you solve Ax = b? (2.1.1) In numerical linear
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationMath 126 Final Exam Solutions
Math 126 Final Exam Solutions 1. (a) Give an example of a linear homogeneous PE, a linear inhomogeneous PE, and a nonlinear PE. [3 points] Solution. Poisson s equation u = f is linear homogeneous when
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationFast Local Laplacian Filters: Theory and Applications
Fast Local Laplacian Filters: Theory and Applications Mathieu Aubry (INRIA, ENPC), Sylvain Paris (Adobe), Sam Hasinoff (Google), Jan Kautz (UCL), and Frédo Durand (MIT) Input Unsharp Mask, not edge-aware
More informationNon-Local Adaptive Structure Tensors. Application to Anisotropic Diffusion and Shock Filtering
Non-Local Adaptive Structure Tensors. Application to Anisotropic Diffusion and Shock Filtering Vincent Doré, Reza Farrahi Moghaddam, Mohamed Cheriet Figure 1: Line-like structural image and two pixels
More informationSpatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ]
CEE 615: Digital Image Processing Spatial Enhancements 1 Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ] Template (Windowing) Operations Template (window, box, kernel)
More informationImage Segmentation: Definition Importance. Digital Image Processing, 2nd ed. Chapter 10 Image Segmentation.
: Definition Importance Detection of Discontinuities: 9 R = wi z i= 1 i Point Detection: 1. A Mask 2. Thresholding R T Line Detection: A Suitable Mask in desired direction Thresholding Line i : R R, j
More informationLecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008
Lecture 6: Edge Detection CAP 5415: Computer Vision Fall 2008 Announcements PS 2 is available Please read it by Thursday During Thursday lecture, I will be going over it in some detail Monday - Computer
More informationTensor-based Image Diffusions Derived from Generalizations of the Total Variation and Beltrami Functionals
Generalizations of the Total Variation and Beltrami Functionals 29 September 2010 International Conference on Image Processing 2010, Hong Kong Anastasios Roussos and Petros Maragos Computer Vision, Speech
More informationMultivariate Random Variable
Multivariate Random Variable Author: Author: Andrés Hincapié and Linyi Cao This Version: August 7, 2016 Multivariate Random Variable 3 Now we consider models with more than one r.v. These are called multivariate
More informationLocal Enhancement. Local enhancement
Local Enhancement Local Enhancement Median filtering (see notes/slides, 3.5.2) HW4 due next Wednesday Required Reading: Sections 3.3, 3.4, 3.5, 3.6, 3.7 Local Enhancement 1 Local enhancement Sometimes
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationHarris Corner Detector
Multimedia Computing: Algorithms, Systems, and Applications: Feature Extraction By Dr. Yu Cao Department of Computer Science The University of Massachusetts Lowell Lowell, MA 01854, USA Part of the slides
More informationLoG Blob Finding and Scale. Scale Selection. Blobs (and scale selection) Achieving scale covariance. Blob detection in 2D. Blob detection in 2D
Achieving scale covariance Blobs (and scale selection) Goal: independently detect corresponding regions in scaled versions of the same image Need scale selection mechanism for finding characteristic region
More informationEdge Detection. CS 650: Computer Vision
CS 650: Computer Vision Edges and Gradients Edge: local indication of an object transition Edge detection: local operators that find edges (usually involves convolution) Local intensity transitions are
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationConjugate gradient acceleration of non-linear smoothing filters
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Conjugate gradient acceleration of non-linear smoothing filters Knyazev, A.; Malyshev, A. TR25-43 December 25 Abstract The most efficient edge-preserving
More informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationImage Compression Based on Visual Saliency at Individual Scales
Image Compression Based on Visual Saliency at Individual Scales Stella X. Yu 1 Dimitri A. Lisin 2 1 Computer Science Department Boston College Chestnut Hill, MA 2 VideoIQ, Inc. Bedford, MA November 30,
More informationComputer Vision Lecture 3
Computer Vision Lecture 3 Linear Filters 03.11.2015 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Demo Haribo Classification Code available on the class website...
More informationNONLINEAR diffusion denoising is a powerful denoising
IEEE TRANSACTION ON IMAGE PROCESSING, EDICS: TEC-PDE 1 An Edge Adapting Laplacian Kernel For Nonlinear Diffusion Filters Mohammad Reza Hajiaboli, Member, IEEE, M. Omair Ahmad, Fellow, IEEE, Chunyan Wang,
More informationMath 180, Final Exam, Fall 2007 Problem 1 Solution
Problem Solution. Differentiate with respect to x. Write your answers showing the use of the appropriate techniques. Do not simplify. (a) x 27 x 2/3 (b) (x 2 2x + 2)e x (c) ln(x 2 + 4) (a) Use the Power
More informationComputational Physics
Interpolation, Extrapolation & Polynomial Approximation Lectures based on course notes by Pablo Laguna and Kostas Kokkotas revamped by Deirdre Shoemaker Spring 2014 Introduction In many cases, a function
More informationLaplacian Eigenmaps for Dimensionality Reduction and Data Representation
Introduction and Data Representation Mikhail Belkin & Partha Niyogi Department of Electrical Engieering University of Minnesota Mar 21, 2017 1/22 Outline Introduction 1 Introduction 2 3 4 Connections to
More informationLecture 1: Graphs, Adjacency Matrices, Graph Laplacian
Lecture 1: Graphs, Adjacency Matrices, Graph Laplacian Radu Balan January 31, 2017 G = (V, E) An undirected graph G is given by two pieces of information: a set of vertices V and a set of edges E, G =
More informationLocal enhancement. Local Enhancement. Local histogram equalized. Histogram equalized. Local Contrast Enhancement. Fig 3.23: Another example
Local enhancement Local Enhancement Median filtering Local Enhancement Sometimes Local Enhancement is Preferred. Malab: BlkProc operation for block processing. Left: original tire image. 0/07/00 Local
More informationSpace-Variant Computer Vision: A Graph Theoretic Approach
p.1/65 Space-Variant Computer Vision: A Graph Theoretic Approach Leo Grady Cognitive and Neural Systems Boston University p.2/65 Outline of talk Space-variant vision - Why and how of graph theory Anisotropic
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationMachine vision. Summary # 4. The mask for Laplacian is given
1 Machine vision Summary # 4 The mask for Laplacian is given L = 0 1 0 1 4 1 (6) 0 1 0 Another Laplacian mask that gives more importance to the center element is L = 1 1 1 1 8 1 (7) 1 1 1 Note that the
More informationBlobs & Scale Invariance
Blobs & Scale Invariance Prof. Didier Stricker Doz. Gabriele Bleser Computer Vision: Object and People Tracking With slides from Bebis, S. Lazebnik & S. Seitz, D. Lowe, A. Efros 1 Apertizer: some videos
More informationScale and Statistics in Variable Conductance Diffusion
Scale and Statistics in Variable Conductance Diffusion Terry S. Yoo Department of Computer Science University of North Carolina Chapel Hill, NC 7599-3175 Introduction Non-linear diffusion or Variable Conductance
More informationScale Space Analysis by Stabilized Inverse Diffusion Equations
Scale Space Analysis by Stabilized Inverse Diffusion Equations Ilya Pollak, Alan S. Willsky and Hamid Krim Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, 77 Massachusetts
More informationA Partial Differential Equation Approach to Image Zoom
A Partial Differential Equation Approach to Image Zoom Abdelmounim Belahmidi and Frédéric Guichard January 2004 Abstract We propose a new model for zooming digital image. This model, driven by a partial
More informationEdge Detection. Introduction to Computer Vision. Useful Mathematics Funcs. The bad news
Edge Detection Introduction to Computer Vision CS / ECE 8B Thursday, April, 004 Edge detection (HO #5) Edge detection is a local area operator that seeks to find significant, meaningful changes in image
More informationModelling Non-linear and Non-stationary Time Series
Modelling Non-linear and Non-stationary Time Series Chapter 2: Non-parametric methods Henrik Madsen Advanced Time Series Analysis September 206 Henrik Madsen (02427 Adv. TS Analysis) Lecture Notes September
More informationEdges and Scale. Image Features. Detecting edges. Origin of Edges. Solution: smooth first. Effects of noise
Edges and Scale Image Features From Sandlot Science Slides revised from S. Seitz, R. Szeliski, S. Lazebnik, etc. Origin of Edges surface normal discontinuity depth discontinuity surface color discontinuity
More informationLecture 3: Compressive Classification
Lecture 3: Compressive Classification Richard Baraniuk Rice University dsp.rice.edu/cs Compressive Sampling Signal Sparsity wideband signal samples large Gabor (TF) coefficients Fourier matrix Compressive
More informationEnergy-Based Image Simplification with Nonlocal Data and Smoothness Terms
Energy-Based Image Simplification with Nonlocal Data and Smoothness Terms Stephan Didas 1, Pavel Mrázek 2, and Joachim Weickert 1 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer
More informationDigital Image Processing. Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009
Digital Image Processing Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 009 Outline Image Enhancement in Spatial Domain Spatial Filtering Smoothing Filters Median Filter
More information