当前位置:文档之家› CT图像重建(X线成像系统)

CT图像重建(X线成像系统)



DE=D×QF
DE单位:希沃特(Sv),D:戈瑞(Gy),QF:无量纲

DE常用于辐射防护的定量计算。
CT剂量
当量剂量和有效剂量

综合以人类为对象的研究结果,专家们对QF进行 了适当的调整,确定了辐射权重因子Wr (radiation weighting factor)的值,并提出了 当量剂量ED( equivalent dose)的概念。
CT图像重建

p即为CT扫描过程中采集到的投影数据

P是输入射线与输出射线强度比值的对数,在数值上等于沿射线 方向上物质的衰减系数的线积分。
CT重建问题:已知物质的X线衰减系数的线积分,如何求解它的 线性衰减系数分布?


数据预处理
投影数据测量误差:

假设X线束是单能的
散射线问题:假定所有到达探测器的都是初级X线光子

CT图像重建
直接反投影法
CT图像重建

反投影法的基本思想
在对某一层面一个方向的扫描完成后,用得到的投影值沿着 扫描路径回抹到体素对应的像素上。改变方向后的多次扫描 形成多次回抹,同一像素上多次回抹的灰度累加即完成图像 重建。

反投影法示例
CT图像重建

第一次运算:水平照射后,将射线和放入图像单元


CT图像重建

直接反投影法的局限:
容易产生星形伪影
产生原因:反投影法把取 自有限物体空间的投影均 匀地回抹(反投影)到了 射线所及的无限空间的各 个像素上,包括原来像素 值为0的点。
CT图像重建
中心切片定理
CT图像重建

中心切片定理:某断层(或它对应的图像)f(x,y)在视角为 时 得到的平行投影(函数)的一维傅里叶变换,等于f(x,y)二维傅 里叶变换F(w1,w2 )过原点的一个垂直切片,且切片与轴w1相交成
CT图像重建
概述 直接矩阵变换法 直接反投影法 中心切片定理 滤波反投影法
CT图像重建
概述
CT图像重建


投影重建图像理论的提出
1917年,奥地利数学家Radon提出了投影重建图像的理论 Radon的观点解决了从函数的线积分求解原函数的问题,即由 物体的一组横断面的投影来重建其横断面图像。


CT图像重建


工程可实现的滤波反投影法的重建步骤:
在某个投射方向上(设角度为 )对物体进行平行束扫描, p(t , ) 获得离散的投影数据 对投影数据 p(t , ) 填补足够多的0以避免“周期间”的干
扰,得到 p ' (t , ) 对补零后的投影 p' (t , ) 进行傅里叶变换,得到 P' (, ) 给频域的投影 P' (, ) 乘以一个滤波函数 H ' ( ) ,得到滤 波后投影的傅里叶变换 G ( , )

CT图像重建

反投影法图解

假设位于扫描范围内只有 一个钉子,则经过钉子进 行一个方向扫描的投影是 一个脉冲函数。 将测得的脉冲信号反投影 到矩阵中去,就得到第一 次反投影。

CT二次线性扫描,将测到的 脉冲信号再反投影到矩阵中去, 得到第二次反投影。 系统每旋转一次,便构成一次 新的反投影。 所有反投影的叠加形成了一个 带有星形伪影的图像,即一个 扩散的圆形区域。

ED=D×Wr×Wt

ED(有效剂量)考虑了人体具体组织器官承受辐射的能力和 危险系数,能更准确地反映不同类型的辐射对人体造成的危 害。
CT剂量
CT剂量分布曲线
CT剂量


CT剂量和普通放射剂量的区别
普通放射过程中剂量接受面积大,剂量一般集中在皮肤表面,而 CT的射线源在不停地旋转,剂量分布比较均匀。 CT采用窄束X线,普通放射检查采用宽束射线。在同样照射条件 下,宽束线散射线多。 CT检查的射线能量高,线质硬,穿透性强,被人体吸收少。 CT的探测器转换效率高,射线利用率高。 CT辐射能量不仅照射了所选择的层面,在层面附近也被能量照 射。 CT剂量不能用常规X线机的病人入射表面剂量(ESD)来表示。



CT剂量
轴向剂量分布曲线
在理想状态下,轴向剂 量分布曲线和灵敏度剖 面线SSP(Section sensitivity profile)重 叠,均呈高斯形状,其 全值半高宽(FWHM) 等于层厚。


CT剂量

在实际中,SSP比剂量分布曲线略窄,尤 其是在探测器侧使用准直器时。




CT剂量
X线的剂量基础 CT剂量分布曲线
CT剂量表示方法
剂量的影响因素
剂量的测量
CT剂量
X线的剂量基础
CT剂量


照射量
用来描述传递给患者的辐射总量。 单位为:伦琴(R) 1R=2.58×10-4C/kg
定义:在1kg空气中产生2.58×10-4C(库仑)静电荷所需要 的X线或 射线的能量。

