时间:2021-05-22
目标:利用python读取dicom文件,并进行处理生成info.txt和raw文件
实现:通过pydicom读取dicom文件
代码:
import numpyimport pydicomimport os# dicom文件所在的文件夹目录PathDicom = '/home/lk/testdata/1.3.6.1.4.1.9328.50.1.42697596859477567872763647333745089432/'# 筛选出文件夹目录下所有的dicom文件lstFilesDCM = []for dirName, subdirList, fileList in os.walk(PathDicom): for filename in fileList: if '.dcm' in filename.lower(): lstFilesDCM.append(os.path.join(dirName, filename))# Get ref fileRefDs = pydicom.read_file(lstFilesDCM[0])# Load dimensions based on the number of rows, columns, and slices (along the Z axis)ConstPixelDims = (int(RefDs.Rows), int(RefDs.Columns), len(lstFilesDCM))# Load spacing values (in mm)ConstPixelSpacing = (float(RefDs.PixelSpacing[0]), float(RefDs.PixelSpacing[1]), float(RefDs.SliceThickness))# save info.txtinfo = ConstPixelDims + ConstPixelSpacingf = open('/home/lk/testdata/1.3.6.1.4.1.9328.50.1.42697596859477567872763647333745089432/info.txt', 'w')for n in info: f.write(str(n)+' ')f.close()# According to location sortinglocation = []for i in range(len(lstFilesDCM)): ds = pydicom.read_file(lstFilesDCM[i]) location.append(ds.SliceLocation)location.sort()# The array is sized based on 'ConstPixelDims'ArrayDicom = numpy.zeros((len(lstFilesDCM), RefDs.Rows, RefDs.Columns), dtype=RefDs.pixel_array.dtype)# loop through all the DICOM filesfor filenameDCM in lstFilesDCM: # read the file ds = pydicom.read_file(filenameDCM) # store the raw image data ArrayDicom[location.index(ds.SliceLocation), :, :] = ds.pixel_array# save rawds = ArrayDicom.tostring()f = open('/home/lk/testdata/1.3.6.1.4.1.9328.50.1.42697596859477567872763647333745089432/1.raw', 'wb')f.write(ds)f.close()代码编写过程遇到的问题及解决方法:
Problem one: pydicom版本问题。
pydicom1.x中读取dicom文件调用pydicom.read_file(filename);
pydicom0.9中读取dicom文件调用dicom.read_file(filename);
Problem two:python中IO操作
(1) f = open(filename, mode)
其中filename为文件的路径, mode为操作标识符:‘r' 表示读, ‘w'表示写,‘a'表示既可读又可写,‘b'表示二进制文件。
(2) f.write(value)
其中参数value必须是字符串类型的。
当然还有一些其他的问题,在这里就不细说了,多入坑才能学的多,切不可烦躁,代码就是要多敲才能得心应手,共勉。
以上这篇python 读取dicom文件,生成info.txt和raw文件的方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
持久化文件读写:f=open('info.txt','a+')f.seek(0)str1=f.read()iflen(str1)==0:f1=open('inf
android之文件操作——读取assets和raw文件下的内容1.分别创建assets文件夹和res/raw文件夹:(要注意的raw文件是在res下new,然
读取txt的数据和把数据保存到txt中是经常要用到的,下面我就总结一下。读txt文件python常用的读取文件函数有三种read()、readline()、re
本文实例讲述了android从资源文件中读取文件流并显示的方法。分享给大家供大家参考。具体如下:在android中,假如有的文本文件,比如TXT放在raw下,要
python读取txt文件:(思路:先打开文件,读取文件,最后用for循环输出内容)1、读取1.1基于pythoncsv库#3.读取csv至字典x,yimpor