最近一直在做實驗，所以好久沒有更新了，先把上周做的一些小的實驗貼出來供大家分享。

在醫生診斷時，是會將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張處理一張存儲到新的數據集里，希望對大家有幫助！

總結

以上是生活随笔為你收集整理的matlab设置固定的窗宽窗位,python实现CT窗宽窗位的调整（即指定HU值保存图像）...的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。