gruvi Edge Aware Anisotropic Diffusion for 3D Scalar Data graphics + usability + visualization Simon Fraser University Burnaby, BC Canada.
|
|
- Oliver Maxwell
- 5 years ago
- Views:
Transcription
1 gruvi graphics + usability + visualization Edge Aware Anisotropic Diffusion for 3D Scalar Data Zahid Hossain Torsten Möller Simon Fraser University Burnaby, BC Canada.
2 Motivation Original Sheep s heart data page 2 of 32
3 Motivation Existing gradient based anisotropic diffusion Original Sheep s heart data page 2 of 32
4 Motivation Existing gradient based anisotropic diffusion Gradient magnitude map page 2 of 32
5 Motivation Existing gradient based anisotropic diffusion Gradient magnitude map Proposed method page 2 of 32
6 Motivation Existing gradient based anisotropic diffusion Note Gradient magnitude map Proposed method Our method yields a consistent smoothing over an iso-surface without undesirable artifacts. page 2 of 32
7 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): = div(h( f ) f ) f t page 3 of 32
8 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): = div(h( f ) f ) f t Carmona et al.(tip, 1998) generalized the PM model in 2D: SF f t = {}}{ h ( ) af nn + bf vv, SF = Stopping Function where v is the orthonormal direction of the normal n. page 3 of 32
9 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): = div(h( f ) f ) f t Carmona et al.(tip, 1998) generalized the PM model in 2D: SF f t = {}}{ h ( ) af nn + bf vv, SF = Stopping Function where v is the orthonormal direction of the normal n. Extension in 3D by Gerig et al.(tmi, 1992): remains isotropic on the tangent plane of the gradient. page 3 of 32
10 Previous Work: Trouble in 3D Original page 4 of 32
11 Previous Work: Trouble in 3D PM model in 3D Original page 4 of 32
12 Previous Work: Trouble in 3D PM model in 3D Original Curvatures taken into account page 4 of 32
13 Previous Work: Trouble in 3D PM model in 3D Original Curvatures taken into account Therefore Weickert classified PM and higher order diffusion processes as non-linear rather than anisotropic. page 4 of 32
14 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). page 5 of 32
15 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. page 5 of 32
16 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. Level Set: The definition of edge is based on the curvature that is measured on the surface of the level set. page 5 of 32
17 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. Level Set: The definition of edge is based on the curvature that is measured on the surface of the level set. Our definition is based on the second derivative measured across the level sets. page 5 of 32
18 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. Level Set: The definition of edge is based on the curvature that is measured on the surface of the level set. Our definition is based on the second derivative measured across the level sets. De-noising: Bilateral filtering, mean-shift filtering and non-linear diffusion are equivalent: Barash et al., (IVC, 2004). page 5 of 32
19 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. Level Set: The definition of edge is based on the curvature that is measured on the surface of the level set. Our definition is based on the second derivative measured across the level sets. De-noising: Bilateral filtering, mean-shift filtering and non-linear diffusion are equivalent: Barash et al., (IVC, 2004). A recent methods namely SRNRAD and ORNRAD (TIP, 2009). page 5 of 32
20 Previous Work: In 3D Krissian et al.(scale-space, 1997) developed a true anistropic diffusion, (KM model), taking curvatures into account as proposed by Whitaker (VBC, 1994). The stopping function had gradient magnitude based parameter. Level Set: The definition of edge is based on the curvature that is measured on the surface of the level set. Our definition is based on the second derivative measured across the level sets. De-noising: Bilateral filtering, mean-shift filtering and non-linear diffusion are equivalent: Barash et al., (IVC, 2004). A recent methods namely SRNRAD and ORNRAD (TIP, 2009). We will compare our method with these two. page 5 of 32
21 Contributions 1 Removal of the gradient magnitude based parameter. page 6 of 32
22 Contributions 1 Removal of the gradient magnitude based parameter. 2 Dynamic adaptation of anisotropy based on local curvatures. page 6 of 32
23 Contributions 1 Removal of the gradient magnitude based parameter. 2 Dynamic adaptation of anisotropy based on local curvatures. 3 Application in de-noising and visualization. page 6 of 32
24 The Proposed Method page 7 of 32
25 Proposed Method: Edge The red solid curve is the error function while the dashed black curve is the second derivative of it. page 8 of 32
26 Proposed Method: Edge The red solid curve is the error function while the dashed black curve is the second derivative of it. A common technique to detect edges in 2D image processing. Basis of 2D transfer functions, Kindlmann and Durkin (VIS 1998). page 8 of 32
27 Four Objectives of Our Diffusion Objectives 1 No diffusion will be performed along the gradient direction. Note This one is different from the classical methods. This prevents blurring across an edge. page 9 of 32
28 Four Objectives of Our Diffusion Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. Note Check for f nn = 0. Note it does not create any problem in constant homogeneous regions. page 9 of 32
29 Four Objectives of Our Diffusion Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. Note Similar to Krissian et al.(1997), i.e the KM model. page 9 of 32
30 Four Objectives of Our Diffusion Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. 4 Diffuse isotropically on the tangent plane of the normal n in regions where the local iso-surface has similar principal curvatures. Note For example, on the surface of a sphere or a slab. page 9 of 32
31 Modeling the PDE Consider the 1D heat, also known as diffusion, equation: f t = cf xx page 10 of 32
32 Modeling the PDE Consider the 1D heat, also known as diffusion, equation: f t = cf xx A 3D extension of the above: f t = hf r 1 r 1 + gf r2 r 2 + wf nn Where the orthonormal bases [r 1,r 2,n] are min curvature, max curvature and normal directions respectively. page 10 of 32
33 Modeling the PDE Consider the 1D heat, also known as diffusion, equation: f t = cf xx A 3D extension of the above: f t = hf r 1 r 1 + gf r2 r 2 + wf nn Where the orthonormal bases [r 1,r 2,n] are min curvature, max curvature and normal directions respectively. Setting g = τh and w = ηh we can simplify the above f ( t = h f r1 r 1 + τf r2 r 2 + ηf nn ), τ,η [0,1] page 10 of 32
34 Modeling the PDE: Objective 1 Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. 4 Diffuse isotropically on the tangent plane of the normal n in regions where the local iso-surface has similar principal curvatures. f ( t = h f r1 r 1 + τf r2 r 2 + ηf nn ), τ,η [0,1] page 11 of 32
35 Modeling the PDE: Objective 1 Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. 4 Diffuse isotropically on the tangent plane of the normal n in regions where the local iso-surface has similar principal curvatures. Set η = 0 f ( t = h f r1 r 1 + τf r2 r 2 + ηf nn ), τ,η [0,1] f ( t = h f r1 r 1 + τf r2 r 2 ), τ [0,1] page 11 of 32
36 Modeling the PDE: Objective 2 Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. 4 Diffuse isotropically on the tangent plane of the normal n in regions where the local iso-surface has similar principal curvatures. page 12 of 32
37 Modeling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) α Plot of h(α) page 13 of 32
38 Modeling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) h can take, as argument, the directional second derivative f nn along the normal n α Plot of h(α) page 13 of 32
39 Modeling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) h can take, as argument, the directional second derivative f nn along the normal n. ) So far: f t = h(f nn) (f r1 r 1 + τf r2 r α Plot of h(α) page 13 of 32
40 Modeling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) h can take, as argument, the directional second derivative f nn along the normal n. ) So far: f t = h(f nn) (f r1 r 1 + τf r2 r α Plot of h(α) Note Diffusion stops when f nn 0, i.e. around the edge locations. Constant homogeneous regions, where f nn 0, would not be affected by diffusion anyway. page 13 of 32
41 Modeling the PDE: Objectives 3 and 4 Objectives 1 No diffusion will be performed along the gradient direction. 2 Diffusion will be stopped around the edge locations. 3 Diffusion will be performed anisotropically along the direction of the minimum curvature. 4 Diffuse isotropically on the tangent plane of the normal n in regions where the local iso-surface has similar principal curvatures. page 14 of 32
42 Modeling the PDE: Objectives 3 and 4 f ) t = h(f nn) (f r1 r 1 + τf r2 r 2 page 15 of 32
43 Modeling the PDE: Objectives 3 and 4 f ) t = h(f nn) (f r1 r 1 + τf r2 r 2 ( ) 2λ κ1 τ = κ 2 where κ 2 > 0,λ Z 1 κ 2 = 0 where, κ 1 : minimum curvature κ 2 : maximum curvature such that, κ 1 κ 2. page 15 of 32
44 Modeling the PDE: Objectives 3 and 4 where, such that, κ 1 κ 2. Note f ) t = h(f nn) (f r1 r 1 + τf r2 r 2 ( ) 2λ κ1 τ = κ 2 where κ 2 > 0,λ Z 1 κ 2 = 0 κ 1 : κ 2 : minimum curvature maximum curvature Both the curvatures are measured from a smoothed version of the data f ρ with a Gaussian filter having a small variance ρ 2. So we will use the notation τ ρ for the rest of the presentation to depict this and call it isotropy function. page 15 of 32
45 Simplification So far we have f ( IF t = h(f {}}{ nn) f r1 r 1 + τ ρ f r2 r 2 ), IF=Isotropy Function page 16 of 32
46 Simplification So far we have f ( IF t = h(f {}}{ nn) f r1 r 1 + τ ρ f r2 r 2 ), IF=Isotropy Function Which can be simplied to. f t = h(f nn) f (κ 1 + τ ρ κ 2 ) page 16 of 32
47 Simplification So far we have f ( IF t = h(f {}}{ nn) f r1 r 1 + τ ρ f r2 r 2 ), IF=Isotropy Function Which can be simplied to. f t = h(f nn) f (κ 1 + τ ρ κ 2 ) Note Has a similar form to Mean Curvature Motion (MCM), yet essentially different. f MCM: t = f (κ 1 + κ 2 ) page 16 of 32
48 Results: Significance of σ Original h( ) = 1 σ = 40 Our: ( ) fnn 2 f t = h(f nn) f (κ 1 + τ ρ κ 2 ), h(f nn ) = 1 (0.9) σ, σ R MCM: f t = f (κ 1 + κ 2 ) page 17 of 32
49 Results page 18 of 32
50 Results: Consistent Smoothing KM method: k = 40 Original Our method: σ = 1 ( ) fnn 2 f t = h(f nn) f (κ 1 + τ ρ κ 2 ), h(f nn ) = 1 (0.9) σ, σ R page 19 of 32
51 Results: Consistent Smoothing KM method: k = 60 Original Our method: σ = 10 ( ) fnn 2 f t = h(f nn) f (κ 1 + τ ρ κ 2 ), h(f nn ) = 1 (0.9) σ, σ R page 19 of 32
52 De-noising Properties page 20 of 32
53 De-noising Properties Note We used the same parameter settings for all our de-noising experiments across all datasets and noise types. With a unit grid spacing in all dimensions we empirically found that t 0.4 is stable for most practical purposes. page 20 of 32
54 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR=26.12 page 21 of 32
55 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile page 21 of 32
56 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile page 21 of 32
57 Results: De-noising (Salt and Pepper) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile page 22 of 32
58 De-noising Comparison page 23 of 32
59 De-noising Comparison Note We compared our anisotropic diffusion with two recent de-noising methods: SRNRAD and ORNRAD proposed by Krissian and Aja-Fernández (2009). We did not change any parameter settings. page 23 of 32
60 Quantitative Measures Mean Squared Error (MSE). Structural Similarity Index (SSIM). Quality Index Based on Local Variance (QILV). Original Blurred Gaussian noise SSIM: (Good) SSIM: (Bad) QILV: (Bad) QILV: (Good) page 24 of 32
61 De-noising: Comparison Gaussian noise Original MSE SSIM QILV Noise SRNRAD ORNRAD Our Noisy Our SRNRAD ORNRAD page 25 of 32
62 De-noising: Comparison Rician noise Original MSE SSIM QILV Noise SRNRAD ORNRAD Our Noisy Our SRNRAD ORNRAD page 25 of 32
63 De-noising: Comparison Note Our method performs similarly if not better in many cases. page 25 of 32
64 De-noising: Comparison Note Our method performs similarly if not better in many cases. And our Matlab implementation is 14 times faster than a multi-threaded C/C++ based implementation of ORNRAD. page 25 of 32
65 Impact on Visualization page 26 of 32
66 2D Transfer Function (With Noise) Gaussian noise Original Noisy Diffused page 27 of 32
67 2D Transfer Function (With Noise) Speckle noise Original Noisy Diffused page 27 of 32
68 2D Transfer Function (Without Noise) Tooth Dataset Original Diffused with our method page 28 of 32
69 2D Transfer Function Sheep s Heart Dataset KM method Original Our method page 29 of 32
70 Conclusion We developed an anistropic diffusion that: Smoothes consistently over an iso-surface Requires much less effort to choose a parameter value. Performs similarly, if not better, to a recent de-noising method. Is easy to implement and fast. Enhances 2D transfer functions. page 30 of 32
71 Acknowledgements We would like to thank: Dr. Karl Krissian and Dr. Santiago Aja-Fernández for providing supports in various occasions. Dr. Philippe Thévenaz for the Sphere phantom data. Natural Science and Engineering Research Council of Canada (NSERC) for funding this project. page 31 of 32
72 Thank you for your attention :) Contact Information Zahid Hossain Dr. Torsten Möller page 32 of 32
Edge Aware Anisotropic Diffusion for 3D Scalar Data
Edge Aware Anisotropic Diffusion for 3D Scalar Data Zahid Hossain and Torsten Möller, Member, IEEE (a) (b) (c) (d) (e) (f) Fig. 1: The left half of the figure demonstrates the consistency in smoothing
More informationFiltering and Edge Detection. Computer Vision I. CSE252A Lecture 10. Announcement
Filtering and Edge Detection CSE252A Lecture 10 Announcement HW1graded, will be released later today HW2 assigned, due Wed. Nov. 7 1 Image formation: Color Channel k " $ $ # $ I r I g I b % " ' $ ' = (
More informationLevel set based Volumetric Anisotropic Diffusion for 3D Image Filtering
Level set based Volumetric Anisotropic Diffusion for 3D Image Filtering Chandrajit L. Bajaj Department of Computer Science and Institute of Computational and Engineering Sciences University of Texas, Austin,
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 informationThe Pre-Image Problem in Kernel Methods
The Pre-Image Problem in Kernel Methods James Kwok Ivor Tsang Department of Computer Science Hong Kong University of Science and Technology Hong Kong The Pre-Image Problem in Kernel Methods ICML-2003 1
More informationFast 3D Mean Shift Filter for CT Images
Fast 3D Mean Shift Filter for CT Images Gustavo Fernández Domínguez, Horst Bischof, and Reinhard Beichel Institute for Computer Graphics and Vision, Graz University of Technology Inffeldgasse 16/2, A-8010,
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 informationQuality Guided Image Denoising for Low-Cost Fundus Imaging
Quality Guided Image Denoising for Low-Cost Fundus Imaging Thomas Köhler1,2, Joachim Hornegger1,2, Markus Mayer1,2, Georg Michelson2,3 20.03.2012 1 Pattern Recognition Lab, Ophthalmic Imaging Group 2 Erlangen
More informationNumerical Methods on the Image Processing Problems
Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006 Objective Develop efficient PDE (partial differential equations)
More informationArtistic Stylization of Images and Video Part III Anisotropy and Filtering Eurographics 2011
Artistic Stylization of Images and Video Part III Anisotropy and Filtering Eurographics 2011 Hasso-Plattner-Institut, University of Potsdam, Germany Image/Video Abstraction Stylized Augmented Reality for
More informationGaussian and Mean Curvature Planar points: Zero Gaussian curvature and zero mean curvature Tangent plane intersects surface at infinity points Gauss C
Outline Shape Analysis Basics COS 526, Fall 21 Curvature PCA Distance Features Some slides from Rusinkiewicz Curvature Curvature Curvature κof a curve is reciprocal of radius of circle that best approximates
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 informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationImage Processing with Nonparametric Neighborhood Statistics
Image Processing with Nonparametric Neighborhood Statistics Ross T. Whitaker Scientific Computing and Imaging Institute School of Computing University of Utah PhD Suyash P. Awate University of Pennsylvania,
More informationNOISE in magnitude magnetic resonance (MR) images is
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 17, NO. 8, AUGUST 2008 1383 Noise and Signal Estimation in Magnitude MRI and Rician Distributed Images: A LMMSE Approach Santiago Aja-Fernández, Carlos Alberola-López,
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 informationMEDICAL IMAGE NOISE REDUCTION AND REGION CONTRAST ENHANCEMENT USING PARTIAL DIFFERENTIAL EQUATIONS
MEDICAL IMAGE NOISE REDUCTION AND REGION CONTRAST ENHANCEMENT USING PARTIAL DIFFERENTIAL EQUATIONS Miguel Alemán-Flores, Luis Álvarez-León Departamento de Informática y Sistemas, Universidad de Las Palmas
More informationA Survey on Magnetic Resonance Image Denoising Methods
A Survey on Magnetic Resonance Image Denoising Methods Snehal More 1, V.V.Hanchate 2 1M.E. Student, PES s MCOE Pune, Maharashtra, India 2Dept. of E&TC Engineering, PES s MCOE Pune, Maharashtra, India ---------------------------------------------------------------------***---------------------------------------------------------------------
More informationGeometry From Imaging Data
Zhang et al. 2004 Geometry From Imaging Data Part II: Geometry Reconstruction Computational Biomechanics Summer Term 2016 Lecture 5/12 Frank Niemeyer Geometry Reconstruction Starting Point Note: intentionally
More informationComputer Vision I. Announcement. Corners. Edges. Numerical Derivatives f(x) Edge and Corner Detection. CSE252A Lecture 11
Announcement Edge and Corner Detection Slides are posted HW due Friday CSE5A Lecture 11 Edges Corners Edge is Where Change Occurs: 1-D Change is measured by derivative in 1D Numerical Derivatives f(x)
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 informationAnnouncements. Edge Detection. An Isotropic Gaussian. Filters are templates. Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3.
Announcements Edge Detection Introduction to Computer Vision CSE 152 Lecture 9 Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3. Reading from textbook An Isotropic Gaussian The picture
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 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 informationImage Enhancement and De-noising by Diffusion based Singular Value Decomposition
Image Enhancement and De-noising by Diffusion based Singular Value Decomposition Nafis uddin Khan ABV-Indian Institute of Information Technology and Management, Gwalior-474010, India K. V. Arya ABV-Indian
More information3D Shape Recovery of Smooth Surfaces: Dropping the Fixed Viewpoint Assumption
IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, VOL., NO., 1 3D Shape Recovery of Smooth Surfaces: Dropping the Fixed Viewpoint Assumption Yael Moses Member, IEEE and Ilan Shimshoni Member,
More informationEvaluation of Diffusion Techniques for Improved Vessel Visualization and Quantification in Three-Dimensional Rotational Angiography
Evaluation of Diffusion Techniques for Improved Vessel Visualization and Quantification in Three-Dimensional Rotational Angiography Erik Meijering 1, Wiro Niessen 1, Joachim Weickert 2, and Max Viergever
More information(Discrete) Differential Geometry
(Discrete) Differential Geometry Motivation Understand the structure of the surface Properties: smoothness, curviness, important directions How to modify the surface to change these properties What properties
More informationIntroduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (3) April 27, 2017 Kenshi Takayama Solid modeling 2 Solid models Thin shapes represented by single polygons Unorientable Clear definition of inside & outside
More informationEdge-Preserving Denoising for Segmentation in CT-Images
Edge-Preserving Denoising for Segmentation in CT-Images Eva Eibenberger, Anja Borsdorf, Andreas Wimmer, Joachim Hornegger Lehrstuhl für Mustererkennung, Friedrich-Alexander-Universität Erlangen-Nürnberg
More information8. Tensor Field Visualization
8. Tensor Field Visualization Tensor: extension of concept of scalar and vector Tensor data for a tensor of level k is given by t i1,i2,,ik (x 1,,x n ) Second-order tensor often represented by matrix Examples:
More informationFeature Preserving Smoothing of 3D Surface Scans. Thouis Raymond Jones
Feature Preserving Smoothing of 3D Surface Scans by Thouis Raymond Jones Submitted to the Department of Electrical Engineering and Computer Science in partial fulfillment of the requirements for the degree
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 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 informationHow does bilateral filter relates with other methods?
A Gentle Introduction to Bilateral Filtering and its Applications How does bilateral filter relates with other methods? Fredo Durand (MIT CSAIL) Slides by Pierre Kornprobst (INRIA) 0:35 Many people worked
More informationEdge detection. Gradient-based edge operators
Edge detection Gradient-based edge operators Prewitt Sobel Roberts Laplacian zero-crossings Canny edge detector Hough transform for detection of straight lines Circle Hough Transform Digital Image Processing:
More informationMultiscale 3D Feature Extraction and Matching
Multiscale 3D Feature Extraction and Matching Hadi Fadaifard and George Wolberg Graduate Center, City University of New York May 18, 2011 Hadi Fadaifard and George Wolberg Multiscale 3D Feature Extraction
More informationEdge Patch Based Image Denoising Using Modified NLM Approach
Edge Patch Based Image Denoising Using Modified NLM Approach Rahul Kumar Dongardive 1, Ritu Shukla 2, Dr. Y K Jain 3 1 Research Scholar of SATI Vidisha, M.P., India 2 Asst. Prof., Computer Science and
More informationNormalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data
Normalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data Tuan Q. Pham, Lucas J. van Vliet Pattern Recognition Group,
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 information3. The three points (2, 4, 1), (1, 2, 2) and (5, 2, 2) determine a plane. Which of the following points is in that plane?
Math 4 Practice Problems for Midterm. A unit vector that is perpendicular to both V =, 3, and W = 4,, is (a) V W (b) V W (c) 5 6 V W (d) 3 6 V W (e) 7 6 V W. In three dimensions, the graph of the equation
More informationSPECKLE REDUCTION IN ECHOCARDIOGRAPHIC IMAGES
4th European Signal Processing Conference (EUSIPCO 006), Florence, Ital September 4-8, 006, copyright by EURASIP SPECKLE REDUCTION IN ECHOCARDIOGRAPHIC IMAGES Nadia SOUAG Images Processing Laborator Faculty
More informationA non-local algorithm for image denoising
A non-local algorithm for image denoising Antoni Buades, Bartomeu Coll Dpt. Matemàtiques i Informàtica, UIB Ctra. Valldemossa Km. 7.5, 07122 Palma de Mallorca, Spain vdmiabc4@uib.es, tomeu.coll@uib.es
More informationA MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space
A MAP Algorithm for AVO Seismic Inversion Based on the Mixed (L 2, non-l 2 ) Norms to Separate Primary and Multiple Signals in Slowness Space Harold Ivan Angulo Bustos Rio Grande do Norte State University
More informationCS-465 Computer Vision
CS-465 Computer Vision Nazar Khan PUCIT 9. Optic Flow Optic Flow Nazar Khan Computer Vision 2 / 25 Optic Flow Nazar Khan Computer Vision 3 / 25 Optic Flow Where does pixel (x, y) in frame z move to in
More informationSEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH
SEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH Ignazio Gallo, Elisabetta Binaghi and Mario Raspanti Universitá degli Studi dell Insubria Varese, Italy email: ignazio.gallo@uninsubria.it ABSTRACT
More informationAn Efficient, Geometric Multigrid Solver for the Anisotropic Diffusion Equation in Two and Three Dimensions
1 n Efficient, Geometric Multigrid Solver for the nisotropic Diffusion Equation in Two and Three Dimensions Tolga Tasdizen, Ross Whitaker UUSCI-2004-002 Scientific Computing and Imaging Institute University
More informationPatch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques
Patch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques Syed Gilani Pasha Assistant Professor, Dept. of ECE, School of Engineering, Central University of Karnataka, Gulbarga,
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 informationThe SIFT (Scale Invariant Feature
The SIFT (Scale Invariant Feature Transform) Detector and Descriptor developed by David Lowe University of British Columbia Initial paper ICCV 1999 Newer journal paper IJCV 2004 Review: Matt Brown s Canonical
More informationImproved Non-Local Means Algorithm Based on Dimensionality Reduction
Improved Non-Local Means Algorithm Based on Dimensionality Reduction Golam M. Maruf and Mahmoud R. El-Sakka (&) Department of Computer Science, University of Western Ontario, London, Ontario, Canada {gmaruf,melsakka}@uwo.ca
More informationEfficient Iterative Semi-supervised Classification on Manifold
. Efficient Iterative Semi-supervised Classification on Manifold... M. Farajtabar, H. R. Rabiee, A. Shaban, A. Soltani-Farani Sharif University of Technology, Tehran, Iran. Presented by Pooria Joulani
More informationLocal Image preprocessing (cont d)
Local Image preprocessing (cont d) 1 Outline - Edge detectors - Corner detectors - Reading: textbook 5.3.1-5.3.5 and 5.3.10 2 What are edges? Edges correspond to relevant features in the image. An edge
More informationComparative Analysis in Medical Imaging
1 International Journal of Computer Applications (975 8887) Comparative Analysis in Medical Imaging Ashish Verma DCS, Punjabi University 1, Patiala, India Bharti Sharma DCS, Punjabi University 1, Patiala,
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 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 informationBayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis
Bayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis Xavier Le Faucheur a, Brani Vidakovic b and Allen Tannenbaum a a School of Electrical and Computer Engineering, b Department of Biomedical
More informationREVIEW PAPER ON IMAGE EDGE DETECTION ALGORITHMS FOR SEGMENTATION
REVIEW PAPER ON IMAGE EDGE DETECTION ALGORITHMS FOR SEGMENTATION Parvita Taya Department of CSE, AIMT, Karnal, Haryana, India Email- parvitataya@yahoo.co.in Abstract Computer vision is the rapid expanding
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationTHE preceding chapters were all devoted to the analysis of images and signals which
Chapter 5 Segmentation of Color, Texture, and Orientation Images THE preceding chapters were all devoted to the analysis of images and signals which take values in IR. It is often necessary, however, to
More informationProjection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT
Projection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT Benedikt Lorch 1, Martin Berger 1,2, Joachim Hornegger 1,2, Andreas Maier 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
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 informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationGeneration of Triangle Meshes from Time-of-Flight Data for Surface Registration
Generation of Triangle Meshes from Time-of-Flight Data for Surface Registration Thomas Kilgus, Thiago R. dos Santos, Alexander Seitel, Kwong Yung, Alfred M. Franz, Anja Groch, Ivo Wolf, Hans-Peter Meinzer,
More informationLocally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling
Locally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling Moritz Baecher May 15, 29 1 Introduction Edge-preserving smoothing and super-resolution are classic and important
More informationON THE GEODESIC TORSION OF A TANGENTIAL INTERSECTION CURVE OF TWO SURFACES IN R Introduction
ON THE GEODESIC TORSION OF A TANGENTIAL INTERSECTION CURVE OF TWO SURFACES IN R 3 B. UYAR DÜLDÜL and M. ÇALIŞKAN Abstract. In this paper, we find the unit tangent vector and the geodesic torsion of the
More informationA Simple Algorithm for Image Denoising Based on MS Segmentation
A Simple Algorithm for Image Denoising Based on MS Segmentation G.Vijaya 1 Dr.V.Vasudevan 2 Senior Lecturer, Dept. of CSE, Kalasalingam University, Krishnankoil, Tamilnadu, India. Senior Prof. & Head Dept.
More informationAdaptive Multiple-Frame Image Super- Resolution Based on U-Curve
Adaptive Multiple-Frame Image Super- Resolution Based on U-Curve IEEE Transaction on Image Processing, Vol. 19, No. 12, 2010 Qiangqiang Yuan, Liangpei Zhang, Huanfeng Shen, and Pingxiang Li Presented by
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 informationON THE GEODESIC TORSION OF A TANGENTIAL INTERSECTION CURVE OF TWO SURFACES IN R Introduction
Acta Math. Univ. Comenianae Vol. LXXXII, (3), pp. 77 89 77 ON THE GEODESIC TORSION OF A TANGENTIAL INTERSECTION CURVE OF TWO SURFACES IN R 3 B. UYAR DÜLDÜL and M. ÇALIŞKAN Abstract. In this paper, we find
More informationLaser sensors. Transmitter. Receiver. Basilio Bona ROBOTICA 03CFIOR
Mobile & Service Robotics Sensors for Robotics 3 Laser sensors Rays are transmitted and received coaxially The target is illuminated by collimated rays The receiver measures the time of flight (back and
More informationSkeleton Extraction via Anisotropic Heat Flow
DIREKOGLU ET AL.: SKELETON EXTRACTION VIA ANISOTROPIC HEAT FLOW 1 Skeleton Extraction via Anisotropic Heat Flow Cem Direkoglu cem.direkoglu@scss.tcd.ie Rozenn Dahyot rozenn.dahyot@scss.tcd.ie Michael Manzke
More informationPart 3: Image Processing
Part 3: Image Processing Moving Window Transform Georgy Gimel farb COMPSCI 373 Computer Graphics and Image Processing 1 / 62 1 Examples of linear / non-linear filtering 2 Moving window transform 3 Gaussian
More informationEE368 Project Report CD Cover Recognition Using Modified SIFT Algorithm
EE368 Project Report CD Cover Recognition Using Modified SIFT Algorithm Group 1: Mina A. Makar Stanford University mamakar@stanford.edu Abstract In this report, we investigate the application of the Scale-Invariant
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 informationChap.12 Kernel methods [Book, Chap.7]
Chap.12 Kernel methods [Book, Chap.7] Neural network methods became popular in the mid to late 1980s, but by the mid to late 1990s, kernel methods have also become popular in machine learning. The first
More informationFinal Review CMSC 733 Fall 2014
Final Review CMSC 733 Fall 2014 We have covered a lot of material in this course. One way to organize this material is around a set of key equations and algorithms. You should be familiar with all of these,
More informationLocal adaptivity to variable smoothness for exemplarbased image regularization and representation
Local adaptivity to variable smoothness for exemplarbased image regularization and representation Charles KERVRANN & Jérôme BOULANGER IRISA / INRIA Rennes, Campus Universitaire de Beaulieu, F-35042 Rennes
More informationCORRELATION VISUALIZATION FOR STRUCTURAL UNCERTAINTY ANALYSIS
International Journal for Uncertainty Quantification, 1(1):xxx xxx, 2012 CORRELATION VISUALIZATION FOR STRUCTURAL UNCERTAINTY ANALYSIS Tobias Pfaffelmoser & Rüdiger Westermann Technische Universität München,
More informationAnnée 2005 numéro 1. David Tschumperlé. Fast Anisotropic Smoothing of Multi-Valued Images. using Curvature-Preserving PDE s
Les cahiers du GREYC Année 2005 numéro 1 David Tschumperlé Fast Anisotropic Smoothing of Multi-Valued Images using Curvature-Preserving PDE s (revised version, May 2006 Groupe de Recherche en Informatique,
More informationMultiscale Properties of Weighted Total Variation Flow with Applications to Denoising and Registration
Multiscale Properties of Weighted Total Variation Flow with Applications to Denoising and Registration Prashant Athavale a,b,, Robert Xu c,d, Perry Radau c,1, Adrian Nachman a,e,f, Graham A. Wright c,d
More informationResearch Article Adaptive Image Enhancement Algorithm Combining Kernel Regression and Local Homogeneity
Mathematical Problems in Engineering Volume 010, Article ID 69353, 14 pages doi:10.1155/010/69353 Research Article Adaptive Image Enhancement Algorithm Combining Kernel Regression and Local Homogeneity
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 informationVehicular shape-based objects classification using Fourier descriptor technique
Journal of Scientific & Industrial Research 484 Vol. 68, June 2009, pp. 484-495 J SCI IND RES VOL 68 JUNE 2009 Vehicular shape-based objects classification using Fourier descriptor technique Raj Bahadur
More information3D Object Recognition using Multiclass SVM-KNN
3D Object Recognition using Multiclass SVM-KNN R. Muralidharan, C. Chandradekar April 29, 2014 Presented by: Tasadduk Chowdhury Problem We address the problem of recognizing 3D objects based on various
More informationRobotics Programming Laboratory
Chair of Software Engineering Robotics Programming Laboratory Bertrand Meyer Jiwon Shin Lecture 8: Robot Perception Perception http://pascallin.ecs.soton.ac.uk/challenges/voc/databases.html#caltech car
More informationCorner Detection. Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology
Corner Detection Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu April 11, 2006 Abstract Corners and edges are two of the most important geometrical
More informationCS 523: Computer Graphics, Spring Differential Geometry of Surfaces
CS 523: Computer Graphics, Spring 2009 Shape Modeling Differential Geometry of Surfaces Andrew Nealen, Rutgers, 2009 3/4/2009 Recap Differential Geometry of Curves Andrew Nealen, Rutgers, 2009 3/4/2009
More informationDenoising the Spectral Information of Non Stationary Image using DWT
Denoising the Spectral Information of Non Stationary Image using DWT Dr.DolaSanjayS 1, P. Geetha Lavanya 2, P.Jagapathi Raju 3, M.Sai Kishore 4, T.N.V.Krishna Priya 5 1 Principal, Ramachandra College of
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 informationEDGE-AWARE IMAGE PROCESSING WITH A LAPLACIAN PYRAMID BY USING CASCADE PIECEWISE LINEAR PROCESSING
EDGE-AWARE IMAGE PROCESSING WITH A LAPLACIAN PYRAMID BY USING CASCADE PIECEWISE LINEAR PROCESSING 1 Chien-Ming Lu ( 呂建明 ), 1 Sheng-Jie Yang ( 楊勝傑 ), 1 Chiou-Shann Fuh ( 傅楸善 ) Graduate Institute of Computer
More informationENGI Parametric & Polar Curves Page 2-01
ENGI 3425 2. Parametric & Polar Curves Page 2-01 2. Parametric and Polar Curves Contents: 2.1 Parametric Vector Functions 2.2 Parametric Curve Sketching 2.3 Polar Coordinates r f 2.4 Polar Curve Sketching
More information1 Introduction The theory of linear scale-space [6, 12] enables the detection and localization of edges while eliminating noise by tracking features a
Local Scale Controlled Anisotropic Diusion with Local Noise Estimate for Image Smoothing and Edge Detection Ping Liang Y. F. Wang College of Engineering Department of Computer Science University of California
More informationGEOMETRY AND GRAPHICS EXAMPLES AND EXERCISES. Preface
GEOMETRY AND GRAPHICS EXAMPLES AND EXERCISES Preface This textbook is intended for students of subjects Constructive Geometry and Computer Graphics at the Faculty of Mechanical Engineering at Czech Technical
More informationGeneralized Scale-Based Image Filtering
Generalized Scale-Based Image Filtering Andre Souza a, Jayaram K. Udupa a, and Anant Madabhushi a a Medical Image Processing Group, Department of Radiology University of Pennsylvania, Philadelphia, PA
More informationFiltering, reconstruction and registration of 3D Ultrasound Images
Filtering, reconstruction and registration of 3D Ultrasound Images Hervé Delingette Projet Epidaure Herve.Delingette Delingette@inria.fr Acknowledgement Johan Montagnat Alexis Roche Karl Krissian Jérémie
More information1 Scalar Transport and Diffusion
ME 543 Scalar Transport and Diffusion February 8 Scalar Transport and Diffusion As we will see, scalar transport is most naturally addressed in Lagrangian form, although it is often used in Eulerian form.
More informationBM3D Image Denoising using Learning-based Adaptive Hard Thresholding
BM3D Image Denoising using Learning-based Adaptive Hard Thresholding Farhan Bashar and Mahmoud R. El-Sakka Computer Science Department, The University of Western Ontario, London, Ontario, Canada Keywords:
More informationLOCAL-GLOBAL OPTICAL FLOW FOR IMAGE REGISTRATION
LOCAL-GLOBAL OPTICAL FLOW FOR IMAGE REGISTRATION Ammar Zayouna Richard Comley Daming Shi Middlesex University School of Engineering and Information Sciences Middlesex University, London NW4 4BT, UK A.Zayouna@mdx.ac.uk
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 information