[论文解读] Fast Finite Shearlet Transform

S. Häuser, Gabriele Steidl|arXiv (Cornell University)|Feb 8, 2012

Mathematical Analysis and Transform Methods参考文献 21被引用 57

一句话总结

本文提出快速有限束剪切变换（FFST），一种基于FFT的平移不变离散束剪切变换，可高效计算图像处理中的束剪切系数。通过使用有限离散束剪切构建Parseval框架，并对傅里叶谱应用保持对称性的修改，该方法即使在偶数尺寸图像上也能确保实值系数，实现对卡通类图像的最优稀疏逼近，满足 $\|f - f_N\|_{L_2}^2 \leq C N^{-2} (\log N)^3$。该实现以GPL许可证免费提供，并通过FFT和适当的索引优化性能。

ABSTRACT

In recent years it has turned out that shearlets have the potential to retrieve directional information so that they became interesting for many applications. Moreover the continuous shearlet transform has the outstanding property to stem from a square integrable group representation. However, to use shearlets and the shearlet transform for reasonable applications one needs fast algorithms to compute a discrete shearlet transform. In this tutorial we present the steps towards an implementation of a fast and finite shearlet transform that is only based on the FFT. Using band-limited shearlets we construct a Parseval frame that provides a simple and straightforward inverse shearlet transform. We provide all proofs and discuss several aspects of our implementation.

研究动机与目标

• 开发一种快速、有限且平移不变的离散束剪切变换，以实现高效的图像处理。
• 通过修改傅里叶谱以保持共轭对称性，确保偶数尺寸图像的实值束剪切系数。
• 提供一种实用的、开源的MATLAB实现（FFST），通过FFT实现性能优化计算。
• 建立一个离散束剪切框架，使其在有限维欧氏空间中构成Parseval框架，从而可通过伴随变换实现简单重构。
• 支持实数与复数束剪切系数，以实现相位分析，并提供清晰的构造指南。

提出的方法

• 该方法通过在全网格上使用均匀、可平方可积的群表示构造离散束剪切，确保数学一致性。
• 采用快速傅里叶变换（FFT）加速所有尺度和方向上束剪切系数的计算。
• 为在偶数尺寸图像上保持傅里叶谱的共轭对称性，将束剪切谱的第一个列和行进行镜像处理，并乘以 $1/\sqrt{2}$，从而在逆FFT后确保系数为实数。
• 离散束剪切在有限维 $L_2$ 空间中构成Parseval框架，可通过伴随变换实现精确重构。
• 通过使用单边频谱支持构造复数束剪切，实现在保持框架性质的同时支持相位分析。
• 实现包含优化的索引、谱计算和模块化文档，便于安装和使用。

实验结果

研究问题

• RQ1 如何利用FFT高效计算有限、平移不变的束剪切变换，同时保持束剪切的数学性质？
• RQ2 为确保偶数尺寸图像的实值束剪切系数，傅里叶谱需要进行何种修改？
• RQ3 所提出的离散束剪切系统如何保持框架性质，并通过伴随变换实现简单重构？
• RQ4 谱对称性对实际图像处理中束剪切系数的数值稳定性和准确性有何影响？
• RQ5 该方法能否同时支持实数与复数束剪切系数，以适应相位分析和几何分离等多样化应用？

主要发现

• FFST在 $N \to \infty$ 时实现对卡通类图像的最优稀疏逼近，误差衰减满足 $\|f - f_N\|_{L_2}^2 \leq C N^{-2} (\log N)^3$，与束剪切系统的理论边界一致。
• 通过镜像并缩放束剪切谱的第一个列和行，该方法确保偶数尺寸图像的束剪切系数为实数，将虚部降至可忽略水平。
• 离散束剪切在有限维 $L_2$ 空间中构成Parseval框架，可通过伴随变换实现精确重构，计算开销极低。
• 性能基准测试表明，基于FFT的实现显著快于非FFT替代方案，且具有高效的索引和内存访问模式。
• 通过单边频谱支持实现复数束剪切，可在保持框架性质和计算效率的同时支持相位分析。
• 开源MATLAB软件包以GPL许可证免费提供，附有清晰的文档和安装说明，便于研究人员和实践者使用。