1 Medical Imaging, SS-2014 Dr. Mohammad Dawood Medical Image Analysis Dr. Mohammad Dawood Department of Computer Science University of Münster Germany
Apr 01, 2015
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