强算子拓扑下基于射线和像素离散化框架的收敛性
Convergence of Ray- and Pixel-Driven Discretization Frameworks in the Strong Operator Topology
摘要 Abstract
断层扫描是医学应用中的核心工具,使医生能够研究患者内部特征。在平行束CT中,Radon变换(二维情况下)常用于模拟测量过程。Radon变换及其伴随变换(称为反投影)的适当离散化至关重要。最常用的离散化方法结合了基于射线的Radon变换和基于像素的反投影,经验报告显示这种方法具有最佳的逼近性能。然而,关于由此产生的逼近误差的严格理解却很少。这些方法涉及三个离散化参数:空间分辨率、探测器分辨率和角度分辨率。通常采用平衡分辨率,即使用相同(或相似)的空间分辨率和探测器分辨率。我们提出了一种新的基于射线和像素离散化的“卷积方法”解释。这使得可以进行结构化的分析,从而解释观察到的行为。特别是,我们在平衡分辨率下证明了基于射线的Radon变换和基于像素的反投影在强算子拓扑下的收敛性,从而从理论上证明了这种方法的合理性。特别地,使用足够高的分辨率可以任意好地逼近Radon变换。
Tomography is a central tool in medical applications, allowing doctors to investigate patients' interior features. The Radon transform (in two dimensions) is commonly used to model the measurement process in parallel-beam CT. Suitable discretization of the Radon transform and its adjoint (called the backprojection) is crucial. The most commonly used discretization approach combines the ray-driven Radon transform with the pixel-driven backprojection, as anecdotal reports describe these as showing the best approximation performance. However, there is little rigorous understanding of induced approximation errors. These methods involve three discretization parameters: the spatial-, detector-, and angular resolutions. Most commonly, balanced resolutions are used, i.e., the same (or similar) spatial- and detector resolutions are employed. We present a novel interpretation of ray- and pixel-driven discretizations as `convolutional methods'. This allows for a structured analysis that can explain observed behavior. In particular, we prove convergence in the strong operator topology of the ray-driven Radon transform and the pixel-driven backprojection under balanced resolutions, thus theoretically justifying this approach. In particular, with high enough resolutions one can approximate the Radon transform arbitrarily well.