matlabbatch{1}.spm.util.import.dicom.data = data_subdir3 ;
matlabbatch{1}.spm.util.import.dicom.outdir = {data_dir2} ;
spm_jobman('run', matlabbatch) ;
存在的问题是:单个病人work,多个病人循环处理经常会出现保存错乱的情况。
3. 使用python代码进行转换,代码如下:
series_file_names = sitk.ImageSeriesReader.GetGDCMSeriesFileNames(path)
series_reader = sitk.ImageSeriesReader()
series_reader.SetFileNames(series_file_names)
img_dt = series_reader.Execute()
#img_dt = sitk.ReadImage(path) #如果是单个dicom文件,请使用这个,注释前面几行
#读取完dicom 直接保存就可
origin =img_dt.GetOrigin()
direction = img_dt.GetDirection()
space = img_dt.GetSpacing()
img_dt.SetOrigin(origin)
img_dt.SetDirection(direction)
img_dt.SetSpacing(space)
sitk.WriteImage(img_dt, savedir+'xxx.nii')
存在的问题:保存得到的图像方向有时候是错的,虽然最终set了direction,并且头文件没有得到保存,且处理起来慢。
4. 使用一些软件手动处理:①SimpleITK:读入dicom,再保存为nii
②DICOM Converter (别用,依托答辩,数据一多就卡住,还收费)
③MRIcroGL(类似于mricron),和mricron一样,GUI程序不能批处理。
那就是MRIcroGL里的dcm2niix结合python,直接上代码:
import os
import sys
sys.path.append(r'C:\Users\Administrator\Desktop\MRIcroGL_windows\MRIcroGL\Resources')
base_path = r'L:\tougu_data\xx'
for i in os.listdir(base_path):
path1 = base_path + '/' + i
for ii in os.listdir(path1):
dist_dir = path1 + '/' + ii
new_name = dist_dir
os.system('chcp 65001')
os.system(r'C:\Users\Administrator\Desktop\MRIcroGL_windows\MRIcroGL\Resources\dcm2niix {0}'.format(dist_dir))
其实就是通过python调用windows shell,然后借助python可以传入变化参数达到批处理。这里有几点需要说明:①os.system(‘chcp 65001’)不加可能会报错编码错误。②dcm2niix 的传入参数有很多,可以自行去shell里输入dcm2niix来查看其功能,比如输出文件名设置等。
题外话
如果你是零基础小白,想快速入门Python是可以考虑的。
一方面是学习时间相对较短,学习内容更全面更集中。
二方面是可以找到适合自己的学习方案
包括:Python安装包+激活码、Python web开发,Python爬虫,Python数据分析,人工智能、机器学习等教程。带你从零基础系统性的学好Python!
零基础Python学习资源介绍
Python学习路线汇总
Python所有方向的技术点做的整理,形成各个领域的知识点汇总,它的用处就在于,你可以按照上面的知识点去找对应的学习资源,保证自己学得较为全面。(学习教程文末领取哈)
Python必备开发工具
温馨提示:篇幅有限,已打包文件夹,获取方式在:文末
Python学习视频600合集
观看零基础学习视频,看视频学习是最快捷也是最有效果的方式,跟着视频中老师的思路,从基础到深入,还是很容易入门的。
实战案例
光学理论是没用的,要学会跟着一起敲,要动手实操,才能将自己的所学运用到实际当中去,这时候可以搞点实战案例来学习。
100道Python练习题
检查学习结果。
面试刷题
上述这份完整版的Python全套学习资料已经上传CSDN官方,朋友们如果需要可以微信扫描下方CSDN官方认证二维码输入“领取资料” 即可领取。