medical image analysis dr. mohammad dawood department of computer science university of münster...
TRANSCRIPT
Medical Image Analysis
Dr. Mohammad Dawood
Department of Computer Science
University of MünsterGermany
2
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Recap
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
4
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Edge detection
5
Medical Imaging, SS-2014
Dr. Mohammad Dawood
What is an “edge”?Discontinuity in Image brightness
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
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
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
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
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+ =
11
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Derivatives of an image
-1 1
1 -2 1
Magnitude of gradient:
Angle:
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
13
Medical Imaging, SS-2014
Dr. Mohammad Dawood
0 1 0
1 -4 1
0 1 0
Laplace operator
H+V Laplace
14
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Cardiac PET
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
* =
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
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+ =
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
19
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Roberts operators0 1
-1 0
1 0
0 -1
Edge detection with spatial operators
+
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
21
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Canny operator th=0.5 th=0.1
22
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Marr-Hildreth operator
Laplacian of the Gaussian (LoG)
23
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Marr Hildreth operator sigma=1 sigma=2
24
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Hough Transform
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
26
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Hough transform for detecting lines
Alternative definition of lines
27
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Hough transform
Similar transforms can be defined for circles, ellipses or other parametric curves
28
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Morphological operations
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
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
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
32
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Closing
Dilation + Erosion
33
Medical Imaging, SS-2014
Dr. Mohammad Dawood
Opening
Erosion + Dilation