傅里叶变换频谱的可视化保存

最新推荐文章于 2024-09-29 15:06:28 发布

zcg1942

最新推荐文章于 2024-09-29 15:06:28 发布

阅读量2.7k

点赞数 1

分类专栏：图像处理

本文链接：https://blog.csdn.net/zcg1942/article/details/81033508

版权

图像处理专栏收录该内容

55 篇文章 5 订阅

订阅专栏

在上一篇关于傅里叶变换的博客中，知道了imshow的一个小trick：对normalize得到的0~1之间的浮点数构成的矩阵会进行放大255的操作，得到可视化的灰度图。即便是在python中也是如此操作的，只不过python中的函数封装得更加严密，我们看不到填充、归一化等操作，当把显示出的图像保存到本地时再打开查看，图像的像素依然在0~1之间，不利于我们之后对于频谱图的处理。所以这里就着重研究一下如何将imshow出的图像原样地保存下来。

链接1提到该函数分为范围归一化与数据值归一化。（Normalizes the norm or value range of an array.），由第三个参数alpha决定，1表示函数用来规范值，2表示规范范围，并且和第四个参数分别规定了范围的下限和上限。由此我们可以看到数值归一化和范围归一化的区别：

数值归一化是利用现有的数据通过某种占比方式重新计算每一个数值的占比作为新的值。

根据范数类型norm_type的不同，分母可以是数值之和，对应L1范数；可以是最大值，对应INF范数；可以是数据的平方之和再开根号，对应L2范数。所以数据值归一化也叫范数归一化。

范围的归一化也叫线性变换（线性拉伸或压缩），因为范围变换引入了范围的上限和下限，根据几何知识可以进行计算。

特别地，当待归一化的数组最小值是0，且范围上下限分别是0、1时，范围归一化的结果和范数INF归一化的结果相同。可以按照如下代码验证：

#include<opencv2/opencv.hpp>
#include<iostream>

using namespace std;

int main()
{
	vector<double>a = { 10, 11, 234, 45, 65, 456};
	cv::normalize(a, a, 2, 300, cv::NORM_MINMAX);
	for (int i = 0; i < a.size(); i++)
	{
		cout << a[i] << endl;
	}
	getchar();
	return 0;
}

为了保证频谱图像本身是可视化图像（灰度值分布在0~255），使用normalize函数：

normalize(mag, mag, 0, 255, CV_MINMAX);//范围归一化在0~255，但是显示出来的是全白图像。

于是决定手动对数组乘255.实现方式是用OpenCV的函数：

void multiply(InputArray src1, InputArraysrc2,OutputArray dst, doublescale=1, int dtype=-1);

//Mat gainMat(mag.rows, mag.cols, CV_32F, Scalar::all(255));
	//multiply(gainMat, mag, mag);//multiply要求两个矩阵类型一样

但结果是显示出的频谱图依然是全白图像。注意到浮点数才想起来imshow对浮点数还会进行一次255的乘法操作，所以关键其实在把数组由float类型转换为可视图像的uchar（8为无符号类型）。同样利用OpenCV函数：convertTo函数（src.convertTo(dst, type, scale, shift）