matlab欠采样,科学网—傅立叶欠采样算子的Matlab代码 - 屈小波的博文

可能由于Compressed sensing MRI的流行，很多人问到傅立叶欠采样算子的Matlab代码写法。请到如下网址下载开源的代码(Demo_Undersampling_FourierDomain)：

根据线性的傅立叶欠采样模型：y=U*F*x，

其中x是图像，y是随机采样到的傅立叶系数(MRI中称为k空间)，F是傅立叶变换，U是下采样算子。U的行数小于列数，导致傅立叶欠才算子Phi=U*F也是行数小于列数的。

但有的人在编程中把数学上的定义和程序实现混淆了 ，导致欠采样之后的傅立叶系数y变成了填零的傅立叶系数(没有采样到的傅立叶系数置0,与x等维度)，这就不满足数学定义式。

合理的利用Phi=U*F和它的adjoint operator的性质，可以得到一些快速计算方法。在我们发表的论文中有利用到：

Patch-based nonlocal operator:Xiaobo Qu, Yingkun Hou, Fan Lam, Di Guo, Jianhui Zhong, Zhong Chen. Magnetic resonance image reconstruction from undersampled measurements using a patch-based nonlocal operator, Medical Image Analysis, vol. 18, no. 6, pp. 843-856, 2014.

Download pdf from: Elsevier or here.

Dataset: T2 weighted brain, water phantom (Both of them are in complex)

Patch-based directional wavelets:Xiaobo Qu, Di Guo, Bende Ning,