光学CT中扩散方程光子传输模型的Monte+Carlo验证.pdfVIP

维普资讯 第33卷第6期 光电工程 Vo1．33，No．6 2006年 6月 Opto—ElectronicEngineering June，2006 文章编号：1003—501X(2006)06—0028—05 光学CT中扩散方程光子传输模型的 MonteCarlo验证 金永明，王加俊 (苏州大学 电子信息学院，江苏 苏州 215021) 摘要：扩散方程模型被广泛应用于光学CT图像的重建中，其精确性直接关系到所成图像的质量。 为了研究扩散方程模型对近红外光子输运过程描述的精确性及其适用范围，本文应用MonteCarlo (Mc)模拟方法，给出了二维组织体光学CT正向问题的数值模拟结果，并与其扩散方程模型的有 限元解进行了分析比较。实验结果表明，基于扩散方程的光子传输模型能较为准确地模拟光子在 强散射介质中的输运过程，适用于该种光学参数条件下的图像重建研究。 关键词：光学CT；正向问题；扩散方程；有限元 中图分类号：TP391．4 文献标识码：A Validationofdiffusionequationphotontransportmodelinoptical tomographyusingM Cmethod JINYong—ming，WANGJia-jun (SchoolofElectronicsandInformationEngineering，SoochowUniversity，Suzhou215021，China) Abstract：Diffusion equ~ion modelhasbeen usedforopticaltomography widely．However,the precisionofthenumericalresultsbasedondiffusionequ~ionisdirectlyrelatedtoimagingquality．So， thevalidiyt scopeofthismodelisneededtostudy．Inthisarticle，numericalresultsofrscatteringtissue basedonMonteCralo(MC)algorithm ispresentedalongwithsomediscussions．Experimentalresults showtheconsistencybetweentheFiniteElementMethod(FEM)solutionofthediffusionequationand theMCsimulation，whichprovesthevalidityofthediuff sionequ~ionmodelofrtheopticaltomography reconstruction． Keywords：Opticaltomography；Forwradproblem；Diffusionequation；Finiteelementmethod 引 言 求解光学CT问题的一个有效的方法是基于正向模型解的迭代图像重建算法 I【J，该方法通过建立给定 模型解与测量值之间的误差范数，并使之极小化而求得满足边界测量值要求的组织体光学参数分布。其具 体实现算法中涉及大量的正向计算问题，故此正向模型的快速有效的求解对于光学CT的实现至关重要。 描述组织体中光子传输行为的数学模型 目前主要有两种：连续模型(基于扩散方程) 和离散模型(基于 MonteCarlo模型) 。扩散方程模型是对辐射传输方程的一阶球谐展开近似，由于它最终表示为相对简单 的椭圆形偏微分方程形式，特别适用诸如有限差分法或有限元法来求解。因而，该模型可模拟任意几何形