投影仿真1.0-二维图像矩阵投影

2D CT投影仿真

一、报告背景及要求

利用给定Shepp-logan头模型图像（256*256的矩阵），利用离散投影方法，编程计算得到如下几何关系下的扇束CT投影：

1. 输入：Shepp-logan图像256×256像素
   plt.imshow(initial_data, cmap=plt.cm.gray, vmin=0.98, vmax=1.05) % 画图方法
2. 输出：扇束投影数据矩阵360*300。（因为探测器是360度旋转扫描，每1度一次角度投影，一次可投影到300个探测器点上）

二、核心算法及实现代码（完成报告用到的核心数学方法及实现该核心算法的代码）

1、核心算法（数学方法）
1.每一根射线下的密度积分计算：
把一根连续的直线进行离散化，并取x轴方向步长step为0.01cm(因为图像边长为20cm，图像的矩阵shape为256×256，所以一个索引的长度为0.078125cm，一个长度单位cm有12.8个索引（元素边长）)，得到n个点，由每个点的坐标位置得到该点在矩阵的什么位置（矩阵索引），由此可得矩阵的元素值。由直线x方向的步长和斜率可以用勾股定理得到每一步长下的射线穿行线长