The point spread function (PSF) of a translation invariant imaging system is its impulse response, which cannot always be measured directly. This is the case in high energy X-ray radiography, and it must be estimated from images of calibration objects indirectly related to the impulse response. When the PSF is assumed to have radial symmetry, it can be estimated from an image of an opaque straight edge. We use a non-parametric Bayesian approach, where the prior probability density for the PSF is modeled as a Gaussian Markov random field and radial symmetry is incorporated in a novel way. Markov Chain Monte Carlo posterior estimation is carried out by adapting a recently developed improvement to the Gibbs sampling algorithm, referred to as partially collapsed Gibbs sampling. Moreover, the algorithm we present is proven to satisfy invariance with respect to the target density. Finally, we demonstrate the efficacy of these methods on radiographic data obtained from a high-energy X-ray diagnostic system at the U.S. Department of Energy's Nevada National Security Site.
翻译:翻译变异成像系统的点分布函数(PSF)是其脉冲反应,不能总是直接测量。在高能X射线射电中,情况就是这样,必须从与脉冲反应间接相关的校准对象的图像中估算出来。当假设PSF具有辐射对称性时,可以从不透明的直边缘的图像中估算出来。我们使用非参数的巴耶西亚方法,即PSF先前的概率密度建模为Gaussian Markov随机场和辐射对称制,以新颖的方式纳入。Markov Caincle Monte Carlo possior的估计是通过对Gibbs取样算法(被称为部分崩溃的Gibbs取样法)进行最新开发的改进来进行的。此外,我们提出的算法被证明符合目标密度的变异性。最后,我们展示了这些方法对从美国能源部内华内地国家安全站高能X射线诊断系统获得的辐射数据的有效性。