CT图像滤波反投影重建算法的研究
- 格式:pdf
- 大小:1.81 MB
- 文档页数:77
反投影重建算法范文反投影重建算法的原理是基于CT扫描的物理特性,即通过测量X射线在物体内的吸收特性来重建图像。
CT扫描使用X射线源从不同角度射入物体内,并在探测器上测量射入和穿过物体的X射线的强度。
这些测量数据被称为投影数据,并被用来重建图像。
反投影重建算法的基本思想是将投影数据进行反投影运算,将每一个探测器测量的X射线强度加到相应的图像位置上。
这样,通过不同角度的反投影运算,可以得到一系列的反投影图像。
然后,将这些反投影图像进行加和处理,得到最终的重建图像。
具体而言,反投影重建算法的步骤如下:1.从CT扫描设备中获取投影数据,包括不同角度下的X射线测量值。
2.对每一个角度的投影数据进行反投影运算。
即将每一个探测器测量的X射线强度加到相应的图像位置上。
3.对所有角度的反投影图像进行加和处理。
即将所有反投影图像的对应位置像素值相加。
4.对加和后的图像进行滤波处理,以降低噪声。
5.得到最终的重建图像。
反投影重建算法的优点是简单易实现,并且可以用于各种不同的CT扫描设备和投影模型。
然而,该算法也存在一些缺点,例如由于对雷达扫描角度的限制,可能无法获取到完整的投影数据;同时,由于探测器的离散性,重建图像可能存在伪影等问题。
为了克服反投影重建算法的缺点,近年来还出现了许多改进和优化的算法。
其中,滤波反投影算法是最常用的一种改进算法之一、滤波反投影算法通过在投影数据中引入滤波操作,来减小伪影和噪声的影响,从而提高重建图像的质量。
另外,统计模型和机器学习方法也广泛应用于图像重建,以提高图像的准确性和鲁棒性。
总之,反投影重建算法是一种广泛应用于医学图像重建的技术。
虽然该算法存在一些缺点,但通过改进和优化,可以得到高质量的重建图像,为医学诊断和临床应用提供有力支持。
基于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]。
以后的工作中也应对抽样间距进行研究。
医学影像学中的图像重建算法讨论一、介绍医学影像学是一门研究如何获取和解释人体内部结构与功能信息的学科。
在医学影像中,图像重建算法起到至关重要的作用。
图像重建算法是利用被记录下来的原始数据生成高质量图像的数学过程。
本文将讨论医学影像学中常见的图像重建算法,并探讨其优缺点以及未来发展方向。
二、传统图像重建算法1. 过滤回投影(Filtered Backprojection)算法过滤回投影算法是在X射线计算机断层扫描中广泛使用的一种图像重建方法。
它基于Radon变换理论,通过测量和记录X射线通过物体不同角度的衰减程度,然后根据这些数据进行反投影并进行必要的滤波处理,最终生成二维切面图像。
这种传统算法实现简单且计算速度快,但在保留细节方面存在一定的限制。
由于过滤回投影方法假设了直线路径上相应点之间具有线性关系,因此在某些情况下可能会出现伪影或模糊问题。
2. 逆问题正则化(Inverse Problem Regularization)方法逆问题正则化方法是通过建立数学模型来描述图像重建问题,并引入正则化项来约束解的平滑性。
这种方法在医学影像中得到广泛应用,例如在正电子发射断层扫描(Positron Emission Tomography,PET)和单光子发射计算机断层扫描(Single-Photon Emission Computed Tomography,SPECT)等领域。
逆问题正则化方法相比于过滤回投影算法能够更好地保留细节信息,并且对于噪声具有一定的稳健性。
然而,在选择合适的正则化参数和处理噪声方面仍存在挑战。
三、基于深度学习的图像重建算法随着深度学习技术的快速发展,越来越多的研究者开始探索将其应用于医学影像图像重建领域。
下面将介绍几种基于深度学习的图像重建算法。
1. 卷积神经网络(Convolutional Neural Networks,CNN)卷积神经网络是当前深度学习中最常用的模型之一。
它通过多个卷积层和汇聚层来提取特征,并通过反卷积层将特征映射恢复到原始图像大小。
反投影重建算法介绍反投影重建算法是一种常用于医学图像重建的技术,它通过将多个投影图像反向投影到空间中,来重建出原始物体的三维形态。
该算法基于X射线摄影和计算机图形学的理论,被广泛应用于医学成像领域,尤其是计算机断层扫描(CT)。
基本原理反投影重建算法的基本原理是通过测量某物体在不同角度下的投影数据,然后将这些数据反向投影到空间中,得到物体的三维重建图像。
具体步骤如下:1.采集投影数据:使用X射线或其他成像设备,获得物体在不同角度下的投影数据。
这些投影数据可以通过测量X射线的衰减程度来获取。
2.投影数据的逆变换:将投影数据的方向进行反转,得到反投影数据。
这一步是算法中最重要的一步,它将每个测量点的投影数据映射回三维空间。
3.反投影重建:将反投影数据从二维空间转换为三维空间。
这可以通过将反投影数据相加来实现。
在这个过程中,需要注意对数据进行插值,以确保重建图像的精确度。
算法细节反投影重建算法的具体细节有多种实现方式,下面列举几种常见的算法:1. 直接反投影算法直接反投影算法是一种最简单的算法,它将每个投影像素直接投影回三维空间。
具体步骤如下:1.对每个投影像素点,计算它在三维空间中的坐标。
这可以通过测量点在投影平面上的位置和投影矩阵来实现。
2.将每个投影像素点的坐标值累加到三维空间的相应位置上。
这一步是反投影的核心过程,通过将每个像素点的坐标值相加,最终可以得到物体的三维重建图像。
3.可选的后处理步骤:根据需要,可以对重建图像进行滤波、增强或其他处理,以提高图像的质量。
2. 迭代反投影算法迭代反投影算法是一种更复杂但更精确的算法,它通过多次迭代求解,逐渐优化重建图像的质量。
具体步骤如下:1.初始化重建图像:将重建图像的像素初始化为某一固定值。
2.迭代求解:重复以下步骤,直到算法收敛或达到最大迭代次数:–计算投影数据的理论值:根据当前重建图像,计算预测的投影数据。
–计算误差:比较预测的投影数据与实际测量数据之间的差异,得到误差。
滤波反投影法:
滤波反投影法根据附件三所给接收信息,采用先修正、后投影重建图像的做法,
可得到原始图像的吸收率信息。
其原理为:在得到某一角度下的投影函数(一维函
数)后,对此函数做滤波处理,得一修正后的滤波函数,再对修正后的滤波函数做
反投影运算,得待检测介质吸收率在正方形托盘中的每一点的分布密度函数。
图1给出了滤波反投影法重建原始图像的流程图。
图1滤波反投影法流程图
反投影法重建原始图像的步骤:
(1)在对应于投影函数的角度下对投影函数做一维Fourier变换;
(2)对(1)得到的变换结果乘以权重因子;
(3)对(2)加权后得到的结果做一维傅立叶;
(4)对(3)所得函数做直接反投影;
(5)改变投影角度,得到180个不同的投影角度,对每一角度,重复上述步骤(1)~(4)。
R-L(Ram-Lak)滤波函数:
此函数的基本条件是二维图像函数的频率是有界的,显然,此题所得附件五的所有数
据满足此条件。
故频域中的滤波函数可表示为:
,
,其它
其函数图像如图 1.
图1R-L滤波函数图像
连续的R-L卷积函数所得结果为:
离散的R-L卷积函数所得结果为:
,
,为偶数
,为奇数
根据上述滤波原理,在本题中,对附件五中数据的具体滤波过程可用Matlab内置的Ram-Lak 命令实现。
西北工业大学硕士学位论文CT图像滤波反投影重建算法的研究姓名:范慧赟申请学位级别:硕士专业:生物医学工程指导教师:卢健康20070301摘要CT,即计算机断层成像术,是用来获取观测目标断层图像的一门技术,它广泛地应用于医学诊断、射电天文学、电子显微和雷达等许多科学领域,在医学诊断中的应用尤其引人注目.cT是由在多个观测角度获得的有关目标的一系列投影数据,通过图像重建技术来得到目标的断层图像的。变换方法,尤其是滤波(或卷积)反投影类型的算法,是商用cT系统中广泛使用的图像重建方法。本文以此为背景,对反投影重建算法展开研究。由于Fouficr变换理论在计算上本身所固有的复数运算特点,由它产生的重建算法会带来如下的问题:增加数学处理的复杂性、需要较大的存储空间和影响计算速度。而在CT图像重建应用的大量场合,投影数据和目标函数都属实数空间范畴。H删ey变换是与Fourier变换平行的一种正弦类正交变换,具有有效的
快速算法存在。与Four{ef变换相比较,Haftlcy变换以其实数运算、正反变换公式的对称性和经济的存储空间利用等为特点,在对实信号或实数据进行处理时可避免复数运算,从而减少计算量,提高运算效率。本文将Hanlev变换引入到cT图像的滤波反投影重建算法中,这~工作表明:在图像重建领域,用Hanlev变换去替代Fo埘cr变换的作用是可行的。在保持精度不变的前提下,由Haftley变换产生的重建算法较FouIi盯方法可节省大约一半的计算量和存储空问。平行射束常规反投影的多方位同时反投影(MsBP)方法是基于重建图像象素和投影射线之问在不同的投影方位上所存在的几何关系而提出的。本文将此方法应用于反投影重建算法中,将象素定位操作中的计算量减少到传统反投影算法的1/8。本文还在导师提出的基于D‘玎的线性定常系统高精度仿真算法的基础上,
研究、推导了基于Hanley变换的逆变换法的公式并给出其误差分析,通过Matlab程序验明了其可行性和正确性。
关键词:cT图像重建,滤波反投影算法,Hanley变换,MsBP方法,逆变换法AbstractComputeIized’Ibmography(CDisatcc岫ologythatob乜Iinsthct蛐。争锄ofthcobscrw;dobjects.niswidelyused缸m姐y∞iendfic
doma缸ssuch鹤medical
dia印∞is、radioa蜘n哪y、clec咖咀mic∞lo盱锄dndaLltsappli龃ti彻sinmcdical
dia舯osisa他a盯estingcspccially.CEwhichobtainsas谢嚣ofprojccti咖dataof
州cctsconoemed‰seVeral“ewing锄gles’can
gctthctomo{即ms
of
thc嘶ec略
thIoughthc僦岫logyofinlagc嫩衄stnlcti∞.11letra璐fo皿methods,髂pcciatly
thefnte-cd0fcoⅡvoIuti衄back蜊ccti∞a190ritlImsa砖widclylIscdinthc
∞删mrcial(了syst锄s.11listh嚣is懿plof器柚d他se砌髓tl嵋algorithmsofthc
舭lc】删ccti∞岫d盯this他∞眦hb∞kground.
IIccau∞ofthcillhe啪t啪plexn岫bcrol'盯撕帆chamcteIisti岱oftheFburi盯
tmsfo皿aIgofithmin00mpu协ti佃,thc
fec0鹏tnIdi蛐al鲥thmbascd咖地wiU
bring∞mepfoblemssuch镐i玳黔鹤ingthe啪plexityofthemath锄aticpm∞豁in舀谢U∞cdmOmm锄oryspa∞柚dsl删ingdown
the伽pu协ti∞spccd.110啊cvcr'
thep叫ecIi∞data锄dtheobjed如ndio嬲theyaubcl伽gtothc∞al姗证bcra坨ainapplic缸io璐ofCrimage咒cod咖cti锄.H枷eytransfo珊is
asillusoidal仃an甜bnn
similartoF0uricrtfansfo咖.FastHanlcy咖sfb珊啪beachjeVcdbyusingthcs咖cI呲ssilnil缸t0thatinF豳tFbuficr
t瑚sf0蚰.Hartleytfjmsfo珊h舔thr∞
adV柚tagesoV盯Fburi盯仃ansf咖.nisarcalnumb盯缸趾sfbnn卸dh硒the
s锄c
f0研ard柚dinver靶仃柚甜b咖f研mul弱.B龉idcsthat,itn∞dsle豁intcmal
stomge
th觚Fbu断仃孤sfo衄.Hcn。c,usiIlgH捌ey缸msfo埘jn
manyrealsi印aI柚d此al
dalapmcc鲻ingc姐avoidtheu∞ofcomplcxaritllnletic∞勰torcduccthe肿mb懿
0f“tllmcticope删。瞄锄d
ilnprovethccompu诅ti衄ale衔cicn哆This
thesis删izes
thefilt锄cdback州edional鲥thmsofCT
imagcreconshⅥcti∞by憾ingHanlcy
t瑚sfo彻.nshowsthatwecanu∞Hartleyt砌sfbnninsteadOfFbllrier缸卸瞎fom
in
fields0ftheimagc
De咖strudi伽.ncnumbc璐0faritllmeticopcrati∞s锄dthe
intemalstoragenecdcdjnimage地OO璐tnldionaIg耐thlnsb弱edonHanley
胁sfb珊a碥almosthalf勰much弱thatinFouricr的nsfo珊whileretainillgthe
prccisi叫.Based∞the嚣卿metrical陀lati衄shipsbetw∞nimagcpi】【clsandp加jection
mys
Ⅱatdif:fe”nt、,icwing
angles,
a
Multiorienta“衄
simunaneouscomputation
of
backpmjecti∞fortomo伊apllicimagcrcconstm甜咖(MsBP)app帕achispfopDsed.
Byu咖gtllismethodinthe
fihercdbackprojectioⅡalgorithIns
ofCr
imagc
rec0嬲tmction,thenmnber0flo∞ti伽Operatio璐0fpixelsfbrinterpolati彻isfcducedbya白ctorof罐.Anerre∞arcllingonthchiglllypr∽i∞simula胁nal鲥thmofu1
systemb猫cd
onDisc”tecosincTransfo】衄(DcDpmposcdbymysupcrvi∞r,thisthcsisp姊o∞s翘inv哪et啪sfbnnb弱cd∞Hartleytr柚sfom.n百V船thc
particular
f0珈ul码蛆d
thcefror耻alysisofthisalg耐tllm.Finallythefe猫ibnity柚dco盯cd∞鹞m
ValidatcdbyusingMatlab芦Dgrams.
Keywords:CTi咖gerec伽st珊c60n'HlteredBac印roj睇6佣AI印rimm,
HanleyT啪sfom,MSBPMe恤od,InVe鼬T瑚sfo哪
ⅡI西北工业大学学位论文知识产权声明书本人完全了解学校有关保护知识产权的规定,即:研究生在校攻读学位期间论文工作的知识产权单位属于西北工业大学。学校有权保留并向国家有关部门或机构送交论文的复印件和电子版。本人允许论文被查阅和借阅。学校可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。同时本人保证,毕业后结合学位论文研究课题再撰写的文章一律注明作者单位为西北工业大学。保密论文待解密后适用本声明。学位论文作者签名::《塑L指导教师签名
。>况订年弓月多D日聊年岁月歹口日
西北工业大学学位论文原创性声明秉承学校严谨的学风和优良的科学道德,本人郑重声明:所呈交的学位论文,是本人在导师的指导下进行研究工作所取得的成果。尽我所知,除文中已经注明引用的内容和致谢的地方外,本论文不包含任何其他个人或集体己经公开发表或撰写过的研究成果,不包含本人或其他已申请学位或其他用途使用过的成果。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式表明。本人学位论文与资料若有不实,愿意承担一切相关的法律责任。学位论文作者签名网年;月如日