数字图像处理实验五.docVIP

数字图像处理 实验五实验目的 学会反投影重建算法的基本概念、 投影定理； 学习卷积反投影重建的概念及意义，图像 实验原理 1、投影重建算法的基本内容： 断层平面中某一点的可看作这一平面内所有经过该点的射线投影之和（的平均值）。 2重建的算法如下所示： (1) (2) 第一种方法中的二维滤波器较难实现。第二种方法交换了滤波与反投影的顺序，此时投影数据是一维的，易于实现。此方法的理论基础就是我们即将提到的投影定理（中心切片定理） 3、卷积反投影重建的步骤： 1）把在固定视角下测得的投影经过滤波，得到滤波后投影； 2）对每一个，把反投射于满足的射线上的所有各点； 3) 将步骤2)中的反投影值对所有进行累加，得重建后的图像。 计算机实现滤波反投影的主要步骤： 1）卷积计算 实验中采用滤波函数，用滤波函数重建的图像其振荡响应减小，对含噪声数据的重建质量也较滤波函数情况为好，但在低频段不及滤波函数的重建质量高。 实验中卷积时是采用的时域方法，利用每个座标值的相互影响的叠加而进行的。 2）射束计算与内插 由于是空间中任意一点的像素坐标，故按照算得的值并不正好为的整数倍，故需要进行内插。 确定后，其可能在和之间，取定，有以下线性内插公式： 3）反投影重建 fp= fopen(atten.rad,r); I=fread(fp,16384,float32); fclose(fp); for