# Level set methods for imaging and application to MRI segmentation - PowerPoint PPT Presentation  Download Presentation Level set methods for imaging and application to MRI segmentation

Level set methods for imaging and application to MRI segmentation
Download Presentation ## Level set methods for imaging and application to MRI segmentation

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
##### Presentation Transcript

1. MRI Segmentation Acknowledgements • Based on results by Denis Neiter (Ecole Polytechnique) during internship 2007, partly using results by Simon Huffeteau (Ecole Polytechnique), internship 2006 • Using Carsten Wolters‘ MR Data

2. MRI Segmentation Problem Setting • Given MR Image(s), find (in an automated way): • the borders between different head compartments (segmentation) • an appropriate map of the normal directions, in particular of the brain surface (classification)- a representation useful for further finite element modelling

3. MRI Segmentation Mathematical Issues • Segmentation needs • to discriminate noise and textures (small scale structures) • to incorporate prior knowledge • - to be flexible with respect to complicated shapes (or even topology) • First two issues treated via regularization, third via level set methods

4. MRI Segmentation Object-Based Segmentation • Classical object based segmentation computes curves (2D) or surfaces (3D) marking the object boundary (contour) • Traditional approach: start with curve and let it evolve towards the contour by some criteria • Velocity of evolving curve determined by two counteracting parts: image-driven part and regularization

5. MRI Segmentation Active Contours - Snakes • Image-driven force related to gradient of the image (local gray-value difference) • Regularization force is (mean) curvature

6. MRI Segmentation Active Contours - Snakes • Popular, but various shortcomings: • needs preprocessing of the image (noise removal, intensity map so that edges are in valleys) • local minima • issues with narrow structures: big trouble in brain images

7. MRI Segmentation Statistical Models • K-Means / C-Means: • Based on optimization, find a 0-1 function (0 in pixels outside, 1 inside) • Optimization goal consists of same parts: image-driven and regularization • No useful boundary representation

8. MRI Segmentation Curve / Surface Representation • Level Set Methods yield boundary representation with appropriate curvatures and subpixel resolution

9. MRI Segmentation Level Set Methods • Osher & Sethian, JCP 1987, Sethian, Cambridge Univ. Press 1999, Osher & Fedkiw, Springer, 2002 • Basic idea: implicit shape representation • with continuous level-set function

10. MRI Segmentation Level Set Methods • Change of front translated to change of function

11. MRI Segmentation Level Set Methods • Implicit representation of dynamic shapes • with time-dependent level set function

12. MRI Segmentation Level Set Methods • Evolution of the shape corresponds to evolution of the level set function (and vice versa) • Movie by • J.Sethian

13. MRI Segmentation Level Set Methods • Topology change is automatic • Movie by • J.Sethian

14. MRI Segmentation Geometric Motion • Start for simplicity with the evolution of a curve • Evolution in a velocity field , each point evolves via ODE

15. MRI Segmentation Geometric Motion • Use any parametric representation • Due to definition of the level set function • Consequently

16. MRI Segmentation Geometric Motion • By the chain rule • Insert ODE for moving points:

17. MRI Segmentation Geometric Motion • For level set function being a solution of • each level set of f is moving with velocity V

18. MRI Segmentation Geometric Motion • In most cases, the full velocity field V is unknown, only normal velocity component v known • Tangential component of the velocity field is not important anyway, it does not change the motion (only change of parametrization)

19. MRI Segmentation Geometric Motion • Normal can be computed from level set function: • By the chain rule

20. MRI Segmentation Geometric Motion • Note that is a tangential direction • Hence, is a normal direction, unit normal is given by

21. MRI Segmentation Geometric Motion • Evolution becomes nonlinear Hamilton-Jacobi equation: • „Level set equation“

22. MRI Segmentation Geometric Motion • Evolution could be anisotropic, i.e. normal velocity depends on the orientation • with one-homogeneous extension H , yields Hamilton-Jacobi equation

23. MRI Segmentation Geometric Motion • Evolution could be of higher order, e.g. normal velocity depends on the mean curvature • Level set equation becomes fully nonlinear second-order parabolic PDE

24. MRI Segmentation Examples • Eikonal equation • Positive velocity field yield monotone advancement of fronts • Arrival time • Solves

25. MRI Segmentation Example: Eikonal Equation

26. MRI Segmentation Examples • Mean curvature flow • Classical example of higher-order geometric motion • Normal velocity equal to curvature of curve (or mean curvature of surface)

27. MRI Segmentation Mean Curvature Flow

28. MRI Segmentation Optimal Geometries • Classical problem for optimal geometry: • Plateau Problem (Minimal Surface Problem) • Minimize area of surface between fixed boundary curves.

29. MRI Segmentation Optimal Geometries • Minimal surface (L.T.Cheng, PhD 2002)

30. MRI Segmentation Optimal Geometries • Wulff-Shapes: Pb in Cu • Surnev et al, J.Vacuum Sci. Tech. A, 1998

31. MRI Segmentation Mumford-Shah • Free discontinuity problems: • find the set of discontinuity from a noisy observation of a function. • Mumford-Shah functional

32. MRI Segmentation Mumford-Shah • Image decomposition

33. MRI Segmentation Mumford-Shah • Limitations

34. MRI Segmentation Improved Model • Decomposition in • 3 parts: smooth, oscillating, edges

35. MRI Segmentation Object-based Mumford-Shah • Chan-Vese: • Approximate smooth component by its mean value inside and outside object • Curve / Surface can be evolved via simple criterion, in each time step mean-value inside and outside need to be computed • Via convex relaxation techniques convergence to global minimum can be ensured

36. MRI Segmentation Level Set Formulation • Level set function y and • Heaviside function H (= 0 negative, = 1 positive)

37. MRI Segmentation Reduced Problem: fixed mean value

38. MRI Segmentation Image Segmentation • Noisy Image

39. MRI Segmentation Image Segmentation • Noise level 10%, l=103

40. MRI Segmentation Regularization • For skull segmentation (smooth) regularization based on length minimization is perfect • For brain structure (sulci) similar issues as for active contours

41. MRI Segmentation MR Results

42. MRI Segmentation MR Results

43. MRI Segmentation Skull Segmentation from MR-PD