灰级窗切片增强--相比较X光更适合于CT,程序如下:
import cv2
import numpy as np
I = cv2.imread('*.png',0)
h, w = I.shape
fa = 40
fb = 160
k = 255/(160-45)
J = np.zeros((h,w))
for i in range(h):
for j in range(w):
if I[i,j] >=0 and I[i,j] < fa:
J[i,j] = 0
elif I[i,j] < fb:
J[i, j] = k*(I[i,j] - fa)
print(I)
print(J)
cv2.imwrite('original.png',I)
cv2.imwrite('new.png', J)
对于fa与fb的选取,主要考虑的是所选择区域的灰度值范围。这个不具有适用性。考虑自适应的选取这两个值是下一步的工作计划。