Summary
在这里，我们介绍了一种结合多平面 3D 重建和彩色多普勒融合的 5D 超声技术，该技术可以同步可视化甲状腺结构和功能信息。通过最大限度地减少盲点，该方法可以快速、精确地定位病变，从而提高诊断准确性，特别是使新手从业者受益。
Abstract
该文提出了一种基于超声数据五维同步重建的新型甲状腺检查技术。原始时间序列被重建为反映解剖结构的 3D 体积数据。实现了三个正交平面的三平面可视化，对整个压盖进行了系统检测。彩色多普勒成像被集成到每个三平面切片中，以绘制血管变化图。这种多模态融合可以在重建的 5D 空间中同步显示结构、功能和血流信息。与传统扫描相比，该技术具有灵活的离线诊断、减少对扫描的依赖、增强的直观解释和全面的多方面评估等优点。通过最大限度地减少监督错误，它可以提高诊断的准确性，特别是对于新手从业者。所提出的5D融合方法可以快速、精确地定位病灶，以便及早发现。未来的工作将探索与生化标志物的整合，以进一步提高诊断精度。该技术在推进甲状腺检查方面具有相当大的临床价值。
Introduction
桥本氏甲状腺炎 （HT） 是最常见的自身免疫性甲状腺疾病 （AITD），是世界上碘充足地区甲状腺功能减退的主要原因1。它的特征是淋巴细胞浸润和针对甲状腺抗原的自身抗体，导致甲状腺结构的破坏和甲状腺功能减退症2。HT 的分期旨在评估严重程度并指导治疗决策。它依赖于促甲状腺激素 （TSH） 和甲状腺自身抗体3 等生化标志物的组合，以及甲状腺超声上可见的超声特征 4,5,6。
超声检查显示特征性发现，包括彩色多普勒 6,7 上的弥漫性回声、异质性回声纹理、微结节和血流量增加。然而，传统的二维 （2D） 灰度超声缺乏定量方法来系统地分析 HT 分期的这些特征8。血管变化的评估也仅限于 2D 模式下的定性目视检查。甲状腺复杂的三维 （3D） 结构进一步阻碍了使用传统 2D 切片进行全面评估 9,10。这些因素会导致影像学盲点和误解，导致敏感性和特异性低，特别是对于经验不足的从业者11,12。
传统的手持式超声扫描集成了实时采集和诊断功能。这种耦合的工作流程依赖性增加了扫描过程中出现疏忽错误的可能性。缺乏空间定位和跟踪也使病变识别和监测不精确12,13。专用的 3D 超声系统已经出现以解决这些限制，并显示出有希望的结果14,15。然而，大多数 3D 超声技术需要复杂的机械扫描机制和专门的换能器，导致成本高昂且采用障碍。
为了克服传统2D和3D超声技术的局限性，本研究提出了一种针对甲状腺检查的新型3D重建和可视化解决方案。使用广泛使用的手持式超声，首先获得多个 2D 扫描以扫描整个甲状腺。然后通过二维序列的空间配准和融合来实现三维体积重建。同时，彩色多普勒帧被共同配准，以创建可视化血流变化的血管图。重建的 3D 灰度体积和彩色血管图最终集成到一个平台中，从而实现同步的多平面可视化和组合结构功能检查。
该三维融合技术从不同方面对复杂的甲状腺形态进行了系统、全面的评估。通过最大限度地减少盲点和实现全局概览，它可以帮助提高诊断准确性并减少监督错误，特别是使新手从业者受益。多模态可视化还有助于快速准确地定位病变，为甲状腺结节和肿瘤的早期诊断和治疗带来希望。此外，该方法还引入了定量 3D 特征分析，该分析以前从未研究过 HT 分期。随着广泛采用，它有可能使目前依赖于经验的超声诊断程序标准化和客观化。通过将手持式 3D 重建、多模态融合、定量特征分析和灵活的可视化协同集成到简化的工作流程中，这种低成本、易于使用的技术代表了传统 2D 超声在推进甲状腺检查方面的诊断强大飞跃。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
本研究已获得北京中医药大学附属顺思淼医院机构评审委员会的批准。该患者是从孙思苗医院甲状腺科招募的。患者接受了甲状腺超声检查，并对该研究给予知情同意。在这项研究中，使用手持设备获取的 4D 超声数据用于重建甲状腺的三平面视图。此外，还实现了实时同步彩色多普勒成像。本研究中使用的软件工具列在 材料表中。
1. 数据收集和准备
- 使用便携式手持式超声设备，将线性阵列换能器横向放置在患者颈部，以在横截面上对甲状腺进行成像。沿着甲状腺的长度缓慢而稳定地滑动探头，同时保持探头的接触和方向。
- 以 33 Hz 的帧速率获取一系列可视化甲状腺形态的横向 B 模式图像。
- 同时，应用彩色多普勒来检测腺体和血管中的血流。从甲状腺上极到下甲状腺两极扫描，覆盖整个腺体。由此产生的动态成像序列包括形成两个 4D 数据集的连续横向切片。
- 加载和浏览 4D B 模式超声数据
- 将所有 DICOM 数据复制到自定义的工作目录。
注意：操作系统和 MATLAB 中的工作目录是相同的。键入每行后按 Enter 键以在 MATLAB 中运行该命令。 - 使用 dicomread 函数将 B 模式 US 数据文件导入 MATLAB，并使用 size 函数查看数据的维度。
- 在计算机上打开 MATLAB。
- 在 “命令 ”窗口中，键入：
VB0 = dicomread（'fname.dcm'）;
其中“fname.dcm”可以替换为 DICOM 数据的实际文件名。这将读取 DICOM 文件并将图像数据存储在变量 VB0 中。 - 若要查看加载数据的大小，请键入：
尺寸（VB0），
注意：此处导入的 4D 数据的尺寸为 768 像素 x 1024 像素 x 3 x 601 层。768 像素 x 1024 像素 x 3 对应于标准 RGB 图像，其中每个像素由三个 24 位深度的通道表示。601 层表示扫描切片的总数。
- 调用 US_B_Show 函数将4D矩阵数据转换为连续的灰度视频序列，连续回放进行详细检查（见 图1）。
- 要使用 DICOM 文件上的 dicomread 函数将步骤 1.4.2.2 中导入的 4D 超声数据矩阵 VB0 转换为连续播放的灰度视频序列，请通过在 MATLAB 命令窗口中键入以下命令来调用 US_B_Show 函数：
US_B_Show（VB0）
其中 VB0 是包含之前导入的超声数据的 4D 矩阵变量。
- 要使用 DICOM 文件上的 dicomread 函数将步骤 1.4.2.2 中导入的 4D 超声数据矩阵 VB0 转换为连续播放的灰度视频序列，请通过在 MATLAB 命令窗口中键入以下命令来调用 US_B_Show 函数：
- 图 1 中的 GUI 显示了用于暂停、前进、快退等的播放按钮。
- 按 播放 按钮开始连续播放帧序列的视频。使用 暂停和播放控制工具图标 可灵活导航任何帧。使用 放大/缩小 按钮在播放过程中动态放大或缩小图像，使用默认 缩放 按钮重置为原始 1 倍视图。
- 单击“ 检查像素值 ”按钮，然后将鼠标移到某个区域上，以将十字准线与像素坐标和强度叠加在一起，以便进行局部分析。
注意：这些交互式控件可以灵活地检查空间和时间上的超声数据特性。
- 将所有 DICOM 数据复制到自定义的工作目录。
- 加载和浏览 4D 彩色多普勒超声数据
- 使用 dicomread 函数将彩色多普勒超声数据文件导入 MATLAB，并使用 size 函数查看数据的尺寸。
注意：此处导入的 4D 数据的尺寸为 768 像素 x 1024 像素 x 3 x 331 层。768 像素 x 1024 像素 x 3 对应于标准 RGB 图像，红色和蓝色代表不同方向的血流。331 层表示扫描切片的总数。 - 使用 US_C_Show 功能将4D矩阵数据转换为连续的彩色视频序列，以便连续播放以进行详细检查（见 图2）。
注意 图 2 中的 GUI 具有与 图 1 的步骤 1.4.4 中描述的相同的交互式控件和操作集。
- 使用 dicomread 函数将彩色多普勒超声数据文件导入 MATLAB，并使用 size 函数查看数据的尺寸。
2. B模和彩色多普勒超声同步观察
注： 图 1 所示的 4D B 模式超声数据和 图 2 所示的 4D 彩色多普勒超声数据在沿时间轴的第四维中包含相同的绝对时间戳。此字段在 DICOM 元数据中记录为 FrameTimeVector。根据该字段中的时间值，可以实时同步 图 1 和 图 2 。
- 使用 dicomread 命令读取两个 4D 文件后，使用两个 4D 矩阵作为输入执行 Synchronize_B_C 函数。
注意： 图 3 显示了仍然可以连续播放的结果视频。现在的不同之处在于，4D B模式超声数据和4D彩色多普勒超声数据在同一视频帧内实时同步。 图 3 中的 GUI 具有与 前面图 1 的步骤 1.4.4 中描述的相同的交互式控件和操作集。
3. 甲状腺同步三平面重建
注意：为了能够更精确地定位和量化病变，本研究根据获取的 4D 超声数据对甲状腺进行了三平面重建，并具有实时交互性。这使临床医生能够快速准确地确定病变，为随后对受影响区域的定量奠定坚实的基础。
- 使用图 1 中的 4D B 模式超声数据作为输入调用 thyroid_triplanar 函数，以推导出三个正交平面（日冕平面、矢状平面和轴向平面），如图 4 所示。
- 图 4 中的十字准线相互作用可以实时检查甲状腺的不同部位。单击并拖动十字准线的中心，对从超声重建的甲状腺解剖结构进行任意 3D 检查。
注意： 图 4 中的 GUI 还可以调整三平面视图的灰度强度范围、对比度和亮度。 - 在图像的任何区域上按住鼠 标左键 并拖动，以实时修改 亮度 和 对比度 级别。松开鼠标按钮以确认并完成调整。
4. 三维血流场同步三平面重建
注意：基于 3D 彩色多普勒超声数据重建 4D 血流场的同步三平面视图对于表征桥本氏甲状腺炎 （HT） 在临床上也很重要。
- 使用图 2 中的 4D C 模式超声数据作为输入调用 thyroid_3D_blood 函数，以推导出三个正交平面（日冕平面、矢状面和轴向平面），如图 5 所示。
- 图 5 中的十字准线相互作用可以实时检查甲状腺的不同部位。单击并拖动十字准线的中心，对从超声重建的甲状腺解剖结构进行任意 3D 检查。
注： 图 5 中的 GUI 还可以调整三平面视图的灰度强度范围、对比度和亮度。 - 在图像的任何区域上按住鼠 标左键 并拖动，以实时修改 亮度 和 对比度 级别。松开鼠标按钮以确认并完成调整。
5. B模三平面视图和彩色多普勒三平面视图同步
注意：基于 图 4 所示的三平面视图，将相应的彩色多普勒血流图像同步到病变位置无疑将有助于桥本甲状腺炎 （HT） 病理进展的诊断和量化。
- 拖动 图 4 中的十字准线交互以定位感兴趣区域，并执行 US_B2C 以获得彩色多普勒三平面视图中的相应位置。
- 拖动 图 5 中的十字准线交互以定位感兴趣区域，并执行 US_C2B 以获得 B 模式三平面视图中的相应位置。
注： 图6 为桥本氏甲状腺炎（HT）病变的精确定位和明确诊断奠定了坚实的超声基础。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
如 图 1 和 图 2 中的图形用户界面 （GUI） 所示，可以连续检查超声扫描序列。 然而，这种二维检查在很大程度上依赖于甲状腺学家的解剖学知识来在精神上重建病变的位置，这对新手来说是一个挑战，并导致缺乏定量一致性。 图 3 将 B 模式灰度与彩色多普勒流图像融合在一起，以便对 2D 扫描进行更明智的检查。
为了便于全面可靠地评估甲状腺疾病， 图 4 显示了 B 模式强度的 3D 空间分布，而 图 5 显示了从彩色多普勒数据重建的 3D 血管图。如GUI所示（图6），实现了沿三个正交平面的结构和功能信息的同步可视化。临床医生可以连续检查甲状腺灰度图像的多平面横截面和相应的血流图像。这种互补方式的无缝整合可以在精确定位和确定甲状腺病理的严重程度方面发挥关键作用。
如果 3D 灰度体积和 3D 彩色多普勒体积本质上是跨越空间和病理维度的 4D 数据，那么将它们在两个同步的三平面可视化中联系起来可以使甲状腺科医生能够在统一的 5D 空间中快速确定病变，并根据关节强度和血流模式提供准确的诊断。
图 1：横向 B 模式超声视频环路。 通过连续横向扫描获得的连续 B 模式帧显示甲状腺形态。 请点击这里查看此图的较大版本.
图 2：横向彩色多普勒超声视频循环。 通过横向扫描获得的连续彩色多普勒帧揭示了甲状腺组织的血流特征。 请点击这里查看此图的较大版本.
图 3：同步 B 模式和多普勒超声。 集成视频同步显示甲状腺结构（灰度）和血流（彩色叠加）。彩色多普勒叠加使用色标描绘流动的方向和速度 - 红色表示流向换能器的流量;蓝色表示远离传感器的流量。 请点击这里查看此图的较大版本.
图 4：从 B 型超声中提取的三平面视图。 使用三平面可视化从 4D B 模式扫描重建的正交日冕、矢状面和轴向平面。 请点击这里查看此图的较大版本.
图 5：从多普勒超声中提取的三平面视图。 从 4D 多普勒扫描重建的正交冠状面、矢状面和轴向平面，以绘制甲状腺组织的血流特征。 请点击这里查看此图的较大版本.
图 6：融合结构和血管数据的同步三平面视图。 融合多平面重建，同步B模和多普勒数据，实现形态学和功能学的综合检查。 请点击这里查看此图的较大版本.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
协议中的关键步骤
虽然 图 1 和 图 2 对检查和诊断具有价值，但从其他角度确定病变位置和视图需要专家经验。对于桥本氏甲状腺炎（HT）的诊断，实时同步 图1 和 图2 也是一个重要而关键的步骤。协议步骤 3.3 是关键步骤之一，如 图 4 所示，主治医师可以交互式检查 3D 甲状腺解剖结构的任意横截面。这对于定位病变和识别异常组织区域至关重要。传统上，手持式超声扫描仅提供 2D 横视图。这不可避免地导致由于盲点而忽视了 3D 病理细节。同样，协议步骤 4.3 生成 3D 血流图，这对于确定病变位置也至关重要。协议步骤 5.1 和 5.2 同步结构和功能甲状腺图像，为临床医生提供更强大的数字智能工具来管理复杂疾病。
修改和故障排除
如果发生重建伪影，则采集扫描范围可能不足。扩大覆盖范围的重复扫描可以克服这个问题。还可以调整切片间距和像素大小等参数。
方法的局限性
手持式超声扫描虽然可以获取各种模式同步的时间戳，但缺乏实时3D探头定位。因此，在甲状腺模型中，只有横向尺寸被精确重建。横向平面的定量测量目前是准确的，而冠状和矢状面的测量有助于病理定位，但目前的定量尺度不可靠。
对现有方法的意义
这种 5D 超声技术通过在全景可视化空间中实现多平面结构检查和血流映射，增强了传统的 2D 扫描。它克服了 2D 超声检查中持续存在的操作员依赖性、盲点和诊断模糊性等局限性。所提出的工作流程为标准化和改变当前甲状腺疾病超声诊断的经验依赖实践奠定了坚实的基础。
潜在应用
该方法可用于精确定位和定量甲状腺结节、肿瘤和炎症性病变，如桥本氏甲状腺炎。它为放射科医生和外科医生提供了用于评估病理学的增强视觉视角。该技术在辅助诊断、治疗计划和手术指导方面具有相当大的潜力。此外，研究团队计划将生化标志物与该5D分析管道相结合，以实现人工智能赋能的甲状腺疾病的精确诊断和量化。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
甲状腺疾病精准定量软件工具，在本研究 材料表 中列为甲状腺疾病精准定量V1.0，是北京智能熵科技有限公司的产品。本软件工具的知识产权归本公司所有。作者没有利益冲突需要声明。
Acknowledgments
本刊得到了陕西省重点研发计划：2023-ZDLSF-56和陕西省“科学家+工程师”团队建设：2022KXJ-019的支持。
Materials
| Name | Company | Catalog Number | Comments |
| MATLAB | MathWorks | 2023B | Computing and visualization |
| Tools for Thyroid Disease Precision Quantification | Intelligent Entropy | Thyroid-3D V1.0 | Beijing Intelligent Entropy Science & Technology Co Ltd. Modeling for Thyroid Disease |
References
- Ragusa, F., et al. Hashimotos' thyroiditis: Epidemiology, pathogenesis, clinic and therapy. Best Pract Res Clin Endocrinol Metab. 33 (6), 101367 (2019).
- Ralli, M., et al. Hashimoto's thyroiditis: An update on pathogenic mechanisms, diagnostic protocols, therapeutic strategies, and potential malignant transformation. Autoimmun Rev. 19 (10), 102649 (2020).
- Soh, S., Aw, T. Laboratory testing in thyroid conditions - pitfalls and clinical utility. Ann Lab Med. 39 (1), 3-13 (2019).
- Cansu, A., et al. Diagnostic value of 3D power Doppler ultrasound in the characterization of thyroid nodules. Turk J Med Sci. 49, 723-729 (2019).
- Haugen, B. R., et al. 2015 American Thyroid Association Management Guidelines for adult patients with thyroid nodules and differentiated thyroid cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid. 26 (1), 1-133 (2016).
- Acharya, U. R., et al. Diagnosis of Hashimoto's thyroiditis in ultrasound using tissue characterization and pixel classification. Proc Inst Mech Eng H. 227 (7), 788-798 (2013).
- Zhang, Q., et al. Deep learning to diagnose Hashimoto's thyroiditis from sonographic images. Nat Commun. 13 (1), 3759 (2022).
- Huang, J., Zhao, J. Quantitative diagnosis progress of ultrasound imaging technology in thyroid diffuse diseases. Diagnostics. 13 (4), 700 (2023).
- Gasic, S., et al. Relationship between low vitamin D levels with Hashimoto thyroiditis. Srp Arh Celok Lek. 151 (5-6), 296-301 (2023).
- Sultan, S. R., et al. Is 3D ultrasound reliable for the evaluation of carotid disease? A systematic review and meta-analysis. Med Ultrason. 25 (2), 216-223 (2023).
- Arsenescu, T., et al. 3D ultrasound reconstructions of the carotid artery and thyroid gland using artificial-intelligence-based automatic segmentation-qualitative and quantitative evaluation of the segmentation results via comparison with CT angiography. Sensors. 23 (5), 2806 (2023).
- Krönke, M., et al. Tracked 3D ultrasound and deep neural network-based thyroid segmentation reduce interobserver variability in thyroid volumetry. PLoS One. 17 (7), e0268550 (2022).
- Hazem, M., et al. Reliability of shear wave elastography in the evaluation of diffuse thyroid diseases in children and adolescents. Eur J Radiol. 143, 109942 (2021).
- Herickhoff, C. D., et al. Low-cost volumetric ultrasound by augmentation of 2D systems: design and prototype. Ultrasound Imaging. 40 (1), 35-48 (2017).
- Seifert, P., et al. Optimization of thyroid volume determination by stitched 3D-ultrasound data sets in patients with structural thyroid disease. Biomedicines. 11 (2), 381 (2023).
