【转】DICOM图像像素值(灰度值)转换为CT值
转自:https://www.cnblogs.com/xuhui24/p/6193032.html
https://zhuanlan.zhihu.com/p/358770379
CT值的单位是Hounsfield,简称为Hu,范围是-1024-3071。用于衡量人体组织对X射线的吸收率,设定水的吸收率为0Hu。
在DICOM图像读取的过程中,我们会发现图像的像素值有可能不是这个范围,通常是0-4096,这是我们常见到的像素值或者灰度值,这就需要我们在图像像素值(灰度值)转换为CT值。
首先,需要读取两个DICOM Tag信息,(0028|1052):rescale intercept和(0028|1053):rescale slope.
然后通过公式:
Hu = pixel * slope + intercept
计算得到CT值。
HU(Hounsfiled Unit)值,是在做医学图像数据预处理不可避免的事情,反映了组织对X射线吸收程度。以水的吸收程度作为参考,即水的HU=0,衰减系数大于水的为正直,小于水的为负值。并以骨皮质和空气的HU值为上限和下限,是在医学影像图像处理技术,处理CT图像不可避免的事情。
HU值的计算公式:(非常重要)
HU = pixel_val * slope + intercept
(如果Slope为1,Intercept为0,则不需要转化的问题,具体Slope和Intercept的值在头文件中都会有说明的)
代码:
import os
import numpy as np
import pydicom as pyd
import matplotlib.pyplot as plt
import cv2
import SimpleITK as sitk
##查看是否需要转化HU值,如果img.RescaleSlope不为1,img.RescaleIntercept不为0,则需要转化的问题
dcm_file='/Pancreas-CT/PANCREAS_0001/11-24-2015-PANCREAS0001-Pancreas-18957/Pancreas-99667/1-001.dcm'
img=pyd.read_file(dcm_file)
img_array=sitk.GetArrayFromImage(sitk.ReadImage(dcm_file))
print(img_array)
print(np.shape(img_array))
print('----')
print(img.RescaleSlope)
print(img.RescaleIntercept)
HU=np.dot(img_array,img.RescaleSlope)+img.RescaleSlope
print(HU)
【转】DICOM图像像素值(灰度值)转换为CT值相关推荐
- DICOM图像像素值(灰度值)转换为CT值
CT值的单位是Hounsfield,简称为Hu,范围是-1024-3071.用于衡量人体组织对X射线的吸收率,设定水的吸收率为0Hu. 在DICOM图像读取的过程中,我们会发现图像的像素值有可能不是这 ...
- 【转】DICOM图像像素值、灰度值与CT值!!!!!
转自:https://blog.csdn.net/syz201558503103/article/details/107238940 做医学图像处理会涉及DICOM图像,初学者应该要对图像像素有一定了 ...
- Dicom 图像--像素值(灰度值)转为CT值
摘要: 1. 概念: 2. 公式: 3. 具体情况: 4. DICOM TAG分类和说明 1. 概念: 数据:Dicom数据的原始像素值(pixel value)是没有实际意义的,所以我们常用的是Di ...
- Dicom 图像像素值(灰度值)转为CT值
CT值的单位是Hounsfield,简称为Hu,范围是-1024-3071.用于衡量人体组织对X射线的吸收率,设定水的吸收率为0Hu. 在DICOM图像读取的过程中,我们会发现图像的像素值有可能不是这 ...
- c语言实现灰度图转换为二值图
转载自:https://www.cnblogs.com/wd1001/p/4571916.html c语言实现灰度图转换为二值图 将上篇得到的灰度图转换为二值图,读取像素数据,低于某一值置0,否则设置 ...
- 小白入门计算机视觉(二) : 图像基本处理----灰度图和二值化
文章目录 解剖图像 图像处理基本原理 RGB模型 灰度 二值化 图像失真问题 从本节开始,我就要正式踏上小白的计算机视觉探索之路,先从图像基础学习吧 解剖图像 要学会图像处理首先就得知道图像的结构,平 ...
- 利用Floyd-Steinberg方法(dithering),将灰度图转换为二值图
灰度图有256级灰度,而二值图只有黑白两色.颜色数目大大降低,直观感觉转换效果不会好.其实人眼类似于一个低通滤波器,你看到的并不是一个一个像素点,而是接受的颜色信息是一个区域内的颜色信息的综合效果. ...
- python批量修改图像像素,修改命名,二值化,划分数据集
目录 修改命名 批量修改文件夹下面所有图片的像素 批量二值化 修改命名 是直接在源文件夹下面修改的 所以记得提前复制下文件再 处理. #-----------------------------批量重 ...
- RGB 转换为灰度图、二值化图
from PIL import Image I = Image.open('lena.png') I.show() L = I.convert('L') #转化为灰度图 L = I.convert(' ...
最新文章
- JSP内置对象-pageContext
- linux tar打包压缩排除某个目录或文件
- python推荐书籍-有哪些 Python 经典书籍?
- unity服务器文件传输,Unity 3D简单C#文件发送到FPT服务器示例脚本?
- android 二次绘制 layout,View的三次measure,两次layout和一次draw
- Linux多线程实践(5) --Posix信号量与互斥量解决生产者消费者问题
- 苹果修复 Gatekeeper 绕过漏洞
- ArcEngine在个人地理数据库下创建要素类
- Ubuntu18.04安装中文字体SimHei
- 将列表类型的数据转换为树形数据
- 国家级版权区块链系统 DCI标准联盟链发布 迅雷区块链提供技术支持
- Linux笔记整理(1)系统的加载和main函数执行准备
- 刷微信点击量的php,PHP一键刷QQ微信支付宝步数代码
- 下沉、重聚、归位背后,是一位学者转身业界带领阿里iDST建立技术「影响力」的四年
- Excel插件《CC办公助手》
- 计算机控制原理跟自动控制原理,上海交大考研自动控制理论和自动控制原理有什么区..._在职考研_帮考网...
- OpenFOAM动网格的壁面边界通量修正
- 大IPD之——解析华为以客户为中心的管理精髓(十)
- Virtual Painter中文版
- 北京科技大学计算机专硕就业,二本考上北京科技大学计算机专硕