《非标准快速傅里叶变换算法综述》.pdfVIP

第19卷 第 3 期 CT理论与应用研究 Vol.19, No.3 2010年9月（33-46） CT Theory and Applications Sep., 2010 薛会, 张丽, 刘以农. 非标准快速傅里叶变换算法综述[J]. CT 理论与应用研究, 2010, 19(3): 33-46. Xue H, Zhang L, Liu YN. Overview of nonuniform fast Fourier transformation[J]. CT Theory and Applications, 2010, 19(3): 33-46. 非标准快速傅里叶变换算法综述 薛会，张丽 ，刘以农 （清华大学工程物理系，北京 100084） 摘要：非标准快速傅里叶变换算法是一种用于处理非标准采样的快速算法，在信号处理领域有 着广泛的应用。本文回顾了 NUFFT 算法的历史和发展，对通用的算法实现形式做了详细描述， 介绍了该算法的最新进展和研究方向。对NUFFT算法在CT 重建中的应用，如傅里叶重建，迭代 重建，稀疏采样重建等也做了介绍。 关键词：非标准离散傅里叶变换；非标准快速傅里叶变换；CT重建 文章编号：1004-4140（2010）03-0033-14 中图分类号：TP301.6 文献标识码：A 快速傅里叶变换算法（Fast Fourier Transform，FFT）的出现，是数字信号处理 发展史上的一个里程碑[1]。Cooley 和 Tukey 提出的 FFT 算法，将N 点数据的离散傅里叶 变换（Discrete Fourier transform，DFT）的计算量由O( N2 ) 减少为 O(N log N) ，使得 快速计算成为可能。它的出现引起了人们极大的兴趣，得到了迅速发展并广泛应用于科 学技术的各个领域。在当今的世界中，FFT 算法几乎无处不在，深刻影响了人类的生活， 由于它对人类社会做出的重大贡献，因而被评价为二十世纪十大算法之一[2]。 遗憾的是，FFT 算法对时域和频域数据的采样具有严格的要求，其数据分布局限于 均匀分布的标准网格上。而实际应用中，往往出现数据采样分布于非均匀网格上的情形， 即非标准离散傅里叶变换（Nonuniform Discrete Fourier Transform，NDFT）问题， 如雷达信号处理[3] [4] [5] [6] 、计算机断层成像 、磁共振成像 、超声成像 等等，传统的 FFT 算 法对此束手无策。为了克服传统 FFT 算法的这一局限性，利用其快速特性，以避免直接 计算 DFT 而付出的沉重代价，非标准快速傅里叶变换（NonUniform Fast Fourier Transform，NUFFT），或称为非均匀快速傅里叶变换（Non-equispaced Fast Fourier Transform）应运而生。 NUFFT 算法最早可以上溯到天文学领域，1975 年，Brouw[7]提出了 gridding 方法， 用于实现极坐标网格采样的离散傅里叶变换，这一方法被 O’Sullivan 引入到 CT 图像 重建领域[8] [9] [10] ，随即在 MRI 图像重建中也得到了应用 。1993 年，Dutt 等 正式研究 了 NDFT 问题的五种基本格式，给出了快速算法，并分析了其误差特性，这可以看做 是 NUFFT 算法的首次出现。从此，NUFFT 算法得到了广泛地关注，新的研究成果不断 涌现[11－31]。 收稿日期：2010-05-27。 基金项目：国家自然科学基金。 34