- //freak.cpp
- #include "stdafx.h"
- #include <cv.hpp>
- #include <highgui.h>
- #include "utils.h"
- #include <iostream>
- using namespace std;
- void freak(char* path1, char* path2, INFO& info, bool show)
- {
- double t1,t2;
- t1=cvGetTickCount();
- initModule_nonfree();
- Mat img1, img2;
- img1=imread(path1,0);
- img2=imread(path2,0);
- if(img1.data==NULL)
- {
- cout<<"The image can not been loaded: "<<path1<<endl;
- system("pause");
- exit(-1);
- }
- if(img2.data==NULL)
- {
- cout<<"The image can not been loaded: "<<path2<<endl;
- system("pause");
- exit(-1);
- }
- vector<KeyPoint> kpts1_freak, kpts2_freak;
- Mat desc1_freak, desc2_freak;
- Ptr<cv::DescriptorMatcher> matcher_l1 = DescriptorMatcher::create("BruteForce-Hamming"); //二进制汉明距离匹配
- vector<vector<DMatch> > dmatches_freak;
- vector<Point2f> matches_freak, inliers_freak;
- SurfFeatureDetector dfreak(200,4);
- dfreak.detect(img1,kpts1_freak);
- dfreak.detect(img2,kpts2_freak);
- info.n1 = kpts1_freak.size();
- info.n2 = kpts2_freak.size();
- FREAK freak;
- freak.compute(img1,kpts1_freak,desc1_freak);
- freak.compute(img2,kpts2_freak,desc2_freak);
- matcher_l1->knnMatch(desc1_freak,desc2_freak,dmatches_freak,2);
- matches2points_nndr(kpts1_freak,kpts2_freak,dmatches_freak,matches_freak,DRATIO);
- info.m=matches_freak.size()/2;
- compute_inliers_ransac(matches_freak,inliers_freak,MIN_H_ERROR,false);
- info.rm=inliers_freak.size()/2;
- t2=cvGetTickCount();
- info.t=(t2-t1)/1000000.0/cvGetTickFrequency();
- Mat img1_rgb_freak = imread(path1,1);
- Mat img2_rgb_freak = imread(path2,1);
- Mat img_com_freak = Mat(Size(img1.cols*2,img1.rows),CV_8UC3);
- if(show == true)
- {
- draw_inliers(img1_rgb_freak,img2_rgb_freak,img_com_freak,inliers_freak,2);
- imshow("freak",img_com_freak);
- waitKey(0);
- }
- return;
- }
- INFO freak_info;
- freak(path1,path2,freak_info,true);
- showInfo(freak_info);
来源: http://www.phpxs.com/code/1004200/