双通道复数Mat的点乘OpenCV实现

1、理论依据:

复数点乘:a=x1+y1i,b=x2+y2;a*b=(x1*x2-y1*y2)+(x1*y2+x2*y1)i

2、代码:

bool  matSplit(Mat &src,Mat &dst_real, Mat &dst_img)
{
    Mat planes[2];
    split(src, planes);
    dst_real = planes[0];//实部
    dst_img =  planes[1];//虚部

    return true;
}

bool matDotMultiply(Mat src1, Mat src2, Mat dst)
{
  
    //src1_real(x1),src1_imag(y1),src2_real(x2),src2_imag(y2)
    Mat src1_real, src1_imag, src2_real, src2_imag;
    
    Mat channel[2];

    //将输入复数Mat数据拆分成虚部与实部
    matSplit(src1, src1_real, src1_imag);
    matSplit(src2, src2_real, src2_imag);

    /*double x1, x2, x3, x4;*/
    //计算src1.*src2
    for (int i = 0; i(i, j) = src1_real.at(i, j)*src2_real.at(i, j)
                - src1_imag.at(i, j)*src2_imag.at(i, j);

            //x1*y2+x2*y1
                channel[1].at(i, j) = src1_real.at(i, j)*src2_imag.at(i, j)
                    + src2_real.at(i, j)*src1_imag.at(i, j);
        }
    }

    //将计算完成的虚部与实部数据合并到输出dst中
    merge(channel, 2, dst);

    return true;
}

转载请注明:Allen101zhang的CSDN博客 https://blog.csdn.net/Allen101zhang/article/details/88918952

你可能感兴趣的:(图像处理)