与彩色多普勒映射相结合的同步三平面重建，可精确快速定位甲状腺病变

摘要

在这里，我们介绍了一种结合多平面 3D 重建和彩色多普勒融合的 5D 超声技术，该技术可以同步可视化甲状腺结构和功能信息。通过最大限度地减少盲点，该方法可以快速、精确地定位病变，从而提高诊断准确性，特别是使新手从业者受益。

摘要

该文提出了一种基于超声数据五维同步重建的新型甲状腺检查技术。原始时间序列被重建为反映解剖结构的 3D 体积数据。实现了三个正交平面的三平面可视化，对整个压盖进行了系统检测。彩色多普勒成像被集成到每个三平面切片中，以绘制血管变化图。这种多模态融合可以在重建的 5D 空间中同步显示结构、功能和血流信息。与传统扫描相比，该技术具有灵活的离线诊断、减少对扫描的依赖、增强的直观解释和全面的多方面评估等优点。通过最大限度地减少监督错误，它可以提高诊断的准确性，特别是对于新手从业者。所提出的5D融合方法可以快速、精确地定位病灶，以便及早发现。未来的工作将探索与生化标志物的整合，以进一步提高诊断精度。该技术在推进甲状腺检查方面具有相当大的临床价值。

引言

桥本氏甲状腺炎 （HT） 是最常见的自身免疫性甲状腺疾病 （AITD），是世界上碘充足地区甲状腺功能减退的主要原因¹。它的特征是淋巴细胞浸润和针对甲状腺抗原的自身抗体，导致甲状腺结构的破坏和甲状腺功能减退^症2。HT 的分期旨在评估严重程度并指导治疗决策。它依赖于促甲状腺激素 （TSH） 和甲状腺自身抗体³ 等生化标志物的组合，以及甲状腺超声上可见的超声特征 ^4,5,6。

超声检查显示特征性发现，包括彩色多普勒 ^6,7 上的弥漫性回声、异质性回声纹理、微结节和血流量增加。然而，传统的二维 （2D） 灰度超声缺乏定量方法来系统地分析 HT 分期的这些特征⁸。血管变化的评估也仅限于 2D 模式下的定性目视检查。甲状腺复杂的三维 （3D） 结构进一步阻碍了使用传统 2D 切片进行全面评估 ^9,10。这些因素会导致影像学盲点和误解，导致敏感性和特异性低，特别是对于经验不足^{的从业者11,12}。

传统的手持式超声扫描集成了实时采集和诊断功能。这种耦合的工作流程依赖性增加了扫描过程中出现疏忽错误的可能性。缺乏空间定位和跟踪也使病变识别和监测不精确^12,13。专用的 3D 超声系统已经出现以解决这些限制，并显示出有希望的结果^14,15。然而，大多数 3D 超声技术需要复杂的机械扫描机制和专门的换能器，导致成本高昂且采用障碍。

为了克服传统2D和3D超声技术的局限性，本研究提出了一种针对甲状腺检查的新型3D重建和可视化解决方案。使用广泛使用的手持式超声，首先获得多个 2D 扫描以扫描整个甲状腺。然后通过二维序列的空间配准和融合来实现三维体积重建。同时，彩色多普勒帧被共同配准，以创建可视化血流变化的血管图。重建的 3D 灰度体积和彩色血管图最终集成到一个平台中，从而实现同步的多平面可视化和组合结构功能检查。

该三维融合技术从不同方面对复杂的甲状腺形态进行了系统、全面的评估。通过最大限度地减少盲点和实现全局概览，它可以帮助提高诊断准确性并减少监督错误，特别是使新手从业者受益。多模态可视化还有助于快速准确地定位病变，为甲状腺结节和肿瘤的早期诊断和治疗带来希望。此外，该方法还引入了定量 3D 特征分析，该分析以前从未研究过 HT 分期。随着广泛采用，它有可能使目前依赖于经验的超声诊断程序标准化和客观化。通过将手持式 3D 重建、多模态融合、定量特征分析和灵活的可视化协同集成到简化的工作流程中，这种低成本、易于使用的技术代表了传统 2D 超声在推进甲状腺检查方面的诊断强大飞跃。

