#include <stdio.h>
#include <stdlib.h>
#include <opencv2\opencv.hpp>
#include <iostream>
#include <math.h>
#include <Eigen\Core>
using namespace cv;
using namespace std;
using namespace Eigen;
int main() {
    Mat image1;
    image1 = imread("Image1.png",0);//圖像大小640 * 480
    if (!image1.data) {
        printf("讀取圖像失敗.\n");
        return -1;
    }
    cout << image1.size() << endl;
    
    /*
    //                                     滑動窗口大小  , block大小     ,block移動步長,   cell大小    ,bins個數(shù)
    HOGDescriptor *hog = new HOGDescriptor(cvSize(64, 48), cvSize(32, 32), cvSize(8, 8), cvSize(16, 16), 9);
    vector<float> descriptors1;
    //           圖像, 特征向量,     window移動步長, padding，用于填充圖片以適應(yīng)大小
    hog->compute(roi1, descriptors1, Size(2, 2), Size(0, 0));
        hog1.compute(image1, descriptors1);
        cout << "descriptors.size = " << descriptors1.size() << endl;//輸出hog特征描繪子的維度
    */
    
    
    Mat roi1(image1, Rect(200,200,300,450));
    Mat roi2(image1, Rect(210, 210, 300, 450));
    Mat roi3(image1, Rect(0, 0, 300, 450));
    //                                     滑動窗口大小  , block大小     ,block移動步長,   cell大小    ,bins個數(shù)
    HOGDescriptor *hog = new HOGDescriptor(cvSize(64, 48), cvSize(16, 16), cvSize(8, 8), cvSize(16, 16), 9);
    vector<float> descriptors1;
    vector<float> descriptors2;
    vector<float> descriptors3;
    
    hog->compute(roi1, descriptors1, Size(64, 28), Size(0, 0));
    
    hog->compute(roi2, descriptors2, Size(64, 28), Size(0, 0));
    hog->compute(roi3, descriptors3, Size(64, 28), Size(0, 0));
    
    VectorXf v1(descriptors1.size());
    VectorXf v2(descriptors1.size());
    VectorXf v3(descriptors1.size());
    for (int i = 0;i < descriptors1.size();i++) {
        v1[i] = descriptors1[i];
        v2[i] = descriptors2[i];
        v3[i] = descriptors3[i];
    }
    float norm1 = (v1-v2).transpose()*(v1-v2);
    float norm2 = (v2-v3).transpose()*(v2-v3);
    cout <<norm1<<":"<< exp(-norm1/250) << endl;
    cout << norm2 << ":"<<exp(-norm2/250) << endl;

    return 0;
    
}