文档CT图像滤波反投影重建算法的研究[宝典]
滤波反投影法重建CT图像实验指导书一、实验目的1.了解傅立叶.docx

滤波反投影法重建CT图像实验指导书一、实验目的1.了解傅立叶变换法、直接反投影法重建CT图像的原理;2.掌握滤波反投影法重建CT图像的原理和基木方法。
二、实验器材装有MATLAB程序的PC机,滤波反投影法图像重建演示软件,投影数据。
三、实验原理CT图像重建问题实际上就是如何从投影数据中解算出成像平面上各像索点的衰减系数。
图像重建的算法有多种,如反投影法、傅立叶变换法、迭代法、滤波反投彩法等。
在介绍算法前,有必要先介绍从投影重建图像的重要依据,即中心切片定理。
1.中心切片定理密度函数/(x,)0在某一•方向上的投影函数g&(/?)的一•维傅立叶变换函数g&(p)是原密度两数f(x,y)的二维傅立叶变换函数F(p,&)在(Q,&)平面上沿同一•方向且过原点的直线上的值。
图1屮心切片定理2.傅立叶变换法如果在不同角度下取得足够多的投影函数数据,并作傅立叶变换,根据小心切片定理,变换后的数据将充满整个(仏叭平而。
一旦频域函数F(w,v)或F(/?,0)的全部值都得到后,将真做傅立叶反变换,就能得到原始的密度函数.f(x,y),即所要重建的图像。
上述图像重建算法称为傅立叶变换法,图2给出了傅立叶变换重建方法的流程图。
图屮指出,对于每次测得的投影数据先作一维傅立叶变换。
根据中心切片定理,可将此变换结果看成二维频率域屮同样角度下过原点的直线上的值。
在不同投影角卜•所得的一维变换函数可在频域中构成完整的二维傅立叶变换函数,将此二维变换函数做一次逆变换,就得到了所要求的空间域屮的密度函数。
为了在二维逆变换屮采用快速傅立叶变换算法,通常在逆变换前要将极坐标形式的频域函数变换成直角坐标形式的数据。
图2傅立叶变换重建图像的过程采用傅立叶变换法重建图像吋,投影函数的一维傅立叶变换在频域中为极处标形式, 把极坐标形式的数据通过插补运算转换为直角坐标形式的数据时,计算工作量较大。
此外, 在极坐标形式的频域数据屮,离原点较远的频率较高的部分数据比较稀疏,当这些位登上的数据转换到肯角坐标下时,需经插补,这将引入一定程度的谋差。
CT系统标定参数与图像重建的探究

CT系统标定参数与图像重建的探究摘要:CT(计算机层析成像)是通过对物体进行多角度的射线投影测量,从而获得物体横断面信息的成像技术。
为将此技术在多个领域进行推广,本文简述了对CT系统的标定参数和图像重建展开研究,关键词:射线光学;图像重建模型;像素矩阵;平行束投影引言自20世纪70年代初第一台CT扫描系统问世以来,其结构和扫描方式不断革新,且性能逐渐强大,至今已发展到第五代电子束CT机。
本文所研究的对象为第一代CT机。
第一代CT机采用平行光束[2]的扫描的方式,按照旋转-平移的方法工作,属于头颅专用机,且工作周期较长,得到的图像较为粗糙。
此外,在CT系统的安装过程中,往往存在误差,其影响图像质量,需引入已知模板来确定CT系统的标定参数,并根据此在图像重建的过程中,对重建后的图像进行修正。
1 模型建立1.1 问题一根据标定模板的几何信息以及其接收数据信息,利用统计学和几何学原理,逐步确定探测器单元之间的距离、CT系统旋转中心在正方形托盘中的位置和CT系统所使用的X射线的180个方向。
1.1.1 探测器单元间距的确定根据统计学知识,对附件2中的数据进行分析,发现规律:标定模板中的小圆在多个方向上的投影中,其所对应的探测器数目恒为29个(极少数数目为28个的情况忽略不计),即其宽度始终不变。
结合题干中提供的标定模板的几何信息,取模板中小圆的直径为8mm。
又因为探测器单元等距排列,从而可在忽略探测器长度的情况下,确定探测器单元的间距。
1.1.2 CT系统旋转中心位置的确定CT系统的旋转中心是必需的标定参数,其在图像重建中起着至关重要的作用,若不对该参数进行精准确定,所得到的重建图像将存在伪影,易造成测量人员判断错误。
由于本文所研究的第一代CT机与第二代在原理上无本质的区别,遂结合标定模板的几何信息以及原始投影数据,借鉴了李堡磊和孟凡勇等人的方法,基于原始平行束投影数据的旋转中心确定法。
1.1.3 系统中X射线的180个方向的确定为简化计算的步骤,利用几何学的原理,以标定模型中的椭圆中点为原点再次建立平面直角坐标系。
反投影重建算法范文

反投影重建算法范文反投影重建算法的原理是基于CT扫描的物理特性,即通过测量X射线在物体内的吸收特性来重建图像。
CT扫描使用X射线源从不同角度射入物体内,并在探测器上测量射入和穿过物体的X射线的强度。
这些测量数据被称为投影数据,并被用来重建图像。
反投影重建算法的基本思想是将投影数据进行反投影运算,将每一个探测器测量的X射线强度加到相应的图像位置上。
这样,通过不同角度的反投影运算,可以得到一系列的反投影图像。
然后,将这些反投影图像进行加和处理,得到最终的重建图像。
具体而言,反投影重建算法的步骤如下:1.从CT扫描设备中获取投影数据,包括不同角度下的X射线测量值。
2.对每一个角度的投影数据进行反投影运算。
即将每一个探测器测量的X射线强度加到相应的图像位置上。
3.对所有角度的反投影图像进行加和处理。
即将所有反投影图像的对应位置像素值相加。
4.对加和后的图像进行滤波处理,以降低噪声。
5.得到最终的重建图像。
反投影重建算法的优点是简单易实现,并且可以用于各种不同的CT扫描设备和投影模型。
然而,该算法也存在一些缺点,例如由于对雷达扫描角度的限制,可能无法获取到完整的投影数据;同时,由于探测器的离散性,重建图像可能存在伪影等问题。
为了克服反投影重建算法的缺点,近年来还出现了许多改进和优化的算法。
其中,滤波反投影算法是最常用的一种改进算法之一、滤波反投影算法通过在投影数据中引入滤波操作,来减小伪影和噪声的影响,从而提高重建图像的质量。
另外,统计模型和机器学习方法也广泛应用于图像重建,以提高图像的准确性和鲁棒性。
总之,反投影重建算法是一种广泛应用于医学图像重建的技术。
虽然该算法存在一些缺点,但通过改进和优化,可以得到高质量的重建图像,为医学诊断和临床应用提供有力支持。
反投影重建算法

反投影重建算法反投影重建算法是一种常用于医学影像领域的图像重建方法,其原理基于X射线的穿透性质和成像原理。
该算法可以从一系列投影图像中恢复出被测物体的内部结构,为医学诊断提供重要的信息。
在了解反投影重建算法之前,我们先来了解一下X射线成像的基本原理。
X射线成像是一种无创检查方法,通过将X射线束通过被测物体后,测量X射线在不同方向上的衰减情况,从而获取被测物体的内部结构信息。
在传统的X射线成像中,使用的是传统的透射成像方法,即将X射线通过被测物体后,直接在胶片或探测器上记录吸收信息。
然而,传统透射成像存在一些局限性,如对于密度相近的组织无法进行有效区分。
反投影重建算法的出现就是为了解决传统透射成像的局限性。
该算法利用了被测物体对X射线的吸收特性,通过在不同方向上测量X 射线的衰减情况,然后将这些投影数据进行反投影操作,即根据投影数据在不同方向上的衰减信息,将其反向投射回原来的空间位置。
这样,就可以通过反投影操作恢复出被测物体的内部结构。
具体来说,反投影重建算法的步骤如下:1. 收集投影数据:通过将X射线束从不同方向照射到被测物体上,并在胶片或探测器上记录下不同方向上的吸收信息,形成一系列投影数据。
2. 反投影操作:将投影数据进行反投影操作,即将各个投影数据反向投射回原来的空间位置。
这个过程可以看作是将每个投影数据沿着其对应的投影方向,从胶片或探测器上反向传播到被测物体内部的过程。
3. 累加重建:将所有反投影后的数据进行累加,得到最终的重建图像。
在反投影操作中,相交的区域会累加,而非相交的区域则不会累加,从而在最终的重建图像中得到被测物体的内部结构信息。
反投影重建算法是一种简单而有效的图像重建方法,在医学影像领域得到了广泛的应用。
与传统的透射成像相比,反投影重建算法可以提供更多的信息,能够更好地区分密度相近的组织,从而提高了诊断的准确性。
此外,反投影重建算法还可以应用于其他领域,如安全检查、材料检测等。
基于Radon变换的滤波反投影重建算法研究

基于Radon变换的滤波反投影重建算法研究作者:张晓瑞来源:《电脑知识与技术》2016年第27期摘要:层析成像技术目前在生物学等领域广泛应用,而Radon变换是投影重建的数学理论基础。
文中简单介绍了Radon变换的基本原理,重点介绍了滤波反投影重建算法,尤其是滤波函数的选取,并对重建结果进行了比较。
关键词:图像重建;Radon变换;滤波反投影中图分类号:TP311 文献标识码:A 文章编号:1009-3044(2016)27-0259-031 引言层析成像技术在医疗、生物等领域具有广泛应用。
图像重建是指通过物体外部测量数据,经过处理从而获得物体的形状信息的技术。
开始主要应用在放射医疗设备中,用于人体各部分图像的显示,即计算机断层摄影技术,简称CT技术,后来逐渐在许多领域获得应用。
透射CT的理论基础是投影重建。
而Radon变换是投影重建的数学基础,它是数学家J.Radon提出来的,被广泛应用于医学、分子生物学等领域,迄今为止,人们已研究出基于Radon变换的多种重建方法。
文中重点介绍的滤波反投影算法也是基于Radon变换的一种变换法重建,目前在CT系统中应用非常广泛。
滤波反投影算法的比较重要的是滤波函数的选取。
4.3 结果分析(1)图像比较:直接反投影算法的重构对像的边缘很不明显,有阴影。
滤波反投影算法重构对象相对来说清晰很多,没有阴影。
(2)重建时间对比:滤波反投影算法的重构时间较长,因为多了卷积、滤波这个步骤,使重构时间加长。
不过在实际应用中,这个时间增加不会有很大影响,但是质量却明显变好,所以,实际应用中一般采用此方法。
5 总结滤波反投影法是重构图像基本常用的算法,也是其他多种算法的基础。
在医学CT 等领域中的应用较为广泛。
但是这种算法的关键是选取的何种滤波函数,会直接影响重建图像的质量。
除了滤波函数对图像质量有着较大的影响外,根据抽样定理,投影数和抽样间距均对重建图像的质量有影响[2]。
以后的工作中也应对抽样间距进行研究。
CT反投影滤波重建算法设计

地理与生物信息学院2011 / 2012 学年第二学期实验报告课程名称:医学图像处理和成像技术实验名称:CT反投影滤波重建算法设计班级学号: B********学生姓名: ****指导教师: ****日期:2012 年 4 月一、实验题目:CT反投影滤波重建算法设计二、实验内容:1.显示图像;2.获得仿真投影数据;3.基于获得的仿真投影数据重建图像。
三、实验要求:1.画出Shepp-Logan头模型,显示尺寸为128×128;2.从头模型中获得投影数据,投影数据格式为180×185;3.基于获得的仿真投影数据重建图像,使用R-L卷积函数,重建尺寸为128×128。
四、实验过程:1.显示图像:①算法实现流程:I. S-L头模型由10个位置、大小、方向、密度各异的椭圆组成,象征一个脑断层图像。
将模型中的椭圆参数写入一个p矩阵中,方便使用其中的数据,并设定所需参数。
II. 使用循环语句给像素赋值:for i=1:10for x….for y…..判断点(x, y)是否在第i个椭圆内;如是,则将第i个椭圆折射指数赋给点(x, y);endendendIII. 显示仿真头模型,使用imshow(f,[])函数显示出图像。
②实验代码:clear all;p=[0 0 0.92 0.69 pi/2 10 -0.0184 0.874 0.6624 pi/2 20.22 0 0.31 0.11 72/180*pi 0-0.22 0 0.41 0.16 108/180*pi 40 0.35 0.25 0.21 pi/2 50 0.1 0.046 0.046 0 60 -0.1 0.046 0.046 0 7-0.08 -0.605 0.046 0.023 0 80 -0.605 0.023 0.023 0 80.06 -0.605 0.046 0.023 pi/2 8];N=256;x=linspace(-1,1,N);y=linspace(-1,1,N);f=zeros(N,N);for i=1:Nfor j=1:Nfor k=1:10A=p(k,3);B=p(k,4);x0=p(k,1);y0=p(k,2);x1=(x(i)-x0)*cos(p(k,5))+(y(j)-y0)*sin(p(k,5));y1=-(x(i)-x0)*sin(p(k,5))+(y(j)-y0)*cos(p(k,5));if((x1*x1)/(A*A)+(y1*y1)/(B*B)<=1) %判断条件f(i,j)=p(k,6);endendendendf=rot90(f);imshow(f,[])③运行结果:2. 获得仿真投影数据:①算法实现流程:I. 取θ∈ [00, 10, ..., 1790], s ∈[-92, -91, ..., 91,92],即表示在[0 1790]范围内角度每隔10度取样,每个角度下有185个探测器。
基于滤波反投影算法的CT系统成像研究

基于滤波反投影算法的CT系统成像研究摘要:CT系统的安装会使得旋转中心发生偏离,从而影响成像质量,因此需要借助于已知结构的样品来标定CT系统的参数,并且利用标定的参数对未知结构的样品进行图像重建。
首先根据直接反投影算法和滤波反投影算法对收集到的数据中的接收信息分别进行图像重建,通过成像图像可知,滤波反投影算法更优;旋转中心可能发生偏移以及CT系统具有初始角度,依次进行旋转、平移、裁剪和残影去除操作,来校正投影图像,从而得到较高质量的图像。
关键词:CT成像原理(影像医学与核医学);滤波反投影法;图像重建;吸收率引言CT(Computed Tomography)是用X线束从多个方向对人体检查部位具有一定厚度的层面进行扫描,由探测器而不用胶片接收透过该层面的X线,转变为可见光后,由光电转换器转变为电信号,再经模拟/数字转换器转为数字,输人计算机处理。
数字矩阵中的每个数字经数字/模拟转换器转为由黑到白不等灰度的小方块,称之为像素,并按原有矩阵顺序排列,即构成CT图像。
所以,CT图像是由一定数目像素组成的灰阶图像,是数字图像,是重建的断层图像。
首先根据直接反投影算法和滤波反投影算法对收集到的数据中的接收信息分别进行图像重建,将图像重建[4-6]的两种结果进行对比,得出效果较好的模型;然后,旋转中心可能发生偏移以及CT系统具有初始角度,通过旋转、平移、裁剪和残影去除等操作来校正投影图像,最后对图像进行标准化调整,从而提高了成像质量。
1 模型的准备与建立1.1 CT成像的数学基础Rand变换如图1所示,直线g是xOy平面内任意一条直线,t是原点到直线g的距离,φ为原点到直线g的垂线与x轴的夹角。
对于xOy平面内任意一条直线可以由(t,φ)唯一确定。
二维平面中函数f(x,y)沿着直线的积分等于其Rand变换。
中心切片定理中心切片定理是CT图像重建算法的基础,在非衍射源情况下,含义是图像在某个视角下平行投影的一维Fourier变换等同于该图像二维Fourier变换的一个中心切片。
滤波反投影重建对CT系统研究

滤波反投影重建对CT系统研究确定CT系统旋转中心的位置、探测器单元的间距及X射线的180个方向[1];利用MATLAB获得模板示意图,对其处理得探测器单元间隔;建立坐标系,让X 射线入射光线与椭圆和圆的交线长度为约束条件,运用最小二乘法求得目标函数,利用遗传算法在MATLAB实现得到X射线的180个方向;重复上述方法旋转中心满足的约束条件和目标函数,得出旋转中心坐标,确定未知介质在正方形托盘的位置、几何形状和吸收率并利用滤波反投影法中的函数实现重建图像求得该介质的吸收率情况和10个位置处的吸收率.标签:遗传算法;滤波反投影;CT系统;最小二乘法0 引言CT系统就是利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像,可以获取样品内部的结构信息。
但是CT系统安装时容易出现误差,所以要进行参数标定,并据此对未知结构的样品进行成像。
在正方形托盘上放置两个均匀固体介质组成的标定模板,确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离及该CT系统使用的X射线的180個方向。
1.求探测器单元之间的距离在正方形托盘上放置两个均匀固体介质组成的标定模板的几何信息,可以利用MATLAB对其进行计算绘制出该模板示意图(见图1)通过对比已知模板示意图可知右侧介质直径为8毫米。
然后对每一点的吸收率进行整理,利用MATLAB中的函数,利用这个函数将矩阵A中的元素数值按大小转化为不同颜色,并在坐标轴对应位置处以这种换颜色染色[2],其中决定坐标范围,是两个二维向量,即的范围内染色,根据图像及处理后的数据发现有84个角度没有穿过椭圆部分,对这84个角度进行统计得到有71个角度有29个探测器单元数据,对统计数据求取平均值估计探测器单元间隔d代入数据可得探测器单元间隔。
2.CT系统使用的X射线的180个方向通过CT成像的基本原理:X射线在穿过物质时,其强度的衰减与物质的厚度及入射辐射强度成正比,比例系数记作k,称为衰减系数,也等同于吸收系数,它代表X射线穿过单位厚度的物质时强度减少的比例,则吸收率W为W=kl根据为每一点的吸收率,则可以知道某次旋转时第个探测器上某点的吸收率和第个探测器上该点的吸收率。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
文档CT图像滤波反投影重建算法的研究[宝典] 西北工业大学学位论文知识产权声明书本人完全了解学校有关保护知识产权的规定,即:研究生在校攻读学位期间论文工作的知识产权单位属于西北工业大学。
学校有权保留并向国家有关部门或机构送交论文的复印件和电子版。
本人允许论文被查阅和借阅。
学校可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。
同时本人保证,毕业后结合学位论文研究课题再撰写的文章一律注明作者单位为西北工业大学。
保密论文待解密后适用本声明。
学位论文作者签名::《~~~塑, 指导教师签名。
,况订年弓月多,日聊年岁月歹口日西北工业大学学位论文原创性声明秉承学校严谨的学风和优良的科学道德,本人郑重声明:所呈交的学位论文,是本人在导师的指导下进行研究工作所取得的成果。
尽我所知,除文中已经注明引用的内容和致谢的地方外,本论文不包含任何其他个人或集体己经公开发表或撰写过的研究成果,不包含本人或其他已申请学位或其他用途使用过的成果。
对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式表明。
本人学位论文与资料若有不实,愿意承担一切相关的法律责任。
学位论文作者签名网年;月如日西北〕,业大学硕士学位论文第一章绪论第一章绪论本章首先概述了图像重建和,,技术的基本概念及图像重建方法的分类,然后介绍了,,的发展及国内外研究现状,最后阐明了本文的研究目的与意义、主要工作和内容安排。
,(,图像重建与,,技术概述由物体的二维截面或断面向该平面内的各个方向作投影,可获得一系列一维投影函数。
由这些一维投影函数来重建该二维截面则称为图像重建。
该技术是随着计算机技术的进步而发展起来并获得广泛应用的图像处理技术,其最典型的应用是医学上的计算机断层成像术(;,,,,,,, ,,,,鲫,,或;;胁,,,,,,,,,(,,,,肿,,,),简称为;,技术。
它用于获取人体头颅、心肺、腹部等内部器官的二维断层图像(故亦称断层摄影技术),对于,射线放射诊断是一个重大突破,具有深远的实际意义,因而受到普遍的重视。
,(,(,图像重建图像重建是图像处理中一个重要研究分支,是指根据对物体的探测获取的数据来重新建立图像,其重要意义在于获取被检测物体内部结构的图像而不对物体造成任何物理上的损伤。
由于具备无损检测技术的显著优点,它在各个不同的应用领域中都显示出独特的重要性。
例如:在医疗放射学、核医学、电子显微、无线电雷达天文学、光显微和全息成像学及理论视觉等领域都多有应用。
在上述的众多领域中,图像重建在医学方面的应用最为显著。
它大大丰富了对于人体内部器官进行无损检测的方法和手段,为疾病的早期正确诊断提供了科学的、准确;的依据。
根据原始数据获取方法及重建原理的不同可分为如下几种:发射断层重建成像(,,,醛,?,,,,,;, ,,,,,,,,,,,,;砷,透射断层重建成像,,,,,,,,,,,,, ;;吼,,,,;, ,,,,掣,,,,,,?,反射断层重建成像(,,,,,;,,?,;粕,,,,,劢衄,,,,,,,,,;,)及核磁共振重建成像(,,印,,,; ,,,?卸;; ,,,,,,,,,,,,,。
西北工业大学硕士学位论文第一章绪论,(,(, ,,技术在各种图像重建算法中,计算机断层成像术即,,技术占有重要的地位。
;,技术的功能是将人体中某一薄层中的组织分布情况,通过射线对该薄层的扫描、检测器对透射信息的采集、计算机对数据的处理,并利用可视化技术在显示器或其他介质上显示出来。
,(,(,图像重建的方法图像重建是,,技术中的一个重要问题。
它的实质是按照采集后的数据,求解图像矩阵中象素,然后重新构造图像的过程;而图像矩阵的求解由计算机完成。
图像重建问题的求解方法根据其特点可分为,大类。
第,类是变换重建方法,也叫解析法呐,其特点是先在连续域解析处理,最后离散化以利用计算机计算。
其中又可分成傅里叶反变换重建法和滤波(或卷积)反投影重建法。
第,类是级数展开重建法,,, ,,,(也叫代数重建法、迭代算法、优化技术等),其特点是从开始就离散化进行分析,从而直接得到数值解。
另外还有将变换法和级数展开法相结合的综合方法,如:迭代变换法、迭代重建重投影、角谐函数重建法和正交多项式展开重建。
变换法的突出优点是实现简单,速度快,对足够精确的投影数据能获得很好的重建质量。
因此目前实用,,系统中,尤其是医用;,系统中广泛采用变换法,特别是滤波反投影类型的算法来进行图像重建。
变换法在技术上有两个主要的限制:(,)噪声特性不好。
因为变换法是基于解析求反公式的闭合形式,要求投影数据是精确的。
对于数据中的噪声,可以通过滤波步骤来适当解决。
如果我们能在投影数据输入给变换法之前,将影响投影数据不精确性的诸多物理因素进行足够的校正,则便可以认为投影数据是相当准确的,从而得到满意的重建。
(乃正因为变换法基于解析求反公式的闭合形式,所以变换法的简单与复杂强烈地依赖于数据采集扫描方式。
如果投影数据不是沿直线的简单线积分,那么可能就得不到解析求反公式的闭合形式,在这样的情况下,变换法就变得无效。
级数展开法能和特定的成像设备及数据采集物理过程的特性相结合,并能利用某些先验信息。
这类方法尤其适合于数据采集过程具有高度统计特性及某些不能获得完整数据集的场合。
并且,级数展开法不像变换法那样强烈依赖于投影数据的采集方式,适合于投影数据不是目标函数沿直线而是沿曲线的积分甚至沿带 ,西北工业大学硕十学位论文第一章绪论状区域的积分。
级数展开法最大缺点是实现起来复杂,速度慢,并且所需的存储空间大。
因此,目前在实用,,系统中多不采用(,(, ,,的发展与研究现状图像重建技术的历史比较悠久。
它的理论起源于,,,,年奥地利数学家,,,?所发表的论文,,,,文中证明了二维或三维物体的密度分布能够通过其无限多个投影来确定,数学上称之为,,,?变换。
但限于当时的技术条件而未能在工程上实现,他的论文一直未被世人重视。
由于科学技术的进步,特别是计算机科学的发展,图像重建问题重新引起了人们的兴趣,相继有不少学者进行了卓有成效的创造性研究。
,,,,年,美国科学家,(,(?,;,首先提出了用,射线投影重建断层图像的计算方法。
,,,,,年,英国肼,公司中央研究所工程师,(,(,叫蚶,,,研制成功诊断头颅用的第一台电子计算机,射线断层摄影装置,并在那一年的英国放射学年会上公布了这一装置及其应用成果,至此,才算真正设计出一个装置用来实现人体断面成像。
在,,,,年,月蒙特利尔召开的第一次国际,,会议上,这一新式的,射线显像技术被正式命名为,;胁,,,,,,(皿,,,,,,,,简称;,。
,,的问世在放射学界引起了爆炸性的轰动,被认为是继伦琴发现,射线后,工程界对放射学诊断的又一划时代的贡献。
,,,,年的诺贝尔生理和医学奖破例地授予两位没有专门医学资历的科学家:,?,,,,,,和,;悯,,;,。
从此,放射诊断学进入了;,时代。
;,装置自问世以来,在其后的十多年内获得了极大发展,已从第一代产品发展到了第五代产品。
除,—,,(用,射线作为辐射源,叫作,射线,,田以外,其他形式的;,也相继问世,如单光子发射;,(,,,;,、正电子发射叫,,,简称,—,,)、核磁共振;,(,,,一叨等均已付诸临床应用。
超声,,、微波;,的研究也取得了极大的进展。
应该指出,各类;,的功能是相互补充而不是相互替代的。
;,作为一种高性能的无损诊断技术显然已在医学成像领域确立了不可动摇的地位。
实际上,临床上使用的,—;,设备的性能近年来又有了明显的提高。
为了进一步提高扫描速度,一些快速扫描装置相继诞生。
螺旋式超高速,,甚至可以对慢运动组织进行商分辨率的检查,并能提供人体组织的三维图像。
,,作为一种技术,既有坚实的数学理论为依托,又有现代微电子技术与计算机技术相支持,它已在工业无损检测、地球物理的研究、农林业和环境保护方 ,西北工业大学硕,学位论文第一章绪论面得到了广泛的应用。
,(,研究目的与意义在图像重建的变换重建方法中,滤波反投影算法(】,,,,,犯,,,;,,,,,;;,,?,,,,)是影像重建中应用最广泛的一种算法,在当代,射线,,系统中几乎都用这种方法构成系统。
滤波反投影算法又称卷积反投影算法,因为频域上的滤波相当于空域上的卷积运算。
滤波反投影法大量实用的原因是它的基本算法很容易用软件和硬件实现,而且在数据质量高的情况下可重建出准确清晰的图像。
,,,,,,变换(,,,,;, ,,姐,,,珊,,,是由著名数学家,(,(,(,,,,,,于,,,,年提出的一种正交变换,,,,它类似于,,,,,,;,变换(,,,,,盱,,柚,,,加,,,,但它定义在实数域内,比,,简单。
离散,变换(,,,;,;,, ,,,,,,,,,,,,,【,,是我国学者王中德提出的实数域上定义的一种正交变换。
,,,,年底,,(,(,,;;,;,,研究了,,,,,,变换的离散形式,称之为离散,,,,,,,变换(,,,;,;,, ,,,,,,,,,,,,叫”,它实际上是,,,的一种特殊形式。
,,,,年,?,,会刊(?,,,黟,,,,,,,, ,,,)出了一期关于,,,,,,,变换的专刊【,,,可见这种变换已在实际应用中得到了广泛重视。
在实序列数据处理中,,,不仅等效于,,,其正反交换又具有相同的形式,而且,,仅需用到实运算,在存储量和复杂性上要比,,更经济,在实数域能够替代,,进行信号和图像的处理。
有学者基于重建图像象素和投影射线之间在不同的投影方位上所存在的几何关系,提出了平行射束常规反投影的多方位同时反投影伽,,,,,,,;,,,,,?,,,,,柚,,,,咖,,,,,,,?,,,,;,,删;;,,?,,,,,)方法,,,。
本文将,,,,,,变换和,,,,方法应用到,,图像的滤波反投影重建方法中,在保持精度不变的前提下,可较大地减少所需的存储空间和计算量。
,(,主要工作和章节安排论文主要完成的工作有: (,) 对图像重建及,,技术的研究和发展情况进行了总结和评述,阐明了进行图像重建的目的和意义。
(,) 对,,图像的反投影重建算法进行了详细的介绍,集中讨论了滤波 ,西北工业大学硕士学位论文第一章绪论反投影重建算法的计算机实现的具体步骤。
(,) 对与,,,,,,,变换平行的一种正弦正交变换—,,,,,,变换进行了探讨,根据用,,,,;,变换导出的滤波反投影重建算法重建,,图像,由此重建算法较,,,,,,,方法将有大约一半计算量和存储空间的节省。
将平行射束常规反投影的多方位同时反投影(,,,,)方法应用于 ;,图像的滤波反投影重建方法中,使象素定位操作中的计算量减少到传统反投影算法的,,,。
(,) 对(刁、(,)中采用的算法分别用,,?曲进行了编程仿真,证明了算法的可行性,并对(,)中的方法与传统反投影方法的计算量和存储空间进行了比较和分析。