C++实现二维图形的傅里叶变换_C 语言

本文实例讲述了C++实现二维图形的傅里叶变换的方法。有一定的借鉴价值。分享给大家供大家参考。

具体代码如下:

// Fourier.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include "stdio.h"
#include "math.h"
#include <cv.h>
#include <highgui.h>
#include "cxcore.h"

int main(int argc, char* argv[])
{
 IplImage *img;
 IplImage *simg;

 CvMat *mat_R;
 CvMat *mat_I;
 CvMat *mat_SRC;
 CvMat *mat_Row;
 CvMat *mat_Col;
 CvMat *dst;
 CvMat *dst_R;
 CvMat *dst_I;
 CvMat *dst_Row;
 CvMat *dst_Col;
 int i,j,k;
 double temp;
 int height,width,step,channels;

 //载入一幅图片
 img=cvLoadImage("c:\\1.bmp",0);

 //mat_R初始化
 mat_R=cvCreateMat(img->height,img->width,CV_64FC1);
 //mat_I初始化
 mat_I=cvCreateMat(img->height,img->width,CV_64FC1);
 //mat_SRC初始化
 mat_SRC=cvCreateMat(img->height,img->width,CV_64FC2);
 //将图片数据存入mat_R(实部)
 cvConvert(img,mat_R);
 //将虚部初始化为零
 cvZero(mat_I);
 //合并实部、虚部
 cvMerge(mat_R,mat_I,NULL,NULL,mat_SRC);
 //创建双通道double类型数组
 dst=cvCreateMat(img->height,img->width,CV_64FC2);
 dst_R=cvCreateMat(img->height,img->width,CV_64FC1);
 dst_I=cvCreateMat(img->height,img->width,CV_64FC1);
 //为循环变量赋值
 height=img->height;
 width=img->width;
 channels=2;
 step=channels*width;
 //局部变量,值为正一或负一
 int check;
 //将输入数据乘以(-1)^(i+j),用于中心化
 for(j=0;j<height;j++)
 {
 for(i=0;i<width;i++)
 {
  check=(i+j)%2>0?1:-1;
  for(k=0;k<channels;k++)
  {
  mat_SRC->data.db[j*step+i*channels+k]=check*mat_SRC->data.db[j*step+i*channels+k];
  }
 }
 }
 //创建一个mat用于临时存储一行数据
 CvMat mat_Header=cvMat(4,4,CV_64FC2);
 mat_Row=cvCreateMat(1,width,CV_64FC2);
 mat_Col=cvCreateMat(1,height,CV_64FC2);
 //创建一个dst用于临时存储一行数据
 dst_Row=cvCreateMat(1,width,CV_64FC2);
 dst_Col=cvCreateMat(height,1,CV_64FC2);
 //为循环变量赋值
 height=img->height;
 width=img->width;
 channels=2;
 step=channels*width;
 //行的傅里叶变换
 for(j=0;j<height;j++)
 {
 //取得第j行数据
 mat_Row=cvGetRow(mat_SRC,&mat_Header,j);
 //正向傅里叶变换
 cvDFT(mat_Row,dst_Row,CV_DXT_FORWARD);
 //执行循环,赋值到dst
 for(i=0;i<width;i++)
 {
  for(k=0;k<channels;k++)
  {
  dst->data.db[j*step+i*channels+k]=dst_Row->data.db[i*channels+k];
  }
 }
 }

 //列的傅里叶变换
 for(i=0;i<width;i++)
 {
 //取得第i列
 mat_Col=cvGetCol(dst,&mat_Header,i);

 //正向傅里叶变换
 cvDFT(mat_Col,dst_Col,CV_DXT_FORWARD);
 //执行循环,赋值到dst
 for(j=0;j<height;j++)
 {
  for(k=0;k<channels;k++)
  {
  dst->data.db[j*step+i*channels+k]=dst_Col->data.db[j*channels+k];
  }
 }
 }

 //分成两个矩阵
 cvSplit(dst,dst_R,dst_I,NULL,NULL);

 //创建临时指针指向dst_R,dst_I
 double *pR,*pI;
 pR=(double *)dst_R->data.ptr;
 pI=(double *)dst_I->data.ptr;
 //创建一张用于显示的图像
 simg=cvCreateImage(cvGetSize(img),8,1);
 //为循环变量赋值
 height=simg->height;
 width=simg->width;
 channels=1;
 step=channels*width;

 for(j=0;j<height;j++)
 {
 for(i=0;i<width;i++)
 {
  for(k=0;k<channels;k++)
  {
  temp=pR[j*step+i*channels+k]*pR[j*step+i*channels+k]+pI[j*step+i*channels+k]*pI[j*step+i*channels+k];
  temp=temp/(height*width);
  simg->imageData[j*step+i*channels+k]=sqrt(temp);
  }
 }
 }

 cvNamedWindow("Mar",CV_WINDOW_AUTOSIZE);
 cvShowImage("Mar",simg);
 cvWaitKey(0);

 cvReleaseMat(&mat_R);
 cvReleaseMat(&mat_I);
 cvReleaseMat(&mat_SRC);
 //cvReleaseMat(&mat_Row);//这里无法正常释放,有待解决
 //cvReleaseMat(&mat_Col);
 cvReleaseMat(&dst);
 cvReleaseMat(&dst_R);
 cvReleaseMat(&dst_I);
 cvReleaseImage(&img);
 cvReleaseImage(&simg);
 return 0;
}

