FBP中的滤波器设计

FBP中的滤波器设计

CT图像重建可采用解析方法重建，速度快，无噪声且数据完备时该方法效果好。
FBP是典型的解析重建方法，适用于2D图像，其中滤波器的设计有多种，基本的滤波器为
$f(\omega )=|\omega |$

离散形式有2种

1. 离散形式的matlab实现代码

filter=2*[0:round(width/2-1), width/2:-1:1]'/width;

其中width表示投影数据的长度。

2. 离散形式的numpy实现为

filt = 2* np.hstack([np.arange(0,np.around(width/2)),np.arange(width/2,1,-1)])/width