python的三维画图(python处理DICOM并计算三维模型体积)
类别:脚本大全 浏览量:630
时间:2022-01-20 00:04:06 python的三维画图
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
您可能感兴趣
- java入坑rabbitmq(Python操作rabbitMQ的示例代码)
- python零基础入门加密与解密(Python玩转加密的技巧推荐)
- 如何用python做一个弹窗(python实现弹窗祝福效果)
- python微信红包代码(php生成微信红包数组的方法)
- python编辑redis(python使用pipeline批量读写redis的方法)
- python分析excel基础数据生成报表(Python实现定制自动化业务流量报表周报功能XlsxWriter模块)
- python开发网站github(使用 Python 玩转 GitHub 的贡献板推荐)
- 零基础学计算机二级python(计算机二级python学习教程1 教大家如何学习python)
- python字符串分析总结(Python 存储字符串时节省空间的方法)
- python排序方法简单(快速排序的四种python实现推荐)
- python创建文件夹(python批量创建指定名称的文件夹)
- python怎样读取txt文件的数据内容(Python实现读取txt文件中的数据并绘制出图形操作示例)
- python画图教程(详解用Python练习画个美队盾牌)
- python 接口测试怎么校验json数据(python接口自动化十七--Json 数据处理---一次爬坑记详解)
- python类的init方法(Python操作配置文件ini的三种方法讲解)
- python后端生成的pdf文件(Python实现截取PDF文件中的几页代码实例)
- 河南尉氏县因地制宜发展果蔬种植 水坡镇绿宝甜瓜变 金瓜(河南尉氏县因地制宜发展果蔬种植)
- 谢广坤,你这么欺负谢腾飞,良心不会痛吗(你这么欺负谢腾飞)
- 乡村爱情15 宋晓峰怀疑自己孩子,腾飞与姜奶奶亲子鉴定出结果(宋晓峰怀疑自己孩子)
- 《乡村爱情13》开播,新版刘能以假乱真,编剧思维进入瓶颈(新版刘能以假乱真)
- 当年的 白洋淀战神 练肌肉 嘎子哥也成为行走的荷尔蒙(当年的白洋淀战神)
- 肌肉小子陈康, 亚洲巨兽 黄哲勋,哪个才是你的菜(肌肉小子陈康亚洲巨兽)
热门推荐
- dedecms怎么调用标签(DEDECMS 留言薄模块的使用方法)
- mysql出现锁表的原因(导致MySQL做全表扫描的几种情况)
- python简易翻译器的运行(Python3.6实现带有简单界面的有道翻译小程序)
- 网站测试包括哪些内容
- MongoDB 日志文件太大
- ubuntu下使用vscode(使用Vscode结合docker进行开发的详细过程)
- apache怎么安装服务(apache后缀名支持 让apache支持apk ipk下载的方法)
- python转换doc到pdf(利用python将图片版PDF转文字版PDF)
- php实现导航地图(php根据地址获取百度地图经纬度的实例方法)
- php常驻进程(PHP 多进程与信号中断实现多任务常驻内存管理实例方法)
排行榜
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9