JoVE Logo
发表
联系我们

登录

需要订阅 JoVE 才能查看此. 登录或开始免费试用。

本文内容

  • 摘要
  • 摘要
  • 引言
  • 研究方案
  • 结果
  • 讨论
  • 披露声明
  • 致谢
  • 材料
  • 参考文献
  • 转载和许可

摘要

我们提出了一个工作流, 以分割和量化的小梁骨2D 和3D 图像基于骨骼的外部边界使用一个 ImageJ 插件。这种方法比当前手工轮廓方法更有效和准确, 并提供了一层层的 quantifications, 在当前的商业软件中是不可用的。

摘要

微计算机断层扫描 (显微 CT) 通常用于评估不同骨质流失情况下小动物骨量和微梁的显微结构特征。然而, 显微 CT 图像小梁分析的标准方法是切片半自动手工轮廓, 这是劳动密集型和易出错的。本文描述的是根据骨的外边界自动分割小梁骨的有效方法, 在适当的情况下, 在合适的情况下, 可自动识别和分割短梁骨。设置了分割参数。为了分析满意的分割参数, 给出了分割结果的图像栈, 将分割参数的所有可能组合按顺序变化, 并通过相关参数的分割结果可以易于目视检查。作为插件的质量控制特性, 模拟标准对象被量化, 在测量量可以与理论值进行比较。该插件报告了小梁性质和小梁厚度的逐层量化, 并能很容易地对所选区域内的这种性质分布进行分析。虽然逐层量化保留了关于小梁骨的更多信息, 并有助于进一步统计分析结构变化, 但此类措施无法从当前商业软件的输出中获得, 只有单一每个示例都报告每个参数的量化值。因此, 所描述的工作流是比较准确、高效地分析小梁骨的较好方法。

引言

