锥束X射线CT投影的蒙特卡罗仿真_精品文档.pdf

锥束X射线CT投影的蒙特卡罗仿真_精品文档.pdf

《锥束X射线CT投影的蒙特卡罗仿真_精品文档.pdf》由会员分享，可在线阅读，更多相关《锥束X射线CT投影的蒙特卡罗仿真_精品文档.pdf（5页珍藏版）》请在冰豆网上搜索。

第42卷第4期2008年4月西安交通大学学报JOURNALOFX17ANJIAOToNGUNIVERSITYV01424Apr2008锥束X射线CT投影的蒙特卡罗仿真闫浩，牟轩沁，罗涛，汤少杰（西安交通大学图像处理与模式识别研究所，710049，西安）摘要：

提出了一种锥束X射线CT（计算机断层扫描）投影的蒙特卡罗计算方法，对X光子逐次抽样并以随机步长跟踪，将探测器区域内的沉积能量转化为投影信号，使用蒙特卡罗计算工具EGSnrc执行后台批处理计算，得到符合成像物理学的仿真投影经过合理加速和参数优化，计算量较通用默认设置降低了2个数量级对生成的锥束投影进行重建，结果显示出预期的射束硬化和散射伪影该方法易于进行二次开发，贴近真实设备，是锥束X射线CT原型机散射校正的基础将仿真投影与真实的物理设备测量数据进行比较，平均误差为118关键词：

蒙特卡罗仿真；锥束投影；计算机断层扫描中图分类号：

TP3914；R8111文献标识码：

A文章编号：

0253987X（2008）04-0414一04MonteCarloSimulationofCone-BeamX-RhyComputerTomographyProjectionYANHao，MOUXuanqin，LUOTao，TANGShaojie（InstituteofImageProcessingandPatternRecognition，XianJiaotongUniversity，Xian710049，China）Abstract：

AMonteCarlosimulationmethodofcone-beamX-raycomputertomography（X-CT）projectionispresentedX-rayphotonsaresampledinseriesandtracedwithrandomdiscretestepsProjectionsignalsareobtainedfromdepositedenergyinthedetectorzoneSimulationsinaccordancewithimagingphysicsaregeneratedbyadoptingtheMonteCarlotoolkitEGSnrc，whichisbatch-programmedforbackgroundcalculationIncontrastwiththecommonlyuseddefaultsettings，accelerationinvestigationandparametersoptimizationleadtoadecreasinginsimulationtimearoundi00timesAsexpected，beam-hardeningandscatteringartifactsarerevealedintheimagesreconstructedfromthesimulatedprojectionThemethodisconvenientforfurtherdevelopingandcanbethebasisofresearchesinthescattercorrectionuponcone-beamX-CTprototypeComparedwiththedatameasuredfromrealmachines，thismethodiSvalidatedwith118errorKeywords：

MonteCarlosimulation；cone-beamprojection；computertomography在锥束X射线CT研究中，开发符合物理实际的投影仿真系统非常必要在投影计算方法中，解析投影计算1-2将X射线等同于数学意义上的直线，无法反映X光子物理层面上的粒子波动性，而蒙特卡罗（MC）投影计算以大数定律为基础，可以有效弥补解析投影计算的不足3本文构建了锥束X射线CT的系统框架，使用Matlab生成各视角输入文件，DOS后台调用蒙特卡罗计算工具EGSnrc完成投影计算与同行的工作相比，本文的研究扩大了系统的通用性啪，同时也解决了软件版权限制51的问题，是开发锥束原型机和散射校正算法的基础1蒙特卡罗方法和EGSnrc工具包蒙特卡罗方法又称随机抽样技巧或统计实验方收稿日期：

2007-0807作者简介：

闫浩（1980一），男，博士生；牟轩沁（联系人），男，教授，博士生导师基金项目：

国家自然科学基金资助项目（60551003，60472004），教育部新世纪优秀人才支持计划资助项目（NCET-05-0828）万方数据第4期囝浩，等：

