IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> 人工智能 -> dicm IMA 医学图像预处理步骤及实现 -> 正文阅读

[人工智能]dicm IMA 医学图像预处理步骤及实现

dcm|dicm是医学图像常见格式(IMA可看作dicm同类型的格式[1]),是直接能从机器上获得的图像格式,但是如果使用深度神经网络进行处理和其他操作需要对原始的图像格式进行处理才能够使用。

首先对于CT图像需要明确的是我们从机器上获得的是X ray穿过物体后投影在探测器上,探测器接收的数据是投影的数值(projection data-->sinogram(projection image) 这个过程也要进行数值转换,如取对数等操作).

dicm应该是对投影后数据前处理完毕进行重建后的图像,但是图像的像素值并不一定以Hu为单位(应该跟厂家有关,比如IMA是西门子CT设备产生的图像,不是HU单位),因此要查看dicm文件的文件头信息确认像素值的单位,如果不是HU需要处理[2]。

对于IMA文件第一步转换为HU格式,第二步进行调窗,第三步可根据情况缩放范围自行选择处理或者不处理送给神经网络或者其他操作。

(1)转换为HU格式:[3,4]中大致讲清楚了转换的方法,但是实现对于我想使用单张图像转换不太好用,因此参考代码有如下改动。

def getCtHU(dicm):
    '''直接传入dicm文件/IMA文件'''
    img=np.array(dicm.pixel_array).astype('int32')
    img[img == -2000.0] = 0
    Hu = np.array(img).astype('float64')

    RescaleIntercept = dicm.RescaleIntercept
    RescaleSlope = dicm.RescaleSlope
    if RescaleSlope != 1:
        Hu = Hu * RescaleSlope
    Hu += RescaleIntercept
    return Hu

(2)获得HU为单位的图像的矩阵后,需要针对HU图像进行调窗,以适应于我们对不同部位的图像进行分析(医学图像中不同部位的适用的HU范围不同,这个范围叫做窗口),在imageJ中窗口的大小叫做window,窗口的中心叫做level,在这几个参考的帖子中有不同的叫法,在此按照ImageJ 的名称[5,6]。

def windowsLevelTransform(Hu, window, level):
    img = Hu
    min = level - float(window) * 0.5;
    max = level + float(window) * 0.5;
    img[img < min] = min
    img[img > max] = max
    norm_ = (img - min) / window
    norm_.astype('float32')
    return norm_

调窗后的数据就是想要查看的数据范围了,在这里对不同部位调窗的数据值设定也不相同[7].


?

调窗前
调窗前
调窗后

?reference

  1. (21条消息) python如何打开.ima文件格式的医学影像(已爬坑,.IMA为西门子CT设备的另一种医学影像文件格式)_巨魔战将-CSDN博客_ima文件

  2. (21条消息) Dicom 图像--像素值(灰度值)转为CT值_小可爱的博客-CSDN博客_ct值转灰度值
  3. 成熟的医学影像“调窗”(window-leveling)算法详解 - 微笑的艾米 - 博客园 (cnblogs.com)Full Preprocessing Tutorial | Kaggle
  4. 处理医疗影像的Python利器:PyDicom - 知乎 (zhihu.com)
  5. (21条消息) 医学图像预处理(三)——windowing(ct对比增强)_normol的博客-CSDN博客_ct图像预处理
  6. CT的窗宽,窗位到底是什么意思? - 知乎 (zhihu.com)
  7. 图像识别 | 计算机断层扫描基础:加窗 - 知乎 (zhihu.com)
  8. (21条消息) 手把手教你用python实现机器学习预测疾病_程序员石磊-CSDN博客_机器学习 疾病预测
  9. 常见医疗扫描图像处理步骤 - 知乎 (zhihu.com)

  人工智能 最新文章
2022吴恩达机器学习课程——第二课(神经网
第十五章 规则学习
FixMatch: Simplifying Semi-Supervised Le
数据挖掘Java——Kmeans算法的实现
大脑皮层的分割方法
【翻译】GPT-3是如何工作的
论文笔记:TEACHTEXT: CrossModal Generaliz
python从零学(六)
详解Python 3.x 导入(import)
【答读者问27】backtrader不支持最新版本的
上一篇文章      下一篇文章      查看所有文章
加:2021-12-10 11:04:02  更:2021-12-10 11:04:17 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/27 1:36:28-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码