OpenCV-复数矩阵点乘ComplexMatrixDotMultiplication

2024-04-23 07:44

本文主要是介绍OpenCV-复数矩阵点乘ComplexMatrixDotMultiplication,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

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

需求说明

       一般用到FFT,就涉及到复数的计算,为了便于调用,我自行封装了一个简单的复数矩阵点乘函数。复数乘法公式:

(a+bi)(c+di)=(ac-bd)+(ad+bc)i

       下面展示代码。

功能函数

// 复数矩阵点乘
cv::Mat ComplexMatrixDotMultiplication(const cv::Mat& A, const cv::Mat& B)
{// 确保A和B都是双通道的矩阵CV_Assert(A.channels() == 2 && B.channels() == 2);// 提取A和B的实部和虚部std::vector<cv::Mat> A_channels, B_channels;cv::split(A, A_channels);cv::split(B, B_channels);cv::Mat A_real = A_channels[0];cv::Mat A_imag = A_channels[1];cv::Mat B_real = B_channels[0];cv::Mat B_imag = B_channels[1];// 计算相乘cv::Mat C_real = A_real.mul(B_real) - A_imag.mul(B_imag);cv::Mat C_imag = A_real.mul(B_imag) + A_imag.mul(B_real);// 合并结果std::vector<cv::Mat> C_channels{ C_real, C_imag };cv::Mat C;cv::merge(C_channels, C);return C;
}

C++&Matlab测试代码

#include <iostream>
#include <time.h>
#include <opencv2/opencv.hpp>using namespace std;
using namespace cv;// 复数矩阵点乘
cv::Mat ComplexMatrixDotMultiplication(const cv::Mat& A, const cv::Mat& B)
{// 确保A和B都是双通道的矩阵CV_Assert(A.channels() == 2 && B.channels() == 2);// 提取A和B的实部和虚部std::vector<cv::Mat> A_channels, B_channels;cv::split(A, A_channels);cv::split(B, B_channels);cv::Mat A_real = A_channels[0];cv::Mat A_imag = A_channels[1];cv::Mat B_real = B_channels[0];cv::Mat B_imag = B_channels[1];// 计算相乘cv::Mat C_real = A_real.mul(B_real) - A_imag.mul(B_imag);cv::Mat C_imag = A_real.mul(B_imag) + A_imag.mul(B_real);// 合并结果std::vector<cv::Mat> C_channels{ C_real, C_imag };cv::Mat C;cv::merge(C_channels, C);return C;
}int main()
{// 定义实际数据的双通道矩阵 A 和 Bcv::Mat A = (cv::Mat_<std::complex<float>>)cv::Mat::zeros(4, 4, CV_32FC2);cv::Mat B = (cv::Mat_<std::complex<float>>)cv::Mat::zeros(4, 4, CV_32FC2);// 自定义 A 和 B 的数据for (int i = 0; i < A.rows; ++i) {for (int j = 0; j < A.cols; ++j) {A.at<std::complex<float>>(i, j) = std::complex<float>(i + 1, j + 1);B.at<std::complex<float>>(i, j) = std::complex<float>(i + 2, j + 2);}}// 调用复数矩阵乘法函数cv::Mat C = GlobalFunc::complexMatrixMultiplication(A, B);// 打印结果std::cout << "Result Matrix C:\n" << C << std::endl;return 0;
}
% 定义 A 和 B 的大小
rows = 4;
cols = 4;% 创建实际数据的双通道矩阵 A 和 B,初始值为零
A = zeros(rows, cols) + 1i * zeros(rows, cols);
B = zeros(rows, cols) + 1i * zeros(rows, cols);% 自定义 A 和 B 的数据
for i = 1:rowsfor j = 1:colsA(i, j) = complex(i, j);B(i, j) = complex(i + 1, j + 1);end
end% 复数矩阵点乘
C = A .* B;% 显示结果
disp('Result Matrix C:');
disp(C);

测试效果

       如果函数有什么可以改进完善的地方,非常欢迎大家指出,一同进步何乐而不为呢~

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

这篇关于OpenCV-复数矩阵点乘ComplexMatrixDotMultiplication的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/928188

相关文章

opencv图像处理之指纹验证的实现

《opencv图像处理之指纹验证的实现》本文主要介绍了opencv图像处理之指纹验证的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学... 目录一、简介二、具体案例实现1. 图像显示函数2. 指纹验证函数3. 主函数4、运行结果三、总结一、

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤

Java中的Opencv简介与开发环境部署方法

《Java中的Opencv简介与开发环境部署方法》OpenCV是一个开源的计算机视觉和图像处理库,提供了丰富的图像处理算法和工具,它支持多种图像处理和计算机视觉算法,可以用于物体识别与跟踪、图像分割与... 目录1.Opencv简介Opencv的应用2.Java使用OpenCV进行图像操作opencv安装j

opencv实现像素统计的示例代码

《opencv实现像素统计的示例代码》本文介绍了OpenCV中统计图像像素信息的常用方法和函数,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 统计像素值的基本信息2. 统计像素值的直方图3. 统计像素值的总和4. 统计非零像素的数量

使用opencv优化图片(画面变清晰)

文章目录 需求影响照片清晰度的因素 实现降噪测试代码 锐化空间锐化Unsharp Masking频率域锐化对比测试 对比度增强常用算法对比测试 需求 对图像进行优化,使其看起来更清晰,同时保持尺寸不变,通常涉及到图像处理技术如锐化、降噪、对比度增强等 影响照片清晰度的因素 影响照片清晰度的因素有很多,主要可以从以下几个方面来分析 1. 拍摄设备 相机传感器:相机传

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

opencv 滚动条

参数介绍:createTrackbar( trackbarname , "hello" , &alpha_slider ,alpha_max ,  on_trackbar )  ;在标签中显示的文字(提示滑动条的用途) TrackbarName创建的滑动条要放置窗体的名字 “hello”滑动条的取值范围从 0 到 alpha_max (最小值只能为 zero).滑动后的值存放在

android-opencv-jni

//------------------start opencv--------------------@Override public void onResume(){ super.onResume(); //通过OpenCV引擎服务加载并初始化OpenCV类库,所谓OpenCV引擎服务即是 //OpenCV_2.4.3.2_Manager_2.4_*.apk程序包,存

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

OpenCV结构分析与形状描述符(11)椭圆拟合函数fitEllipse()的使用

操作系统:ubuntu22.04 OpenCV版本:OpenCV4.9 IDE:Visual Studio Code 编程语言:C++11 算法描述 围绕一组2D点拟合一个椭圆。 该函数计算出一个椭圆,该椭圆在最小二乘意义上最好地拟合一组2D点。它返回一个内切椭圆的旋转矩形。使用了由[90]描述的第一个算法。开发者应该注意,由于数据点靠近包含的 Mat 元素的边界,返回的椭圆/旋转矩形数据