graphics + usability + visualization Edge Aware Anisotropic Diffusion for 3D Scalar Data on Regular Lattices Zahid Hossain MSc.
|
|
- Camilla Shields
- 5 years ago
- Views:
Transcription
1 gruvi graphics + usability + visualization Edge Aware Anisotropic Diffusion for 3D Scalar Data on Regular Lattices Zahid Hossain MSc. Thesis Defense Simon Fraser University Burnaby, BC Canada.
2 Motivation Original Sheep s heart data 2 of 43
3 Motivation Existing gradient based anisotropic diffusion Original Sheep s heart data 2 of 43
4 Motivation Existing gradient based anisotropic diffusion Gradient magnitude map 2 of 43
5 Motivation Existing gradient based anisotropic diffusion Gradient magnitude map Proposed method Our method is consistent and free of undesirable artifacts 2 of 43
6 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): f = div(h( f ) f ) t 3 of 43
7 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): f = div(h( f ) f ) t Carmona et al.(tip, 1998) generalization in 2D: SF f t = {}}{ h ( ) af nn + bf vv, SF = Stopping Function where v = n. 3 of 43
8 Previous Work Perona and Malik (PM) Model (TPAMI, 1990): f = div(h( f ) f ) t Carmona et al.(tip, 1998) generalization in 2D: SF f t = {}}{ h ( ) af nn + bf vv, SF = Stopping Function where v = n. 3D extension by Gerig et al.(tmi, 1992): Isotropic on the tangent plane of the gradient. 3 of 43
9 Previous Work: Trouble in 3D Original 4 of 43
10 Previous Work: Trouble in 3D PM model in 3D Original 4 of 43
11 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. 4 of 43
12 Previous Work: In 3D Whitaker (VBC, 1994) and Krissian et al.(scale-space, 1997) (KM model), took curvatures into account The stopping function had gradient magnitude based parameter. 5 of 43
13 Previous Work: In 3D Whitaker (VBC, 1994) and Krissian et al.(scale-space, 1997) (KM model), took curvatures into account 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. 5 of 43
14 Previous Work: In 3D Whitaker (VBC, 1994) and Krissian et al.(scale-space, 1997) (KM model), took curvatures into account 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. De-noising: Bilateral filtering, mean-shift filtering and non-linear diffusion are equivalent: Barash et al., (IVC, 2004). A recent methods: SRNRAD and ORNRAD (TIP, 2009). 5 of 43
15 Contributions 1 Removal of the gradient magnitude based parameter 6 of 43
16 Contributions 1 Removal of the gradient magnitude based parameter 2 Dynamic adaptation of anisotropy 6 of 43
17 Contributions 1 Removal of the gradient magnitude based parameter 2 Dynamic adaptation of anisotropy 3 Application in de-noising 6 of 43
18 Contributions 1 Removal of the gradient magnitude based parameter 2 Dynamic adaptation of anisotropy 3 Application in de-noising 4 Implementation recipe in arbitrary regular lattices 6 of 43
19 The Proposed Method 7 of 43
20 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 in 2D image processing. 8 of 43
21 Four Objectives of Our Diffusion Objectives 1 No diffusion along the gradient direction. This prevents blurring across an edge. 9 of 43
22 Four Objectives of Our Diffusion Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. Check for f nn = 0. No problem in constant homogeneous regions. 9 of 43
23 Four Objectives of Our Diffusion Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. Similar to Krissian et al.(1997), i.e the KM model. 9 of 43
24 Four Objectives of Our Diffusion Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. 4 Diffuse on the tangent plane where the local iso-surface has similar principal curvatures. For example, on the surface of a sphere or a slab. 9 of 43
25 Modelling the PDE Consider the 1D heat, also known as diffusion, equation: f t = cf xx 10 of 43
26 Modelling 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. 10 of 43
27 Modelling 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] 10 of 43
28 Modelling the PDE: Objective 1 Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. 4 Diffuse on the tangent plane where the local iso-surface has similar principal curvatures. f ( t = h f r1 r 1 + τf r2 r 2 + ηf nn ), τ,η [0,1] 11 of 43
29 Modelling the PDE: Objective 1 Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. 4 Diffuse on the tangent plane 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] 11 of 43
30 Modelling the PDE: Objective 2 Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. 4 Diffuse on the tangent plane where the local iso-surface has similar principal curvatures. 12 of 43
31 Modelling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) α Plot of h(α) 13 of 43
32 Modelling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) Takes the directional second derivative f nn along the normal n as argument α Plot of h(α) 13 of 43
33 Modelling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) Takes the directional second derivative f nn along the normal n as argument. ) So far: f t = h(f nn) (f r1 r 1 + τf r2 r α Plot of h(α) 13 of 43
34 Modelling the PDE: Objective 2 Define h h(α) = 1 (0.9) ( α σ ) 2, σ R h(α) Takes the directional second derivative f nn along the normal n as argument. ) So far: f t = h(f nn) (f r1 r 1 + τf r2 r α Plot of h(α) Diffusion stops when f nn 0, i.e. around the edge locations. Homogeneous regions remains unaffected. 13 of 43
35 Modelling the PDE: Objectives 3 and 4 Objectives 1 No diffusion along the gradient direction. 2 Diffusion stopped around the edge locations. 3 Diffuse along the direction of the minimum curvature. 4 Diffuse on the tangent plane where the local iso-surface has similar principal curvatures. 14 of 43
36 Modelling the PDE: Objectives 3 and 4 f ) t = h(f nn) (f r1 r 1 + τf r2 r 2 15 of 43
37 Modelling the PDE: Objectives 3 and 4 where, 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 15 of 43
38 Modelling the PDE: Objectives 3 and 4 where, 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 Usually λ = 2 works well Curvatures are measured from a blurred version of the data f ρ τ is replaced by τ ρ 15 of 43
39 Simplification So far we have Stopping Function f {}}{( ) t = h(f nn ) f r1 r 1 + τ ρ f r2 r }{{} 2 Isotropy Function 16 of 43
40 Simplification So far we have Which can be simplified to Stopping Function f {}}{( ) t = h(f nn ) f r1 r 1 + τ ρ f r2 r }{{} 2 Isotropy Function. f t = h(f nn) f (κ 1 + τ ρ κ 2 ) 16 of 43
41 Simplification So far we have Which can be simplified to Stopping Function f {}}{( ) t = h(f nn ) f r1 r 1 + τ ρ f r2 r }{{} 2 Isotropy Function. f t = h(f nn) f (κ 1 + τ ρ κ 2 ) Has a similar form to Mean Curvature Motion (MCM), yet essentially different. f MCM: t = f (κ 1 + κ 2 ) 16 of 43
42 Significance of the Stopping Function Original MCM Our method Our: ( ) fnn 2 f t = h(f nn) f (κ 1 + τ ρ κ 2 ), h(f nn ) = 1 (0.9) σ, σ R MCM: f t = f (κ 1 + κ 2 ) 17 of 43
43 Stability Courant-Friedrichs-Lewy (CFL) revealed the following necessary condition for stability: t C d2, d = min{ x, y, z}, C = Courant Number 2 C 0.8 is stable for most practical purposes. 18 of 43
44 Results 19 of 43
45 Results: Consistent Smoothing KM method: k = 40 Original Our method: σ = 1 20 of 43
46 Results: Consistent Smoothing KM method: k = 60 Original Our method: σ = of 43
47 De-noising Properties Note We used the same parameter settings for all our de-noising experiments across all datasets and noise types. 21 of 43
48 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR= of 43
49 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile 22 of 43
50 Results: De-noising (Gaussian) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile 22 of 43
51 Results: De-noising (Salt and Pepper) Noisy, SNR=12.89 Original Diffused, SNR=26.12 Original Noisy Diffused Profile 23 of 43
52 De-noising Comparison Note Comparison with SRNRAD and ORNRAD (Krissian and Aja-Fernández, 2009). We did not change any parameter settings. 24 of 43
53 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) 25 of 43
54 De-noising: Comparison Gaussian noise Original MSE SSIM QILV Noise SRNRAD ORNRAD Our Noisy Our SRNRAD ORNRAD 26 of 43
55 De-noising: Comparison Rician noise Original MSE SSIM QILV Noise SRNRAD ORNRAD Our Noisy Our SRNRAD ORNRAD 26 of 43
56 De-noising: Remarks Our method performs similarly if not better in many cases. 27 of 43
57 De-noising: Remarks Our method performs similarly if not better in many cases. Matlab implementation is 14 faster than a multi-threaded C/C++ based ORNRAD. Our formulation is simpler ORNRAD requires, beside complex statistical measurements, eigen decomposition of structure tensor. 27 of 43
58 Extension Framework for Regular Lattices 28 of 43
59 Lattices in 3D Cartesian Cubic (CC) Body Centric Cubic (BCC) Face Centric Cubic (FCC) A regular lattice can be characterized by a matrix L 29 of 43
60 Motivation For Other Lattices Tai Meng et al. (2011) and Entezari et al. (2008) CC BCC 30 of 43
61 Motivation For Other Lattices (Contd.) Alim et al. (2009) CC BCC 31 of 43
62 Taylor Expansion of Discrete Convolution Discrete convolution of f [k] with a filter [k] can be written as the following: f = fr w [k] = a }{{} n D n f [k] f [k] = f (Lk) n Taylor Coefficient 32 of 43
63 Taylor Expansion of Discrete Convolution Discrete convolution of f [k] with a filter [k] can be written as the following: f = fr w [k] = a }{{} n D n f [k] f [k] = f (Lk) n Taylor Coefficient And a n is a linear function of the filter weights [k] 32 of 43
64 Linear System f = fr w [k] = a n D n f [k], n f [k] = f (Lk) To extract the u-th derivative D u f [k], correct up to a polynomial order of n, the following must be satisfied { a 0, n η [0,n] n = d and n u 1, n = u. 33 of 43
65 Linear System (Contd.) Forms a linear system with the unknown filter weights w i = [k i ] a n 1 a n 2. a n m = of 43
66 Linear System (Contd.) Forms a linear system with the unknown filter weights w i = [k i ] a n 1 a n 2. a n m = g 1,1. g 1, g 1,m. g n,1 g n,2... g n,m w 1 w 2. w m = G W = C 34 of 43
67 Infinite Solution Space The linear system is often not full-rank 35 of 43
68 Infinite Solution Space The linear system is often not full-rank Finding a suitable solution by: Imposing symmetry/anti-symmetry in the filter geometry Minimizing error in the higher polynomial order 35 of 43
69 Infinite Solution Space The linear system is often not full-rank Finding a suitable solution by: Imposing symmetry/anti-symmetry in the filter geometry Minimizing error in the higher polynomial order Usually yields a filter that has small support 35 of 43
70 Results 36 of 43
71 Quality of Normals: Synthetic Data 2-EF 4-EF CC (Mag=3.13, Angle=49.4 ) (2.79, 44.1 ) (2.42, 28.7 ) BCC (2.88, Hossain, Simon Fraser University, Canada ) gruvi graphics + usability + visualization 37 of 43
72 Quality of Normals: Real Data (Iso-Surface) CC (4-EF) BCC (4-EF) 38 of 43
73 Effects of Higher Order Filters (CC) Original 2-EF 4-EF 39 of 43
74 Conclusion We developed... An anisotropic diffusion: Consistent smoothing Less effort to choose a parameter Performs similarly, if not better, to a recent de-noising method Fast and easy to implement 40 of 43
75 Conclusion We developed... An anisotropic diffusion: Consistent smoothing Less effort to choose a parameter Performs similarly, if not better, to a recent de-noising method Fast and easy to implement A framework to design discrete derivative filters: Works on arbitrary regular lattice Any derivative could be estimated Size of the filter could be specified PDE implementation on any regular lattice is made possible 40 of 43
76 Future Work Actual implementation in BCC and FCC De-noising performance in BCC and FCC Different functions for τ Alternate edge detectors 41 of 43
77 Acknowledgements I would like to thank: My supervisor, co-supervisor and the examining committee All the GrUVi members for their relentless support Natural Science and Engineering Research Council of Canada (NSERC) for funding my research 42 of 43
78 Thank you for your attention :) 43 of 43
Image enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationLinear Diffusion. E9 242 STIP- R. Venkatesh Babu IISc
Linear Diffusion Derivation of Heat equation Consider a 2D hot plate with Initial temperature profile I 0 (x, y) Uniform (isotropic) conduction coefficient c Unit thickness (along z) Problem: What is temperature
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationNonlinear Diffusion. 1 Introduction: Motivation for non-standard diffusion
Nonlinear Diffusion These notes summarize the way I present this material, for my benefit. But everything in here is said in more detail, and better, in Weickert s paper. 1 Introduction: Motivation for
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationPartial Differential Equations and Image Processing. Eshed Ohn-Bar
Partial Differential Equations and Image Processing Eshed Ohn-Bar OBJECTIVES In this presentation you will 1) Learn what partial differential equations are and where do they arise 2) Learn how to discretize
More informationConsistent Positive Directional Splitting of Anisotropic Diffusion
In: Boštjan Likar (ed.): Proc. of Computer Vision Winter Workshop, Bled, Slovenia, Feb. 7-9, 2001, pp. 37-48. Consistent Positive Directional Splitting of Anisotropic Diffusion Pavel Mrázek and Mirko Navara
More informationNON-LINEAR DIFFUSION FILTERING
NON-LINEAR DIFFUSION FILTERING Chalmers University of Technology Page 1 Summary Introduction Linear vs Nonlinear Diffusion Non-Linear Diffusion Theory Applications Implementation References Page 2 Introduction
More informationDWI acquisition schemes and Diffusion Tensor estimation
DWI acquisition schemes and Diffusion Tensor estimation A simulation based study Santiago Aja-Fernández, Antonio Tristán-Vega, Pablo Casaseca-de-la-Higuera Laboratory of Image Processing L A B O R A T
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationT H E S I S. Computer Engineering May Smoothing of Matrix-Valued Data. Thomas Brox
D I P L O M A Computer Engineering May 2002 T H E S I S Smoothing of Matrix-Valued Data Thomas Brox Computer Vision, Graphics, and Pattern Recognition Group Department of Mathematics and Computer Science
More informationEstimators for Orientation and Anisotropy in Digitized Images
Estimators for Orientation and Anisotropy in Digitized Images Lucas J. van Vliet and Piet W. Verbeek Pattern Recognition Group of the Faculty of Applied Physics Delft University of Technolo Lorentzweg,
More informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More informationEfficient Beltrami Filtering of Color Images via Vector Extrapolation
Efficient Beltrami Filtering of Color Images via Vector Extrapolation Lorina Dascal, Guy Rosman, and Ron Kimmel Computer Science Department, Technion, Institute of Technology, Haifa 32000, Israel Abstract.
More informationA NEW PARTIAL DIFFERENTIAL EQUATION-BASED APPROACH FOR 3D DATA DENOISING AND EDGE PRESERVING
A NEW PARTIAL DIFFERENTIAL EQUATION-BASED APPROACH FOR 3D DATA DENOISING AND EDGE PRESERVING SORIN POP 1,, OLIVIER LAVIALLE, 3, ROMULUS TEREBES 1, MONICA BORDA 1 Key words: 3-D filtering, Anisotropic diffusion,
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationTensor Visualisation
Tensor Visualisation Computer Animation and Visualisation Lecture 16 Taku Komura tkomura@ed.ac.uk Institute for Perception, Action & Behaviour School of Informatics 1 Tensor Visualisation What is tensor
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say
More informationIntroduction to Nonlinear Image Processing
Introduction to Nonlinear Image Processing 1 IPAM Summer School on Computer Vision July 22, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Mean and median 2 Observations
More informationMedical Image Analysis
Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 23rd January 2006 Outline 1 Recap 2 Edge Enhancement 3 Experimental Results 4 The rest of
More informationTensor Visualisation
Tensor Visualisation Computer Animation and Visualisation Lecture 18 tkomura@ed.ac.uk Institute for Perception, Action & Behaviour School of Informatics Tensors 1 Reminder : Attribute Data Types Scalar
More informationNIH Public Access Author Manuscript Med Image Anal. Author manuscript; available in PMC 2014 July 16.
NIH Public Access Author Manuscript Published in final edited form as: Med Image Anal. 2009 February ; 13(1): 19 35. doi:10.1016/j.media.2008.05.004. Sequential Anisotropic Multichannel Wiener Filtering
More informationRician Noise Removal in Diffusion Tensor MRI
Rician Noise Removal in Diffusion Tensor MRI Saurav Basu, Thomas Fletcher, and Ross Whitaker University of Utah, School of Computing, Salt Lake City, UT 84112, USA Abstract. Rician noise introduces a bias
More informationINTRODUCTION Noise is present in many situations of daily life for ex: Microphones will record noise and speech. Goal: Reconstruct original signal Wie
WIENER FILTERING Presented by N.Srikanth(Y8104060), M.Manikanta PhaniKumar(Y8104031). INDIAN INSTITUTE OF TECHNOLOGY KANPUR Electrical Engineering dept. INTRODUCTION Noise is present in many situations
More informationTime Localised Band Filtering Using Modified S-Transform
Localised Band Filtering Using Modified S-Transform Nithin V George, Sitanshu Sekhar Sahu, L. Mansinha, K. F. Tiampo, G. Panda Department of Electronics and Communication Engineering, National Institute
More informationSeismic Fault Preserving Diffusion
Seismic Fault Preserving Diffusion 1 Olivier Lavialle (1)(), Sorin Pop (1), Christian Germain (1)(), Marc Donias (1), Sebastien Guillon (), Naamen Keskes (), Yannick Berthoumieu (1) 4 5 6 7 8 9 10 11 1
More informationGaussian derivatives
Gaussian derivatives UCU Winter School 2017 James Pritts Czech Tecnical University January 16, 2017 1 Images taken from Noah Snavely s and Robert Collins s course notes Definition An image (grayscale)
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationSOS Boosting of Image Denoising Algorithms
SOS Boosting of Image Denoising Algorithms Yaniv Romano and Michael Elad The Technion Israel Institute of technology Haifa 32000, Israel The research leading to these results has received funding from
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationConfidence and curvature estimation of curvilinear structures in 3-D
Confidence and curvature estimation of curvilinear structures in 3-D P. Bakker, L.J. van Vliet, P.W. Verbeek Pattern Recognition Group, Department of Applied Physics, Delft University of Technology, Lorentzweg
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationGeneralizing the mean intercept length tensor for gray-level images
Generalizing the mean intercept length tensor for gray-level images Rodrigo Moreno, Magnus Borga and Örjan Smedby Linköping University Post Print N.B.: When citing this work, cite the original article.
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationReading. 3. Image processing. Pixel movement. Image processing Y R I G Q
Reading Jain, Kasturi, Schunck, Machine Vision. McGraw-Hill, 1995. Sections 4.-4.4, 4.5(intro), 4.5.5, 4.5.6, 5.1-5.4. 3. Image processing 1 Image processing An image processing operation typically defines
More informationThe Silk Damping Tail of the CMB l. Wayne Hu Oxford, December 2002
The Silk Damping Tail of the CMB 100 T (µk) 10 10 100 1000 l Wayne Hu Oxford, December 2002 Outline Damping tail of temperature power spectrum and its use as a standard ruler Generation of polarization
More informationSpatial smoothing using Gaussian processes
Spatial smoothing using Gaussian processes Chris Paciorek paciorek@hsph.harvard.edu August 5, 2004 1 OUTLINE Spatial smoothing and Gaussian processes Covariance modelling Nonstationary covariance modelling
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationConjugate gradient acceleration of non-linear smoothing filters
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Conjugate gradient acceleration of non-linear smoothing filters Knyazev, A.; Malyshev, A. TR25-43 December 25 Abstract The most efficient edge-preserving
More informationCHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION
59 CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION 4. INTRODUCTION Weighted average-based fusion algorithms are one of the widely used fusion methods for multi-sensor data integration. These methods
More informationRestoration of DWI data using a Rician LMMSE estimator
Restoration of DWI data using a Rician LMMSE estimator Santiago Aja-Fernández, Marc Niethammer, Marek Kubicki, Martha E. Shenton, Carl-Fredrik Westin Abstract This paper introduces and analyzes a Linear
More informationA Tensor Variational Formulation of Gradient Energy Total Variation
A Tensor Variational Formulation of Gradient Energy Total Variation Freddie Åström, George Baravdish and Michael Felsberg Linköping University Post Print N.B.: When citing this work, cite the original
More informationScale Space Analysis by Stabilized Inverse Diffusion Equations
Scale Space Analysis by Stabilized Inverse Diffusion Equations Ilya Pollak, Alan S. Willsky and Hamid Krim Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, 77 Massachusetts
More informationEdge detection and noise removal by use of a partial differential equation with automatic selection of parameters
Volume 24, N. 1, pp. 131 150, 2005 Copyright 2005 SBMAC ISSN 0101-8205 www.scielo.br/cam Edge detection and noise removal by use of a partial differential equation with automatic selection of parameters
More informationNIH Public Access Author Manuscript Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2014 May 19.
NIH Public Access Author Manuscript Published in final edited form as: Med Image Comput Comput Assist Interv. 2009 ; 12(0 1): 919 926. Bias of Least Squares Approaches for Diffusion Tensor Estimation from
More informationAzimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen
Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures
More informationRational Bézier Patch Differentiation using the Rational Forward Difference Operator
Rational Bézier Patch Differentiation using the Rational Forward Difference Operator Xianming Chen, Richard F. Riesenfeld, Elaine Cohen School of Computing, University of Utah Abstract This paper introduces
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More informationTensor Visualisation
Tensor Visualisation Computer Animation and Visualisation Lecture 15 Taku Komura tkomura@ed.ac.uk Institute for Perception, Action & Behaviour School of Informatics 1 Overview Tensor Visualisation What
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More informationEnergy-Based Image Simplification with Nonlocal Data and Smoothness Terms
Energy-Based Image Simplification with Nonlocal Data and Smoothness Terms Stephan Didas 1, Pavel Mrázek 2, and Joachim Weickert 1 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer
More informationIterative Reweighted Least Squares
Iterative Reweighted Least Squares Sargur. University at Buffalo, State University of ew York USA Topics in Linear Classification using Probabilistic Discriminative Models Generative vs Discriminative
More informationCompare and Contrast electrical resistivity measurements in linear 4 probe and Montgomery
Compare and Contrast electrical resistivity measurements in linear 4 probe and Montgomery Rabindra Bhattarai Department of Physics Iowa State University What is Electrical Resistivity Electrical resistivity
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More informationWhere is the bark, the tree and the forest, mathematical foundations of multi-scale analysis
Where is the bark, the tree and the forest, mathematical foundations of multi-scale analysis What we observe is not nature itself, but nature exposed to our method of questioning. -Werner Heisenberg M.
More informationRANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor
RANDOM FIELDS AND GEOMETRY from the book of the same name by Robert Adler and Jonathan Taylor IE&M, Technion, Israel, Statistics, Stanford, US. ie.technion.ac.il/adler.phtml www-stat.stanford.edu/ jtaylor
More informationM. A. Gusmão IF-UFRGS
M. A. Gusmão IF-UFRGS - 217 1 FIP164-217/2 Text 9 Mean-field approximation - II Heisenberg Hamiltonian in wave-vector space As we saw in Text 8, the uniform susceptibility does not diverge in the case
More informationA BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT
A BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT Aleksandra Pižurica, Iris Vanhamel, Hichem Sahli, Wilfried Philips and Antonis Katartzis Ghent University,
More informationTensor fields. Tensor fields: Outline. Chantal Oberson Ausoni
Tensor fields Chantal Oberson Ausoni 7.8.2014 ICS Summer school Roscoff - Visualization at the interfaces 28.7-8.8, 2014 1 Tensor fields: Outline 1. TENSOR FIELDS: DEFINITION 2. PROPERTIES OF SECOND-ORDER
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationDigital Signal Processing
JID:YDSPR AID:53 /FLA [m5g; v1.9; Prn:17/07/2017; 11:02] P.1 (1-13) Digital Signal Processing ( ) 1 Contents lists available at ScienceDirect 2 68 3 69 Digital Signal Processing 6 72 7 www.elsevier.com/locate/dsp
More informationRobust Tensor Splines for Approximation of Diffusion Tensor MRI Data
Robust Tensor Splines for Approximation of Diffusion Tensor MRI Data Angelos Barmpoutis, Baba C. Vemuri, John R. Forder University of Florida Gainesville, FL 32611, USA {abarmpou, vemuri}@cise.ufl.edu,
More informationDiffusion Tensor Imaging quality control : artifacts assessment and correction. A. Coste, S. Gouttard, C. Vachet, G. Gerig. Medical Imaging Seminar
Diffusion Tensor Imaging quality control : artifacts assessment and correction A. Coste, S. Gouttard, C. Vachet, G. Gerig Medical Imaging Seminar Overview Introduction DWI DTI Artifact Assessment Artifact
More informationImage Degradation Model (Linear/Additive)
Image Degradation Model (Linear/Additive),,,,,,,, g x y h x y f x y x y G uv H uv F uv N uv 1 Source of noise Image acquisition (digitization) Image transmission Spatial properties of noise Statistical
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationLeast Squares with Examples in Signal Processing 1. 2 Overdetermined equations. 1 Notation. The sum of squares of x is denoted by x 2 2, i.e.
Least Squares with Eamples in Signal Processing Ivan Selesnick March 7, 3 NYU-Poly These notes address (approimate) solutions to linear equations by least squares We deal with the easy case wherein the
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationTensor-based Image Diffusions Derived from Generalizations of the Total Variation and Beltrami Functionals
Generalizations of the Total Variation and Beltrami Functionals 29 September 2010 International Conference on Image Processing 2010, Hong Kong Anastasios Roussos and Petros Maragos Computer Vision, Speech
More informationTHREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL
THREE-BODY INTERACTIONS DRIVE THE TRANSITION TO POLAR ORDER IN A SIMPLE FLOCKING MODEL Purba Chatterjee and Nigel Goldenfeld Department of Physics University of Illinois at Urbana-Champaign Flocking in
More informationSpace curves, vector fields and strange surfaces. Simon Salamon
Space curves, vector fields and strange surfaces Simon Salamon Lezione Lagrangiana, 26 May 2016 Curves in space 1 A curve is the path of a particle moving continuously in space. Its position at time t
More informationTraveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis
Geophysical Prospecting, 2013, 61, 495 503 doi: 10.1111/j.1365-2478.2012.01067.x Traveltime approximations for inhomogeneous transversely isotropic media with a horizontal symmetry axis Tariq Alkhalifah
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationReview Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.
Spatial Filtering Dr. Praveen Sankaran Department of ECE NIT Calicut January 7, 203 Outline 2 Linear Nonlinear 3 Spatial Domain Refers to the image plane itself. Direct manipulation of image pixels. Figure:
More informationGradient and smoothing stencils with isotropic discretization error. Abstract
Gradient and smoothing stencils with isotropic discretization error Tommy Anderberg (Dated: April 30, 2012) Abstract Gradient stencils with isotropic O(h 2 ) and O(h 4 ) discretization error are constructed
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More informationSurface x(u, v) and curve α(t) on it given by u(t) & v(t). Math 4140/5530: Differential Geometry
Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt ( ds dt )2 Surface x(u, v)
More informationLesson 04. KAZE, Non-linear diffusion filtering, ORB, MSER. Ing. Marek Hrúz, Ph.D.
Lesson 04 KAZE, Non-linear diffusion filtering, ORB, MSER Ing. Marek Hrúz, Ph.D. Katedra Kybernetiky Fakulta aplikovaných věd Západočeská univerzita v Plzni Lesson 04 KAZE ORB: an efficient alternative
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationLecture 11: Differential Geometry
Lecture 11: Differential Geometry c Bryan S. Morse, Brigham Young University, 1998 2000 Last modified on February 28, 2000 at 8:45 PM Contents 11.1 Introduction..............................................
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationQuantitative Understanding in Biology Principal Components Analysis
Quantitative Understanding in Biology Principal Components Analysis Introduction Throughout this course we have seen examples of complex mathematical phenomena being represented as linear combinations
More informationENGINEERINGMATHEMATICS-I. Hrs/Week:04 Exam Hrs: 03 Total Hrs:50 Exam Marks :100
ENGINEERINGMATHEMATICS-I CODE: 14MAT11 IA Marks:25 Hrs/Week:04 Exam Hrs: 03 Total Hrs:50 Exam Marks :100 UNIT I Differential Calculus -1 Determination of n th order derivatives of Standard functions -
More informationCosmological parameters of modified gravity
Cosmological parameters of modified gravity Levon Pogosian Simon Fraser University Burnaby, BC, Canada In collaborations with R. Crittenden, A. Hojjati, K. Koyama, A. Silvestri, G.-B. Zhao Two questions
More informationSV6: Polynomial Regression and Neural Networks
Signal and Information Processing Laboratory Institut für Signal- und Informationsverarbeitung Fachpraktikum Signalverarbeitung SV6: Polynomial Regression and Neural Networks 1 Introduction Consider the
More informationColorado School of Mines Image and Multidimensional Signal Processing
Image and Multidimensional Signal Processing Professor William Hoff Department of Electrical Engineering and Computer Science Spatial Filtering Main idea Spatial filtering Define a neighborhood of a pixel
More informationMULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING. Kaitlyn Beaudet and Douglas Cochran
MULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING Kaitlyn Beaudet and Douglas Cochran School of Electrical, Computer and Energy Engineering Arizona State University, Tempe AZ 85287-576 USA ABSTRACT The problem
More informationProfessor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x
Spring, 997 ME85 - Introduction to Continuum Mechanics Midterm Exam II roblem. (+ points) (a) Let ρ be the mass density, v be the velocity vector, be the Cauchy stress tensor, and b be the body force per
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationClass 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space
Class 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space We have seen that diffraction occurs when, in reciprocal space, Let us now plot this information. Let us designate
More informationNOISE ENHANCED ANISOTROPIC DIFFUSION FOR SCALAR IMAGE RESTORATION. 6, avenue du Ponceau Cergy-Pontoise, France.
NOISE ENHANCED ANISOTROPIC DIFFUSION FOR SCALAR IMAGE RESTORATION Aymeric HISTACE 1, David ROUSSEAU 2 1 Equipe en Traitement d'image et du Signal UMR CNRS 8051 6, avenue du Ponceau 95014 Cergy-Pontoise,
More informationGraceful exit from inflation for minimally coupled Bianchi A scalar field models
Graceful exit from inflation for minimally coupled Bianchi A scalar field models Florian Beyer Reference: F.B. and Leon Escobar (2013), CQG, 30(19), p.195020. University of Otago, Dunedin, New Zealand
More informationStability of Recursive Gaussian Filtering for Piecewise Linear Bilateral Filtering
Stability of Recursive Gaussian Filtering for Piecewise Linear Bilateral Filtering Koichiro Watanabe, Yoshihiro Maeda, and Norishige Fukushima Nagoya Institute of Technology, Nagoya, Japan fukushima@nitech.ac.jp
More information