基于matlab环境下图像重建及吸收率分析.docVIP

基于matlab环境下图像重建及吸收率分析 　　摘 要：针对重建未知介质的几何形状及吸收率，综合运用滤波反投影算法进行图像重建，首先先引入中心切片定理，采集不同视角下的投影并求解其傅里叶变换，汇集成图像的傅里叶变换再经傅里叶反变换，其次利用matlab iradon函数重建得未知介质的几何形状，再次利用坐标变化和标定参数得出未知介质在正方形托盘中的位置，最后通过matlab rgb2gra函数得到重建图形的灰度值矩阵，由灰度值与吸收率的关系待测未知介质吸收率。 　　关键词：图像重建；滤波反投影法；matlab灰度值；吸收率 　　中图分类号：TP391.41 文献标志码：A 文章编号：2095-2945（2018）06-0193-02 　　Abstract： In view of the geometric shape and absorptivity of the unknown medium， the filter back-projection algorithm is used to reconstruct the image. Firstly， the central slice theorem is introduced， and the projection in different angle of view is collected and the Fourier transform is solved. The Fourier transform which gathers the image is then transformed by the inverse Fourier transform， and then the geometric shape of the unknown medium is reconstructed by using the matlab iradon function. The position of unknown media in square pallet is obtained by coordinate change and calibration parameters. Finally， the gray value matrix of reconstructed graph is obtained by Matlab rgb2gra function. According to the relationship between gray value and absorptivity， the absorptivity of unknown medium is determined. 　　Keywords： image reconstruction； filtering inverse projection method； MATLAB gray value； absorption rate 　　引言 　　CT可以在不破坏样品的情况下，利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像，由此获取样品内部的结构信息[1]。一种典型的二维CT系统[2]如图1所示，相互平行的X射线垂直射入512个等距排列的探测器内。然后X射线沿旋转中心逆时针旋转180次得到180组的接受信息数据。 　　CT系统的参数即为CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。我们在已知标定参数的条件下，利用上述CT系统得到的某未知介质的接收信息512*180组数据，就可得出未知介质的几何形状及其内部吸收率，从而了解内部结构信息，这对现实中我们探索未知介质内部信息拥有很重要的意义。 　　1 未知介质的几何形状求解 　　1.1 滤波反投影重建模型 　　待建图像为f（x，y），它的傅里叶变换为F（？棕1，？棕2）= 　　■（？籽，？?），由中心切片定理有，■（？籽，？兹）可通过f（x，y）在不同视角？准下的投影p？准（xr）的傅里叶变换求得[4]。 　　待建图像： 　　将其写成空域变量为xr的傅里叶反变换式得到： 　　1.2 模型的算法 　　Step 1：利用matlab的iradon函数重建出介质的几何形状。 　　Step 2：利用第一问标定参数旋转中心在托盘中的位置和坐标轴变换得到待测介质在正方形托盘的位置。 　　Step 3：将重建出来的图形转换为灰度图。 　　其中白色正方形线段为100*100（mm）的托盘的大小，里面的灰色部分则为重建出来的未知介质形状的灰度图，圆圈表示托盘的中心位置。 　　2 内部吸收率求解 　　2.1 灰度值求解原理及算法 　　由灰度值[5]概念可知，灰度