This work considers sequential edge-promoting Bayesian experimental design for (discretized) linear inverse problems, exemplified by X-ray tomography. The process of computing a total variation type reconstruction of the absorption inside the imaged body via lagged diffusivity iteration is interpreted in the Bayesian framework. Assuming a Gaussian additive noise model, this leads to an approximate Gaussian posterior with a covariance structure that contains information on the location of edges in the posterior mean. The next projection geometry is then chosen through A-optimal Bayesian design, which corresponds to minimizing the trace of the updated posterior covariance matrix that accounts for the new projection. Two and three-dimensional numerical examples based on simulated data demonstrate the functionality of the introduced approach.
翻译:这项工作考虑了用于(分解的)线性反问题(以X射线X射线摄影为示例)的连续边缘促进巴耶西亚实验设计。计算通过低低的显性迭代在图像体内重新吸收的全变型过程在巴耶西亚框架中解释。假设高斯添加式噪声模型,这将导致一个近似高斯后部,并有一个包含后方边缘位置信息的共变结构。下一个预测几何随后通过A-最佳巴耶斯设计选择,该设计可最大限度地减少用于计算新预测的更新后部同化矩阵的踪迹。基于模拟数据的二维和三维数字示例显示了所采用方法的功能。