韩漫免费漫画在线观看方法,《好好疼爱里面》免费看,年轻漂亮的女士护士内衣,妈妈醉酒后把我当爸爸电视剧

文章 > Python基础教程 > 怎样用python加载dicom图片

怎样用python加载dicom图片

头像

FXL

2020-08-24 10:34:193398浏览 · 0收藏 · 0评论

用python加载dicom图片的方法:使用pydicom、CV2、numpy、matplotlib等库即可。pydicom库是专门用来处理dicom图像的python专用库。

python读取DICOM图像,需要以下几个库:pydicom、CV2、numpy、matplotlib。pydicom是专门处理dicom图像的python专用包,numpy高效处理科学计算的包,依据数据绘图的库。

(推荐教程:Python入门教程

安装需要的库

pip install matplotlib
pip install opencv-python
pip install pydicom    
pip install numpy

安装好这些库后就可以对dicom文件操作了。

具体代码如下:

#-*-coding:utf-8-*-
import cv2
import numpy
import dicom
from matplotlib import pyplot as plt

dcm = dicom.read_file("AT0001_100225002.DCM")
dcm.image = dcm.pixel_array * dcm.RescaleSlope + dcm.RescaleIntercept

slices = []
slices.append(dcm)
img = slices[ int(len(slices)/2) ].image.copy()
ret,img = cv2.threshold(img, 90,3071, cv2.THRESH_BINARY)
img = numpy.uint8(img)

im2, contours, _ = cv2.findContours(img,cv2.RETR_LIST,cv2.CHAIN_APPROX_SIMPLE)
mask = numpy.zeros(img.shape, numpy.uint8)
for contour in contours:
 cv2.fillPoly(mask, [contour], 255)
img[(mask > 0)] = 255


kernel = cv2.getStructuringElement(cv2.MORPH_ELLIPSE,(2,2))
img = cv2.morphologyEx(img, cv2.MORPH_OPEN, kernel)


img2 = slices[ int(len(slices)/2) ].image.copy()
img2[(img == 0)] = -2000


plt.figure(figsize=(12, 12))
plt.subplot(131)
plt.imshow(slices[int(len(slices) / 2)].image, 'gray')
plt.title('Original')
plt.subplot(132)
plt.imshow(img, 'gray')
plt.title('Mask')
plt.subplot(133)
plt.imshow(img2, 'gray')
plt.title('Result')
plt.show()

在DICOM图像里,包含了患者的相关信息的字典,我们可以通过dir查看DICOM文件有什么信息,可以通过字典返回相关的值。

import dicom
import json
def loadFileInformation(filename):
 information = {}
 ds = dicom.read_file(filename)
 information['PatientID'] = ds.PatientID
 information['PatientName'] = ds.PatientName
 information['PatientBirthDate'] = ds.PatientBirthDate
 information['PatientSex'] = ds.PatientSex
 information['StudyID'] = ds.StudyID
 information['StudyDate'] = ds.StudyDate
 information['StudyTime'] = ds.StudyTime
 information['InstitutionName'] = ds.InstitutionName
 information['Manufacturer'] = ds.Manufacturer
 print dir(ds)
 print type(information)
 return information

a=loadFileInformation('AT0001_100225002.DCM')
print a
关注

关注公众号,随时随地在线学习

本教程部分素材来源于网络,版权问题联系站长!

小妹妹爱大棒棒免费观看电视剧一 | 无码| 狗和女人戴避孕套会传染吗| 打扑克的全程视频| 慢脚APP官方版| 爸爸你真棒插曲快来救救我小说| 父亲开了两个女儿包演员表| 人与狗胶配方大全高清视频| 无码人妻久久一区二区| 小叔子从后面抱住我的电视剧| 他的手渐渐地也开始不安分起来| 国产免费观看高清电视剧在线观看| 偷偷在线观看免费的电视剧最新| 双飞娘俩国语免费观看视频| 亚洲一区| 《再来一次好吗》免费观看| 器材室轮-J-(4)| 彼岸岛浸画| 脏脏漫| 星空传媒天美传媒有限公司| 日韩AV| 妈妈你真棒里面的插曲叫什么| 女人吃男人鸡要多久恢复| 榴莲APP官网下载1.0.3| 老旺的大肉蟒进进出出| MISSAV| 办公室深深挺进女警小| 苍老师巅峰之作第37集| 冢田诗织| 彩浸福利| 人间中毒车内吸乳被咬后怎么处理| 蜜源APP| 交换老公弄了几次高潮| 欧美| 妈妈的朋友2| 太太你想丈夫被开除吧| 免费看网站在线观看人数在哪直播| 主动把乳头放进男朋友嘴巴里| 《别告诉我丈夫》播放| 边吃奶边扭身体哼唧怎么回事| 欧美黑人狂躁日本妞