基于图形处理器的X射线锥束成像模拟算法.pdfVIP

科 研 J；定黑 与 学 术 交 琉 基于图形处理器的X射线锥束成像模拟算法 杨 涛 。赵 星 (首都师范大学检测成像工程中心，北京 1OOO48) 摘 要：针对 X射线锥束成像模拟计算量大、速度慢的问题，提 出了一种基于图形处理器 (GPU)的快速成像模拟算法。该算法沿着每条射线累加所经过体素对投影值的贡献量，实现了对 X射线成像的模拟。在计算射线与体素的交线长时，采用分类处理交线的方法，减少了增量Sid— don算法的动态分支计算。为了提高投影图像质量，该算法还用GPU硬件线性插值采样取代Sid— don算法的邻近插值采样。对三维Shepp—logan模型的测试结果表明，该算法的速度比基于GPU 的增量 Siddon算法平均提高了44 ，而且 图像质量明显提高。最后，用实测数据进一步验证 了算 法的有效性 。 关键词 ：X射线成像 ；模拟 ；图形处理器；分类；采样 中图分类号：TGll5．28；TP391 文献标志码 ：A 文章编号：i000—6656(2011)12—0006—08 A GPU—BasedAlgorithm fortheSimulationofX-RayCone-Beam Imaging YANG Tao，ZHAO Xing (EngineeringResearchCenterofComputedTomography，CapitalNormalUniversity，Beijing100048，China) Abstract：ToacceleratethesimulationofX-raycone-beam imaging，aGPU (graphicsprocessingunit)based algorithm isproposedinthispaper．Thealgorithm generatesX-rayimagebyaccumulatingthecontributionofvoxels alongeachX-ray．IntersectionlengthsofthesevoxelswithX-rayarecalculatedbyclassifyingtheintersectiontypes ， whichreducesthetime-consuming dynamicbranchescompared to thefamousincrementalSiddon algorithm． To improveimagequality，sampledvaluesalongX-rayarecomputedbyGPU hardwaresupported linearinterpolation insteadofnearestinterpolationusedbytheincrementalSiddonalgorithm．Theexperimentoftheprojection calculationofShepp-loganphantom showsthatthesimulationspeedisimprovedby44 averagelyascomparedto theGPU-basedincrementalSiddonalgorithm andabetterimagequalityisachieved．Finally，theproposedalgorithm isvalidatedbytheexperimentwithrealmeasureddata． Keywords：X-rayimaging；Simulation；Graphicsprocessingunit(GPU)；Classification；Sampling X射线成像模拟技术 卜是在对射线源、被检 辅助检测人员确定检测工件所需的X光机 的管 电 测物以及探测器建模 的基础上，利用计算机模拟 X