对python读取CT医学图像的实例详解

#coding:utf-8
import SimpleITK as sitk
import cv2
 
#LKDS-00058,-102.655469971,108.188810974,438.759994507,12.2279986879
if __name__ == '__main__':
  filename = "F:/cancer_solution/data/train_subset00/LKDS-00058.mhd"
  ds = sitk.ReadImage(filename)
  img_array = sitk.GetArrayFromImage(ds)
  frame_num, width, height = img_array.shape
 
  outpath = "F:/cancer_solution/out/train/LKDS-00058"
  index = -1
  for img_item in img_array:
    index = index + 1
    cv2.imwrite("%s/%d.png"%(outpath,index),img_item)
 
  print "done!"