FBP中的滤波器设计
CT图像重建可采用解析方法重建,速度快,无噪声且数据完备时该方法效果好。
FBP是典型的解析重建方法,适用于2D图像,其中滤波器的设计有多种,基本的滤波器为
f
(
ω
)
=
∣
ω
∣
f(\omega)=|\omega|
f(ω)=∣ω∣
离散形式有2种
- 离散形式的matlab实现代码
filter=2*[0:round(width/2-1), width/2:-1:1]'/width;
其中width表示投影数据的长度。
- 离散形式的numpy实现为
filt = 2* np.hstack([np.arange(0,np.around(width/2)),np.arange(width/2,1,-1)])/width