本文介绍了一种利用三维U-Net卷积神经网络在负重锥形束CT(CBCT)图像中对承受重力载荷的人类脊柱进行分割的实验方案。该方法为研究负重条件下的脊柱力线提供了有价值的工具,并支持在负荷状态下开展脊柱的三维打印与生物力学研究。
方法文章
* These authors contributed equally
本文介绍了一种利用三维U-Net卷积神经网络在负重锥形束CT(CBCT)图像中对承受重力载荷的人类脊柱进行分割的实验方案。该方法为研究负重条件下的脊柱力线提供了有价值的工具,并支持在负荷状态下开展脊柱的三维打印与生物力学研究。
在重力负荷条件下对脊柱排列进行精确分析,对于诊断和管理脊柱侧弯及椎间盘退行性疾病等肌肉骨骼系统疾病至关重要。传统的影像学技术,如二维X线片以及仰卧位CT或MRI,无法完整捕捉脊柱在负重状态下的三维几何结构。承重状态下锥形束计算机断层扫描(CBCT)可在站立位提供高分辨率的三维图像,从而实现对脊柱排列更准确的评估。本研究提出一种基于深度学习的图像分割方案,采用具有三维卷积层和残差连接的U-Net卷积神经网络(CNN),用于分割椎体、骨盆和股骨头。该方案包括图像采集、人工标注、预处理和模型训练。对解剖结构的成功分割可实现对临床相关脊柱参数(如Cobb角和椎体旋转)的精确测量。此外,分割所得的三维模型可导出为STL格式用于3D打印,有助于手术规划和患者教育。同时,该方案还可推广应用于训练其他需在重力负荷条件下评估的解剖结构(如下肢)的分割模型。本方法为在重力负荷条件下精确分割脊柱及其他解剖结构提供了一种可靠且可扩展的工具,在临床诊断和科研研究中均具有实际应用价值。
在负重条件下对脊柱排列进行精确分析,对于诊断和管理肌肉骨骼系统疾病(包括脊柱侧弯、退行性椎间盘疾病和脊柱畸形)至关重要1,2,3,4。传统的影像学技术(如二维X线片)无法完整捕捉脊柱的三维几何结构。尽管常规CT和MRI能够生成三维图像,但这些成像方式要求患者处于仰卧位。相比之下,负重状态下锥形束计算机断层扫描(CBCT)可在患者直立位时提供高分辨率的三维图像,从而实现对脊柱排列更详细的评估5。
从三维影像中分割骨骼结构对于提取具有临床意义的脊柱参数至关重要,例如科布角(Cobb角)、椎体旋转度和骨盆入射角,这些参数是诊断脊柱畸形和规划手术所必需的6。此外,分割还可用于构建解剖结构的三维模型,有助于患者教育和手术准备7。深度学习技术,特别是卷积神经网络(CNN),已广泛应用于医学影像分析,尤其在从图像中分割复杂解剖结构方面表现突出8,9。基于深度学习的分割方法可减少人工工作量,降低操作者间变异,并提高脊柱影像分析的准确性与一致性。
本文介绍了一种基于深度学习的协议,用于利用负重锥形束CT(CBCT)对承受重力负荷的脊柱进行分割。采用三维U-Net卷积神经网络(CNN),在负重条件下自动分割椎体、骨盆和股骨头。将该方法与常规的负重CBCT结合,可提供一种标准化且可重复的工具,实现更精确的脊柱评估,从而支持更优的临床决策并改善患者预后。该方法具有可复制性和可适应性,可用于训练针对其他需要在重力负荷条件下评估的解剖结构(如下肢)的分割模型。
本研究遵循《赫尔辛基宣言》中的伦理标准,并获得了四川大学华西医院伦理审查委员会的批准(批准号:834 / 2021)。纳入研究的受试者被诊断为退行性脊柱侧凸或青少年特发性脊柱侧凸(AIS)。排除标准包括有脊柱手术史、体内存在脊柱植入物,或无法在影像检查过程中保持站立姿势的个体。所有参与本研究的受试者均在图像采集及研究参与前签署了知情同意书。研究全程严格保护所有受试者的隐私与保密性。本研究中使用的设备和软件列于材料表中。
1. 脊柱成像
2. 使用 ITK-SNAP 进行标注
3. 图像预处理
注意:在本方案中,输入图像及相应目标标签数据的图像预处理使用 Python 集成开发环境(IDE)PyCharm 实现。相关步骤通过 Python 代码执行,代码以伪代码形式提供于补充文件1中,以提高可读性并清晰地概述关键步骤。
4. 模型训练
分割结果
本研究中,对三个解剖结构——脊柱、骨盆和股骨头——进行了分割。纳入骨盆和股骨头是因为某些脊柱参数的测量需要这些结构中的解剖标志点。每个结构的数据集均按表2所述划分为训练集和测试集。
采用三种标准指标评估分割性能,包括 Dice 系数、体积相似性(VS)和召回率:

