在医学图像重建的过程中,如果不使用滤波器,重建出来的图像与使用滤波器重建出来的图像的质量相差甚远,效果如下图所示:
由上图我们可以看出,使用滤波器能极大的提高重建图像的质量,其中最为常见的就是RL滤波器与SL滤波器。
首先介绍RL滤波器,RL滤波器是使用窗函数对斜坡滤波器进行截断产生的。如下图所示:
在使用该滤波器的时候,需要将其离散化,离散化之后的滤波器的函数表达式为:
将该滤波器使用在扇束重建中时,将换成
即可。滤波器设计的代码如下:
RL = zeros(1,N);
for k2 = 1:NRL(k2) = -1/(pi*pi*((k2-N/2-1)*d)^2);if mod(k2-N/2-1,2)==0RL(k2) = 0;end
end
RL(N/2+1) = 1/(4*d^2);
然后介绍SL滤波器,SL滤波器是使用Sinc函数对斜坡滤波器进行截断产生的。如下图所示:
在使用该滤波器的时候,同样需要将其离散化,离散化之后的滤波器的函数表达式为:
将该滤波器使用在扇束重建中时,将换成
即可。滤波器设计的代码如下:
SL = zeros(1,N);
for k1 = 1:NSL(k1) = -2/(pi^2*delta^2*(4*(k1-N/2-1)^2-1));
end