基于动态后验模拟的PET图象重建

正电子发射成像是一种有效的生理功能性成像手段 ,但是由于投影数据中噪声大而给重建带来困难 ,为此提出了利用其他高质量的解剖成像结果构造的分割模板先验来进行完全 Bayesian重建以提高重建效果 .分割模板先验可以表示为包含超验参数的 Markov场形式 ,但是它的非凸性和超验参数的存在使得无法用常规的方法得到最大后验估计 ,为此采用动态后验模拟算法计算后验平均估计 .基于满足条件分布的动态后验模拟法可以同时更新象素的密度和超验参数 ,并且容易得到重建的方差和置信区间 .将这种方法和似然估计、最大后验估计结果进行比较 ,重建的结果无论在空间分辨率和抑制噪声方面都有取得了好的效果 .