在MATLAB中查找iradon中重建图像的大小

Question

请考虑以下代码：

P = Phantom(256);
theta = 0:1:179;
R = radon(P, theta);
I = iradon(R, theta);

iradon.m使用

计算重建图像的大小

N = 2*floor(size(R,1)/(2*sqrt(2)))

但为什么这个公式呢？它给出N大约等于投影数除以2的平方根。但是这将如何给出图像的大小？有没有更好的方法来找到给定R和theta的图像大小？

Answer 1

size(R,1)并未提供投影数量，而是投影尺寸。预测数量为size(R,2)。

radon的输出大小为n x m，其中n =各个投影的大小，m是投影数。投影尺寸大于图像尺寸 - 想象一下投影为45度，您需要投影大约是图像的sqrt(2)倍，以免丢失任何信息。 iradon正在进行反向计算以恢复原始图像大小。

实际上，可能由于MATLAB实现radon的方式，重建图像的大小将略大于原始图像。

Answer 2

我认为它与可以放入圆圈内的正方形的最大尺寸有关。圆的直径是馈入iRadon的2D图像的宽度。