作者：翟天保Steven
版权声明：著作权归作者所有，商业转载请联系作者获得授权，非商业转载请注明出处

函数原型

void magnitude(InputArray x, InputArray y, OutputArray magnitude);

参数说明

InputArray类型的x，表示矢量的浮点型X坐标值，也表示复数中的实部。
InputArray类型的y，表示矢量的浮点型Y坐标值，也表示复数中的虚部。
OutputArray类型的magnitude，输出的幅值，同第一个参数x有着同样的尺寸和类型。

计算公式为：

$magnitude=\sqrt{x^2+y^2}$

测试代码

#include<iostream>
#include<opencv2/opencv.hpp>
#include<ctime>
using namespace std;
using namespace cv;void fftshift(cv::Mat &plane0, cv::Mat &plane1);int main(void)
{Mat test = imread("test.jpg", 0);test.convertTo(test, CV_32FC1);//创建通道，存储dft后的实部与虚部（CV_32F，必须为单通道数）cv::Mat plane[] = { test.clone(), cv::Mat::zeros(test.size() , CV_32FC1) };cv::Mat complexIm;cv::merge(plane, 2, complexIm); // 合并通道 （把两个矩阵合并为一个2通道的Mat类容器）cv::dft(complexIm, complexIm, 0); // 进行傅立叶变换，结果保存在自身// 分离通道（数组分离）cv::split(complexIm, plane);// 以下的操作是频域迁移fftshift(plane[0], plane[1]);// 计算幅值cv::Mat mag,mag_log,mag_nor;cv::magnitude(plane[0], plane[1], mag);// 幅值对数化：log（1+m），便于观察频谱信息mag += Scalar::all(1);cv::log(mag, mag_log);cv::normalize(mag_log, mag_nor, 1,0, NORM_MINMAX);cv::Mat BLUR;// 再次搬移回来进行逆变换fftshift(plane[0], plane[1]);cv::merge(plane, 2, BLUR); // 实部与虚部合并cv::idft(BLUR, BLUR);       // idft结果也为复数BLUR = BLUR / BLUR.rows / BLUR.cols;cv::split(BLUR, plane);//分离通道，主要获取通道imshow("original", test / 255);imshow("result", plane[0] / 255);waitKey(0);system("pause");return 0;
}// fft变换后进行频谱搬移
void fftshift(cv::Mat &plane0, cv::Mat &plane1)
{// 以下的操作是移动图像  (零频移到中心)int cx = plane0.cols / 2;int cy = plane0.rows / 2;cv::Mat part1_r(plane0, cv::Rect(0, 0, cx, cy));  // 元素坐标表示为(cx, cy)cv::Mat part2_r(plane0, cv::Rect(cx, 0, cx, cy));cv::Mat part3_r(plane0, cv::Rect(0, cy, cx, cy));cv::Mat part4_r(plane0, cv::Rect(cx, cy, cx, cy));cv::Mat temp;part1_r.copyTo(temp);  //左上与右下交换位置(实部)part4_r.copyTo(part1_r);temp.copyTo(part4_r);part2_r.copyTo(temp);  //右上与左下交换位置(实部)part3_r.copyTo(part2_r);temp.copyTo(part3_r);cv::Mat part1_i(plane1, cv::Rect(0, 0, cx, cy));  //元素坐标(cx,cy)cv::Mat part2_i(plane1, cv::Rect(cx, 0, cx, cy));cv::Mat part3_i(plane1, cv::Rect(0, cy, cx, cy));cv::Mat part4_i(plane1, cv::Rect(cx, cy, cx, cy));part1_i.copyTo(temp);  //左上与右下交换位置(虚部)part4_i.copyTo(part1_i);temp.copyTo(part4_i);part2_i.copyTo(temp);  //右上与左下交换位置(虚部)part3_i.copyTo(part2_i);temp.copyTo(part3_i);
}

测试效果

从图像中对应的数值来看，同计算公式一致，当需要对两个数组进行平方和开根号操作时，也可以直接调用该函数~

如果文章帮助到你了，可以点个赞让我知道，我会很快乐~加油！

OpenCV-计算二维矢量幅值cv::magnitude相关推荐

opencv 计算二维矢量的幅值—magnitude()函数
计算二维矢量的幅值:magnitude()函数该函数用来计算二维矢量的幅值 void magnitude(InputArray x,InputArray y,OutputArray magnitud ...
使用Python，Opencv进行二维直方图的计算及绘制
使用Python,Opencv进行二维直方图的计算及绘制 1. 效果图 2. 源码参考这篇博客将介绍如何使用Python,Opencv进行二维直方图的计算及绘制(分别用Opencv和Numpy计算 ...
matlab 平面偶极子声场,二维矢量声强的误差分析
摘要:为便于识别和定位平面内噪声源,依据双传声器互谱声强法原理,建立二维矢量声强探头的物理模型,推导了二维声强的计算公式,分别在单极子和偶极子声场条件下,利用该探头测量二维声强及定位误差.结果表明: ...
矢量切片_数据粒度均衡的二维矢量瓦片构建方法
作者信息应申1,2,王子豪1,杜志强3,丁火平4, 李翔翔4 (1. 武汉大学资源与环境科学学院,湖北武汉 430079:2. 自然资源部城市国土资源监测与仿真重点实验室,广东深圳 5 ...
【OpenCV-Python】21.OpenCV的二维直方图
21.OpenCV的二维直方图文章目录前言一.OpenCV中的二维直方图二.Numpy中的二维直方图三.直方图示例 1.使用Numpy函数计算直方图 2.使用OpenCV函数计算直方图四. ...
可完全分离的二维矢量图加密域鲁棒可逆水印算法（一）
摘要随着云制造技术的兴起,加密域可逆水印技术逐渐受到了较多的关注.然而,现有的大部分算法不仅只能应用于图像.视频等冗余性较大的载体,而且难以抵御常见的攻击,甚至只能在单一的域中提取水印.为此,本文针 ...
Opencv之二维码识别---QRCodeDetector
OpenCV在对象检测模块中QRCodeDetector有两个相关API分别实现二维码检测与二维码解析: # 检测二维码 bool cv::QRCodeDetector::detect(InputAr ...
梯度下降计实例计算(二维)
梯度下降计实例计算(二维) 文章目录内容简介示例说明绘制曲线进一步讨论参考资料内容简介梯度下降是在机器学习中重要的计算内容.本文就一个具体的示例,展示如何在梯度下降中. 示例说明设损失 ...
基于OpenCV实现二维码发现与定位
基于OpenCV实现二维码发现与定位在如今流行扫描的年代,应用程序实现二维码扫描检测与识别已经是应用程序的标配.特别是在移动端.如果你的应用程序不能自动发现检测二维码,自动定位二维码你都不好意思跟别 ...

OpenCV-计算二维矢量幅值cv::magnitude

函数原型

参数说明

测试代码

测试效果

OpenCV-计算二维矢量幅值cv::magnitude相关推荐

最新文章

热门文章