<强> 1。用SimpleITK读取dicom序列:
进口SimpleITK sitk
进口numpy np
img_path=' F:胰腺数据集\ \ \ \ \ \ \ \输出原始厚\ \ \ \ 1 '
mask_path=' F:胰腺数据集\ \ \ \ \ \ \ \ \ \ groundtruth厚\ \输出1 '
读者=sitk.ImageSeriesReader ()
img_names=reader.GetGDCMSeriesFileNames (img_path)
reader.SetFileNames (img_names)
形象=reader.Execute ()
image_array=sitk.GetArrayFromImage(图片)# z, y、x
读者=sitk.ImageSeriesReader ()
mask_names=reader.GetGDCMSeriesFileNames (mask_path)
reader.SetFileNames (mask_names)
掩码=reader.Execute ()
mask_array=sitk.GetArrayFromImage(面具)# z, y、x
<强> 2。用dicom读取单张dicom图像并显示:
进口日本
进口pylab
ds=dicom.read_file(胰腺”数据集F: \ \ \ \ \ \ \ \ \ \ groundtruth厚\ \输出1 \ \ FILE0001_seg.dcm”)
pixel_bytes=ds.PixelData
# # CT值组成了一个矩阵
照片=ds.pixel_array
# #读取显示图片
pylab.imshow (ds。pixel_array提出=pylab.cm.bone)
pylab.show ()
如果要对dicom图像中的像素值进行修改,继续执行以下代码:
# #修改图片中的元素,不能直接使用data_array,需要转换成PixelData
n, val列举(ds.pixel_array.flat): #例子:零任何& lt;300
如果val & lt;300:
ds.pixel_array.flat [n]=0
ds。PixelData=https://www.yisu.com/zixun/ds.pixel_array.tostring ()
ds.save_as (“newfilename.dcm”)
之前
<强> 3。此外,用pydicom也可读取dicom图像
以上这篇python读取dicom图像示例(SimpleITK和dicom包实现)就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