研究方案

本研究已获得北京中医药大学附属顺思淼医院机构评审委员会的批准。该患者是从孙思苗医院甲状腺科招募的。患者接受了甲状腺超声检查，并对该研究给予知情同意。在这项研究中，使用手持设备获取的 4D 超声数据用于重建甲状腺的三平面视图。此外，还实现了实时同步彩色多普勒成像。本研究中使用的软件工具列在 材料表中。

1. 数据收集和准备

1. 使用便携式手持式超声设备，将线性阵列换能器横向放置在患者颈部，以在横截面上对甲状腺进行成像。沿着甲状腺的长度缓慢而稳定地滑动探头，同时保持探头的接触和方向。
2. 以 33 Hz 的帧速率获取一系列可视化甲状腺形态的横向 B 模式图像。
3. 同时，应用彩色多普勒来检测腺体和血管中的血流。从甲状腺上极到下甲状腺两极扫描，覆盖整个腺体。由此产生的动态成像序列包括形成两个 4D 数据集的连续横向切片。
4. 加载和浏览 4D B 模式超声数据
   1. 将所有 DICOM 数据复制到自定义的工作目录。
      注意:操作系统和 MATLAB 中的工作目录是相同的。键入每行后按 Enter 键以在 MATLAB 中运行该命令。
   2. 使用 dicomread 函数将 B 模式 US 数据文件导入 MATLAB，并使用 size 函数查看数据的维度。
      1. 在计算机上打开 MATLAB。
      2. 在 "命令 "窗口中，键入:
         VB0 = dicomread（'fname.dcm'）;
         其中"fname.dcm"可以替换为 DICOM 数据的实际文件名。这将读取 DICOM 文件并将图像数据存储在变量 VB0 中。
      3. 若要查看加载数据的大小，请键入:
         尺寸（VB0），
         注意:此处导入的 4D 数据的尺寸为 768 像素 x 1024 像素 x 3 x 601 层。768 像素 x 1024 像素 x 3 对应于标准 RGB 图像，其中每个像素由三个 24 位深度的通道表示。601 层表示扫描切片的总数。
   3. 调用 US_B_Show 函数将4D矩阵数据转换为连续的灰度视频序列，连续回放进行详细检查（见 图1）。
      1. 要使用 DICOM 文件上的 dicomread 函数将步骤 1.4.2.2 中导入的 4D 超声数据矩阵 VB0 转换为连续播放的灰度视频序列，请通过在 MATLAB 命令窗口中键入以下命令来调用 US_B_Show 函数:
         US_B_Show（VB0）
         其中 VB0 是包含之前导入的超声数据的 4D 矩阵变量。
   4. 图 1 中的 GUI 显示了用于暂停、前进、快退等的播放按钮。
      1. 按 播放 按钮开始连续播放帧序列的视频。使用 暂停和播放控制工具图标 可灵活导航任何帧。使用 放大/缩小 按钮在播放过程中动态放大或缩小图像，使用默认 缩放 按钮重置为原始 1 倍视图。
      2. 单击" 检查像素值 "按钮，然后将鼠标移到某个区域上，以将十字准线与像素坐标和强度叠加在一起，以便进行局部分析。
         注意:这些交互式控件可以灵活地检查空间和时间上的超声数据特性。
5. 加载和浏览 4D 彩色多普勒超声数据
   1. 使用 dicomread 函数将彩色多普勒超声数据文件导入 MATLAB，并使用 size 函数查看数据的尺寸。
      注意:此处导入的 4D 数据的尺寸为 768 像素 x 1024 像素 x 3 x 331 层。768 像素 x 1024 像素 x 3 对应于标准 RGB 图像，红色和蓝色代表不同方向的血流。331 层表示扫描切片的总数。
   2. 使用 US_C_Show 功能将4D矩阵数据转换为连续的彩色视频序列，以便连续播放以进行详细检查（见 图2）。
      注意 图 2 中的 GUI 具有与 图 1 的步骤 1.4.4 中描述的相同的交互式控件和操作集。

2. B模和彩色多普勒超声同步观察

