磁共振衍生的心肌劳损评估使用特征跟踪

Summary

准确和实用的方法来衡量心肌组织应变等参数，具有很大的临床价值，因为它已被证明，这株是收缩功能障碍的更为敏感，早于经常使用的参数EF标记。

Abstract

目的：准确和实用的方法来衡量心肌组织应变等参数，是很大的临床价值，因为它已被证明，这株是收缩功能障碍的更为敏感，早于经常使用的参数EF标记。当前的CMR技术费时，而且难以在临床实践中实施。特征跟踪是一种技术，可以导致比比较的方法用更少的时间消耗更多automization和定量分析的医学图像的鲁棒性。

方法：在单相自动或手动输入初始化从该系统开始追踪个别图案代表解剖结构，随着时间的推移的位移。这种方法的专业是不需要任何类似CMR的图像，例如标记的方式，事先操纵图像。

结果：该方法很好地适合跟踪肌肉组织和心肌和定量阐述这使血流量。

结论：这种新方法提供了一个强大的和节省时间的程序，定期对CMR成像序列的位移，速度和变形参数量化心肌组织和血液。因此，它可以在临床实践中实施。

Protocol

1。简介

边界的自动检测图像分析中的一个根本问题。在心脏成像，成像左心室心内膜边界的自动检测的可能性将给予客观测量的心室卷，心肌变形（应变）。这是在超声心动图斑点追踪技术完成。可靠的自动边界检测方法的发展，是一项艰巨的任务，并没有收到普遍可靠的解决方案，在心脏磁共振（CMR）的。事实上，在临床实践中，边框是绘制手动操作员或软件检测到心肌和腔^1,2之间的接口。在当前文章中，我们引入了不同的方法，边界不“检测”，而他们是“跟踪”，即遵循时间，从一个可靠的现有瞬时跟踪，是通常，但不一定由有经验的操作手动绘制起点超过一个单一的框架。其次在时间上的个别点组成这样的第一个可靠的跟踪搜索相同的功能，在其附近的一个点在下面的帧。跟踪功能，可以腔组织边界或解剖的元素，沿不同的组织。他们发现，在两个地区的两帧之间的利益的最大似然方法。

本地帧到帧的位移相当于评估当地的速度（位移和时间间隔之间的比例）。在一个点的速度自动评价是从比较这样的点在两个连续帧图像数据的位移决定。已使用这种方法，在几个不同的配方，在许多研究领域。他们属于普通类被称为光流，在先进的图像分析^3,4。它们通常被称为斑点跟踪回波仪成像时使用这样的速度是遵循生理的议案^5，6，但也适用于任何如CMR其他图像模式，这些方法都提到功能跟踪或边境跟踪。

2。材料和方法

特征跟踪方法

心内膜或心外膜的边界是一个2D的CMR电影手动追溯到一个任意帧（见图1）。以及可以追溯到中秋节心肌功能。这种边界被定义为一个N点的序列，鉴定，由他们协调_对的（X _I，Y I）I = 1，... ，N。边界跟踪，通过跟踪每个单点的收益，这种跟踪是基于一个层次多尺度算法和1D跟踪技术的结合，从而保证更高的精度，和二维跟踪，这是必要的，正确地检测到的二维空间扩展功能。

为了第一次捕捉边界的几何大位移，跟踪的方向正交边界本身腔组织的边界是最好辨认的。沿着这个方向进行跟踪，使用的方法如下（见图3）透削减。绘制一条线，穿过隔离墙，通过点和正交。沿透行所采取的像素被放置在列，每列对应一个图像序列帧。沿切透的演变，可以通过这种方式为代表的一次在一个二维的代表性，其中一轴是沿直线的距离和其他轴是时间（见图2）所有的瞬间。这种表示是被称为一个M -模式在超声心动图中的CMR，它对应的“侦察兵”功能相似。为了提高分析的质量，在低信噪比图像不佳的情况下，，空间时间的表示是建行透削减为5个像素的厚度。边界跟踪，然后沿着时空的形象。

