最近一直在做实验，所以好久没有更新了，先把上周做的一些小的实验贴出来供大家分享。

在医生诊断时，是会将CT图像调整成不同的窗来处理的。比如说肺部CT吧，肺窗(窗宽为2000，窗位为-400)用于看小的肺结节；腹部窗(窗宽400，窗位40)用于看大的肺结节和纵隔淋巴结，还有骨窗。所以在处理CT图像的时候会需要选择不同的窗，那么怎么来使得你的.dcm文件.IMA文件.mhd文件以不同的窗宽窗位显示：# -*- coding=utf-8 -*-

import dicom

import os

import cv2

import numpy

from os.pathimport splitext

def getfile(file):

dcm = dicom.read_file(file)

img2 = dcm.pixel_array * dcm.RescaleSlope + dcm.RescaleIntercept

return img2def get_window_size(window_type):

if window_type =='lung':

center = -500

width =2000

elif window_type =='abdomen':

center =40

width =400

elif window_type =='bone':

center =300

width =2000

return center, width#调整CT图像的窗宽窗位

def setDicomWinWidthWinCenter(img_data, window_type):

img_temp = img_data

rows =len(img_temp)

cols =len(img_temp[0])

center, width = get_window_size(window_type)

img_temp.flags.writeable =True

min = (2 * center - width) /2.0 +0.5

max = (2 * center + width) /2.0 +0.5

dFactor =255.0 / (max - min)

for iin numpy.arange(rows):

for jin numpy.arange(cols):

img_temp[i, j] =int((img_temp[i, j]-min)*dFactor)

min_index = img_temp <0

img_temp[min_index] =0

max_index = img_temp >255

img_temp[max_index] =255

return img_temppath ='G:\lowdosect\CH_image/quater_dose_image\L067/'

for root, dirs, filesin os.walk(path):

for i in range(0,len(files),10):

filename = splitext(files[i])[0] +"." +"bmp"

print(filename)

#print(files[i])

im = getfile(path+files[i])

im1 = setDicomWinWidthWinCenter(im,'bone')

cv2.imwrite("G:\lowdosect\CH_image/bone/quarter/" + filename, im1)

还不太会编辑，所以可能空格会有问题，我是在数据集里每10张处理一张存储到新的数据集里，希望对大家有帮助！