![figure-results-2 体积相似性方程:\[1-|FN−FP|/(2×TP+FP+FN)\];统计分析中的公式。](/files/ftp_upload/67781/67781eq002.jpg)

其中,TP 表示真阳性,FP 表示假阳性,FN 表示假阴性。
经过500个训练周期后,基于各项指标对分割结果进行了评估。表3展示了各个解剖结构的性能表现。这些结果表明,所提出的基于U-Net的分割模型在脊柱、骨盆和股骨头的分割任务中均具有良好的有效性。图2展示了解剖结构的代表性分割结果。
分割后应用
三维临床相关参数的测量
通过自动分割生成的三维模型可用于测量一系列临床重要参数,如Cobb角、椎体旋转度、胸椎后凸角、腰椎前凸角、骨盆入射角、骨盆倾斜角和骶骨倾斜角。此处,我们展示使用自主开发的程序对这些参数进行测量,该程序可处理分割后的数据,并对脊柱及相关结构提供准确且自动化的测量结果(图3)。
3D 打印
通过人工智能自动分割及后续处理获得的 3D 模型可导出为 STL 格式,这是一种通用的 3D 模型文件格式,被大多数 3D 打印软件广泛支持。在本研究中,所使用的 3D 打印机为工业级光固化(SLA)打印机,采用紫外光(UV)固化树脂作为打印材料。在分割模型的各个独立部分之间添加了连接结构,以确保形成一个完整连贯的整体,能够作为单一、统一的模型进行打印。该过程实现了患者解剖结构的实体化复制,可用于手术规划、医患沟通以及教学工具。图 4 展示了基于人工智能分割数据生成的 3D 打印模型。

图1:负重状态下锥形束计算机断层扫描(CBCT)成像的患者体位。 患者站立于扫描平台之上,并用安全带固定,以防止扫描过程中发生移动或跌倒。患者需握住上方支撑装置以保持稳定体位,同时CBCT系统在重力负荷条件下获取脊柱的三维图像。请点击此处查看该图的放大版本。

图 2:脊柱、骨盆和股骨头在轴状、矢状和冠状视图中的分割结果。 第一列为 CBCT 图像,随后为真实标签(ground truth),最后为深度学习模型获得的分割结果。第二列和第三列中的红色覆盖区域表示在三种视图中分割出的椎体、骨盆和股骨头区域。图像显示了该模型在捕捉解剖结构方面的高准确性,与真实标签高度一致。请点击此处查看此图的高清版本。

图3:在分割后的脊柱和骨盆结构上测量临床相关参数。 (A)冠状面上T11与L3椎体之间的Cobb角测量值(8.9°)。 (B)矢状面上胸椎后凸角(T1-T12:21.2°)的测量。 (C)矢状面上腰椎前凸角(L1-L5:26.4°)的测量。 (D)轴向面上L2椎体相对于骨盆的椎体旋转(VR:3.3°)。 (E)骨盆参数,包括骨盆倾斜角(PT:21.0°)、骶骨倾斜角(SS:26.1°)和骨盆入射角(PI:47.1°),在矢状面图像上进行测量。 请点击此处查看该图的放大版本。

