函数 radon 用来对给定的二维矩形数组生成一组平行射线投影。
从而实现前面所说的由投影来重建图像
文章目录
一、概念
函数 radon 用来对给定的二维矩形数组生成一组平行射线投影。
该函数的基本语法:R = radon(I,theta)
,其中,I 是一个二维数组,theta 是角度值的一维数组。投影包含在 R 的列中,生成的投影数等于数组 theta 中的角度数。生成的投影长到足以在射线约束旋转时跨越观察的宽度。当射线垂直于矩形数组的主对角线时会出现这种视图。换句话说,对一个大小为 MxN 的输入数组,投影的最小长度为 [M^2 + N^2]^(1/2)
。
当然,其他角度的投影事实上要短得多,且它们要用零来填充,以便所有投影的长度都相同(如要求的那样,R 应是一个矩形数组)。由函数 radon 返回的实际长度,要稍大于每个像素的单位面积计算的主对角线的长度。
函数 radon 一个更一般的语法:[R,xp] = radon(I,theta)
,其中 xp 包含沿着 x’ 轴的坐标值,xp 中的值用于标注图轴。
在 CT 算法模拟中,用来生成一幅著名图像(称为 Sh