没有合适的资源?快使用搜索试试~ 我知道了~
问题背景X射线是一种能够穿透物体的能量射线,并会产生衰减、折射、散射等物理现象。计算机断层成像用畃畯畭異畵畴略畤 畔畯畭畯畧畲畡異畨畹甬 畃畔甩技术依据畘射线穿
资源详情
资源评论
资源推荐
畃畔系统标定与图像重建
摘 要
为解决畃畔系统参数标定及成像问题,本文由相关理论基础入手,剖析模板形态学特
征,以此为切入点建立并优化了标定模型,并且对根据未知介质进行图像重建,根据图形
特性实现了噪声去除。
针对问题一,我们根据基础理论知识,结合几何计算,得出标定模型在各个照射角度
中的投影强度表达式,尝试通过函数拟合的方法求解。在扫描图中,圆和椭圆的 投影函数
相互耦合,为了更好的将其区分开,我们从图像处理的角度入手,通过高通滤波 和形态学
操作大致提取出椭圆曲线边界,利用其宽度信息粗略求出各次投影的角度。再利用粗测的
角度,计算出圆的轨迹,挖去受圆形投影影响的数据点,以此来消除圆的投影函 数对椭圆
投影函数的影响,用剩下的数据对椭圆投影曲线进行拟合,得到了很好的效果,以较高的
精度精确计算出畘光各次照射角度及其他标定相关的参数。
针对问题二、三,考虑到逆畒畡畤畯畮变化的卷积特性,首先使用线性插值得到新的图像
作为逆畒畡畤畯畮变换的输入,再经由坐标变换,内插缩放、滤波变换等操作得到最终的介质
吸收率分布图。
针对问题四,模型检验部分,我们从误差理论角度进行分析,研究标定模型的敏感性,
并尝试人为引入噪声,验证模型的稳定性。在研究过程中我们发现了原定标模型在特定角
度附近精度大幅下降的问题,针对该问题,尝试设计出“三角形”标定模板等若干种模板,
最终兼顾合理性及稳定性,确定“双椭圆”模板作为新标定模板,从理论上证明其稳定性
优于原模版。
关键词 畃畔系统、畒畡畤畯畮变换、图像处理、最小二乘拟合、插值、滤波
由
更多数学建模资料请关注微店店铺“数学建模学习交流”
https://k.weidian.com/RHO6PSpA
1 问题重述
1.1 问题背景
X射线是一种能够穿透物体的能量射线,并会产生衰减、折射、散射等物理现象。计
算机断层成像用畃畯畭異畵畴略畤 畔畯畭畯畧畲畡異畨畹甬 畃畔甩技术依据畘射线穿过待测物体发生吸收衰减的
原理,可以在不破坏物体内 部结构的情况下对其进行三维可视化成像。畃畔技术作为一种
重要的无损检测技术,广泛应用于医学成像、工业探伤、货运安检以及文物复原等领域。
畃畔技术的发展大致经 历四个阶段:第一代畃畔设备通过平行束平移旋转扫描获得投影
数据;第二代设备使用小角度扇形射线束代替平行束;第三代设备仅包含扇形束的旋转扫
描动作,不再采用平移运动;第四代设备将探测 器固定于申甶田
◦
圆周上,仅旋转畘射线源以
解决环形伪像问题。由于被测物与扫描环境的复杂性,畃畔扫描方式日趋灵活。
畃畔技术快速发展,对畃畔设备的精密度要求随之提升。由于安装过程中存在误差,实
际的畃畔成像系统不满足理想成像关系,会使断层图像的重建质量受到影响。因此,使用畃畔系
统前需要借助已知结构的样品进行标定,修正安装误差。
1.2 相关信息
典型初代二维畃畔系统如图??所示,探测器平面由电由甲个等距接收点组成,畘射线垂直
入射。发射器与探测器相对位置固定不变,探测系统绕固定旋转中心逆时针旋转由甸田次。对
每一个畘射线方向,探测器测量经位置固定不动的二维待检测介质吸收衰减后的射线能量,
并经过增益等处理后得到由甸田组接收信息。
图 由町 畃畔系统示意图
甲
1.3 需解决的问题
本文将题述问题归结为以下三个部分,建立数学模型进行分析研究。
问题一:参数标定
根据标定模板的几 何形状与附件由给出的吸收强度, 以及附件甲给出的探测 器接收信
息,标定畃畔系统的 相关参数,包括旋转中心的位置、探测器单元间距、畘射线的由甸田个投
射角度。
问题二、三:CT成像
利用已标定的畃畔系统参数,结合附件申、附件电给出的未知介质接收信息进行成像,确
定未知介质的位置、几何形状与吸收率等信息,并具体给出附件甴所给位置的吸收率。
问题四:模型分析与改进
分析问题一中参数标定的精度与稳定性,自行设 计新模板、建立对 应标定模型以做出
改进。
2 问题分析
本问题主要研究畃畔系统的参数标定与成像。首先通过投影信息确定畘射线入射角度及
探测器单元间距,通过图像重建确定旋转 中心位置,完成系统参数标定;其次由投影信息
直接重建图像,获取未知介质信息;最后分析参数标定的精度与稳定性,设计新 模板并建
立标定模型。
2.1 问题一
分析问题一,建立投影强度曲线的数学模型。首先分析椭圆投影宽度随畘射线角度的
变化规律,以粗略确定各入射角度及旋转 中心。在各入射角度对应的投影强度曲线中,通
过图像重建去除圆形模板影响的点位,拟合椭圆投影强度曲线,求解增益系数、 探测器单
元间距、入 射角度的关系。其次在原强度曲线中,扣除椭圆投影的影响,得到由甸田组精确
的圆形投影强度曲线。拟合圆形投影强度曲线,通过投影宽度求解探测器单元间距,从而
精确求解增益系数与畘射线入射角度。最后通过图像重建,确定旋转中心位置。
申
2.2 问题二、三
畃畔系统的工作原理是由投影重建图像,针对平行束系统的重建算法包括直接反投影
法、滤波反投影法、卷积反投影法等,相关文献中已有详细论述。针 对问题二应以投影信
息与标定所得系统参数作为输入,通过逆畒畡畤畯畮变换重建图像,并对图像进行滤波去噪处
理,以得到未知介质的信息。
2.3 问题四
分析问题一中参数标定的精度与稳定性,针对确定θ取值时所用的畡畲畣畣畯畳θ函数在极值
点附近 对误差的敏感性进行改进。 提出以斜甴电
◦
椭圆代替圆形的 新模板,建立相应标定模
型,对其精度与稳定性进行评估。
3 假设与符号
3.1 模型假设
由甮 入射的畘射线完全平行,忽略相互干涉。
甲甮 探测器单元阵列与二维待检测介质处于同一平面中,忽略载物台在任意角度可能发
生的倾斜。
申甮 附件所给接收信息均为精确值。
3.2 符号说明
4 模型建立与求解
4.1 CT成像的理论基础
X射线与物质的相互作用
当畘射线照射在被检测物体上时,一部分射线能量被物体吸收,使得射线强度发生衰
减。其衰减遵循界畡畭畢略畲畴甭畢略略畲吸收定律畛?畝,呈指数变化。如图??所示,设畘射线初始强度
为I
0
,穿过厚度为x,吸收系数为µ的均匀介质后强度变为I,则有:
I 甽 I
0
e
−µx
用由甩
µx 甽 ln
I
0
I
用甲甩
甴
符号 符号说明 单位
甁d 探测器单元间距 mm
θ
i
第i组畘射线的入射方向 rad
a 椭圆半长轴 mm
b 椭圆半短轴 mm
R 圆半径 mm
D 椭圆投影宽度 mm
f用x, y甩 吸收率函数,表示介质吸收强度的分布
g用s, θ甩 投影函数,表示计算所得的投影强度值,畳为真实投影
坐标
h用l, θ甩 采样函数,表示探测器采样数据拟合所得曲线,畬为像
素坐标,为真实投影长度的
1
/∆l倍
k 增益系数,探测器接收信息与投影强度的比值 mm
−1
k
µ
比例系数,表示由采样信息重建的图像与吸收率数据
的比值
表 由町 符号说明
对于二维平面内的不均匀介质,吸收系数表示为厚度x的函数µ用x甩,此时界畡畭畢略畲畴甭畢略略畲吸
收定律表示为:
Z
µ用x甩dx 甽 ln
I
0
I
用申甩
式??表明,畘射线穿过物体后的射线强度与被检测物体的吸收系数相关,畃畔成像的物
理本质就是通过探测器接收的强度信息计算吸收系数µ的分布。
数据意义的分析解释
相关文献显示, 探测器接收的投影信息应与吸收系数µ用x甩沿投影方向的线积分成正
比,现有图像重建算法均以此为前提。针对本题所给数据,分析附件由可知标定模板吸收
率处处相同;将附件甲中同一入射角度对应的全部接收值求和(如图??所示),可见接收值
总和基本保持不变。可以认为,本题所给数据存在类似关系,即探测器接收信息与吸收率
函数f用x, y甩沿畘射线入射方向的线积分(即投影强度)成正比:
Z
f用x, y甩dl 甽
由
k
h用l, θ甩 用甴甩
式中h用l, θ甩为探测器接收信息,f用x, y甩为介质的吸收率分布,k为增益系数。
电
剩余46页未读,继续阅读
艾斯·歪
- 粉丝: 34
- 资源: 343
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论0