-
Notifications
You must be signed in to change notification settings - Fork 0
/
meningkatkan-kualitas-citra-ctscan.py
48 lines (39 loc) · 1.25 KB
/
meningkatkan-kualitas-citra-ctscan.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
import cv2 as cv
from matplotlib import pyplot as plt
import numpy as np
# baca citra dalam mode grayscale
img = cv.imread('covid_image.png', cv.IMREAD_GRAYSCALE)
mask = cv.imread('mask.png', cv.IMREAD_GRAYSCALE)
# operasi bitwise_and
img_bitwise = cv.bitwise_and(img,img,mask=mask)
# filter
img_noiseRemov = cv.medianBlur(img_bitwise,5)
# plotting histogram
plt.subplot(2, 1, 1)
plt.hist(img_noiseRemov.ravel(),256,[0,256])
plt.xlim([0,256])
plt.show()
# perbaiki contras citra
hist,bins = np.histogram(img_noiseRemov.flatten(),256,[0,256])
cdf = hist.cumsum()
cdf_m = np.ma.masked_equal(cdf,0)
cdf_m = (cdf_m - cdf_m.min())*255/(cdf_m.max()-cdf_m.min())
cdf = np.ma.filled(cdf_m,0).astype('uint8')
img_contras = cdf[img_noiseRemov]
plt.subplot(2, 1, 2)
plt.hist(img_contras.ravel(),256,[0,256])
plt.xlim([0,256])
plt.show()
# threshold
ret,img_th = cv.threshold(img_contras,175,255,cv.THRESH_BINARY)
# menampilkan citra hasil
final = cv.bitwise_and(img_contras,img_contras,mask=img_th)
cv.imshow('Final', final)
# menyimpan citra
cv.imwrite('res_bitwise.jpg',img_bitwise)
cv.imwrite('res_noise_removal.jpg',img_noiseRemov)
cv.imwrite('res_contras.jpg',img_contras)
cv.imwrite('threshold.jpg',img_th)
cv.imwrite('final.jpg',final)
cv.waitKey(0)
cv.destroyAllWindows()