python读取dicom序列_python读取dicom图像（SimpleITK和dicom包实现）

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.d