[论文解读] The Generalized Spherical Radon Transform and Its Application in Texture Analysis

Swanhild Bernstein, Ralf Hielscher+1 more|ArXiv.org|Apr 22, 2005

Medical Image Segmentation Techniques参考文献 17被引用 3

一句话总结

本文引入广义球面Radon变换（GSRT）作为材料科学中纹理分析的几何框架，将$ \mathbb{S}^3 $上的取向概率密度函数$ f $与从衍射实验中获得的角度分布函数$ W $联系起来。通过将GSRT识别为球面平移的球面Radon变换，并利用再生核与径向基函数，作者提出了一种新的方法以从$ W $重建$ f $，阐明了现有方法背后的几何结构，并提高了对实验噪声的鲁棒性。

ABSTRACT

The generalized spherical Radon transform associates the mean values over spherical tori to a function $f$ defined on $\mathbb{S}^3 \subset \mathbb{H}$, where the elements of $\mathbb{S}^3$ are considered as quaternions representing rotations. It is introduced into the analysis of crystallographic preferred orientation and identified with the probability density function corresponding to the angle distribution function $W$. Eventually, this communication suggests a new approach to recover an approximation of $f$ from data sampling $W$. At the same time it provides additional clarification of a recently suggested method applying reproducing kernels and radial basis functions by instructive insight in its involved geometry. The focus is on the correspondence of geometrical and group features but not on the mapping of functions and their spaces.

研究动机与目标

• 为广义球面Radon变换（GSRT）在晶体学纹理分析背景下的几何与群论基础提供理论支撑。
• 阐明$ \mathbb{S}^3 $上的取向概率密度函数$ f $与作为$ f $的GSRT的实验可测量角度分布函数$ W $之间的精确对应关系。
• 利用再生核与径向基函数，为GSRT逆问题提供几何解释。
• 证明对称核的双重Radon变换可恢复其对应的角概率密度函数，从而增强反演稳定性。
• 为未来对比GSRT反演与标准球面Radon变换反演的条件性奠定基础。

提出的方法

• 使用$ \mathbb{S}^3 $中的单位四元数表示$ \mathbb{R}^3 $中的旋转，其中$ T_q(x) = qxq^* $，实现对$ \mathbf{SO}(3) $的双重覆盖。
• 将广义球面Radon变换定义为$ f $在$ \mathbb{S}^3 $中二维球面环面$ T \subset \mathbb{S}^3 $上的平均值，这些环面对应于将固定晶体学方向$ \mathbf{h} \in \mathbb{S}^2 $映射到目标方向$ \mathbf{r} \in \mathbb{S}^2 $的旋转。
• 将GSRT识别为可从旋转样品的平均强度测量中获得的角度分布函数$ W $，从而将实验数据与GSRT直接关联。
• 应用再生核希尔伯特空间理论，结合de la Vallée Poussin核与Abel–Poisson核，构建从采样$ W $中稳定逼近$ f $的方法。
• 推导出核的球面Radon变换的显式表达式，涉及Gegenbauer与勒让德系数，以及Beta与超几何函数。
• 证明对称核$ K $的双重Radon变换可得到其对应的角概率密度函数$ \mathcal{R}\mathcal{R}K $，并以包含$ \cos(\eta/2)^{2\kappa} $与超几何函数的闭式表达式给出。

实验结果

研究问题

• RQ1 GSRT如何通过$ \mathbb{S}^3 $中的四元数与旋转进行几何解释？
• RQ2 取向概率密度函数$ f $与实验可测量的角度分布函数$ W $之间的确切关系是什么？
• RQ3 与标准球面Radon变换相比，GSRT如何改善逆问题的条件性？
• RQ4 对称核的双重Radon变换能否恢复其对应的角概率密度函数？
• RQ5 再生核与径向基函数在稳定从噪声衍射数据中反演GSRT方面起到什么作用？

主要发现

• 广义球面Radon变换被识别为从旋转样品上平均强度测量中获得的角度分布函数$ W $，从而将实验数据与GSRT直接关联。
• GSRT等价于球面平移的球面Radon变换，为基于$ \mathbb{S}^3 $中二维球面环面积分的平均强度数据提供了几何解释。
• 对称核$ K $的双重Radon变换产生对应的角概率密度函数$ \mathcal{R}\mathcal{R}K $，并利用超几何与Beta函数推导出显式表达式。
• 对于Abel–Poisson核，Radon变换为$ \mathcal{R}K = \frac{1 - \kappa^4}{(1 - 2\kappa^2 \cos \eta + \kappa^4)^{3/2}} $，其双重变换涉及椭圆积分。
• 对于de la Vallée Poussin核，双重变换表达为$ \frac{1}{\pi} \frac{2^{1-\kappa}}{\cos \eta_1 \cos \eta_2} \frac{\Gamma(2+\kappa)}{\Gamma(3/2 + \kappa)} \left[ A^{1+\kappa} {}_2F_1(\frac{1}{2}, 1+\kappa, \frac{3}{2}+\kappa, \frac{A}{B}) - B^{1+\kappa} {}_2F_1(\frac{1}{2}, 1+\kappa, \frac{3}{2}+\kappa, \frac{B}{A}) \right] $，其中$ A = 1 + \cos(\eta_1 + \eta_2) $，$ B = 1 + \cos(\eta_1 - \eta_2) $。
• 利用再生核与径向基函数，为从$ W $稳定逼近$ f $提供了具有几何启发性的框架，显著增强了对实验噪声的鲁棒性。