医学影像成像原理

2020/11/4
23
邊緣偵測運算子
2020/11/4
24
邊緣偵測運算子
2020/11/4
25
平滑運算子
2020/11/4
26
閾質分割 (Thresholding)
2020/11/4
27
谢谢
-1 -1 11
2020/11/4
10
銳化濾波器
線性一次空間微分濾波器
3 * 3 mask
-1 0 1 -1 0 1 -1 0 1
-1 -1 -1 000 111
i方向一次微分
j方向一次微分
2020/11/4
11
銳化濾波器
線性一次空間微分濾波器
Sobel filter
-1 0 1 -2 0 2 -1 0 1
m1n1
2020/11/4
19
Low-pass filter
2020/11/4
20
平滑化濾波器
非線性平滑化濾波器
最小值濾波器
g ( i ,j ) m f ( i m i ,j n n ) | 1 { m , n 1 }
最大值濾波器
g ( i ,j ) m f ( i m a ,j n x ) | 1 { m , n 1 }
sum = 0; for(l=0; l<9; l++)
sum += Mask[l] * array[l];
2020/11/4
l_truth = abs(sum);
if (l_truth >= 255)
{
pDoc->cttmp2->m_nVoxel[k][i][j] = 255;
}
else
{
pDoc->cttmp2->m_nVoxel[k][i][j] = l_truth;
mm1nn1
2020/11/4
m1=0 m2=1 n1=0 n2=1
来自百度文库
m1=-1 m2=0 n1=-1 n2=0
m1=-1 m2=1 n1=-1 n2=1
m1=-1
m2=1
n1=0
n2=2
9
銳化濾波器
線性一次空間微分濾波器
2 * 2 mask
i方向一次微分: f i
-1 1 -1 1
j方向一次微分: f j
影像之雜訊去除及影像之平滑化
2020/11/4
7
銳化濾波器
以同一像素近旁之層階值微分或差分 運算
一次空間微分差分濾波器
2 X 2遮罩 3 X 3遮罩 Sobel 濾波器
二次空間微分差分濾波器
2020/11/4
8
銳化濾波器
一次空間微分差分濾波器
m2 n2
g(i, j)
Cm,n•f(im, jn)
}
if (pDoc->cttmp2->m_nVoxel[k][i][j] > max_l)
{
max_l=pDoc->cttmp2->m_nVoxel[k][i][j];
} }
15
拉布拉辛(Laplacian) 濾波器
2020/11/4
16
銳化濾波器
2020/11/4
17
平滑化濾波器
線性平滑化濾波器
-1 -2 -1 000 121
i方向一次微分
j方向一次微分
2020/11/4
12
Sobel filter
2020/11/4
13
銳化濾波器
二次空間微分線性濾波器
拉布拉辛(Laplacian) 濾波器
L
2 i2
2 j2
010 1 -4 1 0 10
1 11 1 -8 1 111
4近旁
8近旁
2020/11/4
Low-pass filter
非線性平滑化濾波器
Median filter
2020/11/4
18
平滑化濾波器
線性平滑化濾波器
Low-pass filter
111 1 11 1 11
11
g(i,j){f(im ,jn)}/9
m 1n1
1 11
121 111
11
g(i, j){
f(im, jn)f(i, j)}/10
医学影像成像原理
什麼是影像處理
幾何處理 影像增強 影像壓縮 影像切割 影像辨識 影像加密 影像描述
2020/11/4
2
影像增強技術
定義 點處理的增強 影像濾波
2020/11/4
3
定義
對影像上之任意像素f(i,j)之新值，
g(i,j)，以像素之周圍層階值局部運算
所得，並經由對整個影像作處理以達 到邊緣強調和雜訊去除的目的。
array[3] = pDoc->cttmp1->m_nVoxel[k][i][j-1]; array[4] = pDoc->cttmp1->m_nVoxel[k][i][j]; array[5] = pDoc->cttmp1->m_nVoxel[k][i][j+1];
array[6] = pDoc->cttmp1->m_nVoxel[k][i+1][j-1]; array[7] = pDoc->cttmp1->m_nVoxel[k][i+1][j]; array[8] = pDoc->cttmp1->m_nVoxel[k][i+1][j+1];
14
int Mask[9] = { 0 , -1, 0 , -1, 4, -1 , 0 , -1, 0 };
for(k=0 ; k<pDoc->cttmp1->m_nFrame ; k++) for(i=1 ; i<pDoc->cttmp1->m_nHeight-1 ; i++) for(j=1 ; j<pDoc->cttmp1->m_nWidth-1 ; j++) { array[0] = pDoc->cttmp1->m_nVoxel[k][i-1][j-1]; array[1] = pDoc->cttmp1->m_nVoxel[k][i-1][j]; array[2] = pDoc->cttmp1->m_nVoxel[k][i-1][j+1];
640
i
f(i,j)
480
2020/11/4
j
4
點處理的增強
灰階反轉
L-1
s
T(r)
2020/11/4
5
0 r L-1
點處理的增強
對比度擴展
L-1
(r2,s2)
s T(r)
2020/11/4
(r1,s1)
6
0
r L-1
影像濾波
g(x,y)=T[f(x,y)]
銳化濾波器
影像之線或邊緣強調
平滑化濾波器
中間值濾波器(Median filter)
g ( i ,j ) m{ f e ( i m d ,j n i ) | a 1 m , n n 1 }
2020/11/4
21
中間值濾波器(Median filter)
2020/11/4
22
影像切割技術
邊緣偵測運算子 平滑運算子 閾質分割 (Thresholding)