使用顺序 3D-3D 注册进行四维 CT 分析

Summary

我们从四维计算机断层扫描数据中分析了关节运动学。顺序 3D-3D 注册方法半自动地提供运动骨骼相对于主体骨骼的运动学，从四维计算机断层扫描数据。

Abstract

四维计算机断层扫描 （4DCT） 提供一系列体积数据，并可视化关节运动。但是，对 4DCT 数据的数值分析仍然很困难，因为在所有体积帧中分割都非常耗时。我们旨在使用顺序 3D-3D 注册技术分析关节运动学，使用 4DCT DICOM 数据和现有软件提供运动骨相对于固定骨骼的运动学。源骨骼的表面数据从 3DCT 重建。修剪的表面数据与 4DCT 中第一帧的表面数据分别匹配。这些修剪曲面按顺序匹配，直到最后一帧。这些过程为 4DCT 的所有帧中的目标骨骼提供位置信息。确定目标骨骼的坐标系后，可以计算任意两个骨骼之间的平移和旋转角度。此 4DCT 分析在复杂结构（如腕骨或牙质骨骼）的运动学分析方面提供了优势。但是，由于运动伪影，无法跟踪快速或大规模运动。

Introduction

联合运动学已经使用许多方法来描述，例如运动捕获传感器、2D-3D注册和尸体研究。每种方法都有特定的优缺点。例如，运动捕捉传感器可以使用红外摄像机测量快速、大规模运动，无论是否对主体^1、2进行传感器。然而，这些方法测量皮肤运动，以推断关节运动学，因此包含皮肤运动^误差3。

卡达韦里克研究已被用来评估运动范围，不稳定性，和接触区域4，5，6。这种方法可以使用CT或光学传感器直接连接到骨头使用针脚或螺钉测量小关节的微小变化。Cadaveric模型可以主要评估被动运动，尽管已使用多个执行器将外力应用于肌腱以模拟动态^运动7。主动关节运动可以通过 2D-3D 配准技术进行测量，将 3DCT 图像与 2D 荧光镜图像相匹配。虽然注册过程的准确性仍然存在争议，但报告的准确性一般足够高，适用于大型关节运动^学8、9。但是，此方法不能应用于狭窄空间中的小骨骼或多个骨骼。

Protocol

这里描述的所有方法都已获得庆应大学医学院机构审查委员会的批准。

注:通过重建移动骨骼围绕固定骨骼的运动来测量关节运动学。对于膝关节运动学，股骨被定义为固定骨骼，而骨质被定义为移动骨骼。

1. CT成像协议

1. 设置 CT 机器。使用 320 探测器行 CT 系统进行 CT 检查，以允许具有 160 mm 颅骨覆盖的多个阶段的 3D 体积数据。例如，在膝关节运动学分析中，图像采集包括 51 次体积扫描，旋转时间为 0.275 s，所有图像均使用半重构重建，因此时间分辨率约为 0.16 s。
2. 使用以下扫描参数:峰值管电压 = 100 kVp;管电流 = 40 mA;扫描覆盖范围 = 160 mm;矩阵大小= 512 x 512 像素;和重建截面厚度和截面间隔 = 0.5 mm。
3. 将学科的目标关节置于 CT 龙门内的 4DCT 考试的起始位置（图1）。
4. 在 CT 考试之前，在规定的考试时间内排练关节从起始位置到结束位置的动作。要求受试者在 10.275 s 扫描时间内移动接头，并获取一系列体积数据。以 DICOM 格式存储顺序卷数据。
5. 对所有目标骨骼执行静态 3DCT，并将数据存储在 DICOM 格式中。

2.....

Results

我们描述了在膝盖延长期间，tibia的运动。膝关节位于CT龙门。三角形枕头用于支撑股骨在起始位置。膝盖在10s的过程中被延长到一个直的位置。除4DCT外，还进行了整个股骨、头骨和骨骨的静态3DCT。重建了整个股骨和头骨的表面数据。骨皮层HU数的阈值设置为250胡，并重建所有51帧的表面数据。

股骨和股骨被修剪成部分表面数据，包括所有4DCT帧通过目视检查4DCT电影数据，这是在预设的4DCT软件中创建。在静态 3DCT 曲面和 4DCT 的第一帧中，绘制了每个段的地标。在股骨上，确定了内侧和侧体和间质切口。在头骨上，关节表面的中端和侧端和骨质管状也被确定为相应的地标。根据这三个地标，股骨和头骨的部分表面数据与4DCT数据的第一帧大致匹配。然后，使用 ICP 算法完全匹配这些曲面。

第一帧股骨和头骨的部分部分与第二帧的整个表面相匹配。.......

Discussion

我们的方法允许可视化和量化整个骨骼的运动，并从 4DCT 数据提供移动骨骼相对于固定骨骼的数值位置数据。已建议使用许多工具测量关节运动学。运动皮肤标记可以分析身体的整体运动很长一段时间。但是，此方法包含皮肤运动错误³。应从相邻骨骼的运动量估算关节运动学。2D-3D 配准方法使用荧光镜和从顺序 2D 图像推断 3D 运动学。尽管分析软件已经发展到可以解释这一点，但转换错误仍然存在。许多尸体研究已经测量关节运动学通过采取CT图像在不同的尸体^位置19。但是，这些代表来自连续静态 3D 图像的被动运动，因此与主动运动在质量上有所不同。

此协议中有几个关键步骤。应精确创建来自 3DCT 的表面数据，因为此质量会影响 4DCT 第一帧的初始曲面配准的准确性。在关节区域周围，骨皮层的阈值可能与骨轴不同。当骨皮层的边界不清楚时，需要进行阈值调整。完成所有帧的表面配准后，应检查重建运动。如果一帧的表面注册失败，则可以通过选取下一帧中的地标并重复协议，从下一帧重新启动自动曲面注册。

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
4DCT 扫描仪	佳能医疗系统（日本栃木）	N/A	4DCT 扫描，静态 3DCT 扫描
AVIZO（9.3.0）*	Thermo Fisher Scientific（美国俄勒冈州）		图像处理软件 根据 CT DICOM 数据和点云数据重建表面。 * Ryan， TM &Walker， A. 类人灵长类动物肱骨和股骨头的小梁骨结构。Anat Rec（霍博肯）。293 （4）， 719-729， doi：10.1002/ar.21139， （2010）.
Meshlab**	ISTI（意大利比萨）	N/A	表面修剪和地标拾取 ** MeshLab：开源网格处理工具。第六届欧洲图形意大利分会会议，第 129-136 页，2008 年 P. Cignoni， M. Callieri， M. Corsini， M. Dellepiane， F. Ganovelli， G. Ranzuglia
VTK（6.3.0）***	Kitware （美国纽约）	N/A	迭代最近点算法。用于 python 语言编程. *** https://vtk.org
Python（3.6.1）	Python Software Foundation	N/A	DICOM 文件处理，以从骨皮层（"dicom.py"模块）中提取点云。 旋转矩阵的计算。（Numpy 模块） 使用 ICP 算法进行顺序图像配准