边界的二维位移，一个标准的2D跟踪（基于光流）在第二个步骤是，每个点独立，在MXM移动窗口总是在以前估计的边界点为中心。二维跟踪两个步骤进行，第一次估计有一半是受雇于第二个跟踪通道中心的移动窗口。窗口，然后从32减少到16个，在另外两个通道。

为了提高精度沿用来估计旋转和扭转的边界的议案，1D跟踪沿着从“水货”的弧形边界（见图3）厚削减内置的时空图像。 ，独立，在每个点的像素采取移动边界沿线，在移动边界点为中心，是摆在列，每列对应一帧图像序列。为了提高分析的质量，最好的捕捉功能，在边境线是5个像素延伸到组织（亚心内膜）。边界跟踪，然后沿着时空的形象与上述相同的步骤进行。为了确保在跟踪边界的空间相干性，3点的中值滤波和3点高斯滤波器（权重0.25，0.5，0.25）是适用于在邻近点在每一步计算的位移。

沿着二维时空图像追踪

本节介绍了一个过程，在二维图像从一个众所周知的立场，在一个瞬间（M型状）一个边界沿着一个方向。

X是定义的水平方向和Y垂直方向。列注明_X I，I = 1 ... M，其中M是图像中的列的数量。跟踪是一个实数_Y I = _Y（X）的离散序列的测定，从一个已知点相应的_列x k Y K开始。

从已知点_Y K _K +1 点的Y位移估计评估之间在 x _k在 x _k +1的整列整列的交叉相关。互相关函数将目前最大，最高的位置给需要两列之间的相似性最大化的垂直_位移值 ，因此Y K +1估计，加入_这样的位移为y k。重复此过程和附近的列所有对之间的结果是整个边境的_估计 Y I，I = 1 ... M.在这里交叉相关计算使用快速傅立叶变换算法，以减少计算时间。

先估计_Y I的进一步细化迭代。为了实现这一目标的图像的一个子集提取，通过采取以上和_低于以前的估计Y I和一个新的形象，其_中心对应的序列y我是生成的，用于校正跟踪几点。这一修改是反复，直到找到不改正的。

改进和更自然的结果，然后由最后的蛇过程[5]实现按照时空的形象，图像的亮度水平， _通过定点 Y K。整个过程使用的时间周期，以确保一个周期性的结果，避免漂移的影响。

特征跟踪技术的限制

边界跟踪技术，像任何斑点追踪的方法，是基于像素的亮度变化进行量化，从一帧到其他。这给出了一个下限需要看到一个斑点，是在一帧的一个像素，移动到相邻的像素在下一帧的速度有关。因此，此限制

公式1
其中Δx是像素的大小和Δt是两帧之间的时间间隔。系数K取决于跟踪算法的质量和能力评估动态子像素的变化。此限制意味着速度远高于此限制是非常准确的估计，这样的精度降低速度值的方法时，低于限额。

这种限制也意味着，在收购一方的帧速率（ΔT减少）增加允许一个大的速度和他们的快速变化（如在等容相）更容易评估。在另一边，帧速率（ΔT减少）增加增加此限制，意味着在速度较低的评价的精度降低，直到它不伴随着类似的空间分辨率的增加（减少Δx的）。

幻影形象准备

一系列人工计算机生成的循环已经准备允许在图像分析过程简单，完全控制的条件下的测试。为此，在一个理想的左心室短轴投影的幻影如下准备。

心内膜和心外膜边界表示半径为R _0（T）和 R（T），分别由两个同心圆。形象准备的纤维环，其中代表中TS两国边界的组织，统一颜色的黑色背景上的灰色。然后，一个8x8顶帽子线性滤波器应用，以避免非物理的不连续性。

心外膜运动是采取以[mm]，为 R _0（T）= 10 +5 COS（2πt/吨），其中 T是心跳期为 T = 1秒 。理论内膜运动是恒定的边界沿线取决于时间，速度只有径向和V0（T）= DR _的0 / dt =-π罪（2πt/吨），厘米/秒。百分比应变，计算相对长的边界，在零时间， 是_ST 0（T）= _{100X（R 0（T）-} R _{0（0））/ 0（0）=} 100（COS（2πt/ T ）-1）/ 3，和应变率如下（1）SR _0（T）= 10 V0 / R _0，在[S ^-1]。心外膜被假定为移动相应的厚度不变，R _1（T）= R _0（T）五毫米 ，或R _1（T）= R _0（0）5毫米 。

