提出了一种形状扩散光学层析成像技术的数值方法,它能同时重建组织器官的边界和器官内的光学参数。假设不同类型的组织器官内具有各异的均匀光学参数,光子传输模型便可由一些耦合的Helmholtz方程描述,并使用边界元法作为正向模型的数值求解方法。区域的复杂光滑边界用傅里叶级数展开参数化方法表示。反演问题采用Levenberg-Marquardt优化算法,并用加入不同噪声水平的模拟数据对算法进行验证。数值模拟结果表明,提出的图像重建算法具有较快的收敛速度和全局收敛性,并能从噪声数据中准确恢复相应区域的形状