【论文笔记】AVSM:结合了仿射配准和vSVF配准的医学图像配准模型
本文是论文《Networks for Joint Affine and Non-parametric Image Registration》的阅读笔记。
文章提出了一个名为AVSM(Affine-vSVF-Mapping)的端到端3D医学图像配准模型,该模型结合了仿射配准和矢量动量参数化静态速度场模型( a vector momentum-parameterized stationary velocity field(vSVF) model)。模型有三个阶段(stage),在第一阶段,通过多步的仿射网络来预测一个仿射变换参数;在第二阶段,使用类UNet网络生成一个动量,从而可以通过平滑计算出一个速度场;在第三阶段,使用一个独立的基于图的vSVF组件来基于当前转换映射做无参的微调。
一、相关工作
当图像具有相似的灰度值分布时,通常采用MSE来评估灰度值的相似性;而在多模配准中,NCC和MI指标会更合适。
微分同胚转换可以通过优化足够平滑的速度场来实现,此时,空间转换可以通过积分来恢复(通过积分得到逆转换)。这样的方法包括LDDMM和Diffeomorphic Demons等。
由于空间转换网络(STN)在配准模型中的使用,实现了模型的端到端的训练。
现有的基于深度学习的配准模型存在诸多局限:
- 这些方法通常会假定图像已经实现预配准过,即进行过刚性或仿射配准。这种预对齐可以通过特定的网络或标准数值优化的方式来进行,在前一种方法中,配准方法不再是端到端的,在后一种方法中,预配准成为了计算的瓶颈;
- 很多方法受限于内存/显存的大小,通常可以通过改为对2D图像会对3D patch进行配准以减少资源开销;
- 这些方法没有进行迭代的改进(refinement)。
本文的贡献有:
- 提出了一个新颖的矢量动量参数化的静态速度场配准模型,矢量动量长允许将平滑和转换参数的预测解构,因此可以得到足够平滑的速度场以及适用于大位移的微分同胚;
- 模型是端到端的,并且结合了仿射配准和vSVF配准;
- 在仿射配准和vSVF配准组件中都采用了多步配准的方法,从而改善了配准的效果;
- 结果只经过一次配准映射得到,避免了多次插值而带来的误差;
- 在仿射配准和vSVF配准组件中都采用了逆一致性损失,从而让模型忽略输入图像的顺序,即从A配准到B和从B配准到A的结果相似。
AVSM的网络架构如上图所示,在仿射配准阶段,采用多步的仿射网络来预测仿射变换参数;在vSVF配准阶段,使用类UNet网络产生一个动量,并由此计算一个速度场,然后初始映射和动量被喂到vSVF组件中来产生最终的转换映射。在以上过程中源图像只做一次插值,由此也避免了图像模糊。
在处理3D图像时,如果图像的大小只有原图的一半,则需要的计算和内存/显存只有原来的18\frac{1}{8}81。
二、方法
1. 多步仿射网络
现有的大多数无参数配准方法由于有正则项惩罚,所以对仿射变换不是不变的,因此在无参配准之前通常会先进行预配准,即仿射配准,以考虑全局的大位移和旋转。
模型应该足够灵活以同时适应小的和大的仿射变形。所以采用多步的方法得到最终的仿射变换,这种策略显著的提高了模型精度和稳定性。
上图是多步仿射配准的网络结构。多步仿射配准网络是一个多次出现的网络,以逐渐的改进仿射变换。为了避免多次三线性插值带来的数值不稳定性和数值消散,文章采用的是直接更新仿射参数而不是在中间步骤对图像进行重采样。用Γ=(Ab)\Gamma=(A \quad b)Γ=(Ab)来表示仿射参数,其中A∈Rd×dA\in\R^{d\times d}A∈Rd×d表示线性转换矩阵,b∈Rdb\in\R^{d}b∈Rd表示位移,ddd是图像的维度。仿射参数更新的规则如下:
A(t)=A~(t)A(t−1),b(t)=A~(t)b(t−1)+b~(t)s.t.A(0)=I,b(0)=0\begin{array}{l} A_{(t)}=\tilde{A}_{(t)} A_{(t-1)}, b_{(t)}=\tilde{A}_{(t)} b_{(t-1)}+\tilde{b}_{(t)} \\ \text {s.t.} \quad A_{(0)}=I, b_{(0)}=0 \end{array} A(t)=A~(t)A(t−1),b(t)=A~(t)b(t−1)+b~(t)s.t.A(0)=I,b(0)=0
其中,A~(t),A(t)\tilde A_{(t)},A_{(t)}A~(t),A(t)表示第ttt步的线性变换矩阵输出和合成结果。相似的,$ \tilde b_{(t)},b_{(t)}表示第表示第表示第t步的仿射变换参数输出和合成结果。仿射映射可以通过下式得到:步的仿射变换参数输出和合成结果。仿射映射可以通过下式得到:步的仿射变换参数输出和合成结果。仿射映射可以通过下式得到:\Phi_{a}^{-1}(x, \Gamma)=A_{\left(t_{\text {last}}\right)} x+b_{\left(t_{\text {last}}\right)}$。
2. 多步仿射网络的损失
损失包括三部分:图像相似性损失La−simL_{a-sim}La−sim,正则项La−regL_{a-reg}La−reg和鼓励变换对称的损失La−symL_{a-sym}La−sym。让I0I_0I0和I1I_1I1分别表示源图像和目标图像,上标st^{st}st和ts^{ts}ts分别表示从I0I_0I0到I1I_1I1的配准和从I1I_1I1到I0I_0I0的配准。
1)图像相似性损失La−sim(I0,I1,Φa−1)L_{a-s i m}\left(I_{0}, I_{1}, \Phi_{a}^{-1}\right)La−sim(I0,I1,Φa−1)
可以是NCC、LNCC(localized NCC)、MSE等。本文使用的是mk-LNCC(multi-kernel LNCC)。让VVV 表示图像的体积,xi,yix_i,y_ixi,yi表示配准后的源图像和目标图像中的第iii个体素,NsN_sNs表示大小为s×s×ss\times s\times ss×s×s的滑动窗口的个数,ζjs\zeta_{j}^{s}ζjs表示中心点在第jjj个体素的窗口,xjˉ,yjˉ\bar{x_j},\bar{y_j}xjˉ,yjˉ表示配准后源图像和目标图像在窗口ζj8\zeta^8_jζj8内的灰度值的平均值。窗口大小为sss的LNCC记为ksk_sks,其定义如下:
κs(x,y)=1Ns∑j∑i∈ζjs(xi−xˉj)(yi−yˉj)∑i∈ζjs(xi−xˉj)2∑i∈ζjs(yi−yˉj)2\kappa_{s}(x, y)=\frac{1}{N_{s}} \sum_{j} \frac{\sum_{i \in \zeta_{j}^{s}}\left(x_{i}-\bar{x}_{j}\right)\left(y_{i}-\bar{y}_{j}\right)}{\sqrt{\sum_{i \in \zeta_{j}^{s}}\left(x_{i}-\bar{x}_{j}\right)^{2} \sum_{i \in \zeta_{j}^{s}}\left(y_{i}-\bar{y}_{j}\right)^{2}}} κs(x,y)=Ns1j∑∑i∈ζjs(xi−xˉj)2∑i∈ζjs(yi−yˉj)2∑i∈ζjs(xi−xˉj)(yi−yˉj)
定义mk-LNCC为具有不同窗口大小的LNCC的加权和。
图像相似性损失为:
La−sim(I0,I1,Γ)=∑iωiκsi(I0∘Φa−1,I1)s.t. Φa−1(x,Γ)=Ax+band ∑iωi=1,wi≥0\begin{array}{l} L_{a-s i m}\left(I_{0}, I_{1}, \Gamma\right)=\sum_{i} \omega_{i} \kappa_{s_{i}}\left(I_{0} \circ \Phi_{a}^{-1}, I_{1}\right) \\ \text { s.t. } \Phi_{a}^{-1}(x, \Gamma)=A x+b \text { and } \sum_{i} \omega_{i}=1, w_{i} \geq 0 \end{array} La−sim(I0,I1,Γ)=∑iωiκsi(I0∘Φa−1,I1) s.t. Φa−1(x,Γ)=Ax+b and ∑iωi=1,wi≥0
2)正则项损失La−reg(Γ)L_{a-reg}(\Gamma)La−reg(Γ)
惩罚合成仿射变换与恒等变换的偏差:
La−reg(Γ)=λar(∥A−I∥F2+∥b∥22)L_{a-r e g}(\Gamma)=\lambda_{a r}\left(\|A-I\|_{F}^{2}+\|b\|_{2}^{2}\right) La−reg(Γ)=λar(∥A−I∥F2+∥b∥22)
其中∥⋅∥F\|\cdot\|_F∥⋅∥F表示Frobenius正则,λar≥0\lambda_{ar}\geq0λar≥0是依赖于epoch的权重系数,在训练开始时较大,以限制大的形变,然后逐渐衰减到0。
3)对称损失La−sym(Γ,Γts)L_{a-sym}(\Gamma,\Gamma^{ts})La−sym(Γ,Γts)
鼓励配准的逆一致性,即鼓励从源图像到目标图像的变换是从目标图像到源图像变换的逆变换,用公式表示为:Ats(Ax+b)+bts=x):La−sym(Γ,Γts)=λas(∥AtsA−I∥F2+∥Atsb+bts∥22)\left.A^{t s}(A x+b)+b^{t s}=x\right):L_{a-s y m}\left(\Gamma, \Gamma^{t s}\right)=\lambda_{a s}\left(\left\|A^{t s} A-I\right\|_{F}^{2}+\left\|A^{t s} b+b^{t s}\right\|_{2}^{2}\right)Ats(Ax+b)+bts=x):La−sym(Γ,Γts)=λas(∥AtsA−I∥F2+∥Atsb+bts∥22),其中λas≥0\lambda_{as}\geq0λas≥0。
总损失为:
La(I0,I1,Γ,Γts)=ℓa(I0,I1,Γ)+ℓa(I1,I0,Γts)+La−sym(Γ,Γts)\begin{aligned} \mathcal{L}_{a}\left(I_{0}, I_{1}, \Gamma, \Gamma^{t s}\right)=& \ell_{a}\left(I_{0}, I_{1}, \Gamma\right)+\ell_{a}\left(I_{1}, I_{0}, \Gamma^{t s}\right) \\ &+L_{a-s y m}\left(\Gamma, \Gamma^{t s}\right) \end{aligned} La(I0,I1,Γ,Γts)=ℓa(I0,I1,Γ)+ℓa(I1,I0,Γts)+La−sym(Γ,Γts)
其中ℓa(I0,I1,Γ)=La−sim(I0,I1,Γ)+La−reg(Γ)\ell_{a}\left(I_{0}, I_{1}, \Gamma\right)=L_{a-s i m}\left(I_{0}, I_{1}, \Gamma\right)+L_{a-r e g}(\Gamma)ℓa(I0,I1,Γ)=La−sim(I0,I1,Γ)+La−reg(Γ)。
3. 矢量动量参数化SVF
1)vSVF方法
为了捕获大的变形并保证微分同胚变换,通常使用流体机制的配准算法。此时,转换映射ϕ3\phi^3ϕ3可以通过速度场v(x,t)v(x,t)v(x,t)对时间积分来得到。控制微分方程为:Φt(x,t)=v(Φ(x,t),t),Φ(x,0)=Φ(0)(x)\Phi_{t}(x, t)=v(\Phi(x, t), t), \Phi(x, 0)=\Phi_{(0)}(x)Φt(x,t)=v(Φ(x,t),t),Φ(x,0)=Φ(0)(x),其中Φ(0)\Phi_{(0)}Φ(0)是初始映射。对于足够平滑的速度场来说,可以获得一个微分同胚变换,通过对不平衡进行惩罚而保证速度场的平滑性,即:
v∗=argminvλvr∫01∥v∥L2dt+sim[I0∘Φ−1(1),I1]s.t. Φt−1+DΦ−1v=0and Φ−1(0)=Φ(0)−1\begin{array}{l} v^{*}=\underset{v}{\operatorname{argmin}} \lambda_{v r} \int_{0}^{1}\|v\|_{L}^{2} \mathrm{d} t+\operatorname{sim}\left[I_{0} \circ \Phi^{-1}(1), I_{1}\right] \\ \text { s.t. } \quad \Phi_{t}^{-1}+D \Phi^{-1} v=0 \quad \text { and } \quad \Phi^{-1}(0)=\Phi_{(0)}^{-1} \end{array} v∗=vargminλvr∫01∥v∥L2dt+sim[I0∘Φ−1(1),I1] s.t. Φt−1+DΦ−1v=0 and Φ−1(0)=Φ(0)−1
其中DDD表示雅克比行列式,∥v∥L2=⟨L†Lv,v⟩\|v\|^2_L=\left\langle L^{\dagger} L v, v\right\rangle∥v∥L2=⟨L†Lv,v⟩是通过微分操作LLL和其共轭L†L^\daggerL†所定义的空间正则。由于动量m=L†Lvm=L^\dagger Lvm=L†Lv,正则又可以表示为:∥v∥L2=⟨m,v⟩\|v\|_{L}^{2}=\langle m, v\rangle∥v∥L2=⟨m,v⟩。SVF配准算法是直接优化速度场vvv,而本文提出一个矢量动量SVF(vSVF)方程,其表达式如下:
m∗=argminm0λvr⟨m0,v0⟩+Sim[I0∘Φ−1(1),I1],s.t. Φt−1+DΦ−1v=0,Φ−1(0)=Φ(∩)−1,v0=(L†L)−1m0\begin{array}{l} m^{*}=\underset{m_{0}}{\operatorname{argmin}} \lambda_{v r}\left\langle m_{0}, v_{0}\right\rangle+\operatorname{Sim}\left[I_{0} \circ \Phi^{-1}(1), I_{1}\right], \text { s.t. } \\ \Phi_{t}^{-1}+D \Phi^{-1} v=0, \Phi^{-1}(0)=\Phi_{(\cap)}^{-1}, v_{0}=\left(L^{\dagger} L\right)^{-1} m_{0} \end{array} m∗=m0argminλvr⟨m0,v0⟩+Sim[I0∘Φ−1(1),I1], s.t. Φt−1+DΦ−1v=0,Φ−1(0)=Φ(∩)−1,v0=(L†L)−1m0
其中,m0m_0m0表示矢量动量,λvr>0\lambda_{vr}>0λvr>0是常数,该方程可以看作是简化的矢量动量参数化的LDDMM方程。该方程的有点事它允许我们通过深度网络来预测动量而显式的控制空间平滑,进而产生平滑的速度场,而不是直接预测速平滑的度场。
2)vSVF配准网络
上图是vSVF配准网络框架,分为两部分:一是动量生成网络,以配准后的源图像和目标图像作为输入,输出低分辨率的动量;二是vSVF配准部分。预测的动量和下采样的初始化映射被输入到vSVF单元,然后输出最后经过上采样得到全分辨率的转换映射。在vSVF单元中,通过对动量进行平滑而得到速度场,然后用它来解对流方程Φ(τ)t−1+DΦ(τ)−1v=0\Phi_{(\tau) t}^{-1}+D \Phi_{(\tau)}^{-1} v=0Φ(τ)t−1+DΦ(τ)−1v=0,然后得到最终的变换映射。初始映射可以是仿射映射,也可以是通过之前的vSVF步骤得到的映射。
动量生成网络
使用四阶Runge-Kutta方法来对时间积分并将所有具有中心差分的空间导数离散化。因此,减少了内存的需要,网络输出的是低分辨率的动量。在实操时,移除了UNet中解码器的最后一层,以产生低分辨率的动量。
vSVF的损失
损失也分为三部分,相似性损失和仿射网络中的相同,正则项损失用来乘法速度场,即:
Lv−reg(m0)=λvr∥v∥L2=λvr⟨m0,v0⟩L_{v-r e g}\left(m_{0}\right)=\lambda_{v r}\|v\|_{L}^{2}=\lambda_{v r}\left\langle m_{0}, v_{0}\right\rangle Lv−reg(m0)=λvr∥v∥L2=λvr⟨m0,v0⟩
其中,v0=(L†L)−1m0v_0=(L^\dagger L)^{-1}m_0v0=(L†L)−1m0,把(L†L)−1(L^\dagger L)^{-1}(L†L)−1看作是具有多高斯核的卷积。
对称性损失
Lv−sym(Φ−1,(Φts)−1)=λvs∥Φ−1∘(Φts)−1−id∥22L_{v-s y m}\left(\Phi^{-1},\left(\Phi^{t s}\right)^{-1}\right)=\lambda_{v s}\left\|\Phi^{-1} \circ\left(\Phi^{t s}\right)^{-1}-i d\right\|_{2}^{2} Lv−sym(Φ−1,(Φts)−1)=λvs∥∥∥Φ−1∘(Φts)−1−id∥∥∥22
其中ididid表示恒等映射,λvs≥0\lambda_{vs}\geq0λvs≥0表示对称权重,(Φts)−1(\Phi^{ts})^{-1}(Φts)−1表示从目标图像到源图像的配准映射,Φ−1\Phi^{-1}Φ−1表示从源图像到目标图像的配准映射。
总损失:
Lv(I0,I1,Φ−1,(Φts)−1,m0,m0ts)=ℓv(I0,I1,Φ−1,m0)+ℓv(I1,I0,(Φts)−1,m0ts)+Lv−sym(Φ−1,(Φts)−1)\begin{aligned} \mathcal{L}_{v}\left(I_{0}, I_{1}, \Phi^{-1},\left(\Phi^{t s}\right)^{-1}, m_{0}, m_{0}^{t s}\right)=\ell_{v}\left(I_{0}, I_{1}, \Phi^{-1}, m_{0}\right) \\ +\ell_{v}\left(I_{1}, I_{0},\left(\Phi^{t s}\right)^{-1}, m_{0}^{t s}\right) \\ +L_{v-s y m}\left(\Phi^{-1},\left(\Phi^{t s}\right)^{-1}\right) \end{aligned} Lv(I0,I1,Φ−1,(Φts)−1,m0,m0ts)=ℓv(I0,I1,Φ−1,m0)+ℓv(I1,I0,(Φts)−1,m0ts)+Lv−sym(Φ−1,(Φts)−1)
其中,ℓv(I0,I1,Φ−1,m0)=Lv−sim(I0,I1,Φ−1)+Lv−reg(m0)\ell_{v}\left(I_{0}, I_{1}, \Phi^{-1}, m_{0}\right)=L_{v-\operatorname{sim}}\left(I_{0}, I_{1}, \Phi^{-1}\right)+L_{v-r e g}\left(m_{0}\right)ℓv(I0,I1,Φ−1,m0)=Lv−sim(I0,I1,Φ−1)+Lv−reg(m0)
具有T步的vSVF模型,其总损失为:
∑τ=1TLv(I0,I1,Φ(τ)−1,Φ(τ)ts,m0(τ),m0(τ)ts)s.t.t. Φ(τ)−1(x,0)=Φ(τ−1)−1(x,1)(Φ(τ)ts)−1(x,0)=(Φ(τ−1)ts)−1(x,1)\begin{array}{c} \sum_{\tau=1}^{T} \mathcal{L}_{v}\left(I_{0}, I_{1}, \Phi_{(\tau)}^{-1}, \Phi_{(\tau)}^{t s}, m_{0}(\tau), m_{0(\tau)}^{t s}\right) \quad \text { s.t.t. } \\ \Phi_{(\tau)}^{-1}(x, 0)=\Phi_{(\tau-1)}^{-1}(x, 1) \\ \left(\Phi_{(\tau)}^{t s}\right)^{-1}(x, 0)=\left(\Phi_{(\tau-1)}^{t s}\right)^{-1}(x, 1) \end{array} ∑τ=1TLv(I0,I1,Φ(τ)−1,Φ(τ)ts,m0(τ),m0(τ)ts) s.t.t. Φ(τ)−1(x,0)=Φ(τ−1)−1(x,1)(Φ(τ)ts)−1(x,0)=(Φ(τ−1)ts)−1(x,1)
三、训练
1. 多步仿射网络的训练
首先训练一个单步的网络,然后用它的参数来初始化多步的网络。在纵向的配准中,训练一个3步的仿射网络,而在测试的时候使用七步的网络;在跨物体配准中,训练一个5步的仿射网络,在测试的时候使用七步的网络。仿射对称因子λas\lambda_{as}λas设为10。仿射正则因子是依赖于epoch的,其定义为:
λar:=CarKarKar+en/Kar\lambda_{a r}:=\frac{C_{a r} K_{a r}}{K_{a r}+e^{n / K_{a r}}} λar:=Kar+en/KarCarKar
其中,CarC_{ar}Car是常数,KarK_{ar}Kar用来控制衰减率,nnn是epoch数,KarK_{ar}Kar设为4,CarC_{ar}Car设为10。
2. 动量生成网络的训练
使用10步和多高斯核(标准差为{0.05, 0.1, 0.15, 0.2, 0.25})以及相关权重{j0.067, 0.133, 0.2, 0.267, 0.333}。在训练时都是训练两步,正则化因子λvr\lambda_{vr}λvr设为10,对称因子λvs\lambda_{vs}λvs设为1e−41e^{-4}1e−4。每个实验有200个epoch,每个epoch有400个batch,一个batch包括一对图像。学习率为5e−45e^{-4}5e−4,没60个epoch衰减为原来的0.5。使用mk-LNCC作为图像相似性指标,(w,w)={(0,3,S/4),(0.7,S/2)}(w,w)=\{(0,3, S/4),(0.7, S/2)\}(w,w)={(0,3,S/4),(0.7,S/2)},其中SSS表示最小的图像维度,滑动窗口的步长设为S/4S/4S/4,采用因子为2的空洞卷积。
四、实验
1. 实验设置
数据集选用的是OAI(Osteoarthritis Initiative)MR数据集,预处理时将灰度值归一化,然后去除小于0.1%和大于99.9%的灰度值。所有图像重采样到192×192×80192\times192\times80192×192×80大小。将无标签的数据按7:3划分为训练集和验证集。使用Dice值来评价配准的结果。选用SyN,Demons,NiftyReg和VoxelMorph作为baseline。
2. 基于优化的多尺度仿射配准
开始时,对低分辨率的图像预测一个粗略的仿射参数,然后使用它作为更高分辨率的初始化参数。使用学习率为1e−41e^{-4}1e−4的随机梯度下降优化器,三个图像尺寸分别为原图大小的k∈{0.25,0.5,1.0}k\in\{0.25, 0.5, 1.0\}k∈{0.25,0.5,1.0},每个迭代200,200,50次。使用mk-LNCC作为相似性准则。在尺寸为1.0时,参数设为(w,s)={(0.3,Sk/4),(0.7,Sk/2)}(w,s)=\{(0.3, S_k/4),(0.7, S_k/2)\}(w,s)={(0.3,Sk/4),(0.7,Sk/2)},在尺寸为0.5和0.25时,参数为(w,s)={(1.0,Sk/2)}(w,s)=\{(1.0, S_k/2)\}(w,s)={(1.0,Sk/2)}。
3. 基于优化的多尺度vSVF配准
采用仿射映射作为初始映射,采用和多尺度仿射配准相同的多尺度策略。采用LBGFS优化器,尺度为{0.25,0.5,1.0}\{0.25, 0.5, 1.0\}{0.25,0.5,1.0},每个尺度迭代60次,使用mk-LNCC作为相似性指标。
4. 实验结果
上图是AVSM的结果图,前五行分别为源图像、目标图像、配准后的图像、具有变形网格的配准后图像和多步仿射配准的图像。后三行为源标签、目标标签、AVSM配准后的标签。
上图是不同配准方法的效果对比表,Affine-opt和vSVF-opt表示基于优化的多尺度仿射和vSVF配准。实验结果显示本文的AVSM模型在跨物体配准中的表现超过了其他方法,在纵向配准中,AVSM的效果也足够好,但是略差于基于优化的方法。这可能是因为形变是微妙的,或者源图像和目标图像非常相似,所以数值优化方法可能更好的对齐。
上图是更直观的表示,越高的越好。
当没有初始的仿射对齐时,VoxelMorph的效果不佳,只有进行了仿射对齐时,VoxelMorph才会达到很好的效果。
为了衡量变化映射的平滑性,使用雅克比行列式来估计:Jϕ(x):=∣Dϕ−1(x)∣J_{\phi}(x):=\mid D \phi^{-1}(x)|Jϕ(x):=∣Dϕ−1(x)∣,当{x:Jϕ(x)<0}\{x:J_\phi(x)<0\}{x:Jϕ(x)<0}时,表示重叠。
上图是消融实验的结果表,在消融实验部分,发现引入多步和逆一致性提升了仿射配准的性能,相比于使用NCC,本文使用的mk-LNCC效果更好。
使用ln(1∣V∣∥Φ−1∘(Φts)−1−id∥22\ln(\frac{1}{|V|}\|\Phi^{-1}\circ(\Phi^{ts})^{-1}-id\|^2_2ln(∣V∣1∥Φ−1∘(Φts)−1−id∥22来衡量对称性,其中VVV表示图像体积大小,Φ\PhiΦ表示通过仿射变换和可变形变换合成的映射。上图显示了有对称损失和无对称损失时的配准结果。
由于不同的方法对边界的处理不同,所以本文只对离边界10体素远的图像进行衡量。
具有2步训练的多步vSVF配准结果,从左图可以看出效果随着step的增加而提升,右图显示同时重叠也增加。
【论文笔记】AVSM:结合了仿射配准和vSVF配准的医学图像配准模型相关推荐
- 【论文笔记】用循环一致性避免形变场重叠的医学图像配准网络
本文是论文<Cycle-Consistent Training for Reducing Negative Jacobian Determinant in Deep Registration N ...
- 【论文笔记09】Differentially Private Hypothesis Transfer Learning 差分隐私迁移学习模型, ECMLPKDD 2018
目录导引 系列传送 Differentially Private Hypothesis Transfer Learning 1 Abstract 2 Bg & Rw 3 Setting &am ...
- 【论文笔记:场景】边缘计算中分布式协作的动态卸载模型:森林火灾管理的一个用例
[摘要]:随着物联网(IoT)的发展,数据量不断增长并变得更加多样化.将数据传输到云端时存在几个问题,例如网络带宽和延迟的限制.这引起了人们对边缘计算研究的极大兴趣,边缘计算处理和分析数据导致的网络终 ...
- 【论文笔记】当Bert炼丹不是玄学而是哲学:Mengzi模型
论文标题:Mengzi: Towards Lightweight yet Ingenious Pre-trained Models for Chinese 论文链接:https://arxiv.org ...
- 学习笔记:医学图像配准简介—附voxelmorph模型
VoxelMorph官方代码:GitHub - voxelmorph/voxelmorph: Unsupervised Learning for Image Registrationb 官方代码讲解: ...
- 论文笔记:CycleMorph: Cycle Consistent UnsupervisedDeformable Image Registration
CycleMorph: Cycle Consistent Unsupervised Deformable Image Registration 针对本篇论文个人总结: 文章最重要的点在于施加循环一致性 ...
- 医学图像配准综述学习
医学图像配准综述学习 目前针对医学图像配准的综述较少,笔者只找到了两篇: 一篇发表在<Machine Vision and Applications>-- Deep Learning in ...
- ICCV 2019 | 微软开源无监督学习的医学图像配准方法:递归级联网络
点击我爱计算机视觉标星,更快获取CVML新技术 编者按:目前,深度学习正广泛应用于医学图像配准领域.无监督机器学习方法能够广泛利用临床中产生的大量原始.无标注医学图像,然而现有算法对于变形大.变化复杂 ...
- 【论文笔记】CIRNet:基于CycleGAN的无监督循环配准模型
本文是论文<UNSUPERVISED THREE-DIMENSIONAL IMAGE REGISTRATION USING A CYCLE CONVOLUTIONAL NEURAL NETWOR ...
最新文章
- C# 数组 二维数组
- 领域驱动设计系列文章汇总
- m1mac安装linux,M1 Mac 能安装 Ubuntu 和 Linux 了 ??
- 年度重磅!中国《营销自动化应用基准报告 2021》正式发布!
- kafka架构:分区机制详解
- java jks 转pfx_JKS和PFX文件相互转换方法
- redis(11)--AOF持久化
- 【CF1199 D,E, F】Welfare State // Matching vs Independent Set // Rectangle Painting 1
- 结构设计模式 - 代理设计模式
- Scrum立会报告+燃尽图(十月二十三日总第十四次)
- 2021年中国物流科技行业投融资现状及未来发展前景分析[图]
- Socket协议脚本编写
- 我的USB DAC改装(MOD)流程(上)
- 2022年,渲染农场都怎么收费,比较便宜的渲染农场测评
- html常用长度度量单位,度量长度的单位是什么 毫米以下的五个计量单位分别是什么?...
- 椭圆检测 Arc Adjacency Matrix-Based Fast Ellipse Detection 代码和数据集使用教程
- sql中的left join 查询
- NOIP2009 pj
- MATLAB小波变换图像处理简单示例
- 任你和QQ陌生人聊天
热门文章
- git使用及上传代码到github
- 计算机状态oxcoooooe9,解决win7 system文件损坏修复的办法
- linux安装activeMQ 5.16 以及集群环境
- anaconda安装后打不开,卡在loading application
- 75道经典AI面试题,我就想把你们安排的明明白白的!(含答案)
- 基于TCP协议的游戏代理接口测试工具<一>:设计初衷与工具构想
- Python 正则表达式+字符串分割(数字/字母/汉字/特殊字符)
- Java中Number转为百分比
- Deep Clustering via Joint Convolutional Autoencoder Embedding and RelativeEntropy Minimization论文笔记
- 透视HTTP协议(一) —— HTTP是什么