HU(Hounsfiled Unit)值,反映了组织对X射线吸收程度。以水的吸收程度作为参考,即水的HU=0,衰减系数大于水的为正直,小于水的为负值。并以骨皮质和空气的HU值为上限和下限。

最近遇到了一个问题,为了站在医生的角度看问题,所以必须将dcm图像的值转换为HU值。在网上也搜了一些帖子,大家都提到了一个公式:
HU = pixel_val * slope + intercept

import dicom
import numpy as np
import SimpleITK as sitkdcm_path = '图片路径'
img = dicom.read_file(dcm_path)
img_array = sitk.GetArrayFromImage(sitk.ReadImage(dcm_path))
HU = np.dot(img_array, img.RescaleSlope) + img.RescaleIntercept

但是,这里会有一个问题,有的图像就已经是CT值(HU值),该图像的Solpe=1,Intercept=0,如果按照上面的公式计算,就还是图像的值;而有的图像则是像素值,所以需要转化,该图像的Sople=1, Intercept=-1024(这是我用LUNA数据得出的结果)。

如有错误,欢迎指正

Dicom图像值与HU值相关推荐

  1. CT医学影像的窗高窗位、CT值(Hu值)

    CT医学影像的窗高窗位.CT值(Hu值) 1.CT医学影像的窗高窗位.CT值(Hu值)学习 CT扫描层面 冠状位. 矢状位.横断位的基本定义及部分正常人体组织的CT值(Hu)[CT扫描层面 冠状位. ...

  2. 医学图像——CT值(Hu值)

    1.CT值定义 CT值的是测定人体某一局部组织或器官密度大小的一种[计量单位],通常称亨氏单位(hounsfield unit ,HU).CT值又叫HU值,HU(Hounsfiled Unit)值,反 ...

  3. (二)CT医学影像的窗高窗位、CT值(Hu值)

    一.简介 1.1. 本文介绍医学影像CT的DICOM文件的窗宽(Window Width,简写WW)窗位(Window Level,简写WL)与CT值(HU)的关系. 比如 WW = 400,WL = ...

  4. CT图像之Hu值变换与窗宽窗位调整

    今天通过分析实验结果,完善了之前的预处理代码,借此博文分享给大家,另外 点击此处可以查看预处理的完整代码,欢迎大家一起来完善. 最重要的是,欢迎大家的批评指正,您的建议和意见将会是我成长源泉. 接下来 ...

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

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

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

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

  7. 利用靶区勾画的RT struct文件分析相关联的Dicom图像中肿瘤部分的CT(HU)值分布(升级版)

    利用靶区勾画的RT struct文件分析相关联的Dicom图像中肿瘤部分的CT(HU)值分布(升级版) % ============================================== ...

  8. 医学图像的CT值与像素值总结及转换代码

    目录 一.CT图像的调窗 1.Window width 2.Window level/center 二.DICOM文件中窗宽窗位对应字段 三.CT值与像素值转换(线性映射) 1.itk-snap软件和 ...

  9. python读取dicom序列_python读取dicom图像(SimpleITK和dicom包实现)

    1. 用SimpleITK读取dicom序列: import SimpleITK as sitk import numpy as np img_path='F:\\dataset\\pancreas\ ...

最新文章

  1. PHP 取302跳转后真实 URL 的两种方法
  2. pyrealsense2 sensor.get_option_description 报错 RuntimeError: object doesn't support option #14
  3. [Swift]LeetCode463. 岛屿的周长 | Island Perimeter
  4. 详解数字电视机顶盒的功能技术与应用3
  5. css sprite 介绍和实例
  6. SpaceX上天容易,Tesla自动驾驶难 | 凌云时刻
  7. ffmpeg 推流命令记载
  8. 物业服务的品质、成本与四保一服的数字化方法
  9. 2013级C++第15周项目——二维数组
  10. 数学杂谈:高维空间向量夹角小记
  11. mysql捕获1300的错误_pt-osc 变更时遇到 “MySQL error 1300” 报错问题解决
  12. 搭建一个专属于两个人的爱情网站,记录生活中的点点滴滴
  13. sql统计各分数段人数示例
  14. 中国网络安全现状:“狼”已经来了
  15. 零基础入门学习Python--永久存储:腌制一缸美味的泡菜
  16. java中的画图公式_GitHub - oneSaber/Compiler-java: 函数绘图语言java实现
  17. 在Linux上yum安装snmp,centos7配置安装snmp
  18. 重庆大学计算机学院课题组,【计算机】计算机学院关于智能计算的大规模优化学术报告圆满结束...
  19. python基础试题(1)
  20. 什么是概念模型?概念模型的作用是什么?

热门文章

  1. 《率土之滨》高流水的背后秘诀:商业模式与系统分析
  2. Python初学者-双人对战小游戏
  3. java 时间格式化 上午下午 am pm,本地windows显示的格式和linux服务器显示格式,风格不一致问题
  4. jackson如何创建json对象和jackson工具类
  5. 网上赚钱的门路方法什么比较靠谱网上做兼职是真的吗?亲身经历告诉你真相!
  6. 自然资源的分类、特点、与经济发展的关系
  7. 美国可以实现几乎对全球无死角的监控,这是可能的吗?
  8. openGL-梁友栋-Barsky算法
  9. java教程pdf文本文档版
  10. 一份标准的STM32工程模板都需要哪些文件?(B站江科大自化协)