感兴趣的朋友可以调试运行一下本文实例,程序美中不足的是会有内存泄漏,主要是mat_Row,mat_Col,dst_Row,dst_Col,有能力的读者可以对此进行修改与完善。相信会有新的收获。

以上是小编为您精心准备的的内容,在的博客、问答、公众号、人物、课程等栏目也有的相关内容,欢迎继续使用右上角搜索按钮进行搜索c++
, 傅里叶变换
二维图形
二维傅里叶变换 c语言、c语言实现傅里叶变换、c语言实现图形化界面、c语言实现图形界面、汇编语言实现图形平移,以便于您获取更多的相关知识。

时间: 2024-10-26 08:17:44

C++实现二维图形的傅里叶变换_C 语言的相关文章

C语言 指针与二维数组详解_C 语言

二维数组在概念上是二维的,有行和列,但在内存中所有的数组元素都是连续排列的,它们之间没有"缝隙".以下面的二维数组 a 为例: int a[3][4] = { {0, 1, 2, 3}, {4, 5, 6, 7}, {8, 9, 10, 11} }; 从概念上理解,a 的分布像一个矩阵: 0   1   2   3 4   5   6   7 8   9  10  11 但在内存中,a 的分布是一维线性的,整个数组占用一块连续的内存: C语言中的二维数组是按行排列的,也就是先存放 a[

C++指针数组、数组指针、数组名及二维数组技巧汇总_C 语言

本文较为详细的分析了关于理解C++指针数组,数组指针,数组名,二维数组的一些技巧.是比较重要的概念,相信对于大家的C++程序设计有一定的帮助作用. 一.关于数组名 假设有数组: int a[3] = {1, 2, 3} 1.数组名代表数组第一个元素的地址,注意,不是数组地址(虽然值相等),是数组第一个元素地址,a 等同于 &a[0]; a+1是第二个元素的地址.比第一个元素地址a(或者&a[0])超出了一个整型指针的大小,在这里是4个字节(byte) cout << a <

C语言数组入门之数组的声明与二维数组的模拟_C 语言

语言中指针与数组这两个概念之间的联系是密不可分的,以至于如果不能理解一个概念,就无法彻底理解另一个概念. C语言中的数组值得注意的地方有以下两点: C语言中只有一维数组,而且数组的大小必须在编译期就作为一个常数确定下来.然而,C语言中数组的元素可以是任何类型的对象,当然也可以是另外一个数组.这样,要"仿真"出一个多维数组就不是一件难事. 对于一个数组,我们只能够做两件事:确定该数组的大小,以及获得指向该数组下标为0的元素的指针.其他有关数组的操作,哪怕它们看上去是以数组下标进行运算的,

直观理解C语言中指向一位数组与二维数组的指针_C 语言

一维数组和指针:对于一位数组和指针是很好理解的: 一维数组名: 对于这样的一维数组:int a[5];  a作为数组名就是我们数组的首地址, a是一个地址常量 . 首先说说常量和变量的关系, 对于变量来说, 用箱子去比喻再好不过了, 声明一个变量就声明一个箱子,比如我们开辟出一个苹果类型的箱子, 给这个变量赋值就是把盛放苹果的箱子中放入一个实实在在的苹果, 这就是变量的赋值.  而对于数组来说, 就是一组类型相同的箱子中,一组苹果箱子, 可以放入不同的苹果. 一维数组空间: 变量被声明后, 我们

Android简明开发教程八:引路蜂二维图形绘制实例功能定义

有了前面对Android平台的介绍,基本上可以开始编写Android应用了,这里将以绘制二维图形为例,对Android开发的一般方 法做过介绍,其中涉及到自定义Application类,扩展View,Intent定义,发送消息,Data Binding(Adapter),和基本UI设计 .示例没有使用Android平台自带的二维图形API,而是调用了引路蜂二维图形库,引路蜂二维图形库Graphics 2D API实现了移 动平台上图形引擎,它能够以一种统一的方式处理各种基本图形(Shape),路

matlab中绘制二维图形

  matlab中有时为了看清方程,更加系统的了解方程,常常需要画图表达,那么为了绘制二维图形就需要了解matlab中常用的图形功能,下面介绍一下matlab中的二维图形功能. plot(x) ------------根据向量绘制二维曲线,例如: x=[1,2,3,5,4,3.5,2,0]; plot(x) plot(x,y)------------根据向量x,y绘制二维图形; 在matlab中利用plot(x,y)绘制y=2x^2(0<=x<=10)d的图形 x=linspace(0,10,

跪求-图像增强关于二维离散快速傅里叶变换

问题描述 图像增强关于二维离散快速傅里叶变换 matlab7环境下自己编程实现fft2函数,及ifft2函数,不能调用工具箱函数,求大神们帮忙解决,分享一下 解决方案 http://wenku.baidu.com/link?url=N3wlVVKKz4uoRNzyV8BlnXi4QiV72zCiNE9oMHyis0z4NyZsP4MfC-AXVVCrgnHoWgDNi7Zw8SVqh8wipdmqZNuFUOG65N2K1rck0xXnm3i

二维图形的矩阵变换(一)——基本概念

原文:二维图形的矩阵变换(一)--基本概念 基本的二维变换可包括旋转.缩放.扭曲,和平移四种,                      而这些几何运算则可以转换为一些基本的矩阵运算:      这几个变换都是线性的,但平移运算不是线性的,不能通过2*2矩阵运算完成.若要将点 (2, 1)在 x 方向将其平移 3 个单位,在 y 方向将其平移 4 个单位. 可通过先使用矩阵乘法再使用矩阵加法来完成此操作.      综合这几种基本运算,数学家们将其统一为一个3*3矩阵,存储形式如下:     

二维图形的矩阵变换(三)——在WPF中的应用矩阵变换

原文:二维图形的矩阵变换(三)--在WPF中的应用矩阵变换 UIElement和RenderTransform 首先,我们来看看什么样的对象可以进行变换.在WPF中,用于呈现给用户的对象的基类为Visual类,但是Visual对象并不具有变换功能,具有变换功能的是它的子类UIElement.这个类也是非常底层的类了,几乎我们所有的常用控件都是继承自它,也就是说,基本上所有的UI对象都是可以应用变换的. 然后,我们在再来看看UIElement中变换种类.UIElement支持两种变换:Render