若用P5=u2+u4…⑤ 式代替④式联立方程组,
∵ ⑤可由①+ ② -③得到 ∴由⑤、①、 ② 、③组成的方程组只有三个独立方 程,方程数少于未知数,方程组无唯一解。

假如把物质的扫描面分成N×N矩阵,只要投影数 据即方程数量足够多,同样可解得每一像素的X线 衰减系数。 在求解方程组时有多种方法,其中之一是直接矩阵 变换法。
CT图像重建

直接矩阵变换法
CT图像重建
假定某物质在扫描面上由4个均匀的部分组成,且衰 减系数分别为u1、 u2、 u3、 u4,并已知它们在水平、 竖直和对角方向的积分。
选择其中四个方程组成独立方 程组: P1=u1+u2 ①
P2=u3+u4
P3=u1+u3 P4=u1+u4

③ ④
CT图像重建

第二次运算:垂直照射后,再将射线和加到图像单 元中
CT图像重建

第三次运算:取对角线方向(右上)的照射,将射 线和再加到图像单元中

第四次运算:取左上方向的对角线照射,并将射线 和加到前面的图像单元中
CT图像重建

最后的运算:从每个图像单元中减去背景值(背景 强度等于某投射角情况下各投影值之和),再将各 吸收系数除以最大公约数,得到最后结果。 其他资料的最后运算:在求出累加值后,再给累加 得到的各个像素除以反投影的次数,也就是除以经 过像素的射线数。


CT剂量
CT剂量
CT剂量
CT剂量

英国1989年CT检查仅占X线诊断的2%,而导致 的国民集体剂量约占总剂量的20%。 1998年,英国CT检查占X线诊断的4%,检查所 致国民集体剂量已上升至40%。 常规X线正位胸片的有效剂量约0.03mSv,而做 CT胸部扫描时,有效剂量最高可达6~7mSv, 是X线正位胸片的200倍。 CT辐射剂量问题越来越受到关注。


CT剂量指数并不直接表征各种扫描所致受检者的剂量, 但与受检者所受剂量密切相关。 与吸收剂量的单位相同:戈瑞(Gy)或mGy

CT剂量


CT剂量指数100(CTDI100)
X线CT旋转一周,将平行于旋转轴Z轴(垂直于横断面)的剂 量分布D(z)沿Z轴从-50mm到+50mm积分,除以层厚T 的商 CTDI100=
ED=D×Wr

当辐射有多个种类和能量时,在一个组织或器官的 当量剂量就是各个辐射所致的当量剂量的和。
CT剂量

ED给出了不同辐射条件下人体发生生物效应的定量描述,但 人体不同组织或器官对辐射的敏感性和给人体造成的危害是 不同的。 因此又提出组织权重因子Wt(tissue weighting factor) 对上述的ED进行进一步修正。由Wr 和Wt两个因子修正后的 吸收剂量称为有效剂量(effective dose,ED)


CT图像重建

对 G ( , ) 进行傅里叶反变换,得到被滤波的时 域投影 g (t , ) 将时域投影 g (t , )反投影并累加到二维离散图像 矩阵。为了提高空间分辨率,经常在反投影前进 行预插值 改变一个投射方向对物体进行扫描,对所得投影 重复上述过程,直到所有投射角度扫描完毕
CT图像重建

探测器和数据采集系统的非线性:探测器的输出 有赖于其之前受到的辐射情况 与扫描物体有关:患者在扫描过程中并非完全静 止不动 焦点外的X线辐射、扫描物质中金属的存在、机 架未对准、球管焦点漂移、机械稳定性、球管转 子颤动等。 数据的预处理和后处理与重建算法同样重要。



CT图像重建
散射辐射使剂量曲线的尾部拉得很长。
轴向剂量分布曲线在很大程度上由焦点尺 寸、CT的几何参数、准直器宽度决定。
CT剂量
CT剂量表示方法
CT剂量


CT剂量指数(CTDI)
CTDI(CT dose index)是CT设备辐射剂量特征的实用 表征量。 迄今得到公认使用的CTDI: CT剂量指数100(CTDI100) 加权CT剂量指数(CTDIw) 容积CT剂量指数(CTDIvol)

CT图像重建

1967年CT研发时所采用的图像重建方法即联立方程组法。


联立方程组法的局限性:
当方程组的规模越来越大时,即便在计算机上编程实现, 其工作量也很大。 需采集远远多于N2个投影数据,因为许多方程是相关的。 当方程的数量超过未知数数量时,方程组的解未必收敛, 因为投影值的测量存在误差。


CT图像重建需解决的问题
就是希望用检测到的投影数据,设法求出该断层的每个像素的线 性衰减系数或其等效值(如密度)后,再以其灰度值绘出图像。 我们求解的是:一个二维分布函数,该函数表示物质的X线衰减 系数。
相关主题