每幅图像大小48毫米，为本组织环广场，和一项决议NxN的。例如图像，如图4所示，板A和B的应变和应变率时间曲线显示在图4中，板C和D准备通过不同的分辨率氮，帧速率FR，和心外膜的运动型循环。
心内膜的跟踪方法应用到的图像，在第一帧上均匀地沿着圆形内膜_间隔一个数N点P。

3。代表性的成果

幻影研究

幻象计算机生成的图像的图像分析方法的应用是在这里进行了分析。计算最终错误的一个全球性的措施是由根均方的百分比差异。平均内膜应变平方米，平均和最大的错误被定义为根

公式（2）
ST _0（T）的精确值，ST（T）是由图像分析计算值，并求和的所有帧延长_N F = FRxT。使用相同的定义是为半径，速度和应变率。跟踪沿心内膜的位置大约是独立的，不同点之间的差异，远低于1％。

结果表一归纳为15个不同空间分辨率的幽灵，帧速率，并外膜边界的运动类型;用来跟踪的心内膜边界点的数目不同的效果也是所示。

在所有情况下，错误非常小的不可分割的数量（半径和应变）和差量（速度和应变率）都涉及到前衍生略大一些。这是预料中的，因为衍生运营商放大错误。结果的质量是退化的分辨率降低时，事实上，相关的像素大小代表（在一个松散的意义）从一帧到其他可读的最小排量的准确性。时间分辨率不显着影响的结果，直到足够的帧速率非常高帧速率的结果，不改善，因为帧帧位移低于像素大小。这表明，帧速率增加时，它不是由伴随着在空间分辨率的增加很少或根本没有效用。

然而，简单的正弦运动，在这里考虑，并不需要一个极端的时间分辨率。同样，使用尽可能少的8点就足够了，按照简单，圆形，心内膜形状。内膜的结果没有明显影响运动型外膜经历。我们也验证的结果是没有显着影响所采用的图像滤波。

图4中的视觉呈现的结果是计算两个瞬间的心内膜边界是在幻影图像（板A和B）的报告。 （板，C和D）的情况下，＃1和＃8小分辨率的情况下，应变和应变率的报道。应变和应变率＃1（正方形）的情况下提出与理论值，平均误差分别为0.6％和3％，等于一个很好的协议。该协议是只有一点点的情况下＃8，图像分辨率减半，有错误的应变和应变率的0.9％和4.5％，分别差。

临床验证^1。

我们比较中旬LV全片圆周谐波相位成像（竖琴）和FT技术心肌劳损191 Duchene肌肉萎缩症_{patie（εCC）}NTS按年龄和心功能不全的严重程度和42岁的匹配，控制科目。匹配的标签和SSFP片回顾展，离线分析。对于整个研究人群（N = 233），意味着金融时报_εCC（-13.3 ± 3.8％），高度相关哈弗_εCC（-13.6 ± 3.4％）与Pearson相关系数为0.899 。由哈弗（-12.52 ± 2.69％）和FT确定的DMD患者平均_{εCC（-12.16} ± 3.12％）均无显着差异（P = NS） 。同样，确定哈弗（-18.85 ± 1.86）和FT（-18.81 ± 1.83），对照组的平均_εCC无显着差异（P = NS） 。我们得出的结论是基于FT _-εCC评估与标签在一个大的DMD患者的心功能不全的广泛人口的图像_所得 εCC高度。

