opencv频域相乘

#include <iostream>
#include <fstream>
#include <sstream>
#include <algorithm>
#include <math.h>
#include <dirent.h>
#include <vector>
#include <opencv2/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
int main()
{
    bool LAB = false;
    bool HOG = true;
    bool FIXEDWINDOW = false;
    bool MULTISCALE = true;
    cv::hal;
    int w=10;
    int h=10;
       float cbufi[h*w];
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
            cbufi[i*w+j]= i*j*0.2+1;
        }
    }
    float cbuf[h*w];
 for (size_t i = 0; i < h; i++) {
     for (size_t j = 0; j < w; j++) {
            cbuf[i*w+j]= i*j*0.1+1;
     }
 }
    float cbuf2i[h*w];
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
                cbuf2i[i*w+j]= i*j*0.9+2;
        }
    }
    float cbuf2[h*w];
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
                cbuf2[i*w+j]= i*j*0.5+2;
        }
    }
    cv::Mat r1(h, w, CV_32FC1, cbuf);
      cv::Mat r2(h, w, CV_32FC1, cbufi);
       std::vector<cv::Mat> pa;
       pa.push_back(r1);
        pa.push_back(r2);
         cv::Mat com;
         cv::merge(pa,com);
//     cv::Mat com1=  complexMultiplication(com,r2);
     cv::Mat r3(h, w, CV_32FC1, cbuf2);
       cv::Mat r4(h, w, CV_32FC1, cbuf2i);
       std::vector<cv::Mat> pb;
       pb.push_back(r3);
        pb.push_back(r4);
         cv::Mat com2;
         cv::merge(pb,com2);
// cv::Mat com2=  complexMultiplication(r3,r4);
      KCFTracker tracker(HOG, FIXEDWINDOW, MULTISCALE, LAB);
      cv::Mat caux;
          cv::mulSpectrums(com,com2, caux, 0, true);
          std::vector<cv::Mat> resv;
          cv::split(caux, resv);
         cv::Mat realMat=resv[0];
          cv::Mat imaginaryMat=resv[1];
          for (int i=0;i<caux.rows;i++)
             {
                 for (int j=0;j<caux.cols;j++)
                 {
//                    float num = caux.at<uchar>(i,j);
                   float realf=  realMat.at<uchar>(i,j);
                     float imaginaryf= imaginaryMat.at<uchar>(i,j);
                      printf("%f+%fi ",realf,imaginaryf);
                 }
                  printf("\n");
             }

}