C++实现二维图形的傅里叶变换
更新时间:2020年4月25日 17:41 点击:2164
本文实例讲述了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,有能力的读者可以对此进行修改与完善。相信会有新的收获。
相关文章
- vector是表示可以改变大小的数组的序列容器,本文主要介绍了C++STL标准库std::vector的使用详解,文中通过示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2022-03-06
- 这篇文章主要介绍了C++中取余运算的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2021-02-23
- 这篇文章主要介绍了C++ string常用截取字符串方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本文通过例子,讲述了C++调用C#的DLL程序的方法,作出了以下总结,下面就让我们一起来学习吧。...2020-06-25
- 本篇文章主要介绍了C++中四种加密算法之AES源代码,小编觉得挺不错的,现在分享给大家,也给大家做个参考。...2020-04-25
- 整数拆分,指把一个整数分解成若干个整数的和。本文重点给大家介绍C++ 整数拆分方法详解,非常不错,感兴趣的朋友一起学习吧...2020-04-25
- 这篇文章主要介绍了C++中Sort函数详细解析,sort函数是algorithm库下的一个函数,sort函数是不稳定的,即大小相同的元素在排序后相对顺序可能发生改变...2022-08-18
- 这篇文章主要介绍了C++万能库头文件在vs中的安装步骤(图文),文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2021-02-23
- PS排除重叠形状是什么意思?很多朋友都不是很清楚,其实方法很简单的,下面小编就为大家介绍介绍一下,不会的朋友可以参考本文,来看看吧。 步骤:1、在PS中,选择“矩形工具...2016-12-31
- 这篇文章主要介绍了C++ bitset用法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本篇文章小编并不是为大家讲解string类型的用法,而是讲解我个人比较好奇的问题,就是string 类型占几个字节...2020-04-25
- 这篇文章主要为大家详细介绍了C++ Eigen库计算矩阵特征值及特征向量,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2020-04-25
- 这篇文章主要介绍了C++ pair的用法实例详解的相关资料,需要的朋友可以参考下...2020-04-25
- 这篇文章主要介绍了VSCode C++多文件编译的简单使用方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-03-29
- 虽然C++11引入了智能指针的,但是开发人员在与内存的斗争问题上并没有解放,如果我门实用不当仍然有内存泄漏问题,其中智能指针的循环引用缺陷是最大的问题。下面通过实例代码给大家介绍c++中的循环引用,一起看看吧...2020-04-25
- 这篇文章主要给大家介绍了关于C++随机点名生成器的相关资料,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- map容器是C++ STL中的重要一员,删除map容器中value为指定元素的问题是我们经常与遇到的一个问题,下面这篇文章主要给大家介绍了关于利用C++如何删除map容器中指定值的元素的相关资料,需要的朋友可以参考借鉴,下面来一起看看吧。...2020-04-25
- 这篇文章主要介绍了C++ 约瑟夫环问题案例详解,本篇文章通过简要的案例,讲解了该项技术的了解与使用,以下就是详细内容,需要的朋友可以参考下...2021-08-15
- 这篇文章主要介绍了C++中cin的用法详细,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-04-25
- 本篇文章是对C++中的常见编译错误进行了详细的分析介绍,需要的朋友参考下...2020-04-25