python读取dicom图像示例(SimpleITK和dicom包实现)


Posted in Python onJanuary 16, 2020

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包实现)就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python3爬取淘宝信息代码分析
Feb 10 Python
python清除字符串中间空格的实例讲解
May 11 Python
Python/ArcPy遍历指定目录中的MDB文件方法
Oct 27 Python
python PyAutoGUI 模拟鼠标键盘操作和截屏功能
Aug 04 Python
django将网络中的图片,保存成model中的ImageField的实例
Aug 07 Python
通过selenium抓取某东的TT购买记录并分析趋势过程解析
Aug 15 Python
python使用socket实现的传输demo示例【基于TCP协议】
Sep 24 Python
Python BeautifulReport可视化报告代码实例
Apr 13 Python
python安装读取grib库总结(推荐)
Jun 24 Python
使用npy转image图像并保存的实例
Jul 01 Python
Python根据字符串调用函数过程解析
Nov 05 Python
Python实现生活常识解答机器人
Jun 28 Python
.dcm格式文件软件读取及python处理详解
Jan 16 #Python
用python解压分析jar包实例
Jan 16 #Python
Python3 实现爬取网站下所有URL方式
Jan 16 #Python
python3爬取torrent种子链接实例
Jan 16 #Python
Python 实现自动获取种子磁力链接方式
Jan 16 #Python
关于初始种子自动选取的区域生长实例(python+opencv)
Jan 16 #Python
Python简单实现区域生长方式
Jan 16 #Python
You might like
MYSQL 小技巧 -- LAST_INSERT_ID
2009/11/24 PHP
php下拉选项的批量操作的实现代码
2013/10/14 PHP
PHP通过加锁实现并发情况下抢码功能
2016/08/10 PHP
解决Laravel 使用insert插入数据,字段created_at为0000的问题
2019/10/11 PHP
Laravel validate error处理,ajax,json示例
2019/10/25 PHP
css动画效果之animation的常用样式
2021/03/09 HTML / CSS
30个最佳jQuery Lightbox效果插件分享
2011/04/11 Javascript
读jQuery之二(两种扩展)
2011/06/11 Javascript
如何制作浮动广告 JavaScript制作浮动广告代码
2012/12/30 Javascript
关于JavaScript与HTML的交互事件
2013/04/12 Javascript
js制作支付倒计时页面
2016/10/21 Javascript
Django中使用jquery的ajax进行数据交互的实例代码
2017/10/15 jQuery
微信小程序云开发之模拟后台增删改查
2019/05/16 Javascript
Vue前端项目部署IIS的实现
2020/01/06 Javascript
解决vue 给window添加和移除resize事件遇到的坑
2020/07/21 Javascript
vue界面发送表情的实现代码
2020/09/11 Javascript
vue实现抽屉弹窗效果
2020/11/15 Javascript
[01:37]PWL S2开团时刻DAY1&2——这符有毒
2020/11/20 DOTA
详解使用Python处理文件目录的相关方法
2015/10/16 Python
Python3中正则模块re.compile、re.match及re.search函数用法详解
2018/06/11 Python
Window 64位下python3.6.2环境搭建图文教程
2018/09/19 Python
Python零基础入门学习之输入与输出
2019/04/03 Python
Python 切分数组实例解析
2019/11/07 Python
py-charm延长试用期限实例
2019/12/22 Python
python opencv把一张图片嵌入(叠加)到另一张图片上的实现代码
2020/06/11 Python
Pycharm2020最新激活码|永久激活(附最新激活码和插件的详细教程)
2020/09/29 Python
谈谈对css属性box-sizing的了解
2017/01/04 HTML / CSS
洛杉矶健身中心女性专用运动服饰品牌:Marika
2018/05/09 全球购物
Kathmandu新西兰官网:新西兰户外运动品牌
2019/07/27 全球购物
优秀党员转正的自我评价
2013/10/06 职场文书
老公给老婆的道歉信
2014/01/10 职场文书
思想品德自我评价
2014/02/04 职场文书
护士求职自荐信范文
2014/03/19 职场文书
国家励志奖学金个人先进事迹材料
2014/05/04 职场文书
党的群众路线教育实践活动领导班子对照检查材料
2014/09/25 职场文书
政风行风整改方案
2014/10/25 职场文书