跨多个CT厂商的腕关节动态CT图像中运动量化的准确性：一项体模研究

《MEDICAL PHYSICS》：Accuracy of motion quantification in dynamic CT images of the wrist joint across multiple CT vendors: A phantom study

　　

背景：四维计算机断层扫描（4D CT）能够动态评估腕关节运动学，并为评估韧带不稳提供非侵入性替代方法。然而，时间分辨率不足以及采集或重建方案的差异可能引入运动相关伪影，影响定量分析。由于这些伪影可模拟病理性腕骨运动，必须量化测量误差以区分真实病理与方法学误差。因此，需要对不同CT系统和采集方案下的运动量化进行系统性评估。

目的：评估不同CT系统的采集和重建方案在腕骨4D CT成像中的运动量化误差。

方法：使用一个包含三个3D打印骨（舟骨、月骨、头状骨）的旋转腕部体模，在来自三家制造商的五个CT系统（Aquilion ONE PRISM和ONE Vision，Canon Medical；Revolution Apex，GE HealthCare；单源和双源SOMATOM Force，Siemens Healthcare）上进行扫描。在不同体模旋转速度下，采集一次静态3D扫描和多次4D扫描，每次持续10秒。单源覆盖0.050–0.300转每秒（rps），双源覆盖0.100–0.600 rps。在两个系统（80 kV/40 mA和120 kV/100 mA）上评估剂量依赖性。图像以全重建和部分重建模式重建，分割并使用点对图像配准进行配准。相对于静态扫描，计算舟骨和头状骨相对于月骨的平移和旋转运动量化误差，并报告每个旋转周期。运动量化误差还表示为每个重建帧归一化运动的函数，以将误差与图像采集期间发生的运动量相关联。使用线性混合模型分析CT系统、骨类型、重建方法和旋转速度的影响。统计显著性定义为p < 0.05。

结果：运动量化误差几乎不受剂量影响。对于0.100 rps，最佳单源系统的误差中位数，全重建头状骨平移为0.29 mm（四分位距：0.21–0.36 mm），旋转为1.37°（0.82°–2.18°）；部分重建平移为0.21 mm（0.16–0.25 mm），旋转为0.70°（0.52°–0.95°）；双源系统误差最低（平移0.14 mm（0.12–0.17 mm），旋转0.48°（0.39°–0.62°））。运动量化误差大致随每个重建帧的归一化运动呈线性缩放，约相当于一个重建帧期间发生运动的15%。

结论：运动量化误差主要由体模旋转速度和时间分辨率决定。所报告的值可为定义定量腕关节运动学的检测阈值提供基础，支持区分真实运动学与运动伪影引起的表观位移。

论文解读：多CT厂商腕关节动态CT图像运动量化准确性体模研究

腕关节损伤如韧带断裂需要早期准确诊断，延迟或漏诊可导致关节不稳定并发展为骨关节炎。目前临床评估依靠临床检查、关节镜和静态影像技术，如X线、计算机断层扫描（CT）和磁共振成像（MRI）。MRI关节造影虽能直接显示内在韧带，但仍是静态评估，无法判断病变是否伴随异常腕骨运动学。动态MRI虽有无辐射优势，但时间和空间分辨率需权衡。四维CT（4D CT）结合高空间和高时间分辨率，可在主动运动过程中捕获腕骨三维运动，实现定性及定量评估。然而，4D CT数据量大，手动判读费时，需自动化分析。在定义诊断阈值前，必须确定自动化分析的检测阈值，即能够与真实骨运动区分的最小运动差异。运动伪影如模糊和重影会降低分割和配准精度，不同CT系统硬件和采集重建协议也影响图像质量。部分旋转重建和双源CT可提高时间分辨率，但跨扫描仪的系统比较研究有限。因此，本研究旨在评估不同CT系统采集和重建协议下腕骨4D CT图像的运动量化误差。

研究人员使用一个旋转腕部体模，包含三个3D打印骨（舟骨、月骨、头状骨），模拟腕骨运动。在五台CT系统上采集静态3D扫描和不同旋转速度的4D扫描：两台佳能宽角CT（Aquilion ONE PRISM和ONE Vision）、一台GE Revolution Apex，以及西门子SOMATOM Force的单源和双源模式。单源系统速度0.050–0.300转/秒，双源0.100–0.600转/秒，每次持续10秒。剂量影响在两个系统（80 kV/40 mA和120 kV/100 mA）上评估。图像采用全重建和部分重建，分割基于区域生长和水平集，配准采用点对图像刚性配准。以静态3D扫描为参考，计算头状骨和舟骨相对月骨的平移和旋转误差。采用线性混合模型分析CT系统、骨类型、重建方法和旋转速度的影响。

在结果方面，研究分三个部分。3.1 静态体模中全重建和部分重建的运动量化误差：在允许部分重建的三个CT系统中，静态体模下测得的头状骨位移在全重建和部分重建之间无差异，变异性略有增大但保持在0.15 mm和0.5°以内。这表明两种重建均可用于动态实验。3.2 运动量化误差与辐射剂量：对WACT-P和DSCT在不同剂量下的扫描显示，剂量变化对运动量化精度无显著影响（所有p > 0.76），因此后续分析采用最低剂量协议。3.3 运动量化误差与时间分辨率：对于单源系统，旋转速度对误差影响最大（平移部分η²=0.69，旋转0.79），骨类型次之（0.56/0.58），重建方法较小（0.06/0.17）。头状骨误差显著高于舟骨（平移0.46 mm vs 0.14 mm，旋转1.93° vs 0.82°，p<0.001）。全重建误差显著高于部分重建（平移0.34 mm vs 0.27 mm，旋转1.60° vs 1.14°，p<0.001）。DSCT数值上最低，但与最佳单源部分重建相比差异未达显著（平移p=0.12，旋转p=0.07）。归一化分析显示，运动量化误差随每帧运动量近似线性增加，平移误差E_trans=0.15x?0.03（R²=0.60），旋转误差E_rotation=0.16x+0.11（R²=0.49）。

在讨论中，研究人员指出，旋转速度即每帧运动量是运动量化误差的主要决定因素，采集协议影响有限。部分重建虽提高时间分辨率，但可能引入额外的重建伪影，因此临床实施前应在静态体模上评估各扫描仪。通过归一化运动参数，可将误差与运动量直接关联，该误差约为每帧运动的15%。骨的位置影响误差，距旋转轴较远的头状骨误差更大。剂量不影响精度，因此可采用低剂量协议。DSCT因时间分辨率最高（66 ms）而精度最高，但需心电图触发，采集时间限制在2秒以内，且要求较高旋转速度，实际临床优势受限。研究局限性包括使用体模而非人体标本、运动模式简单、厂商数量有限等。研究结论为本研究表明，在所研究的条件范围内，腕部体模4D CT成像中的运动量化精度主要由每个重建帧内发生的运动量决定。运动量化误差随每帧归一化运动近似线性增加，与CT系统或重建协议无关。总体而言，单源4D CT协议可为腕部运动定量分析提供足够精度。选择腕部旋转速度时，需权衡减少运动伪影与保留生理相关运动。这些发现可通过定义运动依赖性检测阈值