Non-Rigid Registration I
|
|
- Bethanie Waters
- 6 years ago
- Views:
Transcription
1 Non-Rigid Registration I CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (CS6240) Non-Rigid Registration I 1 / 67
2 Introduction Introduction Non-rigid registration dues with spatial alignment of non-rigid or deformable objects. Can include global rigid alignment, e.g., scaling, rotation, translation. Can be performed with or without known correspondence. Deformable models Can change shape over time. So, transformation is non-linear. Can be represented by equations or connected points. Leow Wee Kheng (CS6240) Non-Rigid Registration I 2 / 67
3 Registration with Known Correspondence Curve Warping Curve Warping Warp a source curve C into a target curve C in 2D space p i C p i C p i = (x i, y i ) denote the points on C. p i = (x i, y i ) denote the corresponding points on C. Warping is nonlinear because the displacement of points cannot be described by a linear transformation. Leow Wee Kheng (CS6240) Non-Rigid Registration I 3 / 67
4 Registration with Known Correspondence Curve Warping A simple representation of a 2D nonlinear transformation is a polynomial function of order m. x i = a pq x p i yq i, p q 0 p + q m, (1) y i = b pq x p i yq i, p q 0 p + q m. (2) Example: 2nd-order Polynomial Function. x i = a 20 x 2 i + a 11 x i y i + a 02 yi 2 + a 10 x + a 01 y + a 00 (3) y i = b 20 x 2 i + b 11 x i y i + b 02 yi 2 + b 10 x + b 01 y + b 00. (4) Eq. 1 and 2 are independent. So, they can be solved individually using model fitting method. Exercise: Curved surface warping by polynomial transformation. Leow Wee Kheng (CS6240) Non-Rigid Registration I 4 / 67
5 Registration with Known Correspondence Image Warping Image Warping Similar method can be used for image warping. p i = (x i, y i ): control points in source image I. p i = (x i, y i ): target locations of control points. Leow Wee Kheng (CS6240) Non-Rigid Registration I 5 / 67
6 Registration with Known Correspondence Image Warping General idea of algorithm: 1 Use known correspondence to compute the nonlinear transformation between p i and p i. 2 For each pixel location p in the target warped image, Compute the real-valued location p is the source image. Use bilinear interpolation to compute the color c of p. Fill pixel p with color c. Notes: In practice, a global transformation applied to the entire source image is not practical may require very high-order polynomial. So, divide the source image into triangular regions using triangulation algorithm (CS5237 Computational Geometry and Applications). Perform a different local transformation within each triangle. Ensure transformations across triangle s edges are smooth. For details, see [Wol90]. Leow Wee Kheng (CS6240) Non-Rigid Registration I 6 / 67
7 Registration with Known Correspondence Image Warping (c) Source image (d) Warped image. Leow Wee Kheng (CS6240) Non-Rigid Registration I 7 / 67
8 Registration with Known Correspondence Further Readings Further Readings 3D volume warping [BMTB02]. Details of image warping [Wol90]. Image morphing (CS5245 lecture notes). Leow Wee Kheng (CS6240) Non-Rigid Registration I 8 / 67
9 Registration with Unknown Correspondence Registration with Unknown Correspondence Registration of a deformable model with unknown correspondence is an ill-posed problem [Had23]. Tikhonov and Arsenin (1977): A function f that maps inputs in X to outputs in Y is well-posed if these conditions are satisfied: 1 Existence. x X, y Y such that y = f(x). 2 Uniqueness. For any pair x, t X, f(x) = f(v) if and only if x = v. 3 Continuity. Let ρ x denote the distance measured in X. For any ǫ > 0, δ = δ(ǫ) such that ρ x (x,v) < δ implies ρ y (f(x), f(v)) < ǫ. Leow Wee Kheng (CS6240) Non-Rigid Registration I 9 / 67
10 Registration with Unknown Correspondence t f f ( t) x δ f ( x) ε If any of the above is not satisfied, then the problem is ill-posed. Leow Wee Kheng (CS6240) Non-Rigid Registration I 10 / 67
11 Registration with Unknown Correspondence Regularization Theory Regularization Theory A general approach to solving ill-posed problem is based on the regularization theory [TA77]. Include information or constraints about possible solutions. Most common constraint: smoothness, i.e., similar inputs map to similar outputs. Tikhonov s regularization theory involves two terms: 1 Standard error term: e.g., mean squared error E S (f) = i (y i f(x i )) 2 (5) 2 Regularization term: E C (f) = Df 2 (6) D is a linear differential operator. Leow Wee Kheng (CS6240) Non-Rigid Registration I 11 / 67
12 Registration with Unknown Correspondence Regularization Theory The total error E to be minimized is E(f) = E S (f) + λ E C (f) = (y i f(x i )) 2 + λ Df 2 (7) i λ is the regularization parameter. E(f) is a functional: function of function. Leow Wee Kheng (CS6240) Non-Rigid Registration I 12 / 67
13 Registration with Unknown Correspondence Regularization Theory The function f(x) that minimizes E(f) takes the form f(x) = i w i G(x,x i ) (8) G(x,x k ) is called the Green s function. w i = (y i f(x i ))/λ is the weight. That is, f is a linear combination of Green s functions. Gaussian function is an example of Green s function. Regularization is frequently used to solve various problems. Leow Wee Kheng (CS6240) Non-Rigid Registration I 13 / 67
14 Medial Axis Fitting Medial Axis Fitting How to fit a curve (1-D model) in 3-D space? Equation of curve in 3-D space does not have the convenient analytic form shown in previous sections. Other ways of representing curves are required. Parametric Equation s = 0 v ( s) s = 1 s s: a parameter, typically 0 s 1 v(s) = (x(s), y(s), z(s)): coordinates of a point on the curve v is a function of s. Leow Wee Kheng (CS6240) Non-Rigid Registration I 14 / 67
15 Medial Axis Fitting First, consider a continuous tubular object. s = 0 v ( s) s p i χ s s = 1 Want to fit a curve to the center line or medial axis. Use parametric equation of curve to represent medial axis. Cross-section χ s at s is perpendicular to the medial axis. v(s) is the centroid of the surface points p i on χ s. v(s) = 1 p i (9) N s p i χ s where N s is the number of points on χ s. Leow Wee Kheng (CS6240) Non-Rigid Registration I 15 / 67
16 Medial Axis Fitting Now, consider a discrete mesh of tubular object. Cross-section χ s may not intersect any mesh point. So, Eq. 9 is not so appropriate. Instead, think of χ s as a thin slab at s that contains a sufficient number of mesh points. That is, a point p i χ s has a normal projection on the medial axis that is close to v(s). Then, can redefine v(s) as v(s) = w i (s)p i p i χ s. (10) w i (s) p i χ s Weight w i (s) is inversely related to the distance between v(s) and the normal projection of p i on the medial axis. Leow Wee Kheng (CS6240) Non-Rigid Registration I 16 / 67
17 Medial Axis Fitting The problem of fitting medial axis becomes one of finding the v(s) that minimize the difference between v(s) and the right-hand-side of Eq. 10. where E(v) = L 0 D(v) = v(s) D(v)ds (11) p i χ s w i (s)p i p i χ s w i (s) 2 (12) Here, 0 s L. Eq. 11 alone is insufficient for optimization because it is ill-posed: there is an infinite number of curves that can fit. To make the problem well-posed, need an additional constraint, e.g., smoothness constraint. Leow Wee Kheng (CS6240) Non-Rigid Registration I 17 / 67
18 Medial Axis Fitting The revised cost function E(v) is [KWL06] E(v) = 1 2 L λ is a regularization weight. v (s) = dv/ds. Minimizing v (s) 2 imposes smoothness. 0 ( D(v) + λ v (s) 2) ds. (13) When E is minimum, v satisfies the following equation, which can be obtained using variational calculus (Appendix A): That is, p i χ s w i (s)[v(s) p i ] 1 2 D(v(s)) λv (s) = 0. (14) ( p i χ s w i (s) ) 1 λv (s) = 0. (15) Leow Wee Kheng (CS6240) Non-Rigid Registration I 18 / 67
19 Medial Axis Fitting Eq. 15 is solved iteratively by looking at change of v over time, i.e., consider v = v(s, t). Then, v t + Then, p i χ s w i (s)[v(s, t) p i ] ( p i χ s w i (s) ) 1 λv (s, t) = 0. (16) v(s, t + 1) = v(s, t) η ( ) 1 w i (s)[v(s, t) p i ] w i (s) λv (s, t) (17) p i χ s p i χ s where η is a small constant time step. Note: v changes in the direction of negative gradient of D. Leow Wee Kheng (CS6240) Non-Rigid Registration I 19 / 67
20 Medial Axis Fitting Example: Fit a curve to the medial axis of a mesh model of aorta. (a) initialization (b) fitting result Demo Exercise: Work out the algorithm/equation for fitting the cross-sectional radii along the medial axis, after fitting the medial axis. This process gives rise to a generalized cylinder model. Leow Wee Kheng (CS6240) Non-Rigid Registration I 20 / 67
21 Active Contour Active Contour Active contour [Coh91, KWT87] also called snake, more like rubber band can be deformed to match any shape The snake model represents a contour v parametrically as v(s) = (x(s), y(s)) 0 s 1 (18) ( x (1), y (1)) ( x ( s ), y ( s )) Leow Wee Kheng (CS6240) Non-Rigid Registration I 21 / 67
22 Active Contour The snake is deformed under the influence of three forces: internal forces: constrain stretching and banding of snake image forces: attract the snake (e.g., edges) The energy E i of the internal forces at point v(s) is given by E i = 1 ( α(s) v (s) 2 + β(s) v (s) 2) (19) 2 v (s) = dv/ds and v (s) = d 2 v/ds 2. α(s) v (s) 2 is the energy associated with stretching or tension or elasticity β(s) v (s) 2 is the energy associated with bending or rigidity Leow Wee Kheng (CS6240) Non-Rigid Registration I 22 / 67
23 Active Contour The energy E associated with the image depends on what you want the snake to be attracted to: dark line: E = I(x, y) bright line: E = I(x, y) edge: E = I(x, y) 2 Marr-Hildreth edge: E = (G σ 2 I) 2 G σ is a Gaussian function Total energy E T of the snake is E T = 1 0 [E i (v(s)) + E(v(s))] ds (20) Want energy E T to be minimum when snake has snapped onto the desired image feature. Leow Wee Kheng (CS6240) Non-Rigid Registration I 23 / 67
24 Active Contour When E T is minimum, v satisfies the following equation, which can be obtained using variational calculus method (Appendix B): (α(s)v (s)) + (β(s)v (s)) + E(v(s)) = 0 (21) To solve Eq. 21, consider change of snake over time, i.e., v(s, t). Then, we have v(s, t) t (α(s)v (s, t)) + (β(s)v (s, t)) + E(v(s, t)) = 0. (22) Leow Wee Kheng (CS6240) Non-Rigid Registration I 24 / 67
25 Active Contour For numerical computation, let s discretize the equation by finite differences. Let h be a small step or difference of s, and v i = v(ih) = (x(ih), y(ih)) α i = α(ih)/h 2 β i = β(ih)/h 4 F = (F x, F y ) = E(v(s)) Apply forward and backward difference alternatively (Why? See Exercises): forward difference : f (x) = f(x + 1) f(x) backward difference : f (x) = f(x) f(x 1) Leow Wee Kheng (CS6240) Non-Rigid Registration I 25 / 67
26 Active Contour Then, Eq. 21 can be discretized into α i (v i v i 1 ) α i+1 (v i+1 v i ) +β i 1 (v i 2v i 1 + v i 2 ) 2β i (v i+1 2v i + v i 1 ) (23) +β i+1 (v i+2 2v i+1 + v i ) + F = 0 Eq. 23 can be written in the matrix form Notes: A is a pentadiagonal banded matrix x = (x 1, x 2,...,x n ) y = (y 1, y 2,...,y n ) F x = (F x1, F x2,...,f xn ) F y = (F y1, F y2,...,f yn ) Ax + F x = 0 (24) Ay + F y = 0 (25) Leow Wee Kheng (CS6240) Non-Rigid Registration I 26 / 67
27 Active Contour For a n-point closed contour (i.e., x n+1 = x 1 ): A = c 1 d 1 e a 1 b 1 b 2 c 2 d 2 e a 2 a 3 b 3 c 3 d 3 e e n a n 1 b n 1 c n 1 d n 1 d n e n 0 0 a n b n c n (26) where a i = β i 1 b i = α i 2β i 1 2β i c i = α i + α i+1 + β i 1 + 4β i + β i+1 d i = α i+1 2β i 2β i+1 e i = β i+1 (27) Leow Wee Kheng (CS6240) Non-Rigid Registration I 27 / 67
28 Active Contour Equations 24 and 25 are solved iteratively by looking at change of v over time. Then, where γ is a small time step. Then, rearranging the terms yield Note: γ(x(t) x(t 1)) + Ax + F x = 0 (28) x(t) = (A + γi) 1 (γx(t 1) F x (t 1)) (29) y(t) = (A + γi) 1 (γy(t 1) F y (t 1)). (30) For a fixed A, the inverse (A + γi) 1 does not change over time. Then, only have to compute it once. Leow Wee Kheng (CS6240) Non-Rigid Registration I 28 / 67
29 Active Contour Example Demo 1 Demo 2 (a) Initial snake points (b) Final snake contour. Additional points are added between initial points before running the snake algorithm to get smoother contour. Have problem snapping onto concave parts of the contour. Leow Wee Kheng (CS6240) Non-Rigid Registration I 29 / 67
30 Active Contour Gradient Vector Flow Gradient Vector Flow The original formulation of snake has two main shortcomings: It does not converge well to concave features Its performance is sensitive to initial guess of snake point positions. Gradient Vector Flow (GVF) overcome the shortcomings [XP97, XP98, XP00]. Based on diffusion of gradient vectors of edge map. Leow Wee Kheng (CS6240) Non-Rigid Registration I 30 / 67
31 Active Contour Gradient Vector Flow Let E(x, y) denote edge map derived from image I(x, y). Then, E is the gradient vector which is normal to the edge. E is maximum at the edge. GVF field g(x, y) = [u(x, y), v(x, y)] is a vector field that minimizes the energy functional E = µ(u 2 x + u 2 y + vx 2 + vy) 2 + E 2 g E 2 dx dy (31) u x = u/ x, similarly for u y, v x, v y. µ: constant that is set according to the amount of noise present. If E is small, the first term dominates, yielding a slowly varying field. If E is large, the second term dominates, keeping g nearly equal to E. Leow Wee Kheng (CS6240) Non-Rigid Registration I 31 / 67
32 Active Contour Gradient Vector Flow Using variational calculus, it can be shown that the GVF field can be found by solving the following Euler equations: µ 2 u (u E x )(E 2 x + E 2 y) = 0 (32) where 2 is the Laplacian operator: µ 2 v (v E y )(E 2 x + E 2 y) = 0 (33) 2 u = 2 u x u y 2 (34) Eq. 32 and 33 can be solved by regarding u and v as functions of time and solving u t = µ 2 u (u E x )(E 2 x + E 2 y) (35) v t = µ 2 v (v E y )(E 2 x + E 2 y). (36) Leow Wee Kheng (CS6240) Non-Rigid Registration I 32 / 67
33 Active Contour Gradient Vector Flow As for the snake model, apply finite difference approximation: u t = 1 (u(k + 1) u(k)) t v t = 1 (v(k + 1) v(k)) t 2 u = u i+1,j + u i,j+1 + u i 1,j + u i,j 1 4u i,j 2 v = v i+1,j + v i,j+1 + v i 1,j + v i,j 1 4v i,j Leow Wee Kheng (CS6240) Non-Rigid Registration I 33 / 67
34 Active Contour Gradient Vector Flow Substituting the approximations into Eq. 35 and 36 yields where u i,j (k + 1) = (1 b i,j t)u i,j (k) + c i,j t+ µ t [u i+1,j (k) + u i,j+1 (k) + u i 1,j (k) + u i,j 1 (k) 4u i,j (k)] v i,j (k + 1) = (1 b i,j t)v i,j (k) + d i,j t+ µ t [v i+1,j (k) + v i,j+1 (k) + v i 1,j (k) + v i,j 1 (k) 4v i,j (k)] Notes: For numerical stability, t should be small. Specifically, t < 1 4µ. For more details, refer to [XP97, XP98]. (37) b = f 2 x + f 2 y (38) c = bf x (39) d = bf y. (40) Leow Wee Kheng (CS6240) Non-Rigid Registration I 34 / 67
35 Active Contour Gradient Vector Flow Combining snake with GVF: Notes: Replace E in Eq. 22 by g. That is, Compute the steady state solutions of u and v, i.e., u(k) and v(k) as k. Replace F x and F y by u and v. Run snake algorithm. GVF is a very general method of propagating edge information to increase the basin of attraction of deformable models. It can be combined with many other deformable models. Leow Wee Kheng (CS6240) Non-Rigid Registration I 35 / 67
36 Active Contour Gradient Vector Flow Example 1: From [XP00]. (a) convergence of contour (b) traditional potential force (c) close-up at concavity: no force to attract the snake towards the bottom of the concavity Leow Wee Kheng (CS6240) Non-Rigid Registration I 36 / 67
37 Active Contour Gradient Vector Flow (a) convergence of contour (b) GVF external force (c) close-up at concavity: forces exist to attract the snake towards the bottom of the concavity Leow Wee Kheng (CS6240) Non-Rigid Registration I 37 / 67
38 Active Contour Gradient Vector Flow Example 2: From [Tia02]. (a) Initial snake. (b) Traditional snake. (c) GVF snake. Demo 1 Demo 2 Leow Wee Kheng (CS6240) Non-Rigid Registration I 38 / 67
39 Active Contour Important Notes Important Notes How to determine the initial snake points? Possible methods: Manual: User provide initial snake points. OK if user inputs are acceptable. Automatic: Use image processing techniques to find approximate locations of image features. User inputs are not required. Semi-automatic: User provide simple inputs, which assist the program to find approximate locations of image features. Need some user inputs. Leow Wee Kheng (CS6240) Non-Rigid Registration I 39 / 67
40 Active Contour Exercises Exercises 1 Derive the discrete estimate of f (x) using forward difference only backward difference only alternative forward and backward differences What are the differences of the estimations derived above? 2 Derive the Euler equations for Gradient Vector Flow using variational calculus. Leow Wee Kheng (CS6240) Non-Rigid Registration I 40 / 67
41 Active Contour Further Readings Further Readings Appendix B: Variational Calculus for Snake Model. Dual contour method [GN94]: One contour expands from inside the target region while the other contracts from outside. This method can reject certain local minima in the snake evolution process. Topology-preserving snakes [HXP01, WJL03]. Snakes with geometric constraints [Ee04, FCH03, SHD01]. Leow Wee Kheng (CS6240) Non-Rigid Registration I 41 / 67
42 Level Set Method Level Set Method Snakes cannot handle applications that require topology changes. (a) (b) (c) Leow Wee Kheng (CS6240) Non-Rigid Registration I 42 / 67
43 Level Set Method (a) Two bubbles expand outward. (b) Later in time, they merge into one, thus changing the topology. (c) With point representation, need to decide which points are located inside the merged bubble; they need to be removed. Doing this in 2D is troublesome. Doing in 3D is very tedious. Level set method [Set96, Set99] solves the problem elegantly by doing it in one higher dimension. Let Γ denote a closed curve in 2D. Define a new function φ(x(t), y(t), t) such that φ(x,y,t = 0) is the distance d of the point (x,y) from Γ. d is positive if (x,y) is outside Γ. d is zero if (x,y) is on Γ. d is negative if (x,y) is inside Γ. Leow Wee Kheng (CS6240) Non-Rigid Registration I 43 / 67
44 Level Set Method The 3D surface (red) φ(x, y, t) is the level set function. Γ at time t = 0 is the intersection of φ(x, y, t = 0) with the xy plane (blue), i.e., the solution of the equation φ(x, y, t = 0) = 0. This is called the zero level set. To compute Γ at a later time, e.g., t = 1, just move the level set function up (or down), and then compute the solution φ(x, y, t = 1) = 0. Leow Wee Kheng (CS6240) Non-Rigid Registration I 44 / 67
45 Level Set Method Now, the problem of propagating Γ over time becomes simpler. F: force that gives the speed of Γ in its normal direction. Then, the change of φ over time t, φ t, is given by the equations φ t + F φ = 0, (41) φ(x, y, t = 0) = given. (42) where φ = (φ 2 x + φ 2 y) 1/2 Leow Wee Kheng (CS6240) Non-Rigid Registration I 45 / 67
46 Level Set Method Notes: Although φ(x, y, t) remains a function, the zero level set φ = 0 corresponds to the propagating contour may change topology and form sharp corners. The idea of level set method is very simple. All the complications come from accurate and stable numerical algorithms for computing φ(x, y, t + t) from φ(x, y, t). See [Set96, Set99] for more details. Leow Wee Kheng (CS6240) Non-Rigid Registration I 46 / 67
47 Level Set Method Force Function Force F is usually formulated as F P + F G. F P is the propagation force: F P = αv(x) (43) α is a positive constant v(x) is a speed function F P controls propagation of surface. F P stops the propagation at the desired boundary. F G is a smoothing force: F G = βv(x)κ(x) (44) β is a positive constant κ(x) is the curvature at x F G smooths out high curvature regions. Leow Wee Kheng (CS6240) Non-Rigid Registration I 47 / 67
48 Level Set Method Possible speed functions Inverse proportion v(x) = (1 + I(x) ) 1 (45) Exponential Sigmoid v(x) = v(x) = exp( γ I(x) ) (46) [ ( 1 + exp I(x) b )] 1 (47) a Leow Wee Kheng (CS6240) Non-Rigid Registration I 48 / 67
49 Level Set Method Notes on numerical implementation: Discretize 3D volume (or 2D plane) into voxels. Compute initial value of level set function at each voxel. Update value at each voxel as the surface propagates. Computationally expensive for high-resolution discretization. Narrow Band Method [Set96, Set99] Update only the voxel values in a narrow band around the propagating surface. So, it s more efficient. Leow Wee Kheng (CS6240) Non-Rigid Registration I 49 / 67
50 Level Set Method Applicaton example: Segmentation of arteries in 2D images. Video demo Leow Wee Kheng (CS6240) Non-Rigid Registration I 50 / 67
51 Level Set Method Application example: Segmentation of great arteries in 3D CT images. 2-D views of segmentation results: Leow Wee Kheng (CS6240) Non-Rigid Registration I 51 / 67
52 Level Set Method 3-D view of segmentation results: Great arteries are segmented using 3D level set. Level set segmentation results have implicit surfaces. Explicit surfaces are obtained using marching cubes algorithm [LC87], followed by surface smoothing. Other video demos: brain, leg Leow Wee Kheng (CS6240) Non-Rigid Registration I 52 / 67
53 Level Set Method Further Readings Further Readings For details of level set method, refer to [Set96, Set99]. Level Set Method Library in ITK: [ITK]. Level Set Method Toolbox: [Mit]. Leow Wee Kheng (CS6240) Non-Rigid Registration I 53 / 67
54 Appendix A: Calculus of Variations Appendix A: Calculus of Variations Calculus of variations or variational calculus is a branch of mathematics that deals with the calculus of functional instead of function. A functional is a function of function, i.e., it maps a function to a real number. Calculus of variations is often used to find the extremum (minimum or maximum) of functional of the form E(v) = b a f(v, v, x)dx (48) where v = v(x) has two continuous derivatives, i.e., v and v exist. E(v) has an extremum only if the Euler-Lagrange equation is satisfied: f v d dx f v = 0. (49) Leow Wee Kheng (CS6240) Non-Rigid Registration I 54 / 67
55 Appendix A: Calculus of Variations Derivation Suppose v is changed to v + δv. Then, v is also changed to v + δv. Since the derivative of v + δv is equal to v + δv, we have d ds (v + δv) = v + δv (50) i.e., d ds δv = δv. (51) Assume that δv and δv are small, and δv = 0 at the boundary points where x = a and x = b. Then, applying Taylor s series expansion yields Therefore, f(v + δv, v + δv, x) = f(v, v, x) + f f δv + v v δv. (52) δe = b a ( ) f f δv + v v δv dx. (53) Leow Wee Kheng (CS6240) Non-Rigid Registration I 55 / 67
56 Appendix A: Calculus of Variations Applying integration-by-parts gives: b a ( ) f v δv dx = = b a = [ f v δv ( ) f d v ds δv dx b a ] b a b a ( d f dx v ( d f dx v ) δv dx ) δv dx (54) So, δe = b a ( f v d dx For δe to be zero for any δv, we need f v d dx f v ) δv dx. (55) f v = 0. (56) Leow Wee Kheng (CS6240) Non-Rigid Registration I 56 / 67
57 Appendix A: Calculus of Variations Application to Medial Axis Fitting Application to Medial Axis Fitting The medial axis is the v(s) that minimizes the cost function E(v): where E(v) = 1 2 and v = v(s), 0 s L. L 0 D(v) = v(s) ( D(v) + λ v 2) ds (57) p i χ s w i (s)p i p i χ s w i (s) 2, (58) Here, f(v,v, s) = 1 2 (D(v) + λ v 2 ). For E(v) to have a minimum, the following must be satisfied, f v d f ds v = 0. (59) Leow Wee Kheng (CS6240) Non-Rigid Registration I 57 / 67
58 Appendix A: Calculus of Variations Application to Medial Axis Fitting Evaluating each partial derivatives of f yields Therefore, f v = 1 D 2 v = 1 2 D f v = λv d f ds v = λv (60) 1 2 D λv = 0. (61) return Leow Wee Kheng (CS6240) Non-Rigid Registration I 58 / 67
59 Appendix B: Variational Calculus for Snake Model Appendix B: Variational Calculus for Snake Model The snake model is given by the equation E T (v) = 1 0 f(v,v,v, s)ds (62) where f(v,v,v, s) = (α v 2 + β v 2 )/2 + E(v), v = v(s), α = α(s), β = β(s), and 0 s 1. Following the derivation in Appendix A, we obtain the condition for E T (v) to have an extremum (Exercise): f v d f ds v + d2 f ds 2 v = 0. (63) Leow Wee Kheng (CS6240) Non-Rigid Registration I 59 / 67
60 Appendix B: Variational Calculus for Snake Model Evaluating each partial derivatives of f yields Therefore, f v = f v = αv E v = E d f ds v = (αv ) f v = βv d 2 ds 2 f v = (βv ). (64) (αv ) + (βv ) + E(v) = 0. (65) return Leow Wee Kheng (CS6240) Non-Rigid Registration I 60 / 67
61 Reference References I L. Balmelli, C. J. Morris, G. Taubin, and F. Bernardini. Volume warping for adaptive isosurface extraction. In Proc. IEEE Visualization, L. D. Cohen. On active contour models and balloons. CVGIP: Image Understanding, 53(2): , citeseer.nj.nec.com/cohen91active.html. X. Ee. Shape-Preserving Elastic Snakes. Honours Year Project Report, Dept. of Computer Science, NUS, Leow Wee Kheng (CS6240) Non-Rigid Registration I 61 / 67
62 Reference References II A. Foulonneau, P. Charbonnier, and F. Heitz. Geometric shape priors for region-based active contours. In Proc. Int. Conf. on Image Processing, S. R. Gunn and M. S. Nixon. A dual active contour for head boundary extraction. In Proc. Colloq. on Image Proc. for Biometric Measurement, pages 6/1 4, J. Hadamard. Lectures on the Cauchy Problem in Linear Partial Differential Equations. Yale University Press, Leow Wee Kheng (CS6240) Non-Rigid Registration I 62 / 67
63 Reference References III X. Han, C. Xu, and J. L. Prince. A topology preserving deformable model using level set. In Proc. IEEE Conf. on Computer Vision and Pattern Recognition, volume 2, pages , ITK: Insight Segmentation and Registration Toolkit, K. J. Kirchberg, A. Wimmer, and C. H. Lorenz. Modeling the human aorta for MR-driven real-time virtual endoscopy. In In Proc. of MICCAI, pages , M. Kass, A. Witkin, and D. Terzopoulos. Snakes: Active contour models. Int. Journal of Computer Vision, 1: , Leow Wee Kheng (CS6240) Non-Rigid Registration I 63 / 67
64 Reference References IV W. E. Lorensen and H. E. Cline. Marching cubes: A high-resolution 3D surface reconstruction algorithm. Computer Graphics, 21(4): , I. Mitchell. A toolbox of level set methods. J. A. Sethian. Examples of the application of level set methods. math.berkeley.edu/~sethian/movies/movieartery.html. J. A. Sethian. Level Set Methods. Cambridge University Press, Leow Wee Kheng (CS6240) Non-Rigid Registration I 64 / 67
65 Reference References V J. A. Sethian. Level Set Methods and Fast Marching Methods. Cambridge University Press, D. Shen, E. H. Herskovits, and C. Davatzikos. An adaptive-focus statistical shape model for segmentation and shape modeling of 3-D brain structures. IEEE Trans. on Medical Imaging, 20(4), A. N. Tikhonov and V. A. Arsenin. Solutions of Ill-Posed Problems. Winston and Sons, Washington, D.C., Leow Wee Kheng (CS6240) Non-Rigid Registration I 65 / 67
66 Reference References VI Tai Peng Tian. Detection of femur fractures in x-ray images. Master s thesis, Dept. Computer Science, National University of Singapore, S. Wang, J. X. Ji, and Z.-P. Liang. Landmark-based shape deformation with topology-preserving constraints. In Proc. Int. Conf. on Computer Vision, G. Wolberg. Digital Image Warping. IEEE Computer Society Press, Leow Wee Kheng (CS6240) Non-Rigid Registration I 66 / 67
67 Reference References VII C. Xu and J. L. Prince. Gradient vector flow: A new external force for snakes. In Proc. IEEE Conf. on Computer Vision and Pattern Recognition, C. Xu and J. L. Prince. Snakes, shapes, and gradient vector flow. IEEE Trans. on Image Processing, 7(3): , C. Xu and J. L. Prince. Gradient vector flow deformable models. In I. Bankman, editor, Handbook of Medical Imaging. Academic Press, Leow Wee Kheng (CS6240) Non-Rigid Registration I 67 / 67
Non-Rigid Image Registration III
Non-Rigid Image Registration III CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (CS6240) Non-Rigid Image Registration
More informationAutomatic Extraction of Femur Contours from Hip X-ray Images
Automatic Extraction of Femur Contours from Hip X-ray Images Ying Chen 1, Xianhe Ee 1, Wee Kheng Leow 1, and Tet Sen Howe 2 1 Dept of Computer Science, National University of Singapore, 3 Science Drive
More informationImage Segmentation II Advanced Approaches
Image Segmentation II Advanced Approaches Jorge Jara W. 1,2 1 Department of Computer Science DCC, U. of Chile 2 SCIAN-Lab, BNI Outline 1. Segmentation I Digital image processing Segmentation basics 2.
More informationSnakes operating on Gradient Vector Flow
Snakes operating on Gradient Vector Flow Seminar: Image Segmentation SS 2007 Hui Sheng 1 Outline Introduction Snakes Gradient Vector Flow Implementation Conclusion 2 Introduction Snakes enable us to find
More informationModel Fitting. CS6240 Multimedia Analysis. Leow Wee Kheng. Department of Computer Science School of Computing National University of Singapore
Model Fitting CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore (CS6240) Model Fitting 1 / 34 Introduction Introduction Model
More informationProblem Solving Assignment 1
CS6240 Problem Solving Assignment 1 p. 1/20 Problem Solving Assignment 1 CS6240 Multimedia Analysis Daniel Dahlmeier National University of Singapore CS6240 Problem Solving Assignment 1 p. 2/20 Introduction
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationAccurate Reconstruction by Interpolation
Accurate Reconstruction by Interpolation Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore International Conference on Inverse Problems and Related Topics
More information1. Two double lectures about deformable contours. 4. The transparencies define the exam requirements. 1. Matlab demonstration
Practical information INF 5300 Deformable contours, I An introduction 1. Two double lectures about deformable contours. 2. The lectures are based on articles, references will be given during the course.
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 informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationSnakes reparameterization for noisy images segmentation and targets tracking
Snakes reparameterization for noisy images segmentation and targets tracking Idrissi Sidi Yassine, Samir Belfkih. Lycée Tawfik Elhakim Zawiya de Noaceur, route de Marrakech, Casablanca, maroc. Laboratoire
More informationEdge Detection and Active Contours
Edge Detection and Active Contours Elsa Angelini, Florence Tupin Department TSI, Telecom ParisTech Name.surname@telecom-paristech.fr 2011 Outline Introduction Edge Detection Active Contours Introduction
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 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 informationSEGMENTATION AND RECONSTRUCTION OF 3D ARTERY MODELS FOR SURGICAL PLANNING
SEGMENTATION AND RECONSTRUCTION OF 3D ARTERY MODELS FOR SURGICAL PLANNING QI YINGYI (B.Sc., FUDAN UNIVERSITY, 2005 ) A THESIS SUBMITTED FOR THE DEGREE OF MASTER OF SCIENCE DEPARTMENT OF COMPUTER SCIENCE
More informationGradient Vector Flow: A New External Force for Snakes
66 IEEE Proc. Conf. on Comp. Vis. Patt. Recog. (CVPR'97) Gradient Vector Flow: A New External Force for Snakes Chenyang Xu and Jerry L. Prince Department of Electrical and Computer Engineering The Johns
More informationA Method of Automated Landmark Generation for Automated 3D PDM Construction
A Method of Automated Landmark Generation for Automated 3D PDM Construction A. D. Brett and C. J. Taylor Department of Medical Biophysics University of Manchester Manchester M13 9PT, Uk adb@sv1.smb.man.ac.uk
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationDr. Ulas Bagci
Lecture 9: Deformable Models and Segmentation CAP-Computer Vision Lecture 9-Deformable Models and Segmentation Dr. Ulas Bagci bagci@ucf.edu Lecture 9: Deformable Models and Segmentation Motivation A limitation
More informationA Finite Element Method for Deformable Models
A Finite Element Method for Deformable Models Persephoni Karaolani, G.D. Sullivan, K.D. Baker & M.J. Baines Intelligent Systems Group, Department of Computer Science University of Reading, RG6 2AX, UK,
More informationMulti-Scale Free-Form Surface Description
Multi-Scale Free-Form Surface Description Farzin Mokhtarian, Nasser Khalili and Peter Yuen Centre for Vision Speech and Signal Processing Dept. of Electronic and Electrical Engineering University of Surrey,
More informationAutomated Segmentation Using a Fast Implementation of the Chan-Vese Models
Automated Segmentation Using a Fast Implementation of the Chan-Vese Models Huan Xu, and Xiao-Feng Wang,,3 Intelligent Computation Lab, Hefei Institute of Intelligent Machines, Chinese Academy of Science,
More informationLeow Wee Kheng CS4243 Computer Vision and Pattern Recognition. Motion Tracking. CS4243 Motion Tracking 1
Leow Wee Kheng CS4243 Computer Vision and Pattern Recognition Motion Tracking CS4243 Motion Tracking 1 Changes are everywhere! CS4243 Motion Tracking 2 Illumination change CS4243 Motion Tracking 3 Shape
More informationMedical Image Segmentation using Level Sets
Medical Image Segmentation using Level Sets Technical Report #CS-8-1 Tenn Francis Chen Abstract Segmentation is a vital aspect of medical imaging. It aids in the visualization of medical data and diagnostics
More informationParticle Filtering. CS6240 Multimedia Analysis. Leow Wee Kheng. Department of Computer Science School of Computing National University of Singapore
Particle Filtering CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore (CS6240) Particle Filtering 1 / 28 Introduction Introduction
More informationSegmentation of Volume Images Using Trilinear Interpolation of 2D Slices
Segmentation of Volume Images Using Trilinear Interpolation of 2D Slices Pouyan TaghipourBibalan,Mehdi Mohammad Amini Australian National University ptbibalan@rsise.anu.edu.au,u4422717@anu.edu.au Abstract
More informationChromosome Segmentation and Investigations using Generalized Gradient Vector Flow Active Contours
Published Quarterly Mangalore, South India ISSN 0972-5997 Volume 4, Issue 2; April-June 2005 Original Article Chromosome Segmentation and Investigations using Generalized Gradient Vector Flow Active Contours
More informationCOMPUTER VISION > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE
COMPUTER VISION 2017-2018 > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE OUTLINE Optical flow Lucas-Kanade Horn-Schunck Applications of optical flow Optical flow tracking Histograms of oriented flow Assignment
More informationPattern Recognition Methods for Object Boundary Detection
Pattern Recognition Methods for Object Boundary Detection Arnaldo J. Abrantesy and Jorge S. Marquesz yelectronic and Communication Engineering Instituto Superior Eng. de Lisboa R. Conselheiro Emídio Navarror
More informationMotion. 1 Introduction. 2 Optical Flow. Sohaib A Khan. 2.1 Brightness Constancy Equation
Motion Sohaib A Khan 1 Introduction So far, we have dealing with single images of a static scene taken by a fixed camera. Here we will deal with sequence of images taken at different time intervals. Motion
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationUnit 3 : Image Segmentation
Unit 3 : Image Segmentation K-means Clustering Mean Shift Segmentation Active Contour Models Snakes Normalized Cut Segmentation CS 6550 0 Histogram-based segmentation Goal Break the image into K regions
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 10 130221 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Canny Edge Detector Hough Transform Feature-Based
More informationcoding of various parts showing different features, the possibility of rotation or of hiding covering parts of the object's surface to gain an insight
Three-Dimensional Object Reconstruction from Layered Spatial Data Michael Dangl and Robert Sablatnig Vienna University of Technology, Institute of Computer Aided Automation, Pattern Recognition and Image
More informationMultimodality Imaging for Tumor Volume Definition in Radiation Oncology
81 There are several commercial and academic software tools that support different segmentation algorithms. In general, commercial software packages have better implementation (with a user-friendly interface
More informationMachine Learning for Signal Processing Lecture 4: Optimization
Machine Learning for Signal Processing Lecture 4: Optimization 13 Sep 2015 Instructor: Bhiksha Raj (slides largely by Najim Dehak, JHU) 11-755/18-797 1 Index 1. The problem of optimization 2. Direct optimization
More informationRegistration-Based Segmentation of Medical Images
School of Computing National University of Singapore Graduate Research Paper Registration-Based Segmentation of Medical Images by Li Hao under guidance of A/Prof. Leow Wee Kheng July, 2006 Abstract Medical
More informationMATH 2400, Analytic Geometry and Calculus 3
MATH 2400, Analytic Geometry and Calculus 3 List of important Definitions and Theorems 1 Foundations Definition 1. By a function f one understands a mathematical object consisting of (i) a set X, called
More informationDynamic 3D Shape From Multi-viewpoint Images Using Deformable Mesh Model
Dynamic 3D Shape From Multi-viewpoint Images Using Deformable Mesh Model Shohei Nobuhara Takashi Matsuyama Graduate School of Informatics, Kyoto University Sakyo, Kyoto, 606-8501, Japan {nob, tm}@vision.kuee.kyoto-u.ac.jp
More informationActive contours in Brain tumor segmentation
Active contours in Brain tumor segmentation Ali Elyasi* 1, Mehdi Hosseini 2, Marzieh Esfanyari 2 1. Department of electronic Engineering, Young Researchers Club, Central Tehran Branch, Islamic Azad University,
More informationLecture 12 Level Sets & Parametric Transforms. sec & ch. 11 of Machine Vision by Wesley E. Snyder & Hairong Qi
Lecture 12 Level Sets & Parametric Transforms sec. 8.5.2 & ch. 11 of Machine Vision by Wesley E. Snyder & Hairong Qi Spring 2017 16-725 (CMU RI) : BioE 2630 (Pitt) Dr. John Galeotti The content of these
More informationMidterm Exam Solutions
Midterm Exam Solutions Computer Vision (J. Košecká) October 27, 2009 HONOR SYSTEM: This examination is strictly individual. You are not allowed to talk, discuss, exchange solutions, etc., with other fellow
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationCurvature Estimation on Smoothed 3-D Meshes
Curvature Estimation on Smoothed 3-D Meshes Peter Yuen, Nasser Khalili and Farzin Mokhtarian Centre for Vision, Speech and Signal Processing School of Electronic Engineering, Information Technology and
More informationPoint-Based Geometric Deformable Models for Medical Image Segmentation
Point-Based Geometric Deformable Models for Medical Image Segmentation Hon Pong Ho 1, Yunmei Chen 2, Huafeng Liu 1,3, and Pengcheng Shi 1 1 Dept. of EEE, Hong Kong University of Science & Technology, Hong
More informationProject Updates Short lecture Volumetric Modeling +2 papers
Volumetric Modeling Schedule (tentative) Feb 20 Feb 27 Mar 5 Introduction Lecture: Geometry, Camera Model, Calibration Lecture: Features, Tracking/Matching Mar 12 Mar 19 Mar 26 Apr 2 Apr 9 Apr 16 Apr 23
More informationGeometric Modeling and Processing
Geometric Modeling and Processing Tutorial of 3DIM&PVT 2011 (Hangzhou, China) May 16, 2011 4. Geometric Registration 4.1 Rigid Registration Range Scanning: Reconstruction Set of raw scans Reconstructed
More informationMean Shift Tracking. CS4243 Computer Vision and Pattern Recognition. Leow Wee Kheng
CS4243 Computer Vision and Pattern Recognition Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore (CS4243) Mean Shift Tracking 1 / 28 Mean Shift Mean Shift
More informationOn 3D Shape Synthesis
On 3D Shape Synthesis Heung-yeung Shum Martial Hebert Katsushi Ikeuchi November 1995 CMU-CS-95-213 School of Computer Science Carnegie Mellon University Pittsburgh, Pennsylvania 15213 c 1995 Carnegie Mellon
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 information03 - Reconstruction. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Spring 17 - Daniele Panozzo
3 - Reconstruction Acknowledgements: Olga Sorkine-Hornung Geometry Acquisition Pipeline Scanning: results in range images Registration: bring all range images to one coordinate system Stitching/ reconstruction:
More informationResearch on the Wood Cell Contour Extraction Method Based on Image Texture and Gray-scale Information.
, pp. 65-74 http://dx.doi.org/0.457/ijsip.04.7.6.06 esearch on the Wood Cell Contour Extraction Method Based on Image Texture and Gray-scale Information Zhao Lei, Wang Jianhua and Li Xiaofeng 3 Heilongjiang
More informationCSE 554 Lecture 6: Fairing and Simplification
CSE 554 Lecture 6: Fairing and Simplification Fall 2012 CSE554 Fairing and simplification Slide 1 Review Iso-contours in grayscale images and volumes Piece-wise linear representations Polylines (2D) and
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More informationTheoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization. M2DO Lab 1,2. 1 Cardiff University
Theoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization M2DO Lab 1,2 1 Cardiff University 2 University of California, San Diego November 2017 A brief description of theory
More informationShape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow
Shape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow Abstract. Finding meaningful 1-1 correspondences between hippocampal (HP) surfaces is an important but difficult
More informationImplicit Active Shape Models for 3D Segmentation in MR Imaging
Implicit Active Shape Models for 3D Segmentation in MR Imaging Mikaël Rousson 1, Nikos Paragios 2, and Rachid Deriche 1 1 I.N.R.I.A. Sophia Antipolis, France E-mail: {Mikael.Rousson,Rachid.Deriche}@sophia.inria.fr
More informationEECS 556 Image Processing W 09. Interpolation. Interpolation techniques B splines
EECS 556 Image Processing W 09 Interpolation Interpolation techniques B splines What is image processing? Image processing is the application of 2D signal processing methods to images Image representation
More informationComputer Vision. Recap: Smoothing with a Gaussian. Recap: Effect of σ on derivatives. Computer Science Tripos Part II. Dr Christopher Town
Recap: Smoothing with a Gaussian Computer Vision Computer Science Tripos Part II Dr Christopher Town Recall: parameter σ is the scale / width / spread of the Gaussian kernel, and controls the amount of
More informationDeformable Contour Method: A Constrained Optimization Approach
Deformable Contour Method: A Constrained Optimization Approach Xun Wang, Lei He, Chia Y. Han, William G. Wee Electrical & Computer Engineering and Computer Science Department University of Cincinnati,
More informationCHAPTER 3 Image Segmentation Using Deformable Models
CHAPTER Image Segmentation Using Deformable Models Chenyang Xu The Johns Hopkins University Dzung L. Pham National Institute of Aging Jerry L. Prince The Johns Hopkins University Contents.1 Introduction
More informationLecture 7: Most Common Edge Detectors
#1 Lecture 7: Most Common Edge Detectors Saad Bedros sbedros@umn.edu Edge Detection Goal: Identify sudden changes (discontinuities) in an image Intuitively, most semantic and shape information from the
More informationSegmentation. Namrata Vaswani,
Segmentation Namrata Vaswani, namrata@iastate.edu Read Sections 5.1,5.2,5.3 of [1] Edge detection and filtering : Canny edge detection algorithm to get a contour of the object boundary Hough transform:
More informationLecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More informationOutline. Reconstruction of 3D Meshes from Point Clouds. Motivation. Problem Statement. Applications. Challenges
Reconstruction of 3D Meshes from Point Clouds Ming Zhang Patrick Min cs598b, Geometric Modeling for Computer Graphics Feb. 17, 2000 Outline - problem statement - motivation - applications - challenges
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 informationGeometrical Modeling of the Heart
Geometrical Modeling of the Heart Olivier Rousseau University of Ottawa The Project Goal: Creation of a precise geometrical model of the heart Applications: Numerical calculations Dynamic of the blood
More informationA Modified Image Segmentation Method Using Active Contour Model
nd International Conference on Electrical, Computer Engineering and Electronics (ICECEE 015) A Modified Image Segmentation Method Using Active Contour Model Shiping Zhu 1, a, Ruidong Gao 1, b 1 Department
More informationInverse and Implicit functions
CHAPTER 3 Inverse and Implicit functions. Inverse Functions and Coordinate Changes Let U R d be a domain. Theorem. (Inverse function theorem). If ϕ : U R d is differentiable at a and Dϕ a is invertible,
More informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More information1. Two double lectures about deformable contours. 4. The transparencies define the exam requirements. 1. Matlab demonstration
Practical information INF 5300 Deformable contours, II An introduction 1. Two double lectures about deformable contours. 2. The lectures are based on articles, references will be given during the course.
More informationImplicit Active Contours Driven by Local Binary Fitting Energy
Implicit Active Contours Driven by Local Binary Fitting Energy Chunming Li 1, Chiu-Yen Kao 2, John C. Gore 1, and Zhaohua Ding 1 1 Institute of Imaging Science 2 Department of Mathematics Vanderbilt University
More informationMatching. Compare region of image to region of image. Today, simplest kind of matching. Intensities similar.
Matching Compare region of image to region of image. We talked about this for stereo. Important for motion. Epipolar constraint unknown. But motion small. Recognition Find object in image. Recognize object.
More informationChapter 3. Automated Segmentation of the First Mitotic Spindle in Differential Interference Contrast Microcopy Images of C.
Chapter 3 Automated Segmentation of the First Mitotic Spindle in Differential Interference Contrast Microcopy Images of C. elegans Embryos Abstract Differential interference contrast (DIC) microscopy is
More informationNonrigid Surface Modelling. and Fast Recovery. Department of Computer Science and Engineering. Committee: Prof. Leo J. Jia and Prof. K. H.
Nonrigid Surface Modelling and Fast Recovery Zhu Jianke Supervisor: Prof. Michael R. Lyu Committee: Prof. Leo J. Jia and Prof. K. H. Wong Department of Computer Science and Engineering May 11, 2007 1 2
More informationChemnitz Scientific Computing Preprints
Roman Unger Obstacle Description with Radial Basis Functions for Contact Problems in Elasticity CSC/09-01 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN
More informationGeometric Transformations and Image Warping
Geometric Transformations and Image Warping Ross Whitaker SCI Institute, School of Computing University of Utah Univ of Utah, CS6640 2009 1 Geometric Transformations Greyscale transformations -> operate
More informationSubpixel Corner Detection Using Spatial Moment 1)
Vol.31, No.5 ACTA AUTOMATICA SINICA September, 25 Subpixel Corner Detection Using Spatial Moment 1) WANG She-Yang SONG Shen-Min QIANG Wen-Yi CHEN Xing-Lin (Department of Control Engineering, Harbin Institute
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 informationENGINEERING TRIPOS PART IIA FINITE ELEMENT METHOD
ENGINEERING TRIPOS PART IIA LOCATION: DPO EXPERIMENT 3D7 FINITE ELEMENT METHOD Those who have performed the 3C7 experiment should bring the write-up along to this laboratory Objectives Show that the accuracy
More informationGVF-Based Transfer Functions for Volume Rendering
GVF-Based Transfer Functions for Volume Rendering Shaorong Wang 1,2 and Hua Li 1 1 National Research Center for Intelligent Computing Systems, Institute of Computing Technology, Chinese Academy of Sciences,
More informationIsometries. 1 Identifying Isometries
Isometries 1 Identifying Isometries 1. Modeling isometries as dynamic maps. 2. GeoGebra files: isoguess1.ggb, isoguess2.ggb, isoguess3.ggb, isoguess4.ggb. 3. Guessing isometries. 4. What can you construct
More informationGlobal Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising
Global Minimization of the Active Contour Model with TV-Inpainting and Two-Phase Denoising Shingyu Leung and Stanley Osher Department of Mathematics, UCLA, Los Angeles, CA 90095, USA {syleung, sjo}@math.ucla.edu
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationSnakes, level sets and graphcuts. (Deformable models)
INSTITUTE OF INFORMATION AND COMMUNICATION TECHNOLOGIES BULGARIAN ACADEMY OF SCIENCE Snakes, level sets and graphcuts (Deformable models) Centro de Visión por Computador, Departament de Matemàtica Aplicada
More informationSegmentation of 3D CT Volume Images Using a Single 2D Atlas
Segmentation of 3D CT Volume Images Using a Single 2D Atlas Feng Ding 1, Wee Kheng Leow 1, and Shih-Chang Wang 2 1 Dept. of Computer Science, National University of Singapore, 3 Science Drive 2, Singapore
More informationExtract Object Boundaries in Noisy Images using Level Set. Literature Survey
Extract Object Boundaries in Noisy Images using Level Set by: Quming Zhou Literature Survey Submitted to Professor Brian Evans EE381K Multidimensional Digital Signal Processing March 15, 003 Abstract Finding
More informationSTATISTICS AND ANALYSIS OF SHAPE
Control and Cybernetics vol. 36 (2007) No. 2 Book review: STATISTICS AND ANALYSIS OF SHAPE by H. Krim, A. Yezzi, Jr., eds. There are numerous definitions of a notion of shape of an object. These definitions
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More informationImage Warping. Srikumar Ramalingam School of Computing University of Utah. [Slides borrowed from Ross Whitaker] 1
Image Warping Srikumar Ramalingam School of Computing University of Utah [Slides borrowed from Ross Whitaker] 1 Geom Trans: Distortion From Optics Barrel Distortion Pincushion Distortion Straight lines
More informationBoundary descriptors. Representation REPRESENTATION & DESCRIPTION. Descriptors. Moore boundary tracking
Representation REPRESENTATION & DESCRIPTION After image segmentation the resulting collection of regions is usually represented and described in a form suitable for higher level processing. Most important
More informationREVIEW I MATH 254 Calculus IV. Exam I (Friday, April 29) will cover sections
REVIEW I MATH 254 Calculus IV Exam I (Friday, April 29 will cover sections 14.1-8. 1. Functions of multivariables The definition of multivariable functions is similar to that of functions of one variable.
More informationChapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationLarge Mesh Deformation Using the Volumetric Graph Laplacian
Large Mesh Deformation Using the Volumetric Graph Laplacian Kun Zhou1 Jin Huang2 John Snyder3 Xinguo Liu1 Hujun Bao2 Baining Guo1 Heung-Yeung Shum1 1 Microsoft Research Asia 2 Zhejiang University 3 Microsoft
More informationTopic 6 Representation and Description
Topic 6 Representation and Description Background Segmentation divides the image into regions Each region should be represented and described in a form suitable for further processing/decision-making Representation
More informationCerebral Artery Segmentation with Level Set Methods
H. Ho, P. Bier, G. Sands, P. Hunter, Cerebral Artery Segmentation with Level Set Methods, Proceedings of Image and Vision Computing New Zealand 2007, pp. 300 304, Hamilton, New Zealand, December 2007.
More informationLecture 10 Segmentation, Part II (ch 8) Active Contours (Snakes) ch. 8 of Machine Vision by Wesley E. Snyder & Hairong Qi
Lecture 10 Segmentation, Part II (ch 8) Active Contours (Snakes) ch. 8 of Machine Vision by Wesley E. Snyder & Hairong Qi Spring 2018 16-725 (CMU RI) : BioE 2630 (Pitt) Dr. John Galeotti The content of
More information