【笔记】DeepIGeoS: A Deep Interactive Geodesic Framework for Medical Image Segmentation
1. 简介
DeepIGeoS 是2019年发表在 TPAMI 中的一篇文章,与2016年发表在 CVPR 上的 Deep Interactive Object Selection 文章不同,本文提出的方法是在自动分割的结果上让用户对其进行标注,细化分割。
创新点:与现有的一些交互式分割方法相比,所需要的用户标注的工作量较少;响应的时间更快,每次用户交互后不再需要对全部模型进行重新训练,而是采取预训练+模拟用户交互的过程。
图 1,两阶段的交互式分割框架
DeepIGeoS 是一个两阶段的框架,使用了两个 CNN,如图1所示。Stage 1 中的 P-Net 用于初始分割(粗分割),将初始分割结果与用户的标注相结合,用作 Stage 2 中 R-Net 的输入,用于 refine 初始分割。其中 CRF-Net(f) 是我们提出的使用自由形式二元势的具有反向传播的 CRF,结合了用户交互作为强约束组成 CRF-Net(fu)。具有如下特性:
- 使用更强大的基础学习模型,即具有自动特征学习的 CNN 可以应付大量的训练集;
- 需要较少的用户交互,用户标注从高质量的自动分割图的基础上开始;
- 可以实时响应用户交互,比较高效,可以节省用户时间。
2. 方法
2.1 基于用户交互的测地线距离图
在 stage 1 初始分割图像的基础上,用户进行标记,标记分为前景和背景,相同标记的交互被转换为距离图。
测地距离图比欧式距离可以更好地反应图像中的边缘和上下文信息,提高同质区域的标签一致性
设 Sf,SbS_f, S_bSf,Sb 分别为前景和背景标记的像素集,iii 为图像 III 的一个像素点,则 iii 到标注集 S(S∈{Sf,Sb})S\left(S\in\left\{S_f,\;S_b\right\}\right)S(S∈{Sf,Sb}) 的无符号测地线距离为:
其中 Pi,jP_{i,\;j}Pi,j 为像素 iii 到 jjj 的所有路径的集合, ppp 是一条可行的路径,由 S∈[0,1]S\in\left[0,\;1\right]S∈[0,1] 进行参数化。 u(s)=p′/∥p′(s)∥\boldsymbol u\left(s\right)=p'/\left\|p'\left(s\right)\right\|u(s)=p′/∥p′(s)∥ 是与路径方向相切的单位向量。
如果没有为前景或背景绘制标记,则相应的测地距离图将填充随机数。
图 2,用户标注转换为测地距离图的示例
如图2所示,(a) 中绿色曲线为初始分割结果,红色为前景标记,青涩为背景标记,(d)、(e) 分别为转换后的基于前景和背景交互的测地线距离图像。原始图像 (a)、初始分割 (c)、测地距离 (d)(e) 图像组合在一起,得到一个四通道的图像,用以作为第二个 CNN (R-Net) 的输入。
2.2 使用扩张卷积的保持分辨率的 CNN
R-Net 改编自 VGG16 (13个卷积层、3个全连接层、5个池化层)
图 3,含有 CRF-Net(f) 的用于 2D/3D 分割 CNN 结构
以图 (a) 为例,前13个卷积层被分为5个块。第1、2块分别有两个卷积层,其余每个块都有3个卷积层。卷积层的参数在深蓝色的矩形框中(包括内核尺寸、输出通道、膨胀参数)。在所有的卷积层中,卷积核的大小都固定为3×3。块1到块6是保留分辨率的。
每个卷积层的步长设置为1。每个块中卷积的输出通道数设置为常数C。为了使用多尺度特征,作者将来自不同块的特征连接起来得到一个长度为5C的组合特征。这个特征被输入到一个分类器中,该分类器由两个附加层实现,如图(a)中的块6所示。这两个层使用大小为1×1、扩张参数为0的卷积核。块6为每个像素提供属于前景或背景类的初始分值。
图 (b) 为用于3D分割的网络结构,其与2D结构类似,只不过为了减少3D图像的内存损耗,作者在图像输入时,使用了下采样,并将块1到块5的输出特征通过1×1×1的卷积块压缩了4倍。
为了获得感受野的指数级增长,VGG16在每个块之后都使用最大池化和下采样层,这种实现方式会以指数级降低特征图的分辨率。因此,为了保持分辨率,本文移除了最大池化层和下采样层,并在每个块中使用扩张卷积。
为了在给定标注时获得更空间一致性的分割并添加硬约束,作者在块6输出的基础上应用CRF。CRF由循环神经网络实现,可以与 P-Net 或 R-Net 联合训练。CRF-Net为每个像素提供正则化预测,该预测被输入交叉熵损失函数层。
2D/3D R-Net 使用与 2D/3D P-Net 相同的结构,只是它的输入通道数为 C1+3C_1+3C1+3 (如上一小节所示),并且由于集成了用户交互,所以 CRF-Net(f) 变为了 CRF-Net(fu)。
2.3 具有自由形式的二元势函数和用户约束的反向传播 CRF-Net
CRF-Net 可以通过反向传播进行训练,本文并没有使用高斯函数,而是扩展了 CRF,使得二元势函数是自由形式的函数,称为 CRF-Net(f)。另外在交互式细化上下文中将用户的标注集成到 CRF-Net(f) 中,称为 CRF-Net(fu)。
即:CRF-Net(f) 连接到 P-Net,CRF-Net(fu) 连接到 R-Net。
设 XXX 是图像 III 中标签集 L={0,1,...,L−1}L=\left\{0,\;1,\;...,\;L-1\right\}L={0,1,...,L−1} 的标签映射,吉布斯分布(Gibbs distribution)建模 CRF 中给定 III 的 XXX 的概率:
即为描述单像素所对应的标签的概率,其中 Z(I)Z\left(I\right)Z(I) 是配分函数的归一化因子。
Eq. (6) 被称为吉布斯能量函数(Gibbs energy),由两部分组成。第一部分称为一元势函数,描述单点像素和标签的关系,本文中的一元势是从 P-Net 或 R-Net 中获得,给每个像素一个分类分(classification scores)。第二部分,二元势函数描述像素对之间的关系,相似的像素分配相同的标签,相差较大的分配不同的标签。最小化能量函数可以得到更为精确的分割。CRF 能使图片尽可能在边界处分割。
Eq. (7) 是二元势函数,其中 μ(xi,xj)\mu\left(x_i,\;x_j\right)μ(xi,xj) 表示 i,ji, ji,j 标签之间的兼容性,由大小为 L×LL×LL×L 的矩阵表示。f(⋅)f\left(\cdot\right)f(⋅) 是关于特征差异 f~i,j=fi−fj{\widetilde{\boldsymbol f}}_{i,\;j}={\boldsymbol f}_i\boldsymbol-{\boldsymbol f}_{\mathbf j}fi,j=fi−fj (其中fi,fj{\boldsymbol f_i}, {\boldsymbol f_j}fi,fj 分别表示为 i,ji, ji,j 的特征向量)和欧式距离 di,jd_{i, j}di,j 的函数,本文并没有将其定义为高斯函数或多个高斯函数的组合,而是将其设置为由全连接神经网络(Pairwise-Net)表示的自由函数,如下图4所示,它可以在训练中学习,有两个隐藏层和一个输出层。
图 4,Pairwise-Net 结构
由于二元势规模越大,导致计算量就越大,所以将复杂的后验分布简化为 QQQ ,在最小散度(KL-divergence)D(Q∥p)D\left(Q\parallel p\right)D(Q∥p) 下对 QQQ 迭代更新,得到原始后验概率分布的近似解。
其中 LLL 是标签集,i,ji, ji,j 是像素对。
对于提出的 CRF-Net(fu),在用户提供的标注集 Sfb=Sf∪SbS_{fb}=S_f\cup S_bSfb=Sf∪Sb 的情况下,强制标注集中像素的概率为1或0,每次迭代的更新规则为:
其中 sis_isi 表示用户提供的像素 iii 的标签,它位于 SfbS_{fb}Sfb 中。
本文使用 CRF-Net 中的实现,通过 RNN 中的多阶段平均场(平均近似场通常用于提升 CRF 的 inference 效率,同时允许通过反向传播学习参数 )法更新 QQQ。每个平均场分裂方程(Eq. 8)分为4个步骤,包括消息传递、兼容性转换、添加一元势和归一化。
2.4 实施细节
2.4.1 测地线距离的变换
光栅扫描算法用于计算测地线距离变换,方法是应用前向扫描和后向扫描,2D使用3×3内核,3D使用3×3×3内核。由于访问连续块中的图像内存,所以速度很快。对于本文提出的具有自由形势的二元势的 CRF-Net,没有使用整个图像中的全连接,而是使用基于局部块的像素连接。原因如下:
- 仅当二元势函数为高斯函数时,效率才会高
- 长距离依存(long-distance dependency)处理医学影像难度大。由于医学影像的对比度通常都很低,在这种情况下,长距离依存可能会导致目标像素的标签被大量具有相似外观的像素所影响
因此,为了保持良好的效率并避免长距离损失,本文定义了以该像素为中心的局部块内的一个像素的成对连接。2D图像的块大小设置为7×7,3D图像的块大小为5×5×3。
2.4.2 自由形式的二元势函数
Eq. (7) 中的 μ(⋅)\mu\left(\cdot\right)μ(⋅) 初始化为 μ(xi,xj)=[xi≠xj]\mu\left(x_i,\;x_j\right)=\left[x_i\neq x_j\right]μ(xi,xj)=[xi=xj],其中 [⋅][\cdot][⋅] 为艾佛森括号(Iverson bracket),满足条件则为1,否则为0。
图4中,第1和第2个隐藏层分别有32和16个神经元。该 Pairwise-Net 是由全卷积神经网络实现的,2D为1×1内核,3D为1×1×1内核。在预训练中使用对比度敏感函数(contrast sensitive function)的近似值初始化 Pairwise-Net:
其中 FFF 是特征向量 fif_ifi 和 fjf_jfj 的维度,ω,σ\omega,\sigmaω,σ 是分别控制函数大小和形状的参数。根据经验分别设置为0.5和0.8。
将 fi,fjf_i, f_jfi,fj 设置为 P-Net 的输入通道(即本文中的图像强度),用以简化实现并获得对比度敏感的二元势能(contrast-sensitive pairwise potentials)。
2.4.3 预训练 Pairwise-Net
生成一个包含100k样本的训练集 T′={X′,Y′}T'=\left\{X',\;Y'\right\}T′={X′,Y′}
其中 X′X'X′ 是模拟连接 f~ij{\widetilde{\boldsymbol f}}_{ij}fij 和 dijd_ijdij 的特征集,Y′Y'Y′ 是模拟 Eq. (11) 的一组预测值集。对于 T′T'T′ 中的每个样本 sss ,特征向量 xs′x'_sxs′ 的维度为 F+1F+1F+1 ,其中前 FFF 个维度表示 f~ij{\widetilde{\boldsymbol f}}_{ij}fij 的值,最后一个维度表示 dijd_{ij}dij 。xs′x'_sxs′ 的第 ccc 个通道用随机数 k′k'k′ 填充,k′∼N(0,2)k'\sim N\left(0,\;2\right)k′∼N(0,2) 表示 c⩽Fc\leqslant Fc⩽F ,k′∼U(0,8)k'\sim U\left(0, \;8\right)k′∼U(0,8) 表示 c=F+1c=F+1c=F+1 。
xs′x'_sxs′ 的预测值 ys′y'_sys′ 的 ground truth 从 Eq. (11) 得到。在生成 X′X'X′ 和 Y′Y'Y′ 后用具有二次损失函数的随机梯度下降(SGD)算法来预训练 Pairwise-Net。
2.4.4 预处理
所有图像都通过训练集的平均值和标准差进行归一化。
通过垂直或水平翻转、角度范围为 [−
【笔记】DeepIGeoS: A Deep Interactive Geodesic Framework for Medical Image Segmentation相关推荐
- 【阅读笔记】Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation...
Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation 作者:Lin Yan ...
- Deep Poincare Map for Robost Medical Image Segmentation 论文阅读
Deep Poincare Map for Robost Medical Image Segmentation于17年3月份发表在arXiv上,作者来自帝国理工学院. 这篇文章提出了一种半自动方法(需 ...
- 阅读论文:Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation
阅读论文:Suggestive Annotation: A Deep Active Learning Framework forBiomedical Image Segmentation 本篇论文发表 ...
- 【论文笔记】Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation
一 概述 1.本文研究的问题为: 在有限的标注工作中,哪些样例应当被标注,以此可以实现最好的识别效果. 2.解决方法: 设计了 一种结合了FCN和主动学习的深度主动学习框架 ,通过在最有效的注 ...
- [论文阅读] Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation
论文地址:https://arxiv.org/abs/1706.04737 发表于:MICCAI'17 Abstract 图像分割是生物医学图像分析的一个基本问题.深度学习的最新进展在许多生物医学图像 ...
- Suggestive Annotation: A Deep Active Learning Framework for Biomedical Image Segmentation
作者:LinYang 1 , Yizhe Zhang 1 , Jianxu Chen 1 , Siyuan Zhang 2 , Danny Z. Chen 1 1Departmentof Comput ...
- [深度学习论文笔记]TransBTSV2: Wider Instead of Deeper Transformer for Medical Image Segmentation
TransBTSV2: Wider Instead of Deeper Transformer for Medical Image Segmentation TransBTSV2:用于医学图像分割的宽 ...
- A Reinforcement Learning Framework for Medical Image Segmentation翻译
摘要 本文介绍了一种使用强化学习方案进行医学图像分割的新方法. 我们使用这个新颖的想法作为一种有效的方法来最佳地找到合适的局部阈值和结构元素值,并在超声图像中分割前列腺. 强化学习代理使用超声图像及其 ...
- 【笔记】An explainable deep machine vision framework for plant stress phenotyping
<An explainable deep machine vision framework for plant stress phenotyping> 背景 目前基于视觉症状的植物胁迫,由 ...
最新文章
- 【前端帮帮忙】第7期 关于节流(throttle)和防抖(debounce)的理解
- html制作棋盘,棋盘的制作方法
- jquery and event
- 1.17 项目实例:模仿斗地主洗牌发牌小游戏
- python生成表格文件_python读取和生成excel文件
- sr650服务器cpu型号,至强Gold 联想ThinkSystem SR650评测
- Qt|C++工作笔记-对虚函数的进一步认识(关键是Qt中如何正确使用,达到想要的效果)
- 杭电1710 (已知二叉树前中序 求后序)
- 启动计算机时页面配置出现问题,开机提示“由于启动计算机时出现了页面配置问题…”...
- 2021-2025年中国船用炉灶行业市场供需与战略研究报告
- sqlite3怎么筛选一个月之前的_怎么查询指定日期 和 这个日期 之前一个月前的数据...
- 桌面魔术可以用透明桌子揭密
- MySQL如何使用索引 较为详细的分析和例子
- XP不能访问WIN7资源
- CentOS升级pcre
- Zookeeper Leader选举算法及选举过程
- 如何推动数字经济新时代产业转型升级
- Python爬虫-中华英才网登陆JS加密登陆
- Arduino使用避障传感器实验
- ALLHiC: 辅助组装简单的二倍体基因组