氡转换成图像对角线后的伪像？

Question

我正在使用MATLAB中图像处理工具箱中的内置氡功能。直到今天，我一直在使用一些自定义功能，这些功能给了我预期的结果。特别是，我正在开发一个数学模型，可以在几个方向上检索点扩散函数（PSF）的投影（基线是0/45/90/135度）。

我准备了一个非常简单的例子，它将显示我正在尝试的问题：

I = zeros(1000,1000);
I(250:750, 250:750) = 1;
theta = [0 45 90 135];
[R,xp] = radon(I,theta);
figure;plot(R);legend('0°','45°','90°','135°');

如果您运行该示例，您将看到45/135°（对角线）的图表显示沿曲线形成锯齿的人工制品。起初我以为它与我使用的采样网格（甚至点数）有关。但是，当使用具有奇数个点的网格时，问题仍然存在。我不太明白这个结果，因为氡变换只是几个方向上的累积积分。因此，我不应该得到这种“锯齿”模式。

我对结果感到很困惑。有人试过同样的问题吗？

提前致谢。

Answer 1

当您使用简单的前向投影仪时，它是混叠工件，我相信这是android.txt函数中实现的。要删除此工件，您需要增加采样数（randon()可能使用相同数量的模体采样，您可能希望将该数量增加到至少为幻像采样数量的两倍），或者实现一种更好的前向投影仪，例如GE的CT图像重建软件中使用的驱动驱动投影仪。