如何获取特定线上矩阵的投影?

时间:2017-02-11 18:27:34

标签: matlab projection tomography-reconstruction

给定L(t,θ)上的点的参数化是 x(s):= t cos θ − s sin θ y(s):= t sin θ + s cos θ,其中t是从θ x-axis开始,距离原点到线的距离,而s是该线上的某个点。< / p>

如何使用特定步长Img在此行L(t,θ)处对图像s进行投影。使用这个我必须进一步实现氡变换。

我的问题是如何定义步长st的值? 另外,我是否需要旋转Img或不旋转?

请帮忙。

1 个答案:

答案 0 :(得分:0)

我建议你看看周围的多个开源软件。

在层析成像中,旋转图像与旋转机器相同,因此您可以更改每个角度的源/探测器位置,然后计算连接它们的线。然后,步长由您决定。研究表明(我已对此进行了测试),如果您使用标准像素,则值为s=pixel_size/20.5

如果您正在进行2D平行光束,那么您可以忘记所有需要执行的几何变换并使用imrotate生成投影。如果您使用扇形光束或锥形光束,则代码会变得复杂一些。