图像重建中常用的滤波器的设计

        在医学图像重建的过程中，如果不使用滤波器，重建出来的图像与使用滤波器重建出来的图像的质量相差甚远，效果如下图所示：

         由上图我们可以看出，使用滤波器能极大的提高重建图像的质量，其中最为常见的就是RL滤波器与SL滤波器。

        首先介绍RL滤波器，RL滤波器是使用窗函数对斜坡滤波器进行截断产生的。如下图所示：

                                          

在使用该滤波器的时候，需要将其离散化，离散化之后的滤波器的函数表达式为：

                          

将该滤波器使用在扇束重建中时，将换成即可。滤波器设计的代码如下：

RL = zeros(1,N);
for k2 = 1:N
    RL(k2) = -1/(pi*pi*((k2-N/2-1)*d)^2);
    if mod(k2-N/2-1,2)==0
        RL(k2) = 0;
    end
end
RL(N/2+1) = 1/(4*d^2);

        然后介绍SL滤波器，SL滤波器是使用Sinc函数对斜坡滤波器进行截断产生的。如下图所示：

                               

在使用该滤波器的时候，同样需要将其离散化，离散化之后的滤波器的函数表达式为：

                               

将该滤波器使用在扇束重建中时，将换成即可。滤波器设计的代码如下：

SL = zeros(1,N);
for k1 = 1:N
    SL(k1) = -2/(pi^2*delta^2*(4*(k1-N/2-1)^2-1));
end