图 4:脊柱与骨盆的 STL 模型及 3D 打印实体模型。左侧图像显示了准备用于 3D 打印的脊柱与骨盆的 STL 格式 3D 模型。右侧图像展示了使用立体光刻(SLA)技术打印出的相应实体模型。请点击此处查看该图的高清版本。
| 参数 | 默认设置 |
| 管电压 | 110 kV |
| 管电流 | 6 mA |
| 帧率 | 12 fps |
| 重建视野(FOV) | 350 mm |
| 层厚 | 2 mm |
| 层间距 | 1 mm |
表1:曝光与重建的默认参数。用于全脊柱负重位锥形束计算机断层扫描成像的默认曝光与重建参数。
| 标签 | 训练数据集 | 测试数据集 |
| 脊柱 | 100 | 20 |
| 骨盆 | 57 | 6 |
| 股骨头 | 47 | 5 |
表2:用于CBCT图像模型训练与测试的数据集。该表格显示了用于训练和测试模型的图像数量,以实现对脊柱、骨盆和股骨头结构的分割。
| 标签 | 骰子系数 | 体积相似性 | 召回率 |
| 脊柱 | 0.93 | 0.954 | 0.89 |
| 骨盆 | 0.925 | 0.965 | 0.93 |
| 股骨头 | 0.96 | 0.971 | 0.959 |
表3:深度学习模型对不同解剖结构(脊柱、骨盆和股骨头)的分割性能。性能指标包括 Dice 系数、体积相似性(VS)和召回率,取值范围均为 0 到 1,数值越接近 1 表示分割准确性越高。
补充文件1:图像预处理(步骤3)和模型训练(步骤4)的伪代码。 请点击此处下载该文件。
本方案介绍了一种基于深度学习的方法,用于利用负重锥形束CT(CBCT)对承受重力负荷的脊柱图像进行分割。关键步骤包括图像采集、人工标注、预处理和模型训练。正确的患者体位和曝光参数对于获取高质量图像至关重要。体位不正或对位不准可能导致图像质量不佳,从而使神经网络难以准确分割结构。预处理也是一个关键步骤,本方案中的预处理包括添加噪声、亮度变化、旋转和位移。这些预处理技术旨在使图像适配神经网络的输入要求,并通过使模型接触多种类型的图像畸变来增强其鲁棒性。
手动标注是该方案的另一个关键组成部分。手动标注的准确性直接影响模型学习和分割椎体、骨盆及股骨头等结构的能力。本方案侧重于椎体的分割,这足以满足测量大多数临床相关参数(如Cobb角和椎体旋转度)的需求。对于需要完整椎骨(包括椎板、横突和棘突)的应用,模型需使用额外的标注数据进行训练。此外,手动标注的变异性也可能带来挑战,可能影响结果的可重复性。为减轻这一问题,我们建议由多名经验丰富的标注人员共同参与,并实施观察者间可靠性检验。
具有3D卷积层和残差连接的U-Net架构非常适合处理3D图像12,13,14。多层级损失函数有助于保留深层网络中的细节信息,从而确保分割的准确性。然而,必须仔细调整训练参数(如学习率和数据增强设置),以防止过拟合并提高模型的鲁棒性。尽管近期研究已提出针对特定任务定制的脊柱分割或检测模型,采用如Transformer或YOLO等神经网络架构15,16,但这些方法通常依赖于CT或MRI脊柱图像。相比之下,本方案采用U-Net架构,该架构在像素级分割任务中更为常用17。更重要的是,该方法旨在利用直立位锥形束CT(upright CBCT)对重力负荷条件下的脊柱进行分割,从而实现对站立位或负重状态下脊柱的评估。这一独特优势可更真实地反映脊柱的排列情况和力学应力状态,对于脊柱畸形的诊断及治疗方案的制定具有重要意义。值得注意的是,该模型在脊柱分割任务中表现出优异性能,Dice系数达到0.93,与使用CT和MRI数据进行脊柱分割的其他基于深度学习的模型相比,性能相当甚至更优9,18,19。
该方法存在一些局限性。本研究使用的数据集虽然足以分割椎体、骨盆和股骨头,但若能进一步扩展,将更有利于处理罕见或复杂的脊柱畸形。更大且更多样化的数据集有助于模型推广到更广泛的病例中。另一个局限性在于CBCT成像对金属伪影较为敏感,导致模型难以准确分割带有金属植入物的术后脊柱。此外,该方法的计算需求较高,需要使用高性能GPU,而并非所有临床和研究环境都能具备相应条件。手动标注过程耗时较长,因此我们未能邀请多位标注人员通过共识方式验证标注的准确性。此外,该模型尚未在外部数据集上进行验证,这限制了其在其他临床环境中的可推广性。
准确分割和测量脊柱相关参数的能力有助于提升诊断和治疗规划水平,尤其在脊柱侧弯等病例中具有重要意义。将模型以STL格式导出用于3D打印的功能在术前规划中也具有重要应用价值,可使外科医生预先演练或参考患者特异性解剖模型。这些实体模型是患者教育的有力工具,有助于患者更好地理解自身病情及治疗方案20。除脊柱成像外,本方案还可拓展用于分割其他承重的解剖结构,如下肢及关节,这在膝骨关节炎或髋关节畸形等病例中尤为有价值。此外,本方案在生物力学研究中亦具潜在应用前景,因为在负荷状态下对脊柱或其他结构进行精细的三维建模,对于分析力学机制和应力分布至关重要21,22。这可为畸形进展过程提供有价值的科学见解,并指导新型治疗策略的开发。
综上所述,本方案提供了一种在重力负荷条件下对脊柱图像进行分割和分析的可靠且高效的方法。尽管存在一些局限性,但其在诊断、手术规划、生物力学研究和教育中的潜在应用,使其成为临床和科研场景中一种有价值的工具。
作者声明,不存在与本方案中所述工作相关的利益冲突。
本研究由四川省重点研发项目(2024YFFK0050 和 2023YFG0126)、四川大学华西医院“1·3·5”工程高峰学科发展基金人工智能项目(ZYAI24013)以及高端人才支持计划“1·3·5”项目(ZYGD23024)资助。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 3D 打印机 (G1400) | UnionTech, 上海, 中国 | G1400 | 工业级立体光刻(SLA)3D 打印机。 |
| CUDA 版本 11.4 | NVIDIA | 11.4 | 用于深度学习模型训练的 GPU 加速。 |
| ITK-SNAP | 开源 | 4.0.2 | 用于三维生物医学图像分割与标注的软件。 |
| Numpy | 开源 | 1.24.4 | 用于图像处理的 Python 库,网站:https://numpy.org/about/ |
| PyCharm | JetBrains | 2023.2.1 | 用于 Python 编程的集成开发环境(IDE)。 |
| Python 版本 3.8 | 开源 | 3.8 | 与 PyTorch 结合用于深度学习模型开发。 |
| PyTorch 版本 1.10 | 开源 | 1.10 | 用于训练分割模型的深度学习框架。 |
| SimpleITK | 开源 | 2.3.1 | 用于图像处理的 Python 库,网站:https://simpleitk.org/ |
| 紫外光固化树脂 | UnionTech, 上海, 中国 | UV-001 | 用于 SLA 3D 打印过程的材料。 |
| 负重锥形束 CT 系统 (WR3D) | Angell Technology, 深圳, 中国 | DX6290 | 用于在重力负荷条件下获取全脊柱图像 |
| 搭载 Intel i7-10700 CPU 的 Windows 10 系统 | 定制组装(多家厂商) | Windows 10 专业版, 22H2 | 模型训练的运行环境;32GB 内存,NVIDIA GeForce RTX 3090(24GB)。 |
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可