注: 图 1 所示的 4D B 模式超声数据和 图 2 所示的 4D 彩色多普勒超声数据在沿时间轴的第四维中包含相同的绝对时间戳。此字段在 DICOM 元数据中记录为 FrameTimeVector。根据该字段中的时间值，可以实时同步 图 1 和 图 2 。

1. 使用 dicomread 命令读取两个 4D 文件后，使用两个 4D 矩阵作为输入执行 Synchronize_B_C 函数。
   注意: 图 3 显示了仍然可以连续播放的结果视频。现在的不同之处在于，4D B模式超声数据和4D彩色多普勒超声数据在同一视频帧内实时同步。 图 3 中的 GUI 具有与 前面图 1 的步骤 1.4.4 中描述的相同的交互式控件和操作集。

3. 甲状腺同步三平面重建

注意:为了能够更精确地定位和量化病变，本研究根据获取的 4D 超声数据对甲状腺进行了三平面重建，并具有实时交互性。这使临床医生能够快速准确地确定病变，为随后对受影响区域的定量奠定坚实的基础。

1. 使用图 1 中的 4D B 模式超声数据作为输入调用 thyroid_triplanar 函数，以推导出三个正交平面（日冕平面、矢状平面和轴向平面），如图 4 所示。
2. 图 4 中的十字准线相互作用可以实时检查甲状腺的不同部位。单击并拖动十字准线的中心，对从超声重建的甲状腺解剖结构进行任意 3D 检查。
   注意: 图 4 中的 GUI 还可以调整三平面视图的灰度强度范围、对比度和亮度。
3. 在图像的任何区域上按住鼠 标左键 并拖动，以实时修改 亮度 和 对比度 级别。松开鼠标按钮以确认并完成调整。

4. 三维血流场同步三平面重建

注意:基于 3D 彩色多普勒超声数据重建 4D 血流场的同步三平面视图对于表征桥本氏甲状腺炎 （HT） 在临床上也很重要。

1. 使用图 2 中的 4D C 模式超声数据作为输入调用 thyroid_3D_blood 函数，以推导出三个正交平面（日冕平面、矢状面和轴向平面），如图 5 所示。
2. 图 5 中的十字准线相互作用可以实时检查甲状腺的不同部位。单击并拖动十字准线的中心，对从超声重建的甲状腺解剖结构进行任意 3D 检查。
   注: 图 5 中的 GUI 还可以调整三平面视图的灰度强度范围、对比度和亮度。
3. 在图像的任何区域上按住鼠 标左键 并拖动，以实时修改 亮度 和 对比度 级别。松开鼠标按钮以确认并完成调整。

5. B模三平面视图和彩色多普勒三平面视图同步

注意:基于 图 4 所示的三平面视图，将相应的彩色多普勒血流图像同步到病变位置无疑将有助于桥本甲状腺炎 （HT） 病理进展的诊断和量化。

1. 拖动 图 4 中的十字准线交互以定位感兴趣区域，并执行 US_B2C 以获得彩色多普勒三平面视图中的相应位置。
2. 拖动 图 5 中的十字准线交互以定位感兴趣区域，并执行 US_C2B 以获得 B 模式三平面视图中的相应位置。
   注: 图6 为桥本氏甲状腺炎（HT）病变的精确定位和明确诊断奠定了坚实的超声基础。

结果

如 图 1 和 图 2 中的图形用户界面 （GUI） 所示，可以连续检查超声扫描序列。 然而，这种二维检查在很大程度上依赖于甲状腺学家的解剖学知识来在精神上重建病变的位置，这对新手来说是一个挑战，并导致缺乏定量一致性。 图 3 将 B 模式灰度与彩色多普勒流图像融合在一起，以便对 2D 扫描进行更明智的检查。

...

讨论

协议中的关键步骤
虽然 图 1 和 图 2 对检查和诊断具有价值，但从其他角度确定病变位置和视图需要专家经验。对于桥本氏甲状腺炎（HT）的诊断，实时同步 图1 和 图2 也是一个重要而关键的步骤。协议步骤 3.3 是关键步骤之一，如 图 4 所示，主治医师可以交互式检查 3D...

材料

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