CT成像介绍第六次课件.pptVIP

透射式CT （Transmission CT） ——衰减衬度CT成像（Attenuation contrast CT, 即通常CT，重点讲） —— 相位衬度CT成像（Phase contrast CT, 简介） 发射式CT（Emission CT, 简介） 散射式CT（Scattering CT，简介） 7.ＣＴ的扫描模式 扇束扫描 左图：肺部的X光影像（不能分辨深度） 右图：肺部的CT断层影像（能分辨深度） 传统X光影像与CT断层影像对比 由断层重构的三维影像 虚拟解剖（科研教学） 虚拟内窥镜（诊断） 虚拟手术（手术设计） 由探测器采集射线衰减量 =投影数据(生数据) 扩展水平成像视野的RT扫描模式及直接重建方法 扩展垂直成像视野的螺旋扫描模式及直接重建方法 针对感兴趣区域的扫描模式及直接重建算法 两种工业CT超分辨成像方法 几种CT重建算法的GPU加速实现方法，使图像重建速度达到国际先进水平 CT扫描参数的自动获取方法 X射线能量谱的间接测量方法 针对射线束流不稳定、探测器之间的非一致性、对能量响应的非线性性，提出了数据校正方法 …… 已有的算法需要将数据通过插值重排成平行束或扇形束，因此计算量大、分辨率低 我们算法的特点: 直接重建,不需要重排 = 速度快，分辨率高 具有并行计算的特点 虚拟探测器的超分辨算法 国家自然科学基金委员会重点项目汇报材料 0 * * DBP图像 CT 图像 部分DBP图像 面整列探测器 我们方法重建结果 重排方法重建结果 垂直成像视野拓展—螺旋扫描模式及其断层直接重建的快速方法 3.２拓展垂直方向视野的研究 Before 2D interpolation Phantom After 2D interpolation 3.3 感兴趣区域扫描模式及直接重建算法 转台平移的超分辨扫描模式 探测器相叠的示意图 平移旋转超分辨扫描模式 等价扫描模式 探测器相叠示意图 特点：转台转动精度要求不高，近似算法 特点：要求转台转动精度高，精确算法 适用范围：探测器单元尺寸与焦点尺寸比较大的CT系统 3.5 CT超分辨成像方法 虚拟焦点的超分辨算法 实际焦点和虚拟焦点的强度分布示意图 适用范围：探测器单元尺寸与焦点尺寸比较小的CT系统 线对卡 The image reconstructed by FBP The image reconstructed by the virtual-focus-based method The image reconstructed by ART 3.6扫描参数自动提取方法 （1）实际的扫描机械难以对准； （2）所需参数也无法直接测量。 经典图像重建算法要求的对准关系 不满足经典对准关系的一般情况 几何伪影的形态分析 (a) (b) (c) 如美国Varian公司型号PaxScan2520的面阵探测器的每帧扫描图像为1920*1536，A/D转换为12位，数据输出为16位。按720个角度采样，数据量约为4GB。相应的最大重建图像为1920*1920*1536，数据量约为10GB。即使重建1024*1024*800的三维图像，数据量仍有1.56GB 3.8CT算法的加速问题 CT常用的硬件加速方法： 多CPU的并行机 价格昂贵 SSE（Streaming SIMD Extension）加速幅度小 图形处理器（GPU）发展速度快 高带宽 高并行性 价格便宜 图形处理器（GPU）加速的缺点：灵活性差——它是为图像显示设计的，需要将算法设计成图形学的绘制形式 目前重建速度仍是制约三维锥束CT应用的主要因素之一 传统的方法：主要集中在FDK算法中反投影部分的加速，其中FFT和滤波等较复杂的计算则放在CPU上实现。纽约州立大学的Muller研究小组的工作具有有代表性 . 2007年论文“Real-time 3D computed tomographic reconstruction using commodity graphics hardware，Phys. Med. Biol. 52, 3405–3419 (2007)”中反投影时间是8.9秒的反投影速度 GPU加速锥束FDK算法研究 我们算法特点：在GPU上实现了FDK重建算法(包括滤波、反投影和累加全过程)，并支持重建1024立方的CT图像。用我们的方法由360个512平方的投影数据重建512立方反投影需要7.752秒. 所有重建时间16.816秒. 优于目前的商用锥束CT机的重建速度 滤波的GPU实现：