转自: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值相关推荐

  1. DICOM图像像素值(灰度值)转换为CT值

    CT值的单位是Hounsfield,简称为Hu,范围是-1024-3071.用于衡量人体组织对X射线的吸收率,设定水的吸收率为0Hu. 在DICOM图像读取的过程中,我们会发现图像的像素值有可能不是这 ...

  2. 【转】DICOM图像像素值、灰度值与CT值!!!!!

    转自:https://blog.csdn.net/syz201558503103/article/details/107238940 做医学图像处理会涉及DICOM图像,初学者应该要对图像像素有一定了 ...

  3. Dicom 图像--像素值(灰度值)转为CT值

    摘要: 1. 概念: 2. 公式: 3. 具体情况: 4. DICOM TAG分类和说明 1. 概念: 数据:Dicom数据的原始像素值(pixel value)是没有实际意义的,所以我们常用的是Di ...

  4. Dicom 图像像素值(灰度值)转为CT值

    CT值的单位是Hounsfield,简称为Hu,范围是-1024-3071.用于衡量人体组织对X射线的吸收率,设定水的吸收率为0Hu. 在DICOM图像读取的过程中,我们会发现图像的像素值有可能不是这 ...

  5. c语言实现灰度图转换为二值图

    转载自:https://www.cnblogs.com/wd1001/p/4571916.html c语言实现灰度图转换为二值图 将上篇得到的灰度图转换为二值图,读取像素数据,低于某一值置0,否则设置 ...

  6. 小白入门计算机视觉(二) : 图像基本处理----灰度图和二值化

    文章目录 解剖图像 图像处理基本原理 RGB模型 灰度 二值化 图像失真问题 从本节开始,我就要正式踏上小白的计算机视觉探索之路,先从图像基础学习吧 解剖图像 要学会图像处理首先就得知道图像的结构,平 ...

  7. 利用Floyd-Steinberg方法(dithering),将灰度图转换为二值图

    灰度图有256级灰度,而二值图只有黑白两色.颜色数目大大降低,直观感觉转换效果不会好.其实人眼类似于一个低通滤波器,你看到的并不是一个一个像素点,而是接受的颜色信息是一个区域内的颜色信息的综合效果. ...

  8. python批量修改图像像素,修改命名,二值化,划分数据集

    目录 修改命名 批量修改文件夹下面所有图片的像素 批量二值化 修改命名 是直接在源文件夹下面修改的 所以记得提前复制下文件再 处理. #-----------------------------批量重 ...

  9. RGB 转换为灰度图、二值化图

    from PIL import Image I = Image.open('lena.png') I.show() L = I.convert('L') #转化为灰度图 L = I.convert(' ...

最新文章

  1. JSP内置对象-pageContext
  2. linux tar打包压缩排除某个目录或文件
  3. python推荐书籍-有哪些 Python 经典书籍?
  4. unity服务器文件传输,Unity 3D简单C#文件发送到FPT服务器示例脚本?
  5. android 二次绘制 layout,View的三次measure,两次layout和一次draw
  6. Linux多线程实践(5) --Posix信号量与互斥量解决生产者消费者问题
  7. 苹果修复 Gatekeeper 绕过漏洞
  8. ArcEngine在个人地理数据库下创建要素类
  9. Ubuntu18.04安装中文字体SimHei
  10. 将列表类型的数据转换为树形数据
  11. 国家级版权区块链系统 DCI标准联盟链发布 迅雷区块链提供技术支持
  12. Linux笔记整理(1)系统的加载和main函数执行准备
  13. 刷微信点击量的php,PHP一键刷QQ微信支付宝步数代码
  14. 下沉、重聚、归位背后,是一位学者转身业界带领阿里iDST建立技术「影响力」的四年
  15. Excel插件《CC办公助手》
  16. 计算机控制原理跟自动控制原理,上海交大考研自动控制理论和自动控制原理有什么区..._在职考研_帮考网...
  17. OpenFOAM动网格的壁面边界通量修正
  18. 大IPD之——解析华为以客户为中心的管理精髓(十)
  19. Virtual Painter中文版
  20. 北京科技大学计算机专硕就业,二本考上北京科技大学计算机专硕

热门文章

  1. 在工程中最常用的 vim使用技巧
  2. implicit request ?
  3. 计算获取最小值和最大值
  4. 次时代各制作插件使用方案以及技巧 包括UV 烘焙 减面等
  5. solr的一些查询语法
  6. P1242 新汉诺塔
  7. win10配置JAVA和tomacat环境变量
  8. app 图标规格参考表
  9. Daily Scrum 10.29
  10. 尚未注册 OLE DB 访问接口 SQLNCLI10 7043 错误