ct图像中的金属伪影校正方法

背景

随着现代医学的发展,通过守住在病人体内植入带有金属物质的假体的情况越来越普遍。最常见的是假牙植入、心脏起搏器,以及越来越多的各种关节和假肢等。这些金属物体相比较于人体组织,对光子有更强的吸收能力,会导致复杂的伪影。
产生金属伪影的产生是多重因素共同作用的结果,这些包括:
1.硬化效应
2.容积效应
3.散射效应
但是,其根本原因仍然与金属本身的高衰减特性有关,物质的高衰减会使X射线硬化,同时也会使散射现象加剧。

金属伪影的表现形式

比如如下的金属假牙ct图:

1.暗带

暗带主要是由射束硬化导致的。
已知X射线经过光电效应、康普顿散射以及电子对效应会产生衰减,其衰减的程度由Lambert-Beers定律描述。 I = I 0 e − ∫ u ( x ) d x I = I_{0}e^{-\int u(x)dx} I=I0​e−∫u(x)dx

实际上,X射线具有很宽的能谱,而大多数物体的衰减系数是与能量相关的,通常情况下对于低能光子具有高的吸收系数,对于高能光子具有地的吸收系数。当具有不同能量的光子设想同一物体时,能量降低的光子被吸收的多,设想抵达下一层结构时,高能光子比例增加,射线穿透性变强(“变硬”),这也就是被称作射束硬化的原因。

将能谱归一化:
I = I 0 ∫ E m i n E m a x p ( E ) e − ∫ u ( E , x ) d x d E I = I_{0}\int_{E_{min}}^{E_{max}}p(E)e^{-\int u(E,x)dx}dE I=I0​∫Emin​Emax​​p(E)e−∫u(E,x)dxdE

我的杯状伪影


原图

杯状伪影的效果

2.条状伪迹

解释1

金属物体的吸收系数通常是人体大部分组织吸收系数的几十倍以上,从而导致金属与组织交界位置的投影数据变化很剧烈,投影数据不连续。
FBP中的滤波器作为一个斜坡滤波器,对高频分量有放大的作用,
因此当投影数据不连续时,滤波过程会使得投影的高频谐波放大很多,从而在反投影过程中产生条状伪影。

实例:

运用matlab中的FBP函数重建之后

自写函数重建

我们取图中所标记数据


matlab函数

自写函数

我的理解


我们从sino域到图像域反投影的过程,对于图像域每个点而言是一个累加的过程,有些相邻的点可能在累加的时候,累加的值会有很大的差别,从而导致了值的差异。(可以通过我程序中的动画来理解这一点)

参考文献:
https://wenku.baidu.com/view/a8acaec190c69ec3d5bb75d0?sfr_fb=0
https://wenku.baidu.com/view/32ff4376541810a6f524ccbff121dd36a22dc4dc.html
https://wenku.baidu.com/view/6c8e753665ce05087732137d.ht…

