Numerical Methods on the Image Processing Problems
|
|
- Erin McBride
- 5 years ago
- Views:
Transcription
1 Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006
2 Objective Develop efficient PDE (partial differential equations) based mathematical models and their numerical algorithms for 1 Noise removal Enhance the quality of images 2 Image segmentation Edge (2D) or surface (3D) detection
3 Objective Develop efficient PDE (partial differential equations) based mathematical models and their numerical algorithms for 1 Noise removal Enhance the quality of images 2 Image segmentation Edge (2D) or surface (3D) detection
4 Objective Develop efficient PDE (partial differential equations) based mathematical models and their numerical algorithms for 1 Noise removal Enhance the quality of images 2 Image segmentation Edge (2D) or surface (3D) detection
5 Applications - Noise Removal Image with 20% impulse noise (left) and denoised image (right)
6 Applications - Noise Removal Image with 10% impulse noise (left) and denoised image (right)
7 Applications - Noise Removal Color image with impulse noise (left) and denoised image (right)
8 Applications - Edge detection
9 Applications - Edge detection Medical Imaging
10 Applications - Edge detection Image Analysis in Materials Science
11 Outline 1 History - PDE based Mathematical Image Processing 2 Image Denoising Conventional approach - Total variation (TV) minimization New models and their numerical procedure - Non-convex diffusion model - Kim and Lim ( 05) - Anisotropic diffusion model - Lim and Williams ( 06) 3 Image Segmentation Conventional Approach Method of Background Subtraction (MBS) for Medical Image Segmentation - Kim and Lim ( 05) 4 Conclusions
12 Outline 1 History - PDE based Mathematical Image Processing 2 Image Denoising Conventional approach - Total variation (TV) minimization New models and their numerical procedure - Non-convex diffusion model - Kim and Lim ( 05) - Anisotropic diffusion model - Lim and Williams ( 06) 3 Image Segmentation Conventional Approach Method of Background Subtraction (MBS) for Medical Image Segmentation - Kim and Lim ( 05) 4 Conclusions
13 Outline 1 History - PDE based Mathematical Image Processing 2 Image Denoising Conventional approach - Total variation (TV) minimization New models and their numerical procedure - Non-convex diffusion model - Kim and Lim ( 05) - Anisotropic diffusion model - Lim and Williams ( 06) 3 Image Segmentation Conventional Approach Method of Background Subtraction (MBS) for Medical Image Segmentation - Kim and Lim ( 05) 4 Conclusions
14 Outline 1 History - PDE based Mathematical Image Processing 2 Image Denoising Conventional approach - Total variation (TV) minimization New models and their numerical procedure - Non-convex diffusion model - Kim and Lim ( 05) - Anisotropic diffusion model - Lim and Williams ( 06) 3 Image Segmentation Conventional Approach Method of Background Subtraction (MBS) for Medical Image Segmentation - Kim and Lim ( 05) 4 Conclusions
15 Outline 1 History - PDE based Mathematical Image Processing 2 Image Denoising Conventional approach - Total variation (TV) minimization New models and their numerical procedure - Non-convex diffusion model - Kim and Lim ( 05) - Anisotropic diffusion model - Lim and Williams ( 06) 3 Image Segmentation Conventional Approach Method of Background Subtraction (MBS) for Medical Image Segmentation - Kim and Lim ( 05) 4 Conclusions
16 History of PDE based Mathematical Image Processing Short history, but has strong impact Image denoising, deconvolution (deblurring), segmentation (edge/surface detection) Image denoising - Total variation (TV) minimization (Osher ( 92), Lions ( 97), Chan ( 98), Kim ( 01)) - Weakness: Edges of images can be easily smeared out due to diffusion property of TV minimization original image (left) and smeared image (right)
17 History of PDE based Mathematical Image Processing Short history, but has strong impact Image denoising, deconvolution (deblurring), segmentation (edge/surface detection) Image denoising - Total variation (TV) minimization (Osher ( 92), Lions ( 97), Chan ( 98), Kim ( 01)) - Weakness: Edges of images can be easily smeared out due to diffusion property of TV minimization original image (left) and smeared image (right)
18 History of PDE based Mathematical Image Processing Short history, but has strong impact Image denoising, deconvolution (deblurring), segmentation (edge/surface detection) Image denoising - Total variation (TV) minimization (Osher ( 92), Lions ( 97), Chan ( 98), Kim ( 01)) - Weakness: Edges of images can be easily smeared out due to diffusion property of TV minimization original image (left) and smeared image (right)
19 History of PDE based Mathematical Image Processing Short history, but has strong impact Image denoising, deconvolution (deblurring), segmentation (edge/surface detection) Image denoising - Total variation (TV) minimization (Osher ( 92), Lions ( 97), Chan ( 98), Kim ( 01)) - Weakness: Edges of images can be easily smeared out due to diffusion property of TV minimization original image (left) and smeared image (right)
20 History of PDE based Mathematical Image Processing Color image denoising (Kim ( 02), Osher ( 03)) - Use RGB color component Color image with 15% impulse noise (left) and denoised image (right)
21 Image Denoising Conventional approach TV minimization model u t σ u γ+1 ( ) u = β (u o u) u Efficiently removes noise Image loses sharpness since the frequency of noise and edges are similar Produces a staircasing (locally constant) effect and nonphysical dissipation
22 Image Denoising Conventional approach TV minimization model u t σ u γ+1 ( ) u = β (u o u) u Efficiently removes noise Image loses sharpness since the frequency of noise and edges are similar Produces a staircasing (locally constant) effect and nonphysical dissipation
23 Image Denoising Conventional approach TV minimization model u t σ u γ+1 ( ) u = β (u o u) u Efficiently removes noise Image loses sharpness since the frequency of noise and edges are similar Produces a staircasing (locally constant) effect and nonphysical dissipation
24 Image Denoising Conventional approach TV minimization model u t σ u γ+1 ( ) u = β (u o u) u Efficiently removes noise Image loses sharpness since the frequency of noise and edges are similar Produces a staircasing (locally constant) effect and nonphysical dissipation
25 Image Denoising - NC Model Non-convex diffusion model Control of nonphysical dissipation Consider min F ɛ,p (u), where u F ɛ,p (u) = Ω ɛ u p dx + λ 2 f u 2. Then ( ) u p ɛ u 2 p λ(f u) = 0, where ɛ u = (u 2 x + u 2 y + ɛ 2 ) 1/2. Sharp (left) and blurry image (right)
26 Image Denoising - NC Model Non-convex diffusion model F ɛ,p Sharp image Blurry image F 0, = = 0.5 F 0, = = 1 F 0.1, F 0, = The strictly convex minimization (p > 1) makes images blurrier. The TV model itself (p = 1 and ɛ = 0) may not introduce blur but its regularization (p = 1 and ɛ > 0) does. When p < 1(non-convex), the model can make the image sharper.
27 Image Denoising - NC Model Non-convex diffusion model F ɛ,p Sharp image Blurry image F 0, = = 0.5 F 0, = = 1 F 0.1, F 0, = The strictly convex minimization (p > 1) makes images blurrier. The TV model itself (p = 1 and ɛ = 0) may not introduce blur but its regularization (p = 1 and ɛ > 0) does. When p < 1(non-convex), the model can make the image sharper.
28 Image Denoising - NC Model Non-convex diffusion model F ɛ,p Sharp image Blurry image F 0, = = 0.5 F 0, = = 1 F 0.1, F 0, = The strictly convex minimization (p > 1) makes images blurrier. The TV model itself (p = 1 and ɛ = 0) may not introduce blur but its regularization (p = 1 and ɛ > 0) does. When p < 1(non-convex), the model can make the image sharper.
29 Image Denoising - NC Model Non-convex diffusion model F ɛ,p Sharp image Blurry image F 0, = = 0.5 F 0, = = 1 F 0.1, F 0, = The strictly convex minimization (p > 1) makes images blurrier. The TV model itself (p = 1 and ɛ = 0) may not introduce blur but its regularization (p = 1 and ɛ > 0) does. When p < 1(non-convex), the model can make the image sharper.
30 Image Denoising - NC Model Non-convex diffusion model New non-convex (NC) model ( ) u u t ɛ u 1+ω ɛ u 1+ω = β (f u), ω ( 1, 1), β 0 Numerical procedures Linearized θ- method. Alternating Direction Implicit (ADI) method Theorem (Stability) The θ- method for the new NC model is stable and holds the maximum principle.
31 Image Denoising - NC Model Non-convex diffusion model New non-convex (NC) model ( ) u u t ɛ u 1+ω ɛ u 1+ω = β (f u), ω ( 1, 1), β 0 Numerical procedures Linearized θ- method. Alternating Direction Implicit (ADI) method Theorem (Stability) The θ- method for the new NC model is stable and holds the maximum principle.
32 Image Denoising - NC Model Non-convex diffusion model New non-convex (NC) model ( ) u u t ɛ u 1+ω ɛ u 1+ω = β (f u), ω ( 1, 1), β 0 Numerical procedures Linearized θ- method. Alternating Direction Implicit (ADI) method Theorem (Stability) The θ- method for the new NC model is stable and holds the maximum principle.
33 Numerical Experiments - NC Model Image with 20% mean zero noise (left), conventional method (middle), new NC method (right)
34 Numerical Experiments - NC Model Image with 20% mean zero noise (left), conventional method (middle), new NC method (right)
35 Numerical Experiments - NC Model Original brain image (left) and enhanced image by new NC method (right)
36 Numerical Experiments - NC Model Horizontal line cuts of brain image and its restored images
37 Image Denoising - AD Model Anisotropic diffusion model Speckle noise is multiplicative and it can be modeled as (Krissian et. al. 04, 05): f = u + un, Corresponding time marching equation: u t u2 ( u ) f + u u u = λ u (f u). u2 f +u u/2 makes the diffusion faster in the lighter region (where the image values are high) and slower in the darker region (where the image values are low). unrealistic and ineffective in practice!
38 Image Denoising - AD Model Anisotropic diffusion model Speckle noise is multiplicative and it can be modeled as (Krissian et. al. 04, 05): f = u + un, Corresponding time marching equation: u t u2 ( u ) f + u u u = λ u (f u). u/2 makes the diffusion faster in the lighter region (where the image values are high) and slower in the darker region (where the image values are low). unrealistic and ineffective in practice! u2 f +u
39 Image Denoising - AD Model Anisotropic diffusion model Speckle noise is multiplicative and it can be modeled as (Krissian et. al. 04, 05): f = u + un, Corresponding time marching equation: u t u2 ( u ) f + u u u = λ u (f u). u/2 makes the diffusion faster in the lighter region (where the image values are high) and slower in the darker region (where the image values are low). unrealistic and ineffective in practice! u2 f +u
40 Image Denoising - AD Model Anisotropic diffusion model Speckle noise is multiplicative and it can be modeled as (Krissian et. al. 04, 05): f = u + un, Corresponding time marching equation: u t u2 ( u ) f + u u u = λ u (f u). u/2 makes the diffusion faster in the lighter region (where the image values are high) and slower in the darker region (where the image values are low). unrealistic and ineffective in practice! u2 f +u
41 Image Denoising - AD Model Consider f = u + ( u u s ) n, where u s : smoothed version of the noised image f. New Anisotropic Diffusion (AD) Model u ( u ) t C u u s α u = β (f u), C > 0, 1/2 < α < 2 u On the regions where noise is present, u u s is relatively big. Diffusion is big enough to reduce the noise efficiently. On the regions where noise is not present, u u s is small. Diffusion is relatively slower.
42 Image Denoising - AD Model Consider f = u + ( u u s ) n, where u s : smoothed version of the noised image f. New Anisotropic Diffusion (AD) Model u ( u ) t C u u s α u = β (f u), C > 0, 1/2 < α < 2 u On the regions where noise is present, u u s is relatively big. Diffusion is big enough to reduce the noise efficiently. On the regions where noise is not present, u u s is small. Diffusion is relatively slower.
43 Image Denoising - AD Model Consider f = u + ( u u s ) n, where u s : smoothed version of the noised image f. New Anisotropic Diffusion (AD) Model u ( u ) t C u u s α u = β (f u), C > 0, 1/2 < α < 2 u On the regions where noise is present, u u s is relatively big. Diffusion is big enough to reduce the noise efficiently. On the regions where noise is not present, u u s is small. Diffusion is relatively slower.
44 Image Denoising - AD Model Consider f = u + ( u u s ) n, where u s : smoothed version of the noised image f. New Anisotropic Diffusion (AD) Model u ( u ) t C u u s α u = β (f u), C > 0, 1/2 < α < 2 u On the regions where noise is present, u u s is relatively big. Diffusion is big enough to reduce the noise efficiently. On the regions where noise is not present, u u s is small. Diffusion is relatively slower.
45 Numerical Procedure - AD Model 1 TFR (texture-free residual) parametrization 1 Set β as a constant: β(x, 0) = β 0. 2 For n = 2, 3, Compute the absolute residual and a quantity G n 1 Res : R n 1 = f u n 1, G n 1 Res = max 0, S m(r n 1 ) R n 1, where S m is a smoothing operator and R n 1 denotes the L 2 -average of R n 1. Update: β n = β n 1 + γ n G n 1 Res, where γ n is a scaling factor having the property: γ n 0 as n.
46 Numerical Experiments - AD Model Cuba Missile Crisis: The original (left above) and restored images by using ITV (right above), ITV-TFR (left below), and AD-TFR (right below) model
47 Medical Image Segmentation - MBS Motivation Medical images can involve noise, diverse artifacts, and unclear edges. Conventional segmentation methods show difficulties when applied to medical imagery. When an appropriate background is subtracted from the given image, the residue can be considered as an essentially binary image.
48 Medical Image Segmentation - MBS Procedure of method of background subtraction
49 Medical Image Segmentation - MBS Procedure of the construction of background 1 Select a coarse mesh {Ω ij } for the image domain Ω and choose a coarse image U c on {Ω ij }. Each element Ω ij in the coarse mesh corresponds to m x m y pixels. 2 Smooth U c. 3 Prolongate U c to the original mesh Ω, for U f. 4 Smooth U f. Assign the result for the background Ũ. Strategies for background construction In step I, choose U c on Ω ij as ra ij + (1 r)m ij, 0 r 1, a ij : arithmetic average, m ij : minimum of U 0 on Ω ij. Ũ must contain only background Information, not objects information. Thus select m = m x = m y such that number of blocks in U c corresponding to objects are smaller than the number of smoothing iterations in step II.
50 Medical Image Segmentation - MBS Procedure of the construction of background 1 Select a coarse mesh {Ω ij } for the image domain Ω and choose a coarse image U c on {Ω ij }. Each element Ω ij in the coarse mesh corresponds to m x m y pixels. 2 Smooth U c. 3 Prolongate U c to the original mesh Ω, for U f. 4 Smooth U f. Assign the result for the background Ũ. Strategies for background construction In step I, choose U c on Ω ij as ra ij + (1 r)m ij, 0 r 1, a ij : arithmetic average, m ij : minimum of U 0 on Ω ij. Ũ must contain only background Information, not objects information. Thus select m = m x = m y such that number of blocks in U c corresponding to objects are smaller than the number of smoothing iterations in step II.
51 Numerical Experiments - MBS Heart: Original (left), Conventional Method (middle), Conventional approach with MBS (right)
52 Numerical Experiments - MBS Hand: Original (left), Conventional Method (middle), Conventional approach with MBS (right)
53 Numerical Experiments - MBS Leukemia: Original (left), Conventional Method (middle), Conventional approach with MBS (right)
54 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
55 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
56 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
57 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
58 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
59 Conclusions 1 Regularization of u can be a significant source of nonphysical dissipation in TV-based models. 2 Non-convex (NC) diffusion model has been introduced in order to simultaneously suppress the noise and enhance edges. 3 Anisotropic diffusion (AD) model is more efficient on SAR image denoising than conventional models. 4 Numerical procedures of NC and AD models are stable. 5 The method of background subtraction (MBS) can be efficiently used as a pre-process of various segmentation methods for medical image segmentation.
Method of Background Subtraction for Medical Image Segmentation
Method of Background Subtraction for Medical Image Segmentation Seongjai Kim Department of Mathematics and Statistics, Mississippi State University Mississippi State, MS 39762, USA and Hyeona Lim Department
More informationIterative Refinement by Smooth Curvature Correction for PDE-based Image Restoration
Iterative Refinement by Smooth Curvature Correction for PDE-based Image Restoration Anup Gautam 1, Jihee Kim 2, Doeun Kwak 3, and Seongjai Kim 4 1 Department of Mathematics and Statistics, Mississippi
More informationImage Processing and related PDEs Lecture 1: Introduction to image processing
Image Processing and related PDEs Lecture 1: Introduction to image processing Yves van Gennip School of Mathematical Sciences, University of Nottingham Minicourse on Image Processing and related PDEs University
More informationSimulation of a Hybrid Model for Image Denoising
Simulation of a Hybrid Model for Image Denoising Ricolindo L. Cariño 1, Ioana Banicescu 2,1, Hyeona Lim 3, Neil Williams 3 and Seongjai Kim 4 1 Center for Computational Sciences ERC, Mississippi State
More informationExact discrete minimization for TV+L0 image decomposition models
Exact discrete minimization for TV+L0 image decomposition models Loïc Denis 1, Florence Tupin 2 and Xavier Rondeau 2 1. Observatory of Lyon (CNRS / Univ. Lyon 1 / ENS de Lyon), France 2. Telecom ParisTech
More informationImage Restoration under Significant Additive Noise
1 Image Restoration under Significant Additive Noise Wenyi Zhao 1 and Art Pope Sarnoff Corporation 01 Washington Road, Princeton, NJ 08540, USA email: { wzhao, apope }@ieee.org Tel: 408-53-178 Abstract
More informationGRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS. Andrey Nasonov, and Andrey Krylov
GRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS Andrey Nasonov, and Andrey Krylov Lomonosov Moscow State University, Moscow, Department of Computational Mathematics and Cybernetics, e-mail: nasonov@cs.msu.ru,
More informationImage denoising using TV-Stokes equation with an orientation-matching minimization
Image denoising using TV-Stokes equation with an orientation-matching minimization Xue-Cheng Tai 1,2, Sofia Borok 1, and Jooyoung Hahn 1 1 Division of Mathematical Sciences, School of Physical Mathematical
More informationFilters. Advanced and Special Topics: Filters. Filters
Filters Advanced and Special Topics: Filters Dr. Edmund Lam Department of Electrical and Electronic Engineering The University of Hong Kong ELEC4245: Digital Image Processing (Second Semester, 2016 17)
More informationAn edge-preserving multilevel method for deblurring, denoising, and segmentation
An edge-preserving multilevel method for deblurring, denoising, and segmentation S. Morigi 1, L. Reichel 2, and F. Sgallari 1 1 Dept. of Mathematics-CIRAM, University of Bologna, Bologna, Italy {morigi,sgallari}@dm.unibo.it
More informationEdge Detection Free Postprocessing for Pseudospectral Approximations
Edge Detection Free Postprocessing for Pseudospectral Approximations Scott A. Sarra March 4, 29 Abstract Pseudospectral Methods based on global polynomial approximation yield exponential accuracy when
More informationBilevel Sparse Coding
Adobe Research 345 Park Ave, San Jose, CA Mar 15, 2013 Outline 1 2 The learning model The learning algorithm 3 4 Sparse Modeling Many types of sensory data, e.g., images and audio, are in high-dimensional
More informationResearch Article Image Enhancement under Data-Dependent Multiplicative Gamma Noise
Applied Computational Intelligence and So Computing, Article ID 981932, 8 pages http://dx.doi.org/10.1155/2014/981932 Research Article Image Enhancement under Data-Dependent Multiplicative Gamma Noise
More informationRemoving a mixture of Gaussian and impulsive noise using the total variation functional and split Bregman iterative method
ANZIAM J. 56 (CTAC2014) pp.c52 C67, 2015 C52 Removing a mixture of Gaussian and impulsive noise using the total variation functional and split Bregman iterative method Bishnu P. Lamichhane 1 (Received
More informationIMAGE FUSION WITH SIMULTANEOUS CARTOON AND TEXTURE DECOMPOSITION MAHDI DODANGEH, ISABEL NARRA FIGUEIREDO AND GIL GONÇALVES
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 15 14 IMAGE FUSION WITH SIMULTANEOUS CARTOON AND TEXTURE DECOMPOSITION MAHDI DODANGEH, ISABEL NARRA FIGUEIREDO AND
More informationImage Denoising Algorithms
Image Denoising Algorithms Xiang Hao School of Compting, University of Utah, USA, hao@cs.tah.ed Abstract. This is a report of an assignment of the class Mathematics of Imaging. In this assignment, we first
More informationNonlocal Multiscale Hierarchical Decomposition on Graphs
Nonlocal Multiscale Hierarchical Decomposition on Graphs Moncef Hidane, Olivier Lézoray, Vinh-Thong Ta, and Abderrahim Elmoataz Université de Caen Basse-Normandie, ENSICAEN, CNRS, GREYC Image Team, 6 Boulevard
More informationRunge Kutta Methods Optimized For Advection Problems
Runge Kutta Methods Optimized For Advection Problems Brown University December 9, 2009 Introduction Classical Runge Kutta Methods Overview Standard Runge Kutta Schemes for Advection Equation Low Dissipation
More informationAn Integrated System for Digital Restoration of Prehistoric Theran Wall Paintings
An Integrated System for Digital Restoration of Prehistoric Theran Wall Paintings Nikolaos Karianakis 1 Petros Maragos 2 1 University of California, Los Angeles 2 National Technical University of Athens
More informationGlobal Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising
Global Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising Shingyu Leung and Stanley Osher Department of Mathematics, UCLA, Los Angeles, CA 90095, USA {syleung, sjo}@math.ucla.edu
More informationAn Enhanced Primal Dual Method For Total Variation-Based Wavelet Domain Inpainting
An Enhanced Primal Dual Method For Total Variation-Based Wavelet Domain Inpainting P. Sravan Kumar 1, H.V. Ram Ravi Kishore 2, V. Ravi Kumar 3 1 MTech, Spatial Information Technology, from JNTU, research
More informationImage Smoothing and Segmentation by Graph Regularization
Image Smoothing and Segmentation by Graph Regularization Sébastien Bougleux 1 and Abderrahim Elmoataz 1 GREYC CNRS UMR 6072, Université de Caen Basse-Normandie ENSICAEN 6 BD du Maréchal Juin, 14050 Caen
More informationICCV 2013 Supplementary Material Super-resolution via Transform-invariant Group-sparse Regularization
ICCV 2013 Supplementary Material Super-resolution via Transform-invariant Group-sparse Regularization October 4, 2013 1 Contents This is the supplementary material for the paper Super-resolution via Transform-invariant
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationProjection Technique for Vortex-Free Image Registration
Projection Technique for Vortex-Free Image Registration Patrick Scheibe 1, Ulf-Dietrich Braumann 1,2, Jens-Peer Kuska 2 1 Translational Center for Regenerative Medicine (TRM) 2 Interdisciplinary Center
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationProject Updates Short lecture Volumetric Modeling +2 papers
Volumetric Modeling Schedule (tentative) Feb 20 Feb 27 Mar 5 Introduction Lecture: Geometry, Camera Model, Calibration Lecture: Features, Tracking/Matching Mar 12 Mar 19 Mar 26 Apr 2 Apr 9 Apr 16 Apr 23
More informationMEDICAL IMAGE ENHANCEMENT BASED ON AN EFFICIENT APPROACH FOR ADAPTIVE ANISOTROPIC DIFFUSION
MEDICAL IMAGE ENHANCEMENT BASED ON AN EFFICIENT APPROACH FOR ADAPTIVE ANISOTROPIC DIFFUSION Shokhan Mahmoud Hama 1 and Muzhir Shaban Al-Ani 2 University of Al-Anbar, Collage of Computer, Anbar, Iraq ABSTRACT
More informationImaging of flow in porous media - from optimal transport to prediction
Imaging of flow in porous media - from optimal transport to prediction Eldad Haber Dept of EOS and Mathematics, UBC October 15, 2013 With Rowan Lars Jenn Cocket Ruthotto Fohring Outline Prediction is very
More informationCurvelet-based non-linear adaptive subtraction with sparseness constraints. Felix J Herrmann, Peyman P Moghaddam (EOS-UBC)
Curvelet-based non-linear adaptive subtraction with sparseness constraints Felix J Herrmann, Peyman P Moghaddam (EOS-UBC) Context Spiky/minimal structure deconvolution (Claerbout, Ulrych, Oldenburg, Sacchi,
More informationMULTISCALE IMAGE REGISTRATION. Dana Paquin. Doron Levy. Eduard Schreibmann. Lei Xing. (Communicated by Yang Kuang)
MATHEMATICAL BIOSCIENCES http://www.mbejournal.org/ AND ENGINEERING Volume 3, Number 2, April 2006 pp. 389 418 MULTISCALE IMAGE REGISTRATION Dana Paquin Department of Mathematics, Stanford University Stanford,
More informationImage Denoising Using Mean Curvature of Image Surface. Tony Chan (HKUST)
Image Denoising Using Mean Curvature of Image Surface Tony Chan (HKUST) Joint work with Wei ZHU (U. Alabama) & Xue-Cheng TAI (U. Bergen) In honor of Bob Plemmons 75 birthday CUHK Nov 8, 03 How I Got to
More informationIterative Image Restoration Combining Total Variation Minimization and a Second-Order Functional
International Journal of Computer Vision 661, 5 18, 2006 c 2006 Springer Science + Business Media, Inc. Manufactured in The Netherlands. DOI: 10.1007/s11263-005-3219-7 Iterative Image Restoration Combining
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 informationing and enhancing operations for shapes represented by level sets (isosurfaces) and applied unsharp masking to the level set normals. Their approach t
Shape Deblurring with Unsharp Masking Applied to Mesh Normals Hirokazu Yagou Λ Alexander Belyaev y Daming Wei z Λ y z ; ; Shape Modeling Laboratory, University of Aizu, Aizu-Wakamatsu 965-8580 Japan fm50534,
More informationA Comparative Study & Analysis of Image Restoration by Non Blind Technique
A Comparative Study & Analysis of Image Restoration by Non Blind Technique Saurav Rawat 1, S.N.Tazi 2 M.Tech Student, Assistant Professor, CSE Department, Government Engineering College, Ajmer Abstract:
More informationProf. Feng Liu. Winter /15/2019
Prof. Feng Liu Winter 2019 http://www.cs.pdx.edu/~fliu/courses/cs410/ 01/15/2019 Last Time Filter 2 Today More on Filter Feature Detection 3 Filter Re-cap noisy image naïve denoising Gaussian blur better
More informationApplication of Proximal Algorithms to Three Dimensional Deconvolution Microscopy
Application of Proximal Algorithms to Three Dimensional Deconvolution Microscopy Paroma Varma Stanford University paroma@stanford.edu Abstract In microscopy, shot noise dominates image formation, which
More informationCAAM Technical Report TR A Fast Algorithm for Image Deblurring with Total Variation Regularization
CAAM Technical Report TR07-10 A Fast Algorithm for Image Deblurring with Total Variation Regularization Yilun Wang, Wotao Yin and Yin Zhang Department of Computational and Applied Mathematics Rice University,
More information1.7.1 Laplacian Smoothing
1.7.1 Laplacian Smoothing 320491: Advanced Graphics - Chapter 1 434 Theory Minimize energy functional total curvature estimate by polynomial-fitting non-linear (very slow!) 320491: Advanced Graphics -
More informationFinding a Needle in a Haystack: An Image Processing Approach
Finding a Needle in a Haystack: An Image Processing Approach Emily Beylerian Advisor: Hayden Schaeffer University of California, Los Angeles Department of Mathematics Abstract Image segmentation (also
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 informationImage Restoration. Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction
Image Restoration Image Restoration Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction Diffusion Term Consider only the regularization term E-L equation: (Laplace equation) Steepest
More informationSolving Vision Tasks with variational methods on the GPU
Solving Vision Tasks with variational methods on the GPU Horst Bischof Inst. f. Computer Graphics and Vision Graz University of Technology Joint work with Thomas Pock, Markus Unger, Arnold Irschara and
More informationYunyun Yang, Chunming Li, Chiu-Yen Kao and Stanley Osher. Speaker: Chiu-Yen Kao (Math Department, The Ohio State University) BIRS, Banff, Canada
Yunyun Yang, Chunming Li, Chiu-Yen Kao and Stanley Osher Speaker: Chiu-Yen Kao (Math Department, The Ohio State University) BIRS, Banff, Canada Outline Review of Region-based Active Contour Models Mumford
More informationDigital Image Processing ERRATA. Wilhelm Burger Mark J. Burge. An algorithmic introduction using Java. Second Edition. Springer
Wilhelm Burger Mark J. Burge Digital Image Processing An algorithmic introduction using Java Second Edition ERRATA Springer Berlin Heidelberg NewYork Hong Kong London Milano Paris Tokyo 12.1 RGB Color
More informationNormalized cuts and image segmentation
Normalized cuts and image segmentation Department of EE University of Washington Yeping Su Xiaodan Song Normalized Cuts and Image Segmentation, IEEE Trans. PAMI, August 2000 5/20/2003 1 Outline 1. Image
More informationProblèmes Inverses et Imagerie. Journées de rentrée des Masters, 6-8 Septembre 2017, IHES, Bures-sur-Yvette France
Problèmes Inverses et Imagerie Ce que les images ne nous disent pas. Luca Calatroni Centre de Mathematiqués Appliquées (CMAP), École Polytechnique CNRS Journées de rentrée des Masters, 6-8 Septembre 2017,
More informationImage Deconvolution.
Image Deconvolution. Mathematics of Imaging. HW3 Jihwan Kim Abstract This homework is to implement image deconvolution methods, especially focused on a ExpectationMaximization(EM) algorithm. Most of this
More informationHigher Order Partial Differential Equation Based Method for Image Enhancement
Higher Order Partial Differential Equation Based Method for Image Enhancement S.Balamurugan #1, R.Sowmyalakshmi *2 S.BALAMURUGAN, M.E (M.B.C.B.S), E.C.E, University college of Engineering, BIT Campus,
More informationAnnouncements. Edges. Last Lecture. Gradients: Numerical Derivatives f(x) Edge Detection, Lines. Intro Computer Vision. CSE 152 Lecture 10
Announcements Assignment 2 due Tuesday, May 4. Edge Detection, Lines Midterm: Thursday, May 6. Introduction to Computer Vision CSE 152 Lecture 10 Edges Last Lecture 1. Object boundaries 2. Surface normal
More informationA Total Variation-Morphological Image Edge Detection Approach
A Total Variation-Morphological Image Edge Detection Approach Peter Ndajah, Hisakazu Kikuchi, Shogo Muramatsu, Masahiro Yukawa and Francis Benyah Abstract: We present image edge detection using the total
More informationUsing Subspace Constraints to Improve Feature Tracking Presented by Bryan Poling. Based on work by Bryan Poling, Gilad Lerman, and Arthur Szlam
Presented by Based on work by, Gilad Lerman, and Arthur Szlam What is Tracking? Broad Definition Tracking, or Object tracking, is a general term for following some thing through multiple frames of a video
More informationDenoising an Image by Denoising its Components in a Moving Frame
Denoising an Image by Denoising its Components in a Moving Frame Gabriela Ghimpețeanu 1, Thomas Batard 1, Marcelo Bertalmío 1, and Stacey Levine 2 1 Universitat Pompeu Fabra, Spain 2 Duquesne University,
More informationTwo-Phase Kernel Estimation for Robust Motion Deblurring
Two-Phase Kernel Estimation for Robust Motion Deblurring Li Xu and Jiaya Jia Department of Computer Science and Engineering The Chinese University of Hong Kong {xuli,leojia}@cse.cuhk.edu.hk Abstract. We
More informationAutomatically Extracting Cellular Structures from Images Generated via Electron Microscopy
Automatically Extracting Cellular Structures from Images Generated via Electron Microscopy Alejandro Cantarero Department of Computer Science University of Colorado, Boulder May 9, 2005 1 Introduction
More informationThe p-laplacian on Graphs with Applications in Image Processing and Classification
The p-laplacian on Graphs with Applications in Image Processing and Classification Abderrahim Elmoataz 1,2, Matthieu Toutain 1, Daniel Tenbrinck 3 1 UMR6072 GREYC, Université de Caen Normandie 2 Université
More informationRegularized Energy Minimization Models in Image Processing. Tibor Lukić. Faculty of Technical Sciences, University of Novi Sad, Serbia
20150715 NOVI SAD Regularized Energy Minimization Models in Image Processing Tibor Lukić Faculty of Technical Sciences, University of Novi Sad, Serbia SSIP Szeged, 2015 NOVI SAD NOVI SAD UNIVERSITY OF
More informationGuided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging
Guided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging Florin C. Ghesu 1, Thomas Köhler 1,2, Sven Haase 1, Joachim Hornegger 1,2 04.09.2014 1 Pattern
More informationEfficient MR Image Reconstruction for Compressed MR Imaging
Efficient MR Image Reconstruction for Compressed MR Imaging Junzhou Huang, Shaoting Zhang, and Dimitris Metaxas Division of Computer and Information Sciences, Rutgers University, NJ, USA 08854 Abstract.
More informationCombining Total Variation and Nonlocal Means Regularization for Edge Preserving Image Deconvolution
Electronic Letters on Computer Vision and Image Analysis 10(1):42-53, 2011 Combining Total Variation and Nonlocal Means Regularization for Edge Preserving Image Deconvolution Binbin Hao and Jianguang Zhu
More informationShape optimization under convexity constraint
Shape optimization under convexity constraint Jimmy LAMBOLEY Université Paris-Dauphine ANR GAOS Work with D. Bucur, I. Fragalà, E. Harrell, A. Henrot, M. Pierre, A. Novruzi 03/04/2012, PICOF J. Lamboley
More informationOptimal Denoising of Natural Images and their Multiscale Geometry and Density
Optimal Denoising of Natural Images and their Multiscale Geometry and Density Department of Computer Science and Applied Mathematics Weizmann Institute of Science, Israel. Joint work with Anat Levin (WIS),
More informationMultiscale TV flow with applications to fast denoising and registration
Multiscale TV flow with applications to fast denoising and registration Prashant Athavale a,robertxu a, Perry Radau b, Adrian Nachman a and Graham Wright b a University of Toronto, Toronto, Canada; b Sunnybrook
More informationSelf-Adaptive Blind Super-Resolution Image Reconstruction
2010 3rd International Congress on Image and Signal Processing (CISP2010) Self-Adaptive Blind Super-Resolution Image Reconstruction Yunfei Bai, Jing Hu, Yupin Luo Tsinghua National Laboratory for Information
More informationBlind Image Deblurring Using Dark Channel Prior
Blind Image Deblurring Using Dark Channel Prior Jinshan Pan 1,2,3, Deqing Sun 2,4, Hanspeter Pfister 2, and Ming-Hsuan Yang 3 1 Dalian University of Technology 2 Harvard University 3 UC Merced 4 NVIDIA
More informationNon-Blind Deblurring Using Partial Differential Equation Method
Volume Issue 3, 3-36, 013, ISSN: 319 8656 Non-Blind Deblurring Using Partial Differential Equation Method Devender Sharma CSE Department HCE,Sonepat, Puneet Sharma CSE Department HCE,Sonepat, Ritu Sharma
More informationWhite Paper. Denoising. February 2007 WP _v01
White Paper Denoising February 2007 WP-03020-001_v01 White Paper Document Change History Version Date Responsible Reason for Change _v01 AK, TS Initial release Go to sdkfeedback@nvidia.com to provide feedback
More informationOpen Access Microscopic Images Restoration of Chinese Herbal Medicine Using the TV Model
Send Orders for Reprints to reprints@benthamscience.ae 266 The Open Automation and Control Systems Journal, 2015, 7, 266-274 Open Access Microscopic Images Restoration of Chinese Herbal Medicine Using
More informationTotal Variation Denoising with Overlapping Group Sparsity
1 Total Variation Denoising with Overlapping Group Sparsity Ivan W. Selesnick and Po-Yu Chen Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 2 Abstract This paper describes
More informationLecture 19: Convex Non-Smooth Optimization. April 2, 2007
: Convex Non-Smooth Optimization April 2, 2007 Outline Lecture 19 Convex non-smooth problems Examples Subgradients and subdifferentials Subgradient properties Operations with subgradients and subdifferentials
More informationModeling Multigrid Algorithms for Variational Imaging
Modeling Multigrid Algorithms for Variational Imaging, Harald Koestler, Reinhard German, Ulrich Ruede Computer Networks and Communication Systems -Nürnberg, Germany April 09, 2010 What do we do? Model
More informationA primal-dual framework for mixtures of regularizers
A primal-dual framework for mixtures of regularizers Baran Gözcü baran.goezcue@epfl.ch Laboratory for Information and Inference Systems (LIONS) École Polytechnique Fédérale de Lausanne (EPFL) Switzerland
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 informationApplication of partial differential equations in image processing. Xiaoke Cui 1, a *
3rd International Conference on Education, Management and Computing Technology (ICEMCT 2016) Application of partial differential equations in image processing Xiaoke Cui 1, a * 1 Pingdingshan Industrial
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationEnhancement of Sharpness and Contrast Using Adaptive Parameters
International Journal of Computational Engineering Research Vol, 03 Issue, 10 Enhancement of Sharpness and Contrast Using Adaptive Parameters 1, Allabaksh Shaik, 2, Nandyala Ramanjulu, 1, Department of
More informationSurface Parameterization
Surface Parameterization A Tutorial and Survey Michael Floater and Kai Hormann Presented by Afra Zomorodian CS 468 10/19/5 1 Problem 1-1 mapping from domain to surface Original application: Texture mapping
More informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationCHAPTER 2 TEXTURE CLASSIFICATION METHODS GRAY LEVEL CO-OCCURRENCE MATRIX AND TEXTURE UNIT
CHAPTER 2 TEXTURE CLASSIFICATION METHODS GRAY LEVEL CO-OCCURRENCE MATRIX AND TEXTURE UNIT 2.1 BRIEF OUTLINE The classification of digital imagery is to extract useful thematic information which is one
More informationgruvi Edge Aware Anisotropic Diffusion for 3D Scalar Data graphics + usability + visualization Simon Fraser University Burnaby, BC Canada.
gruvi graphics + usability + visualization Edge Aware Anisotropic Diffusion for 3D Scalar Data Zahid Hossain Torsten Möller Simon Fraser University Burnaby, BC Canada. Motivation Original Sheep s heart
More informationsmooth coefficients H. Köstler, U. Rüde
A robust multigrid solver for the optical flow problem with non- smooth coefficients H. Köstler, U. Rüde Overview Optical Flow Problem Data term and various regularizers A Robust Multigrid Solver Galerkin
More informationA parameter study of a hybrid Laplacian mean-curvature flow denoising model
J Supercomput (2011) 57:339 356 DOI 10.1007/s11227-010-0417-z A parameter study of a hybrid Laplacian mean-curvature flow denoising model Ioana Banicescu Hyeona Lim Ricolindo L. Cariño Seongjai Kim Published
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 informationWhat is an edge? Paint. Depth discontinuity. Material change. Texture boundary
EDGES AND TEXTURES The slides are from several sources through James Hays (Brown); Srinivasa Narasimhan (CMU); Silvio Savarese (U. of Michigan); Bill Freeman and Antonio Torralba (MIT), including their
More informationSuppression of Missing Data Artifacts. for Deblurring Images Corrupted by
Suppression of Missing Data Artifacts for Deblurring Images Corrupted by Random Valued Noise Nam-Yong Lee Department of Applied Mathematics Institute of Basic Sciences, Inje University Gimhae, Gyeongnam
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
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 informationTOTAL VARIATION WAVELET BASED MEDICAL IMAGE DENOISING. 1. Introduction
TOTAL VARIATION WAVELET BASED MEDICAL IMAGE DENOISING YANG WANG AND HAOMIN ZHOU Abstract. We propose a denoising algorthm for medical images based on a combination of the total variation minimization scheme
More informationSmoothing & Fairing. Mario Botsch
Smoothing & Fairing Mario Botsch Motivation Filter out high frequency noise Desbrun, Meyer, Schroeder, Barr: Implicit Fairing of Irregular Meshes using Diffusion and Curvature Flow, SIGGRAPH 99 2 Motivation
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 informationLevel lines based disocclusion
Level lines based disocclusion Simon Masnou Jean-Michel Morel CEREMADE CMLA Université Paris-IX Dauphine Ecole Normale Supérieure de Cachan 75775 Paris Cedex 16, France 94235 Cachan Cedex, France Abstract
More informationFull-field optical methods for mechanical engineering: essential concepts to find one way
Full-field optical methods for mechanical engineering: essential concepts to find one way Yves Surrel Techlab September 2004 1 Contents 1 Introduction 3 2 White light methods 4 2.1 Random encoding............................................
More informationDr. Ulas Bagci
Lecture 9: Deformable Models and Segmentation CAP-Computer Vision Lecture 9-Deformable Models and Segmentation Dr. Ulas Bagci bagci@ucf.edu Lecture 9: Deformable Models and Segmentation Motivation A limitation
More informationChallenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies
Challenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies Rosemary Renaut, Hongbin Guo, Kewei Chen and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center
More informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
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 informationLecture 4: Nonlinear Filters
Lecture 4: Nonlinear Filters Math 49 Prof. Todd Wittman The Citadel Image Noise Physical effects (dust, raindrops) Electrical noise Caused by fluctuations in electric circuit. Thermal noise Impulse noise
More informationAutomatic Parameter Optimization for De-noising MR Data
Automatic Parameter Optimization for De-noising MR Data Joaquín Castellanos 1, Karl Rohr 2, Thomas Tolxdorff 3, and Gudrun Wagenknecht 1 1 Central Institute for Electronics, Research Center Jülich, Germany
More information