[论文解读] Diffusion, Convection and Erosion on SE(3)/({0} imes SO(2)) and their Application to the Enhancement of Crossing Fibers

Remco Duits, Eric Creusen+2 more|arXiv (Cornell University)|Mar 3, 2011

Advanced Neuroimaging Techniques and Applications参考文献 38被引用 5

一句话总结

该论文在群 $\mathbb{R}^3\rtimes S^2$ 上引入了左不变扩散和汉密尔顿-雅可比（腐蚀）过程，用于建模三维位置与方向，以增强高角分辨率扩散成像（HARDI）中的交叉纤维。通过利用 $SE(3)$ 上的子黎曼几何和卡坦联络，该方法借助通过有限差分和卷积求解的不变偏微分方程，实现了交叉纤维的保留，显著提升了脑部 MRI 中的纤维追踪效果。

ABSTRACT

In this article we study both left-invariant (convection-)diffusions and left-invariant Hamilton-Jacobi equations on the space SE(3)/({0} imes SO(2)) of 3D-positions and orientations naturally embedded in the group SE(3) of 3D-rigid body movements. The general motivation for these (convection-)diffusions and erosions is to obtain crossing-preserving fiber enhancement on probability densities defined on the space of positions and orientations. The linear left-invariant (convection-)diffusions are forward Kolmogorov equations of Brownian motions on SE(3)/({0} imesSO(2)) and can be solved by convolution with the corresponding Green's functions or by a finite difference scheme. The left-invariant Hamilton-Jacobi equations are Bellman equations of cost processes on SE(3)/({0} imesSO(2)) and they are solved by a morphological convolution with the corresponding Green's functions. Furthermore, we consider pseudo-linear scale spaces on the space of positions and orientations that combines dilation and diffusion in a single evolution. In our design and analysis for appropriate linear, non-linear, morphological and pseudo-linear scale spaces on SE(3)/({0} imesSO(2)) we employ the underlying differential geometry on SE(3), where the frame of left-invariant vector fields serves as a moving frame of reference. Furthermore, we will present new and simpler finite difference schemes for our diffusions, which are clear improvements of our previous finite difference schemes. We apply our theory to the enhancement of fibres in magnetic resonance imaging (MRI) techniques for imaging water diffusion processes in brain white matter. We provide experiments of our crossing-preserving evolutions on neural images of a human brain containing crossing fibers.

研究动机与目标

• 开发在三维位置与方向空间（$\mathbb{R}^3\rtimes S^2$）上的不变扩散与形态学过程，以提升医学影像中的纤维增强效果。
• 通过利用李群几何建模方向依赖的信号传播，解决扩散 MRI 中保留交叉纤维的挑战。
• 设计稳定且高效的有限差分格式，用于求解 $\mathbb{R}^3\rtimes S^2$ 上的左不变偏微分方程。
• 通过根植于 $SE(3)$ 的统一几何框架，统一线性、非线性和伪线性尺度空间。
• 将该框架应用于真实 HARDI 和 DTI 数据，展示在人体脑部图像中实现更优的纤维追踪效果。

提出的方法

• 该方法采用 $SE(3) = \mathbb{R}^3\rtimes SO(3)$ 上的左不变向量场作为移动框架，在 $\mathbb{R}^3\rtimes S^2$ 上定义内在微分算子。
• 线性对流-扩散过程被建模为 $\mathbb{R}^3\rtimes S^2$ 上的前向柯尔莫哥洛夫方程，通过 $\mathbb{R}^3\rtimes S^2$ 卷积与格林函数或有限差分求解。
• 非线性腐蚀过程通过左不变汉密尔顿-雅可比方程（贝尔曼方程）建模，利用格林函数的形态学卷积求解。
• 伪线性尺度空间通过群流形上的混合框架，将扩散与腐蚀结合为单一演化过程。
• 利用 $SE(3)$ 上的底层卡坦联络与测地线，推导数值格式的解析近似与稳定性边界。
• 有限差分格式通过角向与空间增量矩阵设计，稳定性通过基于群结构导出的步长边界保证。

实验结果

研究问题

• RQ1 如何在 $\mathbb{R}^3\rtimes S^2$ 上定义左不变扩散与腐蚀过程，以在扩散 MRI 中保留交叉纤维？
• RQ2 卡坦联络与子黎曼几何在 $SE(3)$ 上实现不变偏微分方程以增强纤维方面起到何种作用？
• RQ3 $\mathbb{R}^3\rtimes S^2$ 上的有限差分格式在稳定性与精度方面相较于先前实现有何差异？
• RQ4 伪线性尺度空间能否有效结合扩散与形态学操作，以提升纤维追踪效果？
• RQ5 所提出方法在具有复杂交叉结构的 HARDI 图像中，能在多大程度上增强纤维检测？

主要发现

• 所提出的左不变扩散与腐蚀过程在合成与真实 HARDI 数据中成功保留了交叉纤维，优于标准各向同性方法。
• $\mathbb{R}^3\rtimes S^2$ 上的有限差分格式被证明比先前实现更具稳定性和精度，所导出的稳定性边界确保了收敛性。
• 基于 $SE(3)$ 的几何结构能够精确建模测地线路径与方向空间中的曲率效应，这对纤维追踪至关重要。
• 推导出 $\mathbb{R}^3\rtimes S^2$ 上代价过程格林函数的解析估计，支持高效的数值实现。
• 在人体脑部 MRI 上的实验表明，自适应的左不变演化能以更高保真度检测并增强复杂的纤维交叉，优于传统方法。
• 时间积分核的海森堡近似实现了高效的局部分析与核近似，提升了计算可扩展性。