`

(转载)OpenCV基础篇之图像的DFT频域变换

 
阅读更多

转载自:http://xiahouzuoxin.github.io/notes

程序及分析

/*
 * FileName : fft2.cpp
 * Author   : xiahouzuoxin @163.com
 * Version  : v1.0
 * Date     : Wed 30 Jul 2014 09:42:12 PM CST
 * Brief    : 
 * 
 * Copyright (C) MICL,USTB
 */

#include <iostream>
#include <cv.h>
#include <highgui.h>
#include "imgproc/imgproc.hpp"

using namespace std;
using namespace cv;

int main(int argc, char *argv[])
{
    if (argc < 2) {
        cout<<"Usage:./fft2 [image name]"<<endl;
        return -1;
    }

    // Read as grayscale image
    Mat image = imread(argv[1], CV_LOAD_IMAGE_GRAYSCALE);
    if (!image.data) {
        cout << "Read image error"<<endl;
        return -1;
    }

    Mat padded;
    int m = getOptimalDFTSize(image.rows);  // Return size of 2^x that suite for FFT
    int n = getOptimalDFTSize(image.cols);
    // Padding 0, result is @padded
    copyMakeBorder(image, padded, 0, m-image.rows, 0, n-image.cols, BORDER_CONSTANT, Scalar::all(0));

    // Create planes to storage REAL part and IMAGE part, IMAGE part init are 0
    Mat planes[] = {Mat_<float>(padded), Mat::zeros(padded.size(), CV_32F) };
    Mat complexI;
    merge(planes, 2, complexI);

    dft(complexI, complexI);

    // compute the magnitude and switch to logarithmic scale
    split(complexI, planes);
    magnitude(planes[0], planes[0], planes[1]);
    Mat magI = planes[0];

    // => log(1+sqrt(Re(DFT(I))^2+Im(DFT(I))^2))
    magI += Scalar::all(1);
    log(magI, magI);

    // crop the spectrum
    magI = magI(Rect(0, 0, magI.cols & (-2), magI.rows & (-2)));
    Mat _magI = magI.clone();
    normalize(_magI, _magI, 0, 1, CV_MINMAX);

    // rearrange the quadrants of Fourier image so that the origin is at the image center
    int cx = magI.cols/2;
    int cy = magI.rows/2;

    Mat q0(magI, Rect(0,0,cx,cy));    // Top-Left
    Mat q1(magI, Rect(cx,0,cx,cy));   // Top-Right
    Mat q2(magI, Rect(0,cy,cx,cy));   // Bottom-Left
    Mat q3(magI, Rect(cx,cy,cx,cy));  // Bottom-Right

    // exchange Top-Left and Bottom-Right
    Mat tmp;
    q0.copyTo(tmp);
    q3.copyTo(q0);
    tmp.copyTo(q3);

    // exchange Top-Right and Bottom-Left
    q1.copyTo(tmp);
    q2.copyTo(q1);
    tmp.copyTo(q2);

    normalize(magI, magI, 0, 1, CV_MINMAX);

    imshow("Input image", image);
    imshow("Spectrum magnitude before shift frequency", _magI);
    imshow("Spectrum magnitude after shift frequency", magI);
    waitKey();

    return 0;
}

本程序的作用是:将图像从空间域转换到频率域,并绘制频域图像。

  1. 二维图像的DFT(离散傅里叶变换),

    图像的频域表示的是什么含义呢?又有什么用途呢?图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。图像的边缘部分是突变部分,变化较快,因此反应在频域上是高频分量;图像的噪声大部分情况下是高频部分;图像大部分平缓的灰度变化部分则为低频分量。也就是说,傅立叶变换提供另外一个角度来观察图像,可以将图像从灰度分布转化到频率分布上来观察图像的特征。

    频域在图像处理中,就我所知的用途主要在两方面:图像压缩和图像去噪。关于这两点将在下面给出图片DFT的变换结果后说明。

    有关DFT的更多性质请参考胡广书教授的《数字信号处理》教材。

  2. 请注意读图片的函数与之前有所不同:

    Mat image = imread(argv[1], CV_LOAD_IMAGE_GRAYSCALE);

    CV_LOAD_IMAGE_GRAYSCALE参数表示将原图像转换为灰度图后读入,这是因为后面的DFT变换都是基于二维信号的,而彩色图像是三维信号。当然,也可以对RGB每一通道都进行DFT运算。

  3. DFT算法的原理要求输入信号的长度最好为2^n,这样可以使用快速傅里叶变换算法(FFT算法)进行加速。所以程序中使用

    copyMakeBorder(image, padded, 0, m-image.rows, 0, n-image.cols, BORDER_CONSTANT, Scalar::all(0));

    填充0使横纵长度都为2^n。

    对于一维信号,原DFT直接运算的复杂度是O(N^2),而快速傅里叶变换的复杂度降低到O(Nlog2(N)),假设N为512,足足提高了512/9≈57倍。

  4. 由DFT的性质知,输入为实信号(图像)的时候,频域输出为复数,因此将频域信息分为幅值和相位。频域的幅值高的代表高频分量,幅值低的地方代表低频分量,因此程序中使用

    // => log(1+sqrt(Re(DFT(I))^2+Im(DFT(I))^2))
    magI += Scalar::all(1);
    log(magI, magI);
    
    // crop the spectrum
    magI = magI(Rect(0, 0, magI.cols & (-2), magI.rows & (-2)));
    Mat _magI = magI.clone();
    normalize(_magI, _magI, 0, 1, CV_MINMAX);

    进行log幅值计算及归一化幅值(归一化目的主要是方便将频域通过图像的形式进行显示)。

  5. 关于频域中心平移:将图像的高频分量平移到图像的中心,便于观测。

    int cx = magI.cols/2;
    int cy = magI.rows/2;
    
    Mat q0(magI, Rect(0,0,cx,cy));    // Top-Left
    Mat q1(magI, Rect(cx,0,cx,cy));   // Top-Right
    Mat q2(magI, Rect(0,cy,cx,cy));   // Bottom-Left
    Mat q3(magI, Rect(cx,cy,cx,cy));  // Bottom-Right
    
    // exchange Top-Left and Bottom-Right
    Mat tmp;
    q0.copyTo(tmp);
    q3.copyTo(q0);
    tmp.copyTo(q3);
    
    // exchange Top-Right and Bottom-Left
    q1.copyTo(tmp);
    q2.copyTo(q1);
    tmp.copyTo(q2);

    其原理就是将左上角的频域和右下角的互换,右上角和左下角互换。

    请注意:频域点和空域点的坐标没有一一对应的关系,两者的关系只是上面的DFT公式所见到的。

  6. 本程序因为使用到图像处理相关的函数,所以包含了头文件imgproc/imgproc.hpp,该文件位于OpenCV安装目录的include/opencv2/目录下,在编写Makefile时也要增加相关的头文件路径和库,本程序使用的Makefile如下:

    TARG=fft2
    SRC=fft2.cpp
    LIB=-L/usr/local/lib/
    INC=-I/usr/local/include/opencv/ -I/usr/local/include/opencv2
    CFLAGS=
    
    $(TARG):$(SRC)
      g++ -g -o $@ ${CFLAGS} $(LIB) $(INC) \
          -lopencv_core -lopencv_highgui -lopencv_imgproc \
          $^
    
    .PHONY:clean
    
    clean:
      -rm $(TARG) tags -f

    其中Makefile中的\表示换行(反斜杠后不能再有任何字符,包括空格),如上库增加了-lopencv_imgproc,头文件路径增加了-I/usr/local/include/opencv2

效果

dft

  1. 上图从左到右分别是:原始灰度图(我大爱的杨过啊)、频域平移前的频域图像、频域中心平移后的频域图像。

  2. 提到图像频域变换的用途:压缩和去噪。压缩的原理就是在频域中,大部分频域的值为0(或接近0,可以进行有损压缩,如jpeg图像),只要压缩频域中的少数非0值即可达到图片压缩的目的。去噪则是通过频域的滤波实现,因为噪声大部分情况下体现为高频信号,使用低通滤波器即可滤除高频噪声(当然,也会带来损失,那就是边缘会变得模糊(之前说过,边缘也是高频信号))。

 
 
分享到:
评论

相关推荐

    fourier-transform:opencv+vs2013 二维傅里叶变换

    本篇文章将深入探讨傅里叶变换的基本原理、在OpenCV中的实现方法以及如何在VS2013环境下进行编译和运行。 首先,让我们理解什么是傅里叶变换。傅里叶变换是一种数学工具,它可以将一个函数(在图像处理中通常是二维...

    图片加盲水印,没有用opencv

    值得注意的是,虽然这里没有使用OpenCV库,但OpenCV也提供了傅里叶变换相关的函数,如`cv2.dft()`和`cv2.idft()`,对于更复杂的图像处理任务,OpenCV通常会提供更多的功能和优化。然而,对于基础的傅里叶变换和水印...

    vc++图像处理全系列

    "VC图像处理系列之四_傅立叶变换篇.pdf"会讲解如何使用离散傅里叶变换(DFT)和快速傅里叶变换(FFT)对图像进行频域分析。傅里叶变换可以将图像从空间域转换到频率域,这样可以方便地识别和处理图像中的高频和低频...

    数字图像处理(对灰度图像的各种处理)

    1. 离散傅立叶变换:DFT可以揭示图像的高频和低频成分,常用于图像去噪、频域滤波和频谱分析。在VC6.0中,可以使用库函数如fftw3来实现DFT计算。 2. 离散余弦变换:DCT在图像压缩,如JPEG格式中起到关键作用,因为...

    数字图像处理实验系统的设计及与实现

    例如,从RGB到灰度的转换,使用公式R*0.299 + G*0.587 + B*0.114,傅里叶变换则通过DFT计算图像的频域特性。 4. 图像分析:如图像分割、特征提取等,这些是图像识别和分类的基础。例如,可以使用阈值分割法、区域...

    VC图像处理压缩编码

    本篇将深入探讨图像增强、霍夫曼压缩编码、边缘提取以及傅里叶变换等关键知识点,并结合实际的VC(Visual C++)编程环境,阐述如何实现这些技术。 首先,图像增强是对原始图像进行处理以改善其视觉质量或突出某些...

    c++数字图像处理集合

    12. 图像变换:包括傅立叶变换(`dft`)、离散余弦变换(DCT)等,常用于频域分析和图像压缩。 13. 图像配准:通过寻找两个或多个图像之间的对应关系,实现图像对齐。 14. 图像拼接:将多个图像合并成一个大图像,...

    DigitalProcess.rar_QT 傅里叶_qt 灰度_傅里叶_灰度_边缘检测 qt

    QT结合OpenCV库,可以实现离散傅里叶变换(DFT),从而进行频域分析,例如高频部分可能对应图像的细节,而低频部分则代表大范围的色彩变化。 "DigitalProcess"这个软件很好地展现了QT在数字图像处理中的应用,不仅...

    VC实现的数字水印图像处理

    以DFT为例,这个过程包括将图像转换到频域,嵌入水印信息,然后再转换回空间域。嵌入水印时,可以选择在高频区域或低频区域操作,因为这些区域的变化对人眼感知影响较小。提取水印则需要通过同样的变换过程,从处理...

    数字水印与信息隐藏源代码

    在源代码的实现中,我们可能看到包括图像处理库(如OpenCV)和信号处理算法(如离散傅里叶变换DFT)的运用,以确保嵌入的水印不仅难以被察觉,还能在数据经过各种形式的处理后依然可以检测和提取出来。 信息隐藏...

Global site tag (gtag.js) - Google Analytics