Random walks for deformable image registration Dana Cobzas
|
|
- Moses Willis
- 5 years ago
- Views:
Transcription
1 Random walks for deformable image registration Dana Cobzas with Abhishek Sen, Martin Jagersand, Neil Birkbeck, Karteek Popuri Computing Science, University of Alberta, Edmonton
2 Deformable image registration?t I target J source Outline: Continuous formulation of deformable registration Discrete formulation and related work Diffusion and regularization from continuous to discrete Random walker for deformable registration Examples and discussion
3 Deformable image registration?t transformation I target I,J :R, R d T : R d, T =Idu J source Eu= E D I, J T u E R u Data similarity regularization u=argmin u E u u=(u x,u y ) displacement field
4 Variational image registration u=argmin u x I x, J xu x d x Data similarity (point-wise score) x J x, u x d x Adaptable regularization + EXAMPLE: SSD Linear isotropic regularization I x, J xux=i x J xu x 2 J x, u x=w J x 2 u x 2 u y 2 scalar valued diffusivity
5 Euler-Lagrange equations min u E D I, J T x w. u x 2 u y 2 d x E D u x w. u x =0 E D u y w. u y =0 Linear isotropic diffusion Continuous solution to deformable registration: Nonlinear optimization (gradient descent, Gauss-Newton...) Semi-implicit scheme to solve the EL PDEs Regularization can be decoupled from the data term update demon's algorithm Main limitation local minima > need for better frameworks
6 Discrete deformable registration u i D I I i, i=1..n J i, i=1..n J Image = Graph G = (N, E) d k Quantized deformations D = {d 1, d 2,...,d K } Labels L = {u 1, u 2,...,u K } Note: D can be very large Ex: deformations [0-20pix] 41 3 = labels
7 Discrete deformable registration u i D I I i, i=1..n J i, i=1..n J Image = Graph G = (N, E) Registration as labeling : MRF Quantized deformations D = {d 1, d 2,...,d k } Labels L = {u 1, u 2,...,u k } Note: D can be very large Ex: deformations [0-20pix] 41 3 = labels Eu= i N i u i i, j E ij u i, u j data similarity Regularization Interaction
8 Related work Eu= i N i u i i, j E ij u i, u j Global optimum ij metric graph cut optimization [T.W. Tang et al. MICCAI 2007][L. Tang et al. MICCAI 2010] More complex interaction terms linear programming method based on primal-dual principle [Glocker et al. MIA 2008]
9 Related work Eu= i N i u i i, j E ij u i, u j Global optimum ij metric graph cut optimization [T.W. Tang et al. MICCAI 2007][L. Tang et al. MICCAI 2010] More complex interaction terms linear programming method based on primal-dual principle [Glocker et al. MIA 2008] We propose : Gaussian MRF for deformable registration solved using random walker algorithm [Grady PAMI 2006, CVPR 2005] Fast algorithm Can easily incorporate a large number of displacements Robustness to noise Unique global minimum (has to be relaxed)
10 Related work Eu= i N i u i i, j E ij u i, u j ij metric graph cut optimization [T.W. Tang et al. MICCAI 2007][L. Tang et al. MICCAI 2010] More complex interaction terms linear programming method based on primal-dual principle [Glocker et al. MIA 2008] Here: Gaussian MRF for deformable registration solved using random walker algorithm [Grady PAMI 2006, CVPR 2005] Fast algorithm Can easily incorporate a large number of displacements Robustness to noise Unique global minimum (has to be relaxed) Regularization E R u= u 2 dxdy diffusion u=0 NEXT: - diffusion and regularization - diffusion on graphs - RW for image registration - Experiments and discussion
11 Diffusion u t Uniform =div u=u u t Isotropic =div w u u0=u 0 u0=u 0 u(0) w s 2 = 1 s 2 s= u 0 u(t) u(t)... Scalar-valued diffusivity Changes the metric of the space Slows down/blocks diffusion at edges Steady state
12 Diffusion u t Uniform =div u=u u t Isotropic =div w u u0=u 0 u0=u 0 u(0) w s 2 = 1 s 2 s= u 0 u(t) u(t)... Steady state Scalar-valued diffusivity Changes the metric of the space Slows down/blocks diffusion at edges Numerical solution: Parabolic PDEs fully implicit time discretization > elliptic PDEs solve with implicit or semi-implicit FD schemes
13 Regularization and diffusion f u? Energy functional : Euler Lagrange Equation : E u u= 1 2 f u 2 u 2 dxdy u=argmin u E u u u f =u data regularization
14 Regularization and diffusion f u? Energy functional : Euler Lagrange Equation : data E u u= 1 2 f u 2 u 2 dxdy u=argmin u E u u u f regularization Uniform diffusion u =u t =u Fully implicit time discretization of the uniform diffusion with initial value f and time step α
15 Regularization Uniform Isotropic w E R u= u 2 dxdy E R u= w f 2 u 2 dxdy Euler Lagrange u f =u u f = w. u
16 Regularization Uniform Isotropic w E R u= u 2 dxdy E R u= w f 2 u 2 dxdy Euler Lagrange u f =u u f = w. u Deformable registration min u E D I, J T x w. u x 2 u y 2 d x E D u x w. u x =0 E D u y w. u y =0
17 u : u=0 ux=f B, x K B Discrete diffusion Finite differences + implicit scheme Boundary conditions u i-1,j-1 u i-1,j u i-1,j+1 u i,j-1 u i,j u i,j+1 u i+1,j-1 u i+1,j u i+1,j+1 Pixel location u Nx1 L NxN Lu=0 [ L u B B][ u B] u B L f =0 L u u u = B T f B
18 u: u=0 ux=f B, x K B Discrete diffusion Finite differences + implicit scheme Boundary conditions u i-1,j-1 u i-1,j u i-1,j+1 u i,j-1 u i,j u i,j+1 u i+1,j-1 u i+1,j u i+1,j+1 Pixel location Lu=0 [ L u B B][ u B] u B L f =0 L u u u = B T f B Generalize L isotropic diffusion Combinatorial Laplacian L ij = d i = k w ik w ij if i= j if i, j E 0 otherwise Positive semi-definite (PSD)
19 u=0 ux=f B, x K B Discrete diffusion Boundary conditions Combinatorial Dirichlet integral u: Finite differences + implicit scheme 1 u i-1,j-1 u i-1,j u i-1,j u i,j-1 u i,j u i,j+1 1 u i+1,j-1 u i+1,j u i+1,j+1 Pixel location D [u]= u 2 dx D[u ]=u t A T A u=u t L u A (ij)k incidence matrix (combinatorial gradient) Critical points = minima [L-PSD] Lx=0 Lu=0 [ L u B B][ u B] u B L f =0 L u u u = B T f B Generalize L isotropic diffusion Combinatorial Laplacian L ij = d i = k w ik w ij if i= j if i, j E 0 otherwise Positive semi-definite (PSD)
20 Random walker algorithm Lu=0 [ L u B B][ u B] u B L f =0 L u u u = B T f B [Grady PAMI 2006] Combinatorial Laplacian L ij = d i = k w ik if i= j if i, j E 0 otherwise w ij w ij =exp I i I j 2 f B seeds ; u i prob that a random walker from node i first reaches a seed.
21 RW for deformable registration u i D Quantized deformations D = {d 1, d 2,...,d K } Labels L = {u 1, u 2,...,u K } I i, i=1..n Registration as labeling : MRF J i, i=1..n E u= i N i u i i, j E ij u i,u j data similarity potential Regularization Interaction
22 RW for deformable registration u i D Quantized deformations D = {d 1, d 2,...,d K } Labels L = {u 1, u 2,...,u K } I i, i=1..n Registration as labeling : MRF J i, i=1..n E u= i N i u i i, j E ij u i,u j data similarity potential Gaussian MRF : Regularization Interaction E k u k = i N l=1: K,l k i l u i k 2 i k 1 u i k 2 ij E w ij u i k u j k 2 Modifications : 1. u i k = prob of node i having label u k 2. Gaussian MRF with interaction u i k,u j k =w ij u i k u j k 2 3. 'priors' for probability of a displacement i k =exp I i J id k 2
23 Gaussian MRF for deformable registration Displacement probability i k =exp I i J id k 2 1 if displacement d k fits perfectly the similarity score Small - otherwise Data potential i u i k = l =1: K,l k i l u i k 2 i k 1 u i k 2 Displacement l k has high prob. > low prob. for label k Displacement k has high prob. > high prob. for label k
24 Gaussian MRF for deformable registration Displacement probability i k =exp I i J id k 2 1 if displacement d k fits perfectly the similarity score Small - otherwise Data potential Regularization Interaction i u i k = l =1: K,l k i l u i k 2 i k 1 u i k 2 Displacement l k has high prob. > low prob. for label k u i k,u j k =w ij u i k u j k 2 w ij =exp J i J j 2 Displacement k has high prob. > high prob. for label k Recall : regularization term (Dirichlet integral) Encourages two neighbors to have similar labels if the image intensity is similar. D[u]= w J 2 u 2 dx D[ u]=u t Lu
25 RW with priors [following Grady CVPR 2005] E k u k = i N l=1 : K, l k i l u i k 2 i k 1 u i k 2 ij E w ij u i k u j k 2 Combinatorial Laplacian L L ij = d i = k w ik if i= j if i, j E 0 otherwise w ij k =diag k
26 RW with priors [following Grady CVPR 2005] E k u k = i N l=1: K,l k i l u i k 2 i k 1 u i k 2 ij E w ij u i k u j k 2 Combinatorial Laplacian L L ij = d i = k w ik if i= j if i, j E 0 otherwise w ij k =diag k E k u k = l=1 : K,l k u kt l u k 1 u k T k 1 u k u kt Lu k Minimum when K L l=1 l u k = k Note: The combined matrix still positive semi-definite One equation system per label!!
27 Augmented graph K L l=1 l u k = k Combinatorial Laplace equation for an augmented graph w ij =exp J i J j 2 Large (1) if neighbors have similar color (allows diffusion)
28 Augmented graph K L l=1 l u k = k Combinatorial Laplace equation for an augmented graph u u 1l u K i l =exp I i J id l 2 w ij =exp J i J j 2 Large (1) if Matching color (SSD) Large (1) if neighbors have similar color (allows diffusion)
29 Implementation Computational cost Need to solve for each displacement label (actually K-1) Size of K L l=1 l u k = k Multi-resolution pyramid Loose global optimality displacements # range ¼ 60 [15,15] ½ 40 [-10,10] 1 30 [-7,7]
30 Implementation Computational cost Need to solve for each displacement label (actually K-1) Size of K L l=1 l u k = k Multi-resolution pyramid Loose global optimality displacements # range ¼ 60 [15,15] ½ 40 [-10,10] 1 30 [-7,7]
31 Synthetic image + Experiments: Quality of deformations known deformation field Demon's alg. RW 256x256 image 3 levels resol. 160sec (MATLAB) Ang. err. (deg) 1.58 ± ± 0.97 Mag. Err. (deg) 1.12 ± ± 2.89 SSD [0-255]
32 Brain MRI with WM/GM/CSF segmentation source target Reg. Demons SSD=8.82 Reg. RW SSD=3.38 Data from Brain segmentation Repository Dice
33 Dice Abdominal CT with muscle segmentation source target Data from CCI, Edmonton Reg. Demons SSD=15.46 Reg. RW SSD=9.64
34 Discussion CONTRIBUTION: New discrete formulation of deformable image registration based on the RW method Image dependent regularization term Global solution with respect to discretization LIMITATIONS AND EXTENSIONS : Computational cost : Solve a large equation system for each label FFD Requirement of a point-wise similarity score Eu= i N i u i i, j E ij u i, u j Possible solution Lower dim models : Compute displacements only at control points Approximate non-local scores in a neighborhood defined by the control points FFD
Random walks and anisotropic interpolation on graphs. Filip Malmberg
Random walks and anisotropic interpolation on graphs. Filip Malmberg Interpolation of missing data Assume that we have a graph where we have defined some (real) values for a subset of the nodes, and that
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationMotion Estimation (I)
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
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 informationMulti-class DTI Segmentation: A Convex Approach
Multi-class DTI Segmentation: A Convex Approach Yuchen Xie 1 Ting Chen 2 Jeffrey Ho 1 Baba C. Vemuri 1 1 Department of CISE, University of Florida, Gainesville, FL 2 IBM Almaden Research Center, San Jose,
More informationIntroduction to motion correspondence
Introduction to motion correspondence 1 IPAM - UCLA July 24, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Why estimate visual motion? 2 Tracking Segmentation Structure
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationApplied Mathematics 505b January 22, Today Denitions, survey of applications. Denition A PDE is an equation of the form F x 1 ;x 2 ::::;x n ;~u
Applied Mathematics 505b January 22, 1998 1 Applied Mathematics 505b Partial Dierential Equations January 22, 1998 Text: Sobolev, Partial Dierentail Equations of Mathematical Physics available at bookstore
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationSpace-Variant Computer Vision: A Graph Theoretic Approach
p.1/65 Space-Variant Computer Vision: A Graph Theoretic Approach Leo Grady Cognitive and Neural Systems Boston University p.2/65 Outline of talk Space-variant vision - Why and how of graph theory Anisotropic
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More 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 informationFantope Regularization in Metric Learning
Fantope Regularization in Metric Learning CVPR 2014 Marc T. Law (LIP6, UPMC), Nicolas Thome (LIP6 - UPMC Sorbonne Universités), Matthieu Cord (LIP6 - UPMC Sorbonne Universités), Paris, France Introduction
More informationEfficient Inference in Fully Connected CRFs with Gaussian Edge Potentials
Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials by Phillip Krahenbuhl and Vladlen Koltun Presented by Adam Stambler Multi-class image segmentation Assign a class label to each
More informationLearning Spectral Graph Segmentation
Learning Spectral Graph Segmentation AISTATS 2005 Timothée Cour Jianbo Shi Nicolas Gogin Computer and Information Science Department University of Pennsylvania Computer Science Ecole Polytechnique Graph-based
More informationGeneralized Newton-Type Method for Energy Formulations in Image Processing
Generalized Newton-Type Method for Energy Formulations in Image Processing Leah Bar and Guillermo Sapiro Department of Electrical and Computer Engineering University of Minnesota Outline Optimization in
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 informationOptic Flow Computation with High Accuracy
Cognitive Computer Vision Colloquium Prague, January 12 13, 2004 Optic Flow Computation with High ccuracy Joachim Weickert Saarland University Saarbrücken, Germany joint work with Thomas Brox ndrés Bruhn
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationModulation-Feature based Textured Image Segmentation Using Curve Evolution
Modulation-Feature based Textured Image Segmentation Using Curve Evolution Iasonas Kokkinos, Giorgos Evangelopoulos and Petros Maragos Computer Vision, Speech Communication and Signal Processing Group
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationCS598 Machine Learning in Computational Biology (Lecture 5: Matrix - part 2) Professor Jian Peng Teaching Assistant: Rongda Zhu
CS598 Machine Learning in Computational Biology (Lecture 5: Matrix - part 2) Professor Jian Peng Teaching Assistant: Rongda Zhu Feature engineering is hard 1. Extract informative features from domain knowledge
More informationSolving linear systems (6 lectures)
Chapter 2 Solving linear systems (6 lectures) 2.1 Solving linear systems: LU factorization (1 lectures) Reference: [Trefethen, Bau III] Lecture 20, 21 How do you solve Ax = b? (2.1.1) In numerical linear
More informationMarkov Random Fields
Markov Random Fields Umamahesh Srinivas ipal Group Meeting February 25, 2011 Outline 1 Basic graph-theoretic concepts 2 Markov chain 3 Markov random field (MRF) 4 Gauss-Markov random field (GMRF), and
More informationMAP Examples. Sargur Srihari
MAP Examples Sargur srihari@cedar.buffalo.edu 1 Potts Model CRF for OCR Topics Image segmentation based on energy minimization 2 Examples of MAP Many interesting examples of MAP inference are instances
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More 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 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 informationSpectral Clustering. Spectral Clustering? Two Moons Data. Spectral Clustering Algorithm: Bipartioning. Spectral methods
Spectral Clustering Seungjin Choi Department of Computer Science POSTECH, Korea seungjin@postech.ac.kr 1 Spectral methods Spectral Clustering? Methods using eigenvectors of some matrices Involve eigen-decomposition
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationAdaptive diffeomorphic multi-resolution demons and its application to same. modality medical image registration with large deformation 1
Adaptive diffeomorphic multi-resolution demons and its application to same modality medical image registration with large deformation 1 Chang Wang 1 Qiongqiong Ren 1 Xin Qin 1 Yi Yu 1 1 (School of Biomedical
More informationIntroduction to Machine Learning Midterm Exam Solutions
10-701 Introduction to Machine Learning Midterm Exam Solutions Instructors: Eric Xing, Ziv Bar-Joseph 17 November, 2015 There are 11 questions, for a total of 100 points. This exam is open book, open notes,
More informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationNONRIGID IMAGE REGISTRATION WITH TWO-SIDED SPACE-FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS
NONRIGID IMAGE REGISTRATION WITH TWO-SIDED SPACE-FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS Clarissa C. Garvey, Nathan D. Cahill, Andrew Melbourne, Christine Tanner, Sebastien Ourselin, and David J. Hawkes
More informationLearning a Degree-Augmented Distance Metric from a Network. Bert Huang, U of Maryland Blake Shaw, Foursquare Tony Jebara, Columbia U
Learning a Degree-Augmented Distance Metric from a Network Bert Huang, U of Maryland Blake Shaw, Foursquare Tony Jebara, Columbia U Beyond Mahalanobis: Supervised Large-Scale Learning of Similarity NIPS
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
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 informationFacebook Friends! and Matrix Functions
Facebook Friends! and Matrix Functions! Graduate Research Day Joint with David F. Gleich, (Purdue), supported by" NSF CAREER 1149756-CCF Kyle Kloster! Purdue University! Network Analysis Use linear algebra
More informationConvex Optimization Algorithms for Machine Learning in 10 Slides
Convex Optimization Algorithms for Machine Learning in 10 Slides Presenter: Jul. 15. 2015 Outline 1 Quadratic Problem Linear System 2 Smooth Problem Newton-CG 3 Composite Problem Proximal-Newton-CD 4 Non-smooth,
More informationSelf-Tuning Semantic Image Segmentation
Self-Tuning Semantic Image Segmentation Sergey Milyaev 1,2, Olga Barinova 2 1 Voronezh State University sergey.milyaev@gmail.com 2 Lomonosov Moscow State University obarinova@graphics.cs.msu.su Abstract.
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationCovariance Matrix Simplification For Efficient Uncertainty Management
PASEO MaxEnt 2007 Covariance Matrix Simplification For Efficient Uncertainty Management André Jalobeanu, Jorge A. Gutiérrez PASEO Research Group LSIIT (CNRS/ Univ. Strasbourg) - Illkirch, France *part
More informationGrouping with Bias. Stella X. Yu 1,2 Jianbo Shi 1. Robotics Institute 1 Carnegie Mellon University Center for the Neural Basis of Cognition 2
Grouping with Bias Stella X. Yu, Jianbo Shi Robotics Institute Carnegie Mellon University Center for the Neural Basis of Cognition What Is It About? Incorporating prior knowledge into grouping Unitary
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationActive Contours Snakes and Level Set Method
Active Contours Snakes and Level Set Method Niels Chr. Overgaard 2010-09-21 N. Chr. Overgaard Lecture 7 2010-09-21 1 / 26 Outline What is an mage? What is an Active Contour? Deforming Force: The Edge Map
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationRoadmap. Introduction to image analysis (computer vision) Theory of edge detection. Applications
Edge Detection Roadmap Introduction to image analysis (computer vision) Its connection with psychology and neuroscience Why is image analysis difficult? Theory of edge detection Gradient operator Advanced
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
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 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 informationMedical Image Analysis
Medical Image Analysis CS 593 / 791 Computer Science and Electrical Engineering Dept. West Virginia University 20th January 2006 Outline 1 Discretizing the heat equation 2 Outline 1 Discretizing the heat
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationManifold Coarse Graining for Online Semi-supervised Learning
for Online Semi-supervised Learning Mehrdad Farajtabar, Amirreza Shaban, Hamid R. Rabiee, Mohammad H. Rohban Digital Media Lab, Department of Computer Engineering, Sharif University of Technology, Tehran,
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationConnectedness of Random Walk Segmentation
Connectedness of Random Walk Segmentation Ming-Ming Cheng and Guo-Xin Zhang TNList, Tsinghua University, Beijing, China. Abstract Connectedness of random walk segmentation is examined, and novel properties
More informationarxiv: v2 [cs.cv] 17 May 2017
Misdirected Registration Uncertainty Jie Luo 1,2, Karteek Popuri 3, Dana Cobzas 4, Hongyi Ding 5 William M. Wells III 2,6 and Masashi Sugiyama 7,1 arxiv:1704.08121v2 [cs.cv] 17 May 2017 1 Graduate School
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationMarching on the BL equations
Marching on the BL equations Harvey S. H. Lam February 10, 2004 Abstract The assumption is that you are competent in Matlab or Mathematica. White s 4-7 starting on page 275 shows us that all generic boundary
More informationThe University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.
The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment 1 Caramanis/Sanghavi Due: Thursday, Feb. 7, 2013. (Problems 1 and
More informationRandom coordinate descent algorithms for. huge-scale optimization problems. Ion Necoara
Random coordinate descent algorithms for huge-scale optimization problems Ion Necoara Automatic Control and Systems Engineering Depart. 1 Acknowledgement Collaboration with Y. Nesterov, F. Glineur ( Univ.
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationIntelligent Systems:
Intelligent Systems: Undirected Graphical models (Factor Graphs) (2 lectures) Carsten Rother 15/01/2015 Intelligent Systems: Probabilistic Inference in DGM and UGM Roadmap for next two lectures Definition
More informationWeighted Nonlocal Laplacian on Interpolation from Sparse Data
Noname manuscript No. (will be inserted by the editor) Weighted Nonlocal Laplacian on Interpolation from Sparse Data Zuoqiang Shi Stanley Osher Wei Zhu Received: date / Accepted: date Abstract Inspired
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationGeneralization to inequality constrained problem. Maximize
Lecture 11. 26 September 2006 Review of Lecture #10: Second order optimality conditions necessary condition, sufficient condition. If the necessary condition is violated the point cannot be a local minimum
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
More informationVariational Methods in Image Denoising
Variational Methods in Image Denoising Jamylle Carter Postdoctoral Fellow Mathematical Sciences Research Institute (MSRI) MSRI Workshop for Women in Mathematics: Introduction to Image Analysis 22 January
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Part I: basic optimization concepts in R n Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada
More informationComputer Vision Group Prof. Daniel Cremers. 14. Clustering
Group Prof. Daniel Cremers 14. Clustering Motivation Supervised learning is good for interaction with humans, but labels from a supervisor are hard to obtain Clustering is unsupervised learning, i.e. it
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationSupport Vector Machines and Kernel Methods
2018 CS420 Machine Learning, Lecture 3 Hangout from Prof. Andrew Ng. http://cs229.stanford.edu/notes/cs229-notes3.pdf Support Vector Machines and Kernel Methods Weinan Zhang Shanghai Jiao Tong University
More informationKernel metrics on normal cycles for the matching of geometrical structures.
Kernel metrics on normal cycles for the matching of geometrical structures. Pierre Roussillon University Paris Descartes July 18, 2016 1 / 34 Overview Introduction Matching of geometrical structures Kernel
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationLow Rank Matrix Completion Formulation and Algorithm
1 2 Low Rank Matrix Completion and Algorithm Jian Zhang Department of Computer Science, ETH Zurich zhangjianthu@gmail.com March 25, 2014 Movie Rating 1 2 Critic A 5 5 Critic B 6 5 Jian 9 8 Kind Guy B 9
More informationLocally-biased analytics
Locally-biased analytics You have BIG data and want to analyze a small part of it: Solution 1: Cut out small part and use traditional methods Challenge: cutting out may be difficult a priori Solution 2:
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationLinear & nonlinear classifiers
Linear & nonlinear classifiers Machine Learning Hamid Beigy Sharif University of Technology Fall 1396 Hamid Beigy (Sharif University of Technology) Linear & nonlinear classifiers Fall 1396 1 / 44 Table
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationLinear & nonlinear classifiers
Linear & nonlinear classifiers Machine Learning Hamid Beigy Sharif University of Technology Fall 1394 Hamid Beigy (Sharif University of Technology) Linear & nonlinear classifiers Fall 1394 1 / 34 Table
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationDiscriminative Learning and Big Data
AIMS-CDT Michaelmas 2016 Discriminative Learning and Big Data Lecture 2: Other loss functions and ANN Andrew Zisserman Visual Geometry Group University of Oxford http://www.robots.ox.ac.uk/~vgg Lecture
More informationPart 7: Structured Prediction and Energy Minimization (2/2)
Part 7: Structured Prediction and Energy Minimization (2/2) Colorado Springs, 25th June 2011 G: Worst-case Complexity Hard problem Generality Optimality Worst-case complexity Integrality Determinism G:
More informationLaplacian Eigenimages in Discrete Scale Space
Laplacian Eigenimages in Discrete Scale Space Martin Tschirsich and Arjan Kuijper,2 Technische Universität Darmstadt, Germany 2 Fraunhofer IGD, Darmstadt, Germany Abstract. Linear or Gaussian scale space
More informationMarkov Random Fields and Bayesian Image Analysis. Wei Liu Advisor: Tom Fletcher
Markov Random Fields and Bayesian Image Analysis Wei Liu Advisor: Tom Fletcher 1 Markov Random Field: Application Overview Awate and Whitaker 2006 2 Markov Random Field: Application Overview 3 Markov Random
More informationGeneralized Gradient Descent Algorithms
ECE 275AB Lecture 11 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San Diego p. 1/1 Lecture 11 ECE 275A Generalized Gradient Descent Algorithms ECE 275AB Lecture 11 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San
More informationHandout to Wu Characteristic Harvard math table talk Oliver Knill, 3/8/2016
Handout to Wu Characteristic Harvard math table talk Oliver Knill, 3/8/2016 Definitions Let G = (V, E) be a finite simple graph with vertex set V and edge set E. The f-vector v(g) = (v 0 (G), v 1 (G),...,
More informationWhy should you care about the solution strategies?
Optimization Why should you care about the solution strategies? Understanding the optimization approaches behind the algorithms makes you more effectively choose which algorithm to run Understanding the
More informationA hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
More information