锥束x射线CT投影的蒙特卡罗仿真415法，以概率统计理论为基础，能够逼真地描述事物的特点及物理实验过程，可以部分代替物理实验6在医学物理方面，主要应用于剂量计算和辐射防护模拟EGSllrc是用于模拟光子和电子在物质中的传输行为的蒙特卡罗计算工具包7，EGSnrc模拟的能量范围从lkeV到几百GeV，其模拟能力已为众多不同研究领域的研究人员接受和证实8EGSnrc具有多平台移植、源代码开放以及免费的优点，本文使用了EGSnrc的C+界面Egs+2锥束X射线CT投影系统的构建21锥束x射线CT投影系统的结构和参数锥束X射线投影系统采用源点和平板探测器固定，模体运动的模式，通过坐标系变换，等同于商用CT机球管旋转的扫描效果，如图1所示物体在Xoy平面内的转动和Z轴方向的平动等效为X射线源点轨迹，可以是圆、螺旋线、标准马鞍线或其他解析曲线描述的运动模式UDV为平板探测器所在平面的二维局部坐标平面；D为探测器所在平面原点；s为X光源发点；C为物体中心图1锥束X射线CT投影系统结构投影系统的组成和参数如下：

探测器屏幕的物质成分为CsI、Gd202S等闪烁物或气体探测器；光谱为单能谱（X射线源能量为10-150keV）或者多能谱（X射线源管电压为70，80，140k9）；系统光子数为1081010，由探测单元数调整；其他参数默认设置见表l，均具有一定的可选范围22模体的构造Egs+界面版本将基本几何构造从层扩展到圆柱、立方体、球体、椎等，新版本V4一r225增加了椭圆柱针对实际原型机上使用的Phan500模体，设计的数字Phan500模体如图2所示模体分A、B、C3个部分，分别用于考查空间分辨率、信噪表1锥束X射线CT成像系统结构参数X焦点近焦远焦尺寸ramX源点轨迹抽样率源到中心距离roan中心到探测器距离mm探测器行71数探测器行17,I单元尺寸mm比和一致性、密度分辨率组织成分确定参照文献Do的方法，结合EGSnrc中物质断面成分软件PEGSData生成所需的pegs输入文件在密度分辨率的模体中，本底物质的CT值为50Hu，CT值为45Hu、55tIu和60Hu的软组织圆柱体对应生物体中的脑质、硬膜外血肿和眼球组织，用来评价算法的抗噪性设计骨物质圆柱是为观察因射束硬化和散射造成的杯形和条纹伪影，以便进一步评价散射校正算法飘蘸（a）立体空间示意图f一12201硬镪硬礅丝域a014（c）一致性模体截面（d）空间分辨率模体截面图2数字Phan500模体（长度单位：

mm）23EGSllFC模拟计算中的加速策略蒙特卡罗计算时间与光子总数、模体复杂度、不确定性（精度要求）等多种因素相关在锥束X射线CT投影的蒙特卡罗计算中，信噪比要求较高，而且需要多个视角的投影，如果使用常规设置进行计算，需要大量时间，本节将研究计算中的加速策略

（1）通过光子截止能夕。

和电子截止能P。

等光电运输参数的设置，关闭当前光谱能量下对计算结2）几椒叭6355瓢Lk8（堕影筮万方数据西安交通大学学报第42卷果没有影响的物理过程根据成像物理学理论，医用X射线在生物组织中主要发生光电吸收和康普顿散射，两个物理过程都产生能量较低的荧光X射线和具有动能的电子其中具有动能的电子转化为生物热能，仅对生物分子具有破坏作用，却不影响光子在探测器上的响应在投影计算中，理论上可以使P。

等于P。

加光谱的最高能量，关闭电子二次运输的所有模拟，将携能电子沉积在当前位置，节省计算时间。

使用Machlett固定阳极球管，管电压为120kV考察模拟时间和均方误差两项指标，先打开所有相关物理过程的开关，夕。

设为最小值，单独改变e。

，实验结果如图3a所示确定P。

后，进一步改变p。

值，实验结果如图3b所示实验结果说明e。

设置弹性较大，对模拟速度起关键作用以模拟结果具有相同确定性为前提，记为光谱最高能量，P髓。

和P即。

分别是夕。

和e。

合理值的上限加速需满足如下关系P饥t+P黜一s一P瞄t

（1）户cLIt+e鲈。

一s一P

（2）P。

egate（3）PcutP私。

（4）下面以图3为例，说明在实际模拟计算中，如何选取满足上述条件的e。

和乡。

由实验得到e弘。

一0609MeV，Pga。

一o025MeV在模拟时，由于户。

对模拟速度影响不大，可选择p。

=o001MeV，eeut设置的自由度较大（例如在图3a的情