python处理DICOM并计算三维模型体积
作者:qd_tudou 发布时间:2022-07-16 09:24:29
标签:python,DICOM,三维模型
在已知DICOM和三维模型对应掩膜的情况下,计算三维模型的体积。
思路:
1、计算每个体素的体积。每个体素为长方体,x,y为PixelSpacing,z为层间距
使用pydicom.read_file读取DICOM文件,dcm_tag.PixelSpacing获取像素间距,dcm_tag.SliceLocation 获取层间距
2、计算体素的个数
代码如下:
from PIL import Image
import numpy as np
import pydicom
import os
def get_pixels_No(bmp_data_dir):
pixels_No = 0
bmp_files = os.listdir(bmp_data_dir)
for bmp in bmp_files:
bmp_file = os.path.join(bmp_data_dir,bmp)
img = Image.open(bmp_file)
img_array = np.array(img)
# img_array.dtype为布尔类型,需要转换为Int类型,其累加和恰好为体素总和
img_array_int = img_array.astype(int)
pixels_No = pixels_No+img_array_int.sum()
return pixels_No
def get_pixel_info(dcm_data_dir):
pixel_infos = []
dcm_files = os.listdir(dcm_data_dir)
dcm_file_1 = os.path.join(dcm_data_dir,dcm_files[0])
dcm_tag_1 = pydicom.read_file(dcm_file_1)
# 获取像素间距.
spacex, spacey = dcm_tag_1.PixelSpacing
# 获取层间距
# 有些 dcm图像并不是按照InstanceNumber进行排序的,不能直接用最后一张的slicelocation减去第一张,再除以张数
SliceLocations = []
ImagePositon_z = []
for dcm in dcm_files:
dcm_file = os.path.join(dcm_data_dir, dcm)
dcm_tag = pydicom.read_file(dcm_file)
SliceLocations.append(dcm_tag.SliceLocation)
ImagePositon_z.append(dcm_tag.ImagePositionPatient[2])
SliceLocations_max =max(SliceLocations)
SliceLocations_min =min(SliceLocations)
ImagePositon_z_max = max(ImagePositon_z)
ImagePositon_z_min = min(ImagePositon_z)
print(SliceLocations_max)
print(SliceLocations_min)
print(ImagePositon_z_max)
print(ImagePositon_z_min)
if SliceLocations_max - SliceLocations_min < 1e-10:
spacez = abs(ImagePositon_z_max - ImagePositon_z_min)/(len(dcm_files)-1)
else:
spacez = abs(SliceLocations_max - SliceLocations_min)/(len(dcm_files)-1)
pixel_infos = [spacex, spacey, spacez]
return pixel_infos
def get_volume(dcm_data_dir,bmp_data_dir):
pixel_infos = get_pixel_info(dcm_data_dir)
pixels_No = get_pixels_No(bmp_data_dir)
volume=pixel_infos[0]*pixel_infos[1]*pixel_infos[2]*pixels_No/1000
return volume
# dcm = pydicom.read_file(r"E:\20181210090945_LENG HONGYING F-44Y\Venous\0000.dcm")
# print(dcm)
# print(dcm.ImagePositionPatient[2])
# print(dcm[0x0020, 0x0032].keyword,dcm[0x0020, 0x0032].value)
volume=get_volume(r"E:\20181210090945_LENG HONGYING F-44Y\Venous",r"E:\20181210090945_LENG HONGYING F-44Y\Results\LungL")
print("体积为%.1f"%volume)
来源:https://blog.csdn.net/yaoliuwei1426/article/details/85680152
0
投稿
猜你喜欢
- Python 很强大其原因就是因为它庞大的三方库 , 资源是非常的丰富 , 当然也不会缺少关于音频的库关于音频, PyAudio 这个库,
- 目录arrow模块的使用获取arrow对象时间形式转换获取数据修改时间总结Python中有很多时间和日期处理的库,有time、datetim
- 一般情况下,导出超时可能都是以下三种情况:一、sql语句复杂,查询时间过长;二、处理查询后数据逻辑冗余;三、数据量过大导致响应超时。接下来分
- asp使用WScript.Shell获取电脑的网络配置信息Option Explicit Dim WSHShe
- 报错:Uncaught SyntaxError: Unexpected token o in JSON at position 1at JS
- <script> Function.prototype.$bind=function(object) {  
- 人生苦短,快学Python!今天介绍的是psutil模块,它是一个跨平台库 https://github.com/giampaolo/psu
- asp三天学好ADO对象之第一天 今天说一下Recordset 对象的属性1、CursorType 属性AdOpenForwardOnly:
- Pandas库十分强大,但是对于切片操作iloc, loc和ix,很多人对此十分迷惑,因此本篇博客利用例子来说明这3者之一的区别和联系,尤其
- 前言最近做了几个简单的爬虫python程序,于是就想做个窗口看看效果。首先是,窗口的话,以前没怎么接触过,就先考虑用Qt制作简单的ui。这里
- 1 引言在进行图像处理过程中,我们经常会遇到一些和物体轮廓相关的操作,比如求目标轮廓的周长面积等,我们直接使用Opencv的findCont
- 一、为何人工智能(AI)首选Python?读完这篇文章你就知道了。我们看谷歌的TensorFlow基本上所有的代码都是C++和Python,
- 如果你正在运行使用MySQL的Web应用程序,那么你把密码或者其他敏感信息保存在应用程序里的机会就很大。保护这些数据免受黑客或者窥探者的获取
- asp使用session来防止表单多次被提交的方法。formtest.asp' 表单文件<%Randomize&nb
- 这篇文章主要介绍了Python3打包exe代码2种方法实例解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,
- optimizer.param_groups: 是长度为2的list,其中的元素是2个字典;optimizer.param_groups[0
- 本文实例讲述了php7 参数、整形及字符串处理机制修改。分享给大家供大家参考,具体如下:参数处理机制修改一、重复参数命名不再支持。重复的参数
- 1、下载python3.6的安装包: wget https://www.python.org/ftp/p
- 项目编写过程中,总能遇见对字典进行排序什么的,如果要实现多条件排序只需要下面几行代码实现。充分体现了python的好处了。teamitems
- 1 简介在日常运行程序的过程中常常涉及到循环迭代过程,对于执行时间很短的程序来说倒无所谓,但对于运行过程有明显耗时的涉及循环迭代的程序,为其