傅里叶变换学习参考资料

最新推荐文章于 2024-09-10 20:07:15 发布

m1nt

最新推荐文章于 2024-09-10 20:07:15 发布

阅读量2.6k

点赞数

分类专栏：数字图像处理编程语言文章标签：图像处理算法数学

本文链接：https://blog.csdn.net/Tuanz7/article/details/78261988

版权

这篇博客整理了学习傅里叶变换的参考资料，包括离散傅里叶变换（DFT）教程、图像处理书籍及实例，还介绍了如何使用傅里叶变换进行文本校正和滤波操作，提供了修改后的代码实现。

摘要由CSDN通过智能技术生成

以下是我学习傅里叶变换过程中参考的一些书籍、网站等，记录下来，方便以后查阅：

DFT：Discrete Fourier Transform 离散傅里叶变换

我参考（4. DFT实现算法）修改了（5. 官方文档）的代码，主要删除了5中“剪切和重分布幅度图象限”这一步。替换的方法是在进行傅里叶变换之前用(-1)^(x+y)乘以输入的图像函数，就可以将傅里叶变换的原点F(0,0)移到(M/2,N/2)上。效果与官方代码类似。

#include <iostream>
#include <math.h>
#include <opencv2/core.hpp>
#include <opencv2/highgui.hpp>
#include <opencv2/imgproc.hpp>


using namespace std;
using namespace cv;


int main()
{   
    //读入图像灰度图
    Mat I = imread("Spectrum.tif", IMREAD_GRAYSCALE);       

    //判断图像是否加载成功
    if (I.empty())
    {
        cout << "图像加载失败!" << endl;
        return -1;
    }
    else
        cout << "图像加载成功!" << endl << endl;

    //(-1)^(x+y)乘以输入的图像函数，将傅里叶变换的原点F(0,0)移到(M/2,N/2)上。
    Mat padded(I.rows,I.cols,CV_32F);
    for (int i = 0;i < I.rows;i++)
    {
        for (int j = 0;j < I.cols;j++)
        {
            float temp= I.at<uchar>(i, j);
            padded.at<float>(i, j) = temp*pow(-1, i + j);
        }
    }

    //将planes融合合并成一个多通道数组complexI , 初始时实部就是图像
    Mat planes[] = { Mat_<float>(padded), Mat::zeros(padded.size(),CV_32F) };
    Mat complexI;
    merge(planes, 2, complexI);     

    //进行傅里叶变换
    dft(complexI, complexI);

    //计算幅值，planes[0] = Re(DFT(I),planes[1] = Im(DFT(I)) , 即planes[0]为实部,planes[1]为虚部
    split(complexI, planes);        
    magnitude(planes[0], planes[1], planes[