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

yipeiwu_com6年前Python基础

需要安装OpenCV和SimpleItk。

SimpleItk比较简单,直接pip install SimpleItk即可。

代码如下:

#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!"

如上所示,就将CT影像解析成了多个单幅图片。

python读取CT医学图像

python读取CT医学图像

目前还不太理解(x,y,z)坐标是如何对应的。

但祈世间人无病,何愁架上药生尘。癌症,终有一天可以被战胜。

以上这篇对python读取CT医学图像的实例详解就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持【听图阁-专注于Python设计】。

相关文章

详解Numpy中的广播原则/机制

广播的原则 如果两个数组的后缘维度(从末尾开始算起的维度)的轴长度相符或其中一方的长度为1,则认为它们是广播兼容的。广播会在缺失维度和(或)轴长度为1的维度上进行。 在上面的对arr每...

详解python项目实战:模拟登陆CSDN

详解python项目实战:模拟登陆CSDN

前言 今天为大家介绍一个利用Python模拟登陆CSDN的案例,虽然看起来很鸡肋,有时候确会有大用处,在这里就当做是一个案例练习吧,提高自己的代码水平,也了解Python如何做到模拟登陆...

numpy中的meshgrid函数的使用

numpy官方文档meshgrid函数帮助文档https://docs.scipy.org/doc/numpy/reference/generated/numpy.meshgrid.ht...

python实现两个经纬度点之间的距离和方位角的方法

最近做有关GPS轨迹上有关的东西,花费心思较多,对两个常用的函数总结一下,求距离和求方位角,比较精确,欢迎交流! 1. 求两个经纬点的方位角,P0(latA, lonA), P1(la...

使用Python生成200个激活码的实现方法

题目:使用 Python 生成 200 个不重复的激活码 编写思路 # 激活码一般是由26个大写字母和10个数字任意组合而成 # 长度为12位或者16位的居多激活码 # 一个激活码里的字...