CTA血管减影详细教程附练习数据

CTA血管减影方法

概述

  • 任务:通过减影的方法分割脑血管。
  • 成像方式:注射增强药物前扫描的CTABrainBaseline图像和注射增强药物之后扫描的CTABrainContast图像。
  • 扩展模块:需要安装SlicerElastix和SlicerVMTK。
  • 3DSlicer版本:Slicer-4.10或更高版本。

数据下载



说明:CTABrainBaseline为未注射增强药物的CT,CTABrainContast为注射增强药物后的CTA。


方法

  • 使用 General registration (Elastix) 模块对CTABrainBaseline和CTABrainContast进行配准,最大程度地减少由于患者运动引起的伪影。
  • 使用 Subtract scalar volumes 模块从CTABrainContast图像中减去CTABrainBaseline图像。
  • 使用 Vesselness filtering 模块移除减影不全的骨板。
  • Segment editor 模块中,使用 Threshold 效果提取血管。
  • 使用 Islands 效果 Keep largest island 方法去除多余的伪像。
  • (可选)使用 Centerline extraction 模块计算容器中心线。

要点:

  • 如果CTABrainBaseline图像和CTABrainContast图像之间只有刚性的患者运动,则 generic rigid (all) 可以使用配准预设,它会尝试通过应用刚性配准来对齐图像。它可能比使用变形配准更快捷。
  • 通过按Ctrl + F查找模块并开始输入其名称
  • 调整图像的窗宽/窗位以提高二维视窗中血管的可见性:在工具栏中单击“窗宽/窗位”鼠标模式(4.11版本才有),然后单击Ctrl +鼠标左键并拖动(从容器点开始)以设置最佳窗宽/窗位水平。

示例

下图动画显示CTABrainBaseline和CTABrainContast图像之间的对比,明显可见患者增强扫描前后头部有运动。

e600014f2122806e6b68bed5b1d54ada148f012b_2_1035x430

  • 转到 General registration (Elastix) 模块
  • 对于 Fixed volume 选择CTABrainBaseline
  • 用于 Moving volume 选择CTABrainContast
  • Preset 选择 generic (all)
  • 对于 Output volume ,选择 Create new Volume as... 并输入 CTABrainContastAligned
  • 单击 Apply 并等待,直到配准完成(可能需要几分钟)

配准后的CTABrainBaseline和CTABrainContastAligned可见患者运动造成的误差减少。

image-002

  • 转到 Subtract scalar volumes 模块以从CTABrainContastAligned减去CTABrainBaseline
  • Input volume 1 选择 CTABrainContrastAligned
  • Input volume 2 选择 CTABrainBaseline
  • 对于 Output volume ,选择 Create new Volume as... 并输入 CTABrainSubtracted
  • 请点击 Apply

  • 转到 Vesselness filtering 模块
  • Input volume 选择 CTABrainSubtracted
  • 对于 Seed point 选择 Create new MarkupsFiducial 并单击图像中的血管点。提示:选择一个不容易分割的点,例如靠近骨骼表面但未被减影完全抑制的点。
  • 单击开始 Preview ,等待过滤器预览完成。

  • 缩小切片视图并检查结果:骨骼的某些部分被检测为血管。通过打开 Advanced 截面并将 Suppress plates 值增加到35%(在此示例中使用的其他值: Minimum vessel diameter = 1体素, Maximum vessel diameter = 5体素, Vessel contrast = 52, Suppress blobs = 10%)来增加对骨板的抑制。

  • 可选:通过删除当前点,然后在切片查看器中放置一个新点,在几个不同的位置检查过滤结果(单击旁边的垃圾桶图标 Seed point ,单击其旁边的箭头按钮,然后单击图像)。单击 Preview 按钮查看过滤结果。
  • 单击 Start 以对整个图像执行血管过滤。该操作可能需要几分钟。
  • 可选:在“ Volumes”模块中调整血管图像的可视化。将window( W )值设置为1.0,将level( L )值设置为0.5。设置 ThresholdOff

  • 转到 Segment editor 模块。
  • 单击 Add 按钮创建一个新的Segment。
  • 双击 Color 列中的绿色矩形,开始键入 Artery 并单击 Enter 以确认Segment的名称和颜色。
  • 选择“Threshold effect”,将较低的阈值设置为0.7-以显示相关的血管,但仍保持可控制的伪影数量。点击 Apply

  • 以3D形式显示结果:取消选中按钮的 Surface smoothing 下拉菜单, Show 3D 以使3D显示更新更快,并防止血管缩小。然后单击 Show 3D button 以在3D视图中查看血管。

  • 要删除断开的非血管结构,请单击“ Islands 效果”,使用默认 Keep largest island 选项,然后单击“确定” Apply

阈值= 0.07的最终结果:

image

阈值= 0.15的最终结果:

image

曹主任您好,请问Vesselness filtering 模块在哪个版本中可以下载?4.10.2版本扩展模块无法搜索到。

需要在扩展模块管理器中安装SlicerVMTK模块。

谢谢曹主任指导。

教程很实用,谢谢曹主任

曹主任,从这步可否出一个GIF的教程,seed point总是选择不好。

在医院pace系统取得数据分不清CTABRAINCONTAST和CTA…BASELINE

dicom数据载入3DSlicer后自动区分各个序列。

谢谢

为什么我选择“ 转到 Vesselness filtering 模块”
后显示“fialed to load”


老师 您好!为什么我自己提取的数据,无论怎么调shift,都调不出血管

如果独立显卡的电脑,请将渲染方式更改为“GPU”模式。


曹老师:您好! 在Subtract scalar volumes` 模块以从CTABrainContastAligned减去CTABrainBaseline后,4.13版本3d出现周围有一圈黑影,4.11版本3d就根据读不出数据

先进行配准操作,之后再执行subtract操作,一次检查过程中,因患者呼吸影响,增强前后图像有少许移动,配准后再减影操作会减少此情况发生。另外减影Interpolation order默认为1即可。