ct图像中的金属伪影相关推荐

  1. CT图像中肺结节的精准分割 MATLAB算法源代码

    <CT图像中肺结节的精准分割> MATLAB算法源代码 资料:ppt,设计报告,源代码 任务:从给定的原始CT或超声图像数据中先手动圈定或选择包含完整目标的图像子区域,再从子区域中自动分割 ...

  2. 基于ITK-SNAP实现肺部CT图像中肺叶的标注

    时间:2019/06/12-2019/06/15 内容:肺叶图像标注,完成了10个病例的标注 材料:原始的肺部CT图像,以及采用U-net分割的肺实质的mask. 工具:itk-snap, paint ...

  3. CT图像中的窗宽和窗位

    通常我们称CT图像为DICOM(Digital Imaging and Communications in Medicine)图像,但是DICOM不是一个图像或者文件格式,它是包含了传输.存储.检索. ...

  4. 利用深度学习对医学CT图像中(LIDC-IDRI)的肺结节进行良恶性判断2

    在上篇博文中,我们详细介绍了如何分割肺实质,并根据标注信息提取肺结节,本文主要介绍如何利用CNN网络训练分类模型,辅助医生作出判断. 本文网络结构见论文: http://downloads.hinda ...

  5. (二十二:2020.11.09)论文学习之《CT中伪影的识别和规避》

    脚踏实地地解决CT的伪影问题(一)<Artifacts in CT: Recog-nition and Avoidance> 讲在前面 摘要 介绍 物理伪影 一.Beam Hardenin ...

  6. CT伪影:根本原因与规避方法

    CT伪影很常见,且原因复杂.掌握关于伪影的相关知识显得非常重要,因为它们会模仿病理学改变(例如部分容积伪影)或将图像质量降低到无法诊断的水平. 今天我们根据造成伪影的根本原因对CT伪影进行分类,并简要 ...

  7. 图像去除栅格,条纹伪影的办法(附python代码)

    医学图像中,有些时候由于使用栅格作为准束用的手段,这会导致栅格的信息叠加在图像中,形成所谓的栅格伪影,比如DR图像中的栅格伪影,如下所示: 比如像这张图像中的条纹,这就是栅格伪影.栅格伪影对于正常的图 ...

  8. acdsee扫描没有图像_详解CT图像常见伪影成因及解决方法

    CT 图像伪影 (artifact) 是指重建图像上与实际解剖结构不相符的密度异常变化.CT 图像比传统平扫X线更容易出现伪影,这是因为 CT 图像是由成千上万独立的原始测量数据重建而得,而计算机重建 ...

  9. CT图像常见伪影及解决方法

    - 前言 - CT的伪影理论上可被定义为图像中被重建数值与物体真实衰减系数之间的差异,简单来说,对于图像重建过程中不该出现在图像上的影像,可认为其是伪影(antifacts). - 01 伪影的分类  ...

最新文章

  1. HDU——1106排序(istringstream的使用、STLvector练习)
  2. elk系列7之通过grok分析apache日志
  3. struts2.xml详解
  4. Python 技术篇-socket套接字实现服务器客户端消息传递实例演示,UDP实现
  5. 由于启动用户实例的进程时出错,导致无法生成 SQL Server 的用户实例。该连接将关闭...
  6. 《ASP.NET Core 微服务实战》译者序
  7. apple手表android手机,Apple Watch 4发布了,安卓手机用户如何选择呢?
  8. html 地址 点击召唤高德,高德地图api 点聚合+海量点+点击事件(根据地区或坐标进行定位)...
  9. 《java程序设计》结对编程-四则运算整体总结
  10. TENSORFLOW 指定使用GPU跑
  11. sparklyr 1.0发布,有哪些新功能?
  12. Android service Binder用法
  13. logout退出功能是怎么实现的?login登陆功能室怎么实现的
  14. 易语言:实例讲解所有流程控制命令
  15. rpg存档修改网站_Switch上国宝级欧美版RPG游戏推荐
  16. 视频插帧算法之光流法
  17. 手机怎么在线拍照翻译英语?这个方法不能不知道
  18. 你还在用分页?试试 MyBatis 流式查询,真心强大!
  19. json schema php解析,php的json校验json-schema
  20. IntelliJ IDEA详细安装步骤

热门文章

  1. 2020届毕业生秋招 算法岗(偏NLP)总结
  2. Likecloud 吃吃吃
  3. 「3」线性代数(期末复习)
  4. Axmath符号面板放大
  5. Activity理解
  6. 深度学习Week14-yolov5s.yaml文件解读(YOLOv5)
  7. java ftl导出word文档 在wps可以打开在office打不开总结自己的一点心得
  8. matlab做能耗制动,能耗制动matlab仿真.doc
  9. Harmony OS操作系统使用体验二三事
  10. 几个可以画网络拓扑图的开源库