SimpleITK读取医学影像文件

最新推荐文章于 2025-07-08 10:39:37 发布

转载最新推荐文章于 2025-07-08 10:39:37 发布 · 7.6k 阅读

深度学习同时被 2 个专栏收录

5 篇文章

订阅专栏

Tips

5 篇文章

订阅专栏

本文详细介绍了ITK和SimpleITK在图像分析和处理中的应用，包括X射线和CT图像的读取、矩阵信息获取及CT影像的拆分方法，提供了丰富的代码示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原文：https://blog.youkuaiyun.com/songchunxiao1991/article/details/79074604

SimpleITK
ITK是一个开源、跨平台的框架，提供给开发者增强功能的图像分析和处理套件（推荐使用）。

Note:注意SimpleITK不支持中文，即路径中不能有中文

X射线图像对应的读取

import SimpleITK as sitk
filename = '*.dcm'
itk_img = sitk.ReadImage(filename)
img_array = sitk.GetArrayFromImage(itk_img)
print(img_array.shape)

输出：（对应信息：frame_num, width, height）注意：img_array对应的是图像的矩阵信息

(1, 3328, 2560)

推荐用法：

import SimpleITK as sitk
filename = '*.dcm'
itk_img = sitk.ReadImage(filename)
img_array = sitk.GetArrayFromImage(itk_img)[0]
print(img_array.shape)
输出：

(3328, 2560)

CT图像的读取

import SimpleITK as sitk
filename = '*.mhd'
itk_img = sitk.ReadImage(filename)
img_array = sitk.GetArrayFromImage(itk_img)
print(img_array.shape)
输出：（对应信息：frame_num, width, height）帧参数：frame_num，代表CT扫描层数注意：img_array对应的是图像的矩阵信息

(133, 512, 512)

将CT影像拆分成多个单幅图像

方法一：

import SimpleITK as sitk
import cv2
filename = '*.mhd'
itk_img = sitk.ReadImage(filename)
img_array = sitk.GetArrayFromImage(itk_img)
frame_num, width, height = img_array.shape
outpath = 'F:/data/LUNA16'#存放拆分得到的图像路径
index = -1
for img_item in img_array:
    index = index + 1
    cv2.imwrite("%s/%d.png" % (outpath, index), img_item)
print("done!")

方法二：

import SimpleITK as sitk
import os
import cv2
filename = '*.mhd'
path = 'F:/data/LUNA16'   #存放拆分得到图像的路径
itk_img = sitk.ReadImage(filename)
img_array = sitk.GetArrayFromImage(itk_img)
for i, im in enumerate(img_array):
  cv2.imwrite(os.path.join(path, '{}.png'.format(i)), im)

如下图所示是CT影像拆分出来的部分图像。拆分得到的图像大小均为512*512，总数是：frame_num张。