Top Banner
1 Medical Imaging, SS-2014 Dr. Mohammad Dawood Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany
33

Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

Apr 01, 2015

Download

Documents

Rowan Biles
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

Medical Image Analysis

Dr. Mohammad Dawood

Department of Computer Science

University of MünsterGermany

Page 2: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

2

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Recap

Page 3: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

3

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Grayscale transformations1. Linear2. Logarithmic3. Power law

Point operations

Local operators

4. Histogram Equalization5. Adpative/Local Hist Eq6. Color space7. Fourier transform8. Spatial filtering

3 3 3 3 0

3 5 3 3 0

3 3 3 3 0

0 0 0 0 0

0 0 0 0 0

1 1 1

1 1 1

1 1 1

Page 4: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

4

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Edge detection

Page 5: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

5

Medical Imaging, SS-2014

Dr. Mohammad Dawood

What is an “edge”?Discontinuity in Image brightness

Page 6: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

6

Medical Imaging, SS-2014

Dr. Mohammad Dawood

15 15 15 15 15

15 15 15 15 15

15 15 15 15 15

0 0 0 0 0

0 0 0 0 0

1

-1

0 0 0 0 0

0 0 0 0 0

15 15 15 15 15

0 0 0 0 0

0 0 0 0 0

* =

Recognizing the edge

Page 7: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

7

Medical Imaging, SS-2014

Dr. Mohammad Dawood

15 15 15 15 15

15 15 15 15 15

15 15 15 15 15

0 0 0 0 0

0 0 0 0 0

1

0

-1

0 0 0 0 0

0 0 0 0 0

15 15 15 15 15

15 15 15 15 15

0 0 0 0 0

* =

Increasing edge thickness- easier to detect and better connected edges

Page 8: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

8

Medical Imaging, SS-2014

Dr. Mohammad Dawood

1 1 1

0 0 0

-1 -1 -1* =

15 15 15 15 15

15 15 15 15 15

15 15 15 15 15

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

45 45 45 45 45

45 45 45 45 45

0 0 0 0 0

Strengthening the edges

Page 9: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

9

Medical Imaging, SS-2014

Dr. Mohammad Dawood

1 0 -1

1 0 -1

1 0 -1

1 1 1

0 0 0

-1 -1 -1

Edge detection with spatial operators

Prewitt operators

Page 10: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

10

Medical Imaging, SS-2014

Dr. Mohammad Dawood

2 1 0

1 0 -1

0 -1 -2

Adding operators

1 0 -1

1 0 -1

1 0 -1

1 1 1

0 0 0

-1 -1 -1+ =

Page 11: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

11

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Derivatives of an image

-1 1

1 -2 1

Magnitude of gradient:

Angle:

Page 12: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

12

Medical Imaging, SS-2014

Dr. Mohammad Dawood

-1 1

First derivative

Forward difference

Backward difference

Central difference

1 -1

-0.5 0 0.5

MRI Spine fw bw cd bw_i bw+bw_i

Page 13: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

13

Medical Imaging, SS-2014

Dr. Mohammad Dawood

0 1 0

1 -4 1

0 1 0

Laplace operator

H+V Laplace

Page 14: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

14

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Cardiac PET

Page 15: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

15

Medical Imaging, SS-2014

Dr. Mohammad Dawood

1 2 1

0 0 0

-1 -2 -1

15 15 15 15 15

15 15 15 15 15

15 15 15 15 15

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

60 60 60 60 60

60 60 60 60 60

0 0 0 0 0

Gaussian+Gradient

* =

Page 16: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

16

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Sobel operators

1 0 -1

2 0 -2

1 0 -1

1 2 1

0 0 0

-1 -2 -1

Edge detection with spatial operators

Page 17: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

17

Medical Imaging, SS-2014

Dr. Mohammad Dawood

2 2 0

2 0 -2

0 -2 -2

1 0 -1

2 0 -2

1 0 -1

1 2 1

0 0 0

-1 -2 -1+ =

Page 18: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

18

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Scharr operators

3 0 -3

10 0 -10

3 0 -3

3 10 3

0 0 0

-3 -10 -3

Edge detection with spatial operators

Page 19: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

19

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Roberts operators0 1

-1 0

1 0

0 -1

Edge detection with spatial operators

+

Page 20: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

20

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Canny operator

1. Gaussian for noise reduction

2. Calculation of edges (sobel operator)

3. Non-maximum suppression, no neighbor should have a higher gradient except in the same direction

0 : if intensity > the intensities in the N and S directions45 : if intensity > the intensities in the NW and SE directions90 : if intensity > the intensities in the W and E directions135 : if intensity > the intensities in the NE and SW directions

4. Hysteresisdelete edges below threshold 1keep edges above threshold 2keep edges between thresholds, if one neighbor is above threshold 2

Page 21: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

21

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Canny operator th=0.5 th=0.1

Page 22: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

22

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Marr-Hildreth operator

Laplacian of the Gaussian (LoG)

Page 23: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

23

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Marr Hildreth operator sigma=1 sigma=2

Page 24: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

24

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Hough Transform

Page 25: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

25

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Hough transform for detecting lines

A line can be defined as:

Take the edge map of the image I

Look for the neighbors of a pixel and determine m and b

Accumulate the m and b in an accumulator array

Find the maxima of the accumulator array

Transform them back to image space

Page 26: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

26

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Hough transform for detecting lines

Alternative definition of lines

Page 27: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

27

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Hough transform

Similar transforms can be defined for circles, ellipses or other parametric curves

Page 28: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

28

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Morphological operations

Page 29: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

29

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Morphological operators

Operations are based on Set Theory and require a structure element

Basic morphological operations are:1. Erosion2. Dilation3. Opening4. Closing

Page 30: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

30

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Erosion

If A is an image and B is a structure element then

0 0 0 0 0

0 1 1 1 0

0 0 1 1 0

0 0 1 1 0

0 0 0 1 0

0 0 0

0 1 1

0 1 0

0 0 0 0 0

0 0 1 0 0

0 0 1 0 0

0 0 0 0 0

0 0 0 0 0

X

Page 31: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

31

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Dilation

0 0 0 0 0

0 1 1 1 0

0 0 1 1 0

0 0 1 1 0

0 0 0 1 0

0 0 0

0 1 1

0 1 0

0 1 1 1 0

1 1 1 1 0

0 1 1 1 0

0 1 1 1 0

0 0 1 1 0

X

Page 32: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

32

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Closing

Dilation + Erosion

Page 33: Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany.

33

Medical Imaging, SS-2014

Dr. Mohammad Dawood

Opening

Erosion + Dilation