Clear Sky Science · zh

BoneDensityMapping:用于处理和可视化骨密度数据的 R 包

· 返回目录

为何骨强度与我们每个人都息息相关

随着人们寿命延长,即便是轻微跌倒导致的骨折也越来越常见,常常带来疼痛、手术和漫长的康复期。医生知道,骨骼变薄变弱是这一问题的重要原因,但现有测试通常只为整个髋部或脊柱区域提供单一评分。本文介绍了 BoneDensityMapping —— 一个将医学扫描转换为生动三维地图的免费软件工具包,它能准确显示个体骨骼在何处强或弱,以及这些模式在人群间如何不同。

从常规扫描到详细骨密度地图

现代计算机断层扫描(CT)在许多常见医疗检查中已经捕获到丰富的三维骨骼图像。定量计算机断层扫描利用这些影像来估算骨骼中每个微小部位的矿物质含量,这是强度的重要指标。与将信息平均为单一数值的标准骨密度检测不同,这种方法可以揭示骨表面及内部密度的高峰与低谷。挑战在于,此类扫描产生庞大而复杂的数据集,若无专门工具难以处理和解读。

Figure 1. CT 扫描数据如何被转换为彩色三维地图,以显示骨骼何处强或弱。
Figure 1. CT 扫描数据如何被转换为彩色三维地图,以显示骨骼何处强或弱。

为清晰且可重复工作而构建的工具箱

BoneDensityMapping 以 R 统计语言编写,任何人都可以免费下载使用。它提供一组构建模块函数,引导用户遵循一致的工作流程:导入扫描文件和三维骨骼模型,检查二者是否正确对齐,估算骨表面与内部的密度,并将这些数值转换为带颜色编码的图像。该包兼容用于医学影像和三维网格的常见文件格式,并提供简单的检查以确保解剖标志和模型位于扫描体积的正确位置。

同时可见骨的外壳与核心

该软件区分骨骼的坚硬外壳与海绵状内部。对于外壳,它可以直接读取每个表面点在扫描中对应位置的密度,或沿着短线深入骨内以找到该层中最高的密度值。对于内部,它会在骨模型内填充规则间距的点网格,并查找每个点的密度。随后这些数值被转换为颜色并绘制到三维模型或截面上,为单个个体或多人群体提供直观的骨强度分布图。

比较不同年龄和群体的骨骼

为展示该包的功能,作者分析了来自年轻和年长成人的六块腕骨。研究者先选取一块骨作为模板,在其表面布置数千个均匀分布的点。借助解剖标志作为锚点,在其他骨上找到对应点,以便逐一比较相同位置。软件随后在这些对应点估算密度、对整组值进行平均,并为年轻组和年长组生成并排的密度图。最后,它在单一骨模型上创建差异图,突出显示年长骨骼通常低于年轻骨骼密度的区域。

Figure 2. 样本点和 CT 切片如何揭示年轻和年长腕骨之间的详细骨密度差异。
Figure 2. 样本点和 CT 切片如何揭示年轻和年长腕骨之间的详细骨密度差异。

这对未来护理与研究意味着什么

作者认为,使这类分析变得简单、透明且可复现,有助于更多门诊和研究团队更好地利用现有扫描数据。由于工具是开放且基于脚本的,其他科学家可以检查、重用并改编这些方法,甚至将其扩展到测量组织密度的其他成像类型,如磁共振成像。尽管用户仍需使用单独工具轮廓化扫描中的骨骼并具备一定的 R 使用经验,但计划中的功能扩展(例如自动标志定位和动画)旨在降低这些门槛。总体而言,BoneDensityMapping 提供了一种实用方法,将大规模影像数据集转化为清晰的骨强度图景,支持更好的衰老、疾病与治疗研究。

引用: Telfer, S., Lacambra, L. BoneDensityMapping: an R package for processing and visualizing bone density data. Sci Rep 16, 15324 (2026). https://doi.org/10.1038/s41598-026-46582-y

关键词: 骨密度, 骨质疏松, CT 成像, R 软件, 舟骨