白云岛资源网 Design By www.pvray.com
1. 用SimpleITK读取dicom序列:
import SimpleITK as sitk import numpy as np img_path='F:\\dataset\\pancreas\\Output\\thick\\original\\1' mask_path='F:\\dataset\\pancreas\\Output\\thick\\groundtruth\\1' reader = sitk.ImageSeriesReader() img_names = reader.GetGDCMSeriesFileNames(img_path) reader.SetFileNames(img_names) image = reader.Execute() image_array = sitk.GetArrayFromImage(image) # z, y, x reader = sitk.ImageSeriesReader() mask_names = reader.GetGDCMSeriesFileNames(mask_path) reader.SetFileNames(mask_names) mask = reader.Execute() mask_array = sitk.GetArrayFromImage(mask) # z, y, x
2. 用dicom读取单张dicom图像并显示:
import dicom import pylab ds=dicom.read_file("F:\\dataset\\pancreas\\Output\\thick\\groundtruth\\1\\FILE0001_seg.dcm") pixel_bytes = ds.PixelData ##CT值组成了一个矩阵 pix = ds.pixel_array ##读取显示图片 pylab.imshow(ds.pixel_array, cmap=pylab.cm.bone) pylab.show() 如果要对dicom图像中的像素值进行修改,继续执行以下代码: ##修改图片中的元素,不能直接使用data_array,需要转换成PixelData for n,val in enumerate(ds.pixel_array.flat): # example: zero anything < 300 if val < 300: ds.pixel_array.flat[n]=0 ds.PixelData = ds.pixel_array.tostring() ds.save_as("newfilename.dcm")
3. 此外,用pydicom也可读取dicom图像
以上这篇python读取dicom图像示例(SimpleITK和dicom包实现)就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
白云岛资源网 Design By www.pvray.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
白云岛资源网 Design By www.pvray.com
暂无评论...
RTX 5090要首发 性能要翻倍!三星展示GDDR7显存
三星在GTC上展示了专为下一代游戏GPU设计的GDDR7内存。
首次推出的GDDR7内存模块密度为16GB,每个模块容量为2GB。其速度预设为32 Gbps(PAM3),但也可以降至28 Gbps,以提高产量和初始阶段的整体性能和成本效益。
据三星表示,GDDR7内存的能效将提高20%,同时工作电压仅为1.1V,低于标准的1.2V。通过采用更新的封装材料和优化的电路设计,使得在高速运行时的发热量降低,GDDR7的热阻比GDDR6降低了70%。