Mohammad-Djafari A, Demoment G
Appl Opt. 1987 May 1;26(9):1745-54. doi: 10.1364/AO.26.001745.
In diffraction tomography, the generalized Radon theorem relates the Fourier transform (FT) of the diffracted field to the two-dimensional FT of the diffracting object. The relationship stands on algebraic contours, which are semicircles in the case of Born or Rytov first-order linear approximations. But the corresponding data are not sufficient to determine uniquely the solution. We propose a maximum entropy method to reconstruct the object from either the Fourier domain data or directly from the original diffracted field measurements. To do this, we give a new definition for the entropy of an object considered as a function of R(2) to C. To take into account the presence of noise, a chi-squared statistic is added to the entropy measure. The objective function thus obtained is minimized using variational techniques and a conjugate-gradient iterative method. The computational cost and practical implementation of the algorithm are discussed. Some simulated results are given which compare this new method with the classical ones.
在衍射层析成像中,广义拉东定理将衍射场的傅里叶变换(FT)与衍射物体的二维傅里叶变换联系起来。这种关系基于代数轮廓,在玻恩或里托夫一阶线性近似情况下,这些轮廓是半圆。但相应的数据不足以唯一确定解。我们提出一种最大熵方法,可从傅里叶域数据或直接从原始衍射场测量值重建物体。为此,我们给出了一个将物体熵视为从(R(2))到(C)的函数的新定义。为了考虑噪声的存在,在熵度量中添加了卡方统计量。使用变分技术和共轭梯度迭代方法对由此得到的目标函数进行最小化。讨论了该算法的计算成本和实际实现。给出了一些模拟结果,将这种新方法与经典方法进行了比较。