
传统CT重建的FBP算法代码,可直接运行
5星
- 浏览量: 0
- 大小:None
- 文件类型:None
简介:
这是一段用于实现传统计算机断层扫描(CT)图像重建功能的FBP算法代码。用户可以下载并直接在相关环境中运行,适用于科研和教学用途。
MATLAB提供了大量函数来方便地完成FBP算法的实现。
1. FBP算法原理:中心切片定理(CST)指出,原数据投影的一维傅立叶变换等于该原始数据的二维傅立叶变换。
经过以下步骤应该能够得到原始图像:
- 投影
- 一维傅里叶变换
- 滤波
- 二维傅里叶反变换
2. 投影相关知识:
2.1 正投影:对每一根通过像素的射线进行积分,所得值作为该角度下的权值。
对一组数据P执行Radon变换(即做正投影),会得到两个结果[R, xp] = radon(P,theta):
- xp表示的是投影线条的数量
- R则是对应于每个xp在指定的theta角下进行线积分后的数值,也就是权值。
2.2 反投影:反投影过程是利用上面正投影得到的R权值,将这些R值重新投射回xy坐标系中。
如果一个点(x, y)满足 x*cos(theta) + y*sin(theta) = R,则说明该点位于某条特定角度theta下的投影线上。
全部评论 (0)
还没有任何评论哟~


