滤波反投影算法的Matlab实现

课程大作业，正好记录下～

原理部分（分2块）

参考书：《数字图像处理》–刚雷萨斯
《医学成像的基本原理》–黄力宇

1.雷登变换的实现

原理

雷登变换是一种描述投影的数学方法。

其原理示意图如上，就是对 ø = 1 : step : 180 中的每个ø，图片在这个角度上的多条线积分构成一组投影，共有180组，即构成原图的雷登变换。

初步想法如下：
1.既然要求线积分，那么关键在于积分的上下限，那么对于原图在某条线上的线积分，其上下限可以由这条线与图片边界的交点确定。
2.得到上下限后，积分的方法就相对灵活了，我的想法是，可以将原图中由上下限确定的小矩阵取出来，用Matlab的diag（）函数构造其对角阵，再将对角阵中元素相加，即可得到线积分。

这个方法中的积分部分确实有点意思，但是实现起来发现计算上下限的过程却十分蠢，因为交点存在6种情况，且其中的精度控制和异常处理都很复杂，遂放弃。