从三室心脏MRI影像检测主动脉瓣病变
Detecting Aortic Valve Pathology from the 3-Chamber Cine Cardiac MRI View
摘要
背景
- 心脏磁共振(CMR)是量化心脏容量、功能和血流量的金标准。
- 定制的MR脉冲序列定义了对比机制,采集几何形状和定时,可以在CMR期间应用,以实现独特的组织表征。
- 每个病人都有所有可能的获取选择是不切实际的。
在三腔(3-CH) CMR影像中定位主动脉瓣。主动脉瓣有两种主要的异常类型。狭窄:瓣膜变窄导致血液不能充分流出;不全(反流):无法阻止血液回流到左心室。
本文方法
开发并评估了一个深度学习系统,以准确分类主动脉瓣异常,为需要的患者提供进一步的定向成像。
受低级图像处理任务的启发,我们提出了一个多级网络
- 生成热图来定位主动脉瓣小叶的铰链点和主动脉狭窄或反流射流
- 在从三家NHS医院获得的临床CMR研究数据集(n = 1017例患者)上训练和评估了我们的所有模型。我们的结果(平均准确率= 0.93,f1得分= 0.91)表明,专家指导的基于深度学习的特征提取和分类模型为进一步定向成像提供了可行的策略,从而提高了CMR扫描的效率和效用。
方法
本文方法由六个主要步骤组成,如图2所示:
- 3CH CMR图像中链点的自动定位
- 主动脉瓣小叶和狭窄或反流射流引起的病理曲线的热图估计
- 估计热图中的曲线跟踪
- 每帧检测曲线的量化
- 跨帧特征总结
- 患者分类
前两步被定义为回归任务,我们训练相同网络的副本来估计这些任务中的热图。
The Proposed Network for Heat Map Regression
该网络包含三个相同的子网,对前一个子网的输出进行了细化。每个子网络都类似于U-Net,具有编码器、解码器和跳过连接。为了提高子网络对铰链点和瓣膜小叶/病理射流的敏感性,用在低级图像处理任务中表现得更好的密集块替换了前两个和最后两个卷积块
每个密集块包含5个3 × 3卷积层,然后是1 × 1卷积层。
对于铰链点定位,网络生成三个热图:每个铰链点一个热图,两个铰链点一个热图,以约束它们在单个帧中的位置。
对于曲线热图估计,网络产生4张热图;两个是主动脉瓣小叶两个是病理喷流。用损失函数L = L(≈H1, H1) + L(≈H2, H2) + L(≈H3, H3)来优化这两个网络的参数,其中L是均方误差损失。H和~ H分别表示地热图标签及其估计
从热图定位铰链点
我们通过取3CH CMR视图每帧每个铰链点对应的热图的最大值来检测估计热图中铰链点的位置。
由于我们的训练数据集中帧的有限表示,某些帧的预测可能不准确。因此,我们使用帧间铰链点的中间值在中间铰链点位置裁剪CMR图像。裁剪后的框架用于主动脉小叶和曲线状结构的热图回归,这是病理指示
Pathology Classification
使用了一种简单的跟踪方法:从最大值的位置开始跟踪热图中的脊点,直到达到停止阈值。
我们这样做是为了检测预测曲线热图中的任何潜在曲线。当热图的最大值超过初始阈值时,跟踪开始。
然后,从预测曲线中提取特征。尽管任何病理曲线的存在都足以将3CH CMR视图分类为异常,但假曲线的存在-通常是在采集过程中由人工产生的-使得这种分类不太准确。
主动脉瓣从一个框架移动到另一个框架可以获得重要的信息并揭示细微的异常
例如,临床医生会怀疑狭窄的主动脉瓣(在心脏周期内不能完全打开)
因此,我们量化每条曲线——主动脉瓣小叶和病理流——基于:
(i)它们与铰链点的接近程度:从曲线的中点到图像中心的距离,即在预测瓣膜铰链的中点;
(ii)它们相对于连接铰链点的直线的方向
(iii)其长度
(iv)是曲线的概率,能够区分真曲线和假曲线
我们通过将生成的热图视为曲线的概率图,并取跟踪曲线位置的采样概率的平均值来实现这一点(特征(iv))。考虑到每帧有四种曲线类型(两叶和两个流),我们总共分析了16个特征。这些特征(见图3)与临床医生对主动脉瓣异常的3CH CMR的解释相似。
式中,Ff表示框架f的特征,包括c1···c4四条曲线的特征,包括右冠状动脉尖叶、狭窄性血流病理曲线、非冠状动脉尖叶和反流血流病理曲线
L、A、D、P分别表示曲线长度、铰点与曲线的夹角、曲线中间到图像patch中心的距离以及沿曲线取曲线概率的平均值。
对于具有k帧CMR图像的患者,存在一个大小为R(k×16)的特征集,其中k取决于患者数据。在不正常的情况下,一张电影中的多个画面会显示主动脉异常。
随着严重程度的增加,更多的帧将包含病理曲线。因此,我们用¯F = K (F1, F2,···Fk)来总结k帧以减少缺失异常病例的可能性,其中¯F表示患者的总结特征集,通过使用总结技术k,如下所述,在k帧上
本文提出了四种特征总结方法
(i)计算跨帧特征的中位数(median)
(ii)计算跨帧特征的平均值(mean)
(iii)使用具有所有曲线概率最大值的单帧特征(SingleFrameMaxCurveP),无论曲线类型如何
(iv)使用显示每种类型的最大曲线概率的曲线特征(MaxCurveP)。除了第三种方法外,所有技术都使用多个框架来检测主动脉瓣异常患者。前两种技术假设特征是独立的,而后两种技术则保留特征与曲线的隶属关系
Imaging Data and Manual Annotation
使用两个独立的数据集进行热图回归和病理分类。对于前一项任务,我们使用了来自80名患者的1221个独特框架,其中铰链点,主动脉瓣小叶和病理流由三位心脏病专家手动注释。利用σ = 5像素的类高斯核对二值标注进行平滑处理,生成热图。
900帧用于训练,100帧用于验证,221帧用于绩效评估。我们还使用了1000个健康帧来增加曲线估计的训练集的大小。数据分割基于患者选择。在分类任务中,我们获得了1017名患者的队列,他们有二元标签来描述他们是否有正常或异常的主动脉瓣。496/1017例患者主动脉瓣异常。其中主动脉瓣狭窄184例,主动脉瓣反流222例,混合性瓣膜病变90例。
患者的平均帧数为31±15。图像分辨率在1.17 × 1.17 ~ 1.56 × 1.56像素之间
结果
从三室心脏MRI影像检测主动脉瓣病变相关推荐
- 医学影像检测方法(B超、DR、CT、MRI)
医学影像检测方法 当涉及到医学影像学时,B超(超声波检查).DR(数字X射线摄影).CT(计算机断层扫描)和MRI(磁共振成像)是常见的诊断工具.以下是对这四种影像技术的基本概念.优缺点以及它们在不同 ...
- MRI影像学习笔记(三)
欢迎扫码关注微信公众号"医影杂记" 按照解剖部位划分 颅脑 Brain MRI是最常用也是最重要的颅脑病变检查手段.与CT相比,MRI没有骨性伪影,具有更好的软组织分辨能力.并且可 ...
- MRI影像学习笔记(四)
欢迎扫码关注微信公众号"医影杂记" 前列腺Prostate MRI是前列腺增生和前列腺癌诊断的常用检查手段之一,通常使用T1和T2扫描,使用最多的扫描方向为横轴位.此外,也会选用动 ...
- MRI影像学习笔记(一)
*推荐搜索关注微信公众号:医影杂记 MRI影像检查有一个突出特点,就是有着多种多样的成像序列.这些成像序列能够产生各具特点的MRI图像,不仅能够反映人体解剖形态,而且能够反映人体血流和细胞代谢等生理功 ...
- MRI影像学习笔记(五.完结)
欢迎扫码关注微信公众号"医影杂记" 按照阅片增强功能划分 挂片协议与自动布局 Hanging Protocol and Smart Layout 即使是常规的单次MRI检查,也会包 ...
- 谈PS-based架构的线扫描(Line-scan)影像检测系统
机械视觉应用在各种产业的生产制造及品质检测已是行之有年, 利用机械视觉可以提升检测精度或加速生产速度,因此逐渐变成许多生产检测设备必备的一环. 目前市面上的影像检测系统大多采用面扫描(Area-sca ...
- MRI影像学习笔记(二)
*推荐搜索关注微信公众号:医影杂记 广义的fMRI成像包括了弥散加权成像(DWI),灌注加权成像(PWI),磁共振波普成像(MRS)和血氧饱和水平依赖成像(BOLD).而狭义的fMRI则单独指血氧饱和 ...
- [医学影像基础] 脑部MRI影像
[参考文献] 非医学专业,不保证内容准确性,此笔记仅个人学习备忘 通过调节射频脉冲序列,可以得到同一区域不同模态的影像.Repetition Time (TR)是连续射频脉冲在一个slice上作用的时 ...
- 用人工智能算法检测皮肤病变的方式训练医生
十几年来,维也纳医科大学教授的皮肤科医生哈拉尔德·基特勒一直用医学经验教授学生如何诊断皮肤病变,今年下学期开始,他将加入利用人工智能算法诊断皮肤病变的课程.这一算法系统来源于基特勒帮助组织过的一场比赛 ...
最新文章
- Win10 + Python + GPU版MXNet + VS2015 + RTools + R配置
- VTK:图像索贝尔Sobel用法实战
- angularJs自定义模块
- mysql如果数据不存在,则插入新数据,否则更新
- 在linux中查找运行程序句柄,如何查找我的进程在Linux中打开的文件句柄?
- B. File List
- 添加蜂窝注册表及永久存储
- c#.net实体类序列化方法 为什么要序列化
- 获取指定年月的全部日期放入数组
- python中使用selenium模块登录QQ邮箱
- [业余项目]黄金点游戏
- 【学习资料】Gradle Android重要的学习资料
- 倾斜摄影——3维建模软件PhotoScan教程(附安装包+教学视频)
- 64位Win10 Modelsim破解及证书LICENSE.TXT无法生成解决方法
- jsp中空格字符怎么写_jsp多个空格符号怎么打
- 肠道核心菌属——巨单胞菌属(Megamonas)
- 如何使定时器1ms进入中断一次?
- javaOOP -- 笔记
- 关于通信方面的总结(通信协议、通信端口)
- ORAN专题系列-30:5G基站如何升级到O-RAN基站 - FHGW(FrontHaul Gateway)的时钟同步系统
热门文章
- Unity Shader - ddx/ddy偏导函数测试,实现:锐化、高度图、Flat shading应用、高度生成法线
- 计算机网络进阶 ---- 网络类型 ---- 二层封装协议 ---- HDLC ---- PPP ---- pap认证 ---- chap认证 ---- 详解
- 【Windows 安装JDK8】如何安装java JDK8
- 【免杀前置课——Windows编程】五、窗口控件——什么是控件、Windolws 窗口两大类、Windows标准控件/通用控件、控件响应的接收、创建窗口制作不同控件
- 如何使用Logitech Harmony遥控器控制智能家居设备
- 算法提高 ADV-91 素数判断
- js、css引用文件的下载方式
- STC12驱动PCF8575
- RocksDB 架构
- java处理图片与base64编码互相转换