小梁骨显微 CT 分析是跟踪不同骨丢失情况下动物骨骼形态学变化的标准方法1,2,3, 其中几个变量与报告了骨骼的结构4。但是, 此类参数在长骨骼的骺中不均匀分布5, 并且仅报告每个样本的每个结构变量的汇总值或平均值 (由当前商用微 CT 机6,7, 但单个值不能完全表示分析区域中测量参数的特征。一层一层定量的小梁骨不仅保留了对每个变量的更多信息, 而且还能够在分析区域中对此类变量的分布进行分析, 便于随后对结构进行统计分析。在不同条件下的更改5。因此, 该方法的目的是量化微梁骨的显微 ct 扫描在每个切片水平, 这是没有在任何商业可用的微 ct 分析包目前。

为了有效地分割小梁骨切片, 需要自动分割方法。然而, 目前的显微 CT 分析标准技术是基于人工交互式轮廓, 然后是半自动插补, 将小梁骨从皮质隔间分离出来, 这是劳动密集型、易出错和与大量运算符偏差关联8,9,10。自动分割方法11,12报告, 但这种方法只有在良好的分离小梁骨和皮质骨的区域, 但不是在没有明确的分离区域。另外, 不同的样本12需要不同的分割参数, 通过尝试各种参数组合12, 手动选择适用于骨骼样本组的满意的分割参数是很繁琐的,即使在设置所有相关参数时, 分割过程也是自动的。由于骨外边界与扫描背景有最大的对比, 长骨骼的骺皮质壳在所选择的分析区域中表现出很少的变化, 根据长骨骼外边界轮廓分割方法可以可靠准确地将小梁骨从皮质壳中分离出来。这种分割方法的优点是, 分割是基于背景和骨骼的外边界之间的差异, 而不是小梁和皮质骨骼之间的差异6,12, 13, 因此通常很容易找到一个组合的分割参数, 这是令人满意的一组骨样本, 有助于更可靠地分析不同组之间的小梁变化。

在每个切片级别上, 报告了2D 分析的面积、周长和二维 (2D) 厚度, 而在 3D quantifications 中报告了体积、表面和三维 (3D) 厚度。目前的图像分析工具通常不报告此类信息, 表明所报告的程序可应用于需要此类信息的一般图像。

研究方案

涉及动物主题的程序是按照《动物保育和使用指南》 (NIH 出版, 第八版, 2011) 进行的, 并已由武汉机构动物保育和使用委员会审查和批准。大学.

1. 软件安装

  1. 安装 ImageJ 软件。从 https://imagej.nih.gov/ij/下载与64位 Java 捆绑在一起的 ImageJ (版本 1.51p) 软件的 Windows 版本。将下载的软件解压到一个文件夹中, 随后将其称为 "ImageJ 目录"。
    注: 小梁分析插件需要64位 java 运行时间 (版本 1.8) 和64位 windows 操作系统, 最好是64位 windows 7 操作系统。
  2. 安装小梁分析插件。请求小梁分析插件从 http://www.bomomics.com, 并解压的内容到 ImageJ 的插件目录, 这是 "ImageJ 目录/插件"。
    注: 该插件可以作为一个自由版本获得, 其中5指定相邻切片的测量报告, 或作为商业版本, 其中的切片范围可以指定和测量。

2. 为小梁分析准备3D 数据集

  1. 在标准扫描协议5之后, 用微型 CT 扫描大鼠股骨, 然后以可导入 ImageJ 的格式 (例如) tiff 格式保存数据。如果在单个样本管中同时扫描多个骨样本, 首先将数据导入 ImageJ, 然后使用 ImageJ 的图像处理工具裁剪其他样本, 从而分离出每个骨骼。随后, 将生成的图像保存为可在以后导入 ImageJ 的格式。
    注意: 分析中使用的代表性示例图像文件 (辅助文件 1)。
  2. 模拟标准2D 和3D 对象。
    1. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击圆圈按钮, 然后在弹出窗口中输入200作为直径, 然后单击OK以生成直径为200像素的模拟圆 (图 3B)。以 tiff 格式保存生成的圆。
    2. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击正方形按钮, 然后在弹出窗口中输入200作为侧面长度, 然后单击OK以生成具有200像素的边长 (图 3B) 的模拟正方形。以 tiff 格式保存生成的正方形。
    3. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击矩形按钮, 然后在弹出窗口中输入200作为宽度和100作为高度, 然后单击确定生成具有200像素宽度和100像素高度的模拟矩形 (图3 B)。以 tiff 格式保存生成的矩形。
      注: 一个圆 (直径: 200 像素), 正方形 (边长: 200 像素) 和一个矩形 (宽度: 200 像素; 高度: 100 像素) 为后续分析保存。
    4. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击多维数据集按钮, 然后在弹出窗口中输入30作为侧面长度, 然后单击确定生成模拟多维数据集。最后, 单击插件 | 3D |卷查看器查看生成的多维数据集并以 tiff 格式保存它 (图 3C)。
    5. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击长方体按钮, 然后输入80作为长度, 40 作为宽度, 30 作为高度在弹出窗口中, 然后单击确定生成模拟长方体。最后, 单击插件 | 3D |卷查看器查看生成的长方体并以 tiff 格式保存它 (图 3C)。
    6. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击球面按钮, 然后在弹出窗口中输入30作为直径, 然后单击确定生成模拟球体。最后, 单击插件 | 3D |卷查看器查看生成的球体并以 tiff 格式保存它 (图 3C)。
    7. 打开 ImageJ 软件。在插件 |BoMomics |模拟对象菜单, 单击钢瓶按钮, 然后在弹出窗口中输入30作为直径和100作为高度, 然后单击确定生成模拟气缸。最后, 单击插件 | 3D |卷查看器查看生成的气缸并以 tiff 格式保存它 (图 3C)。
      注: 球体 (直径:30 像素), 立方体 (边长:30 像素), 长方体 (长度:80 像素; 宽度:40 像素; 高度:30 像素) 和圆筒 (直径:30 像素; 高度: 100 像素) 为后续分析保存。

3. 性能分析参数

  1. 打开 ImageJ 软件, 打开或导入扫描的图像。
  2. 滑动底部滚动条以选择一个切片, 然后单击图像 |调整 |阈值按钮。在 "弹出式阈值" 窗口中, 通过手动检查调整最小和最大阈值, 以确保骨骼与背景分离良好, 并将最小阈值记录为皮质骨阈值。
    注: 由于不同的微 CT 制造商使用不同的比例因子来存储图像中的 X 射线衰减系数, 所以实际的阈值应该通过手动检查或按照特定制造商的建议.在我们的实践中, 设置最小阈值的 ISQ 文件生产的微型 CT 机到 6,000-7, 000 可以可靠地分离骨骼的扫描背景。
  3. 单击插件 |BoMomics |Trab 参数分析按钮。在弹出窗口中, 将切片索引设置为代表切片的位置, 并设置皮质骨 ("Cort 骨")、范围 ("范围") 和步骤 ("步长") 值, 用于计算用于分析分割参数的一组皮质门限, 其中皮质骨阈值由步骤 3.2 (图 1) 获取。默认范围和步骤值为2000和400为大多数微 CT 图像工作。将切片索引设置为示例数据集的 5, 并保留其他设置的默认值。
    注: "Cort 骨" 是皮质骨的门槛, 用于参数分析的阈值将从最低阈值更改为最高阈值, 并增加步长值, 其中最低阈值是 Cort 的值。骨范围和最高阈值是 Cort 骨 + 范围。如果 Cort 骨骼、范围和步骤的 set 值分别为6000、1000和 500, 则皮质骨骼的最低阈值为6,000-1、000 = 5000, 最高阈值为 6000 + 1000 = 7000, 分析中使用的阈限为 5000,5500、6000、6500和7000。
  4. 设置噪音直径 ("噪音直径"), 步骤 ("步骤"), 和范围 ("范围") 值, 用于指定分析中的一组噪声值, 以及用于计算一组孔值的孔直径 ("孔直径")、步骤 ("步长") 和范围 ("范围") 值。通常, 默认设置适用于大多数微 CT 骨样, 其中噪声、步骤和范围分别为5、5、2、孔径、步长和范围分别为15、5、2。保留示例数据集的默认设置。
    注: "噪音直径" 是噪音抑制过滤器的直径, 而 "孔直径" 是皮质骨内孔的直径。用于参数分析的噪声和孔值可以类似于上面描述的皮质阈值, 使用指定的噪声/孔、范围和步长值。由于在采集图像时扫描参数的变化, 噪声和孔参数应根据图像质量进行实证确定, 因为没有一般的噪声范围和孔参数对所有扫描都有好处。阈值、噪音或孔的最低数值必须大于或等于 0, 并且如果使用提供的参数计算出的最低值为负值, 则将特定的最低值设置为0。
  5. 单击确定以执行参数分析。在 "参数分析结果" 窗口中直观地检查分割结果, 然后选择切片图层, 以便准确地勾勒出骨骼外部边界 (图 1B)。随后, 从参数分析结果表中与所选切片层对应的项中检索分析参数 (表 1)。

4. 小梁分析

  1. 小梁骨的分割
    1. 打开 ImageJ 软件, 然后打开或导入扫描的图像。
    2. 单击插件 |BoMomics |Trab 细分按钮, 并填写适当的分析参数。设置 "开始", "轮廓边界", "Trab"。骨骼 "," 降噪直径 "," 孔填充直径 ", 和" 皮质厚度直径 "分别为5、7200、7000、6、12和25。
      注: "开始" 和 "结束" 指定选定切片范围的小梁骨分割, "轮廓边界" 是对应的异形 "Cort 骨" 参数, "降噪直径" 到 "噪音直径"。参数, 并 "孔填充直径" 到 "孔直径" 参数。"皮质厚度直径" 是排除外皮质骨的指定厚度。trab.骨骼 "是提取小梁骨的门槛 (图 2), 其中设置是使用分析分析参数命令确定的, 如步骤3.5 中所述。
    3. 单击确定执行小梁分割。在Trab 细分结果窗口 (图 2B) 中直观地检查分割结果。以 tiff 格式 (图 2B) 中的分段小梁骨窗口保存所提取的小梁骨, 可由其他软件进一步分析。
  2. 小梁骨分析。
    1. 打开 ImageJ 软件, 然后打开或导入扫描的图像。
    2. 单击插件 |BoMomics |Trab 分析按钮, 并填写适当的分析参数, 如 "开始"、"轮廓边界"、"Trab"。骨骼 "," 降噪直径 "," 孔填充直径 ", 和" 皮质厚度直径 "如上所述 (图 3A), 其中设置是使用" 分析分析参数 "命令确定的, 如步骤3.5 中所述。设置 "开始", "轮廓边界", "Trab"。骨骼 "," 降噪直径 "," 孔填充直径 ", 和" 皮质厚度直径 "分别为5、7200、7000、6、12和25。
      注意: 在插件的自由版本中, 五个相邻切片从指定的 "开始" 切片索引开始选择用于度量, 而在商业版本中, 用户可以指定任意数量的切片。
    3. 在 "结果报告" 部分中, 为要测量的参数选择一个或多个选项, 其中所选区域的小梁骨量 (BV)、总容积 (TV) 和厚度均测量为二维 (2D) 或三维(3D) 可供选择, 通过三个复选框, 即 "BV 电视", "2D", 和 "3D"。选择复选框 "2D" 和 "3D", 然后单击 "确定" 执行小梁分析 (图 3A,表 2)。
      注: 当检查 "bv 电视" 时, 无论 "2D" 和 "3D" 的选择状态如何, 都报告了 bv、电视和强度的原始测量方法, 并将分割后的小梁骨提取并显示在一个新窗口中, 可以保存并进一步分析其他软件.在检查 "2D" 时, 报告了在每个切片级别上用平板模型测量的 BV、电视、强度和厚度的原始测量值。如果检查 "3D", 则在没有任何模型假设的情况下直接计算每个体素的三维厚度, 然后报告每个切片级别取样的 BV、电视、强度和三维厚度的原始度量值。但是, 如果未选择复选框, 则使用上述参数集对小梁骨进行分段, 但不报告任何测量。

5. 量化模拟物体

  1. 打开 ImageJ 软件, 然后打开一个模拟图像。在这里, 打开模拟球体的直径为30像素的例子。
  2. 选择插件 |BoMomics |Trab 分析按钮并填写相应的分析参数, 如前所述。保留默认值为 "开始"、"结束"、"轮廓边界"、"Trab"。骨骼 ", 并设置" 降噪直径 "," 降噪直径 "," 孔填充直径 ", 和" 皮质厚度直径 "到 0 (图 3)。
    注意:对于模拟对象, 没有意外的信号噪声, 也可能没有相应的皮层壳。因此, 应适当地设置这些值的参数 (默认设置为零)。对于 3D quantifications, 在指定的 "开始" 切片之前和之后使用插件的自由版本处理30切片, 而只有由 "开始" 和 "结束" 切片指定的区域中的切片才能通过插件的商业版本进行分析。
  3. 在结果报告部分中, 为要测量的参数选择 "2D" 和 "3D", 然后单击确定对模拟对象执行小梁分析 (表 3)。

6. 小梁测量和数据表示的校准: 剖面分析区域中小梁的分布

  1. 根据微 CT 供应商的说明, 从扫描数据集中获取校准信息。
    注意: 只有原始测量的 BV, 电视, 强度和厚度是由插件报告。为了获得与其他软件生成的报告类似的度量值, 需要对结果进行校准。
  2. 打开 "Microsoft Excel", 然后打开 "报告的结果" 表。根据以下表达式计算新 excel 列中校准的骨量 (bv)、总容积 (TV)、骨矿物质含量 (BMC)、骨量分数 (bv/TV) 和骨密度 (BMD)。
  3. 注: 扫描分辨率 (分辨率、µm)、灰度缩放 (缩放)、密度单位 (毫克 HA/cm3)、密度斜率 (斜率) 和密度截距 (拦截) 可以从扫描的微 CT 图像或元文件 (如 Scanco ISQ 文件) 中提取。因此, 校准的措施计算如下:
    假设 bv、电视和强度都是原始测量, BVc和电视c是校准值、分辨率、缩放、斜率、截距均来自于微 CT 图像的校准文件。
    bvc = BV x 分辨率    3 [微米3]
    电视c = 电视 x 分辨率3 [微米3]
    bv/电视 = bv ÷电视
    BMC     = (强度÷缩放 x 斜率-BV x 拦截) x 分辨率3 x 10-12 [毫克 HA]
    骨密度-BMC ÷电视c x 1012 [毫克 HA/cm3]
  4. 使用 Microsoft Excel 软件 (图 4) 创建针对切片层 (X) 的校准度量值 (Y) 的 XY (散点) 图。

结果

小梁分析插件旨在自动分割和量化小梁骨的准确性。最初, 骨外边界被发现和划定后, 一个孔充填手术, 其中任何孔内的骨外皮质壳填补。然后进行侵蚀手术, 排除外皮质骨, 并获得分割小梁骨。最后, 对分段区小梁骨的测量方法进行量化。

由于微 CT 图像固有的噪音, 分割使用预定义的任意参数往往无法准确识别骨骼的外部边界?...

讨论

本研究介绍了一种用于分析小梁骨的 ImageJ 插件, 它具有自动、高效和人性化的功效。该插件还可用于量化任何2D 或3D 对象, 用于逐层度量区域、卷和厚度。目前, 通过标准的显微 CT 分析, 每个样品的每一个小梁参数只有一个测量值, 不能完全代表所选分析区域中测量实体的特征。所述插件报告每个样本的每一个参数的逐层数量, 在选定的分析区域充分保留测量参数的分布信息, 因此更先进和敏感的?...

披露声明

提出了一种用于量化物体厚度的底层算法的专利。作者与其他人合作, 注册了 bomomics.com 网站, 托管免费版本插件, 在这里提供咨询和量化图像分析服务的要求。

致谢

这项工作得到了赠款 NFSC 81170806 的部分支持。作者感谢武汉大学口腔医学院的微 CT 核心设施, 帮助对大鼠股骨进行扫描和分析。

材料

NameCompanyCatalog NumberComments
ImageJNIHimagejAny version with a java 1.8 run time
trabecular analysis pluginBomomicsbomomicsfree or commercial version
Micro CT scannerScancoμ-50micro CT from any vendor
Computer SystemLenovoany brand
Windows Operating SystemMicrosoftWindows 7 x64any 64-bit Windows operating system 
Office SoftwareMicrosoftOffice 2010any speadsheet software that has xy chart function

参考文献

  1. Ruegsegger, P., Koller, B., Muller, R. A microtomographic system for the nondestructive evaluation of bone architecture. Calcif Tissue Int. 58 (1), 24-29 (1996).
  2. Muller, R., Ruegsegger, P. Micro-tomographic imaging for the nondestructive evaluation of trabecular bone architecture. Stud Health Technol Inform. 40, 61-79 (1997).
  3. Clark, D. P., Badea, C. T. Micro-CT of rodents: state-of-the-art and future perspectives. Phys Med. 30 (6), 619-634 (2014).
  4. Bart, Z., Wallace, J. Microcomputed Tomography Applications in Bone and Mineral Research. Advances in Computed Tomography. 2, 121-127 (2013).
  5. Ji, Y., Ke, Y., Gao, S. Intermittent activation of notch signaling promotes bone formation. Am J Transl Res. 9 (6), 2933-2944 (2017).
  6. Jiang, Y., Zhao, J., White, D. L., Genant, H. K. Micro CT and Micro MR imaging of 3D architecture of animal skeleton. J Musculoskelet Neuronal Interact. 1 (1), 45-51 (2000).
  7. Laib, A., et al. 3D micro-computed tomography of trabecular and cortical bone architecture with application to a rat model of immobilisation osteoporosis. Med Biol Eng Comput. 38 (3), 326-332 (2000).
  8. Cole, H. A., Ichikawa, J., Colvin, D. C., O'Rear, L., Schoenecker, J. G. Quantifying intra-osseous growth of osteosarcoma in a murine model with radiographic analysis. J Orthop Res. 29 (12), 1957-1962 (2011).
  9. Jensen, M. M., Jorgensen, J. T., Binderup, T., Kjaer, A. Tumor volume in subcutaneous mouse xenografts measured by microCT is more accurate and reproducible than determined by 18F-FDG-microPET or external caliper. BMC Med Imaging. 8, 16 (2008).
  10. Soviero, V. M., Leal, S. C., Silva, R. C., Azevedo, R. B. Validity of MicroCT for in vitro detection of proximal carious lesions in primary molars. J Dent. 40 (1), 35-40 (2012).
  11. Kohler, T., Stauber, M., Donahue, L. R., Muller, R. Automated compartmental analysis for high-throughput skeletal phenotyping in femora of genetic mouse models. Bone. 41 (4), 659-667 (2007).
  12. Buie, H. R., Campbell, G. M., Klinck, R. J., MacNeil, J. A., Boyd, S. K. Automatic segmentation of cortical and trabecular compartments based on a dual threshold technique for in vivo micro-CT bone analysis. Bone. 41 (4), 505-515 (2007).
  13. Dougherty, G. Quantitative CT in the measurement of bone quantity and bone quality for assessing osteoporosis. Med Eng Phys. 18 (7), 557-568 (1996).
  14. Doube, M., et al. BoneJ: Free and extensible bone image analysis in ImageJ. Bone. 47 (6), 1076-1079 (2010).
  15. Bouxsein, M. L., et al. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J Bone Miner Res. 25 (7), 1468-1486 (2010).

转载和许可

请求许可使用此 JoVE 文章的文本或图形

请求许可

探索更多文章

133ImageJCT

This article has been published

Video Coming Soon

JoVE Logo

政策

使用条款

隐私

联系我们

向图书馆推荐

JoVE 新闻简报

科研

教育

发表

图书馆员

关于 JoVE

版权所属 © 2025 MyJoVE 公司版权所有，本公司不涉及任何医疗业务和医疗服务。