幻影＃	帧速率	Resol，N	心外膜的议案	N P	εř	εr 最高	εV	εV 最大	ε 简	εSR 最大	ε 圣	ε 圣 最大
1	32	401	与远藤	16	0.57	1.10	3.98	8.26	2.95	6.55	0.59	1.34
2	32	401	没有运动	16	0.27	0.60	3.62	8.53	2.90	5.99	0.75	1.29
3	64	401	与远藤	16	0.35	0.62	1.75	4.92	3.08	10.47	0.34	0.86
4	64	401	没有运动	16	0.39	0.73	1.78	5.48	2.38	6.46	0.32	0.63
5	128	401	与远藤	16	0.22	0.40	1.57	3.66	2.73	8.39	0.35	0.89
6	16	401	与远藤	16	0.41	0.92	12.10	18.27	7.12	15.17	0.88	1.93
7	32	201	与远藤	16	0.43	0.82	3.84	7.03	4.21	7.79	0.83	1.52
8	32	201	没有运动	16	0.49	0.90	4.00	8.54	4.54	10.84	0.93	1.64
9	32	101	与远藤	16	2.58	3.87	5.70	17.44	9.22	18.27	4.44	6.77
10	32	401	与远藤	64	0.32	0.73	3.98	8.09	3.11	8.86	0.48	1.10
11	32	401	没有运动	64	0.33	0.53	3.56	7.83	2.78	7.44	0.70	0.98
12	32	401	与远藤	32	0.44	1.00	3.99	8.16	2.99	6.82	0.63	1.49
13	32	401	与远藤	8	0.20	0.41	3.64	6.84	3.09	9.10	0.43	0.74
14	32	401	没有运动	8	0.22	0.41	3.36	7.42	2.76	5.93	0.58	1.01
15 *	32	401	与远藤	16	1.35	2.42	5.54	9.68	7.02	16.58	2.62	4.46

表1：魅影心内膜边界跟踪分析：根均方和最大percentual的错误[％]为在不同的幻象参数对应的主要数量计算。参数标记大胆表明幻影＃1的变化。被认为是帧速率，分辨率和跟踪点的数量的依赖。心外膜运动类型的影响，被认为是两个极限情况时心外膜的边界不动（运动）或远藤（没有增厚）的动议。最后的幻象（*）是建立无滤波，带亮度突然改变一个像素的基本逐步图像。 10％以上的错误标记为粗体。

图1。CMR图像的左心室长轴（左图）和短轴（右图）与追踪的心内膜边界，绘制上顶部。

图2。时空的代表性，在空间的图像序列，沿切透。透削减为出发点，在图3。时间演化的出发点，自动跟踪，报道。

图3。左心室，透裁员和削减弧形边界平行，长轴图像。

图4。幻影研究。两幅图像（＃2的情况下）最大扩展（板）和收缩（板B），计算心内膜边界点是重叠的。应变（板C）和应变率（板D）计算实际值的比较显示两个不同的幽灵（例＃1和＃8）。

图5（一）全球圆周（黑色曲线）和节段性应变正常患者（彩色曲线）的例子。全球圆周（黑色曲线）和节段性circumferental starain郁闷的左心室功能和左束支传导阻滞（二）患者（彩色曲线）的例子。请注意不同的计时显示左室不同步的存在峰值圆周应变

Discussion

我们描述了一种新方法“功能跟踪”来衡量SSFP图像所得的心肌组织中的变形和位移参数。这种方法最近被确认对全球高峰在一个人口众多的杜氏肌营养不良症（1）患者的周向应变的评估竖琴。我们发现这两种技术之间的一个极好的相关性。特征跟踪技术是在使用相对简单，不需要任何额外的成像序列，如标记，检测或密集。此外临床验证峰值全球ECC的评估，区域应变数据可以得出，以及金融时报“似乎是一个强大的技术能够估计的纵向和径向的LV的应变以及纵向右心室应变。

为了准确地确定有关时间和空间分辨率的图像技术的限制，需要进行更多的研究。 FT允许跨其他方式，如超声心动图的比较，因为其基本原理是相同的。

结论

特征跟踪已在临床实践中，因为它省去了耗时的分析和额外的采集程序需要集成的潜力。 FT可用于评估的左，右心室全球和节段性株，左室不同步。正在进行进一步的临床研究，以评估临床心脏病学的FT有用。

Materials

Name	Company	Catalog Number	Comments
Diogenes® Feature Tracking software	TomTec Imaging Systems, Fulda, Germany