本博客下载链接包含修改的word版本, 可免费下载阅览学习, 也可作为数学建模相关课程作业修改上交:
链接:https://pan.baidu.com/s/1HxzDk3q0p6y2xpuJyxPgvw?pwd=qtnc
提取码:qtnc
本题针对一种二维CT获取样品内部结构信息的工作方式及成像原理,意在通过借助已知结构样品进行参数标定,消除系统误差,而后对未知结构的样品进行成像,得出该未知介质的相关信息。
问题一中,要求通过标定模板的相关,确定CT系统的旋转中心、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。本文利用标定模板几何参数,通过条数、以及探测器单元间距相等等信息,首先计算出探测器单元之间的距离为0.2778mm。根据X射线与标定模板的几何特性,以椭圆短轴所在直线为x轴,长轴所在直线为y轴建立直角坐标系,得旋转中心在所建立的坐标系中的坐标为(-9.0396,5.7060)。最后通过近似确定X射线180个方向的旋转角度具有高度线性相关性,根据部分确定数据拟合出整体旋转角度。前十五组旋转角度为:
29.6493 31.0027 31.5572 32.6476 33.6795 34.6478 35.6483 36.6490 37.6475 38.6486 39.6488 40.6491 41.6487 42.6497 43.6492
问题二中,要求通过已求得的标定参数,确定未知介质在正方形托盘中的位置、几何形状和吸收率等信息,并具体给出所要求十个位置的吸收率。本文依据CT层析成像原理,利用逆拉东变换作出重构图像,并利用Excel中数据分布计算出原位置介质的相关性质。所求十个位置的吸收率为:
序号 | 1 | 2 | 3 | 4 | 5 |
吸收率 | 0.0185 | 0.9701 | -0.0001 | 1.1703 | 1.0287 |
序号 | 6 | 7 | 8 | 9 | 10 |
吸收率 | 1.4425 | 1.2815 | -0.0074 | -0.0025 | 0.0194 |
问题三中,与问题二计算方法相似得十个位置的吸收率为:
序号 | 1 | 2 | 3 | 4 | 5 |
吸收率 | 0.0736 | 2.5305 | 7.1327 | 0.0332 | 0.8015 |
序号 | 6 | 7 | 8 | 9 | 10 |
吸收率 | 3.1487 | 5.5714 | 0.0530 | 8.1362 | 0.0272 |
问题四中,要求对问题一中的标定模型进行改进以减小误差并增加稳定性,本文利用在问题一求解过程中遇到的问题进行思考,首先适当增大模板减小偶然误差,其次做出投射图像后应容易找到极值,并且图像应有一定的对称性;图像扫描后应尽量少地得出重复数据。
关键词:CT层析成像 Radon变换与逆变换 吸收率 MATLAB算法