function [Result] = spectrum2D(I)
% I is a gray image
% 'Input Image should be gray!' A = rgb2gray(I);A = I; % A 要是一个灰度图像
F = fftshift(fft2(A)); % fft2() 是快速傅里叶变换函数% fftshift() 目的是为了将图像横纵划分成四块，对角线互换。% 经过变换以后的图像 F 中心是低频，往外频率增大。% 以下的操作是为获得可以显示的图像。
% 如果直接显示 F ，得到的很可能是一个全白的图。
% 原因是，F内像素的值一般都很大，远大于255.
% 而matlab的imshow函数黑白对应0-255，所以不做归一化处理的图像，显示都是白色的。B = abs(real(F)); % F 有可能是负值，而且数据是中心对称的，中心为最亮点（可以参考C++例程中的说明）
[m, n] = size(B); % 获得行数、列数
R = reshape(B, 1, m*n); % 将B映射成1行m*n列的矩阵
Temp = mapminmax(R, 0, 255); % 将R的像素值归一化到0-255之间
Result = reshape(Temp, m, n); % imshow(Result); % 建议调用的时候使用：imshow(spectrum2D(灰度图像矩阵))% imshow(spectrum2D( rgb2gray(RGB图像矩阵)))end

mex -setup
mbuild -setup
# mcc -W cpplib:名字  -T link:lib matlab文件.m
mcc -W cpplib:spectrum2D -T link:lib spectrum2D.m

// Spectrum2D.cpp : 定义控制台应用程序的入口点。
// 配置手顺：
// 1. vs2015（工程属性设置为64bit）+opencv320，64 bit windows，配置好（不懂得话去网上搜吧，很容易）。
// 2. 加载matlab生产的dll、头文件、lib，不需要cpp。（不懂得话，百度搜“怎么加载dll”）
// 3. 上步的头文件应该会提示 mclmcrrt.h mclcppclass.h找不到，这是matlab安装目录下的文件，加载到（属性->VC++ ->包含目录）
// 4. 加载matlab库文件 mclmcrrt.lib mclmcr.lib
//      （属性->VC++ ->库目录）路径：MATLAB2013a\extern\lib\win64\microsoft
//
// 作者：路边的十元钱硬币
// 时间：20170630
// 最后的话，做注释可以的话，最好还是英文，因为vs平台，汉字注释有可能会造成编译错误。
//

#include "stdafx.h"
#include <opencv2\opencv.hpp>
#include "spectrum2D.h"/// mwArray 是matlab用的数据类型，不论是矩阵还是数字，当需要作为参数传递到matlab生产的函数中时，要转化成这种数据类型
/// 使用方法见代码
mwArray Mat2mwArray(cv::Mat src)
{CV_Assert(src.type() == CV_8UC1);mwArray dst(src.rows, src.cols, mxUINT8_CLASS); /// 初始化，可以理解成矩阵cv::Mat src_t = src.t();dst.SetData(src_t.data, src.rows*src.cols); /// 赋值return dst;
}cv::Mat mwArry2Mat(mwArray src, int rows, int cols)
{if(src.IsEmpty()) /// 是否为空return cv::Mat();cv::Mat dst = cv::Mat::zeros(rows, cols, CV_64FC1);for(int j(0); j<rows; ++j){double* pdata = dst.ptr<double>(j);for(int i(0); i<cols; ++i){pdata[i] = src(j+1,i+1); /// 元素访问（行号，列号）
        }}return dst;
}int _tmain(int argc, _TCHAR* argv[])
{cv::Mat src = cv::imread("1.JPG", CV_LOAD_IMAGE_GRAYSCALE);cv::imshow("src", src);/// matlab mwArray 初始化，必须要做，不然报错if( !spectrum2DInitialize()){std::cout << "Could not initialize!" << std::endl;return -1; }mwArray mlImg = Mat2mwArray(src);mwArray mlResult ;spectrum2D(1,mlResult, mlImg);/**简单说下mwArray的用法：mwArray a_int_varible(1, 1, mxINT8_CLASS); /// 定义mwArray a_double_varible(1, 1, mxDOUBLE_CLASS);a_int_varible(1,1) = 10; /// 赋值a_double_varible(1,1) = 10.0;int a = a_int_varible(1,1);**/cv::Mat result = mwArry2Mat(mlResult, src.rows, src.cols);cv::imshow("FFT2", result);cv::waitKey(0);return 0;
}