2010-09-23 21:44:23 +08:00
|
|
|
/*M///////////////////////////////////////////////////////////////////////////////////////
|
|
|
|
//
|
|
|
|
// IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
|
|
|
|
//
|
|
|
|
// By downloading, copying, installing or using the software you agree to this license.
|
|
|
|
// If you do not agree to this license, do not download, install,
|
|
|
|
// copy or use the software.
|
|
|
|
//
|
|
|
|
//
|
|
|
|
// Intel License Agreement
|
|
|
|
// For Open Source Computer Vision Library
|
|
|
|
//
|
|
|
|
// Copyright (C) 2000, Intel Corporation, all rights reserved.
|
|
|
|
// Third party copyrights are property of their respective owners.
|
|
|
|
//
|
|
|
|
// Redistribution and use in source and binary forms, with or without modification,
|
|
|
|
// are permitted provided that the following conditions are met:
|
|
|
|
//
|
|
|
|
// * Redistribution's of source code must retain the above copyright notice,
|
|
|
|
// this list of conditions and the following disclaimer.
|
|
|
|
//
|
|
|
|
// * Redistribution's in binary form must reproduce the above copyright notice,
|
|
|
|
// this list of conditions and the following disclaimer in the documentation
|
|
|
|
// and/or other materials provided with the distribution.
|
|
|
|
//
|
|
|
|
// * The name of Intel Corporation may not be used to endorse or promote products
|
|
|
|
// derived from this software without specific prior written permission.
|
|
|
|
//
|
|
|
|
// This software is provided by the copyright holders and contributors "as is" and
|
|
|
|
// any express or implied warranties, including, but not limited to, the implied
|
|
|
|
// warranties of merchantability and fitness for a particular purpose are disclaimed.
|
|
|
|
// In no event shall the Intel Corporation or contributors be liable for any direct,
|
|
|
|
// indirect, incidental, special, exemplary, or consequential damages
|
|
|
|
// (including, but not limited to, procurement of substitute goods or services;
|
|
|
|
// loss of use, data, or profits; or business interruption) however caused
|
|
|
|
// and on any theory of liability, whether in contract, strict liability,
|
|
|
|
// or tort (including negligence or otherwise) arising in any way out of
|
|
|
|
// the use of this software, even if advised of the possibility of such damage.
|
|
|
|
//
|
|
|
|
//M*/
|
|
|
|
|
|
|
|
#include "precomp.hpp"
|
2013-05-07 10:45:11 +08:00
|
|
|
#include <limits>
|
2014-08-01 22:11:20 +08:00
|
|
|
#include "opencl_kernels_features2d.hpp"
|
2010-09-23 21:44:23 +08:00
|
|
|
|
2011-05-27 20:15:36 +08:00
|
|
|
#if defined(HAVE_EIGEN) && EIGEN_WORLD_VERSION == 2
|
2010-09-23 21:44:23 +08:00
|
|
|
#include <Eigen/Array>
|
|
|
|
#endif
|
|
|
|
|
|
|
|
namespace cv
|
|
|
|
{
|
|
|
|
|
2014-08-12 03:26:39 +08:00
|
|
|
/////////////////////// ocl functions for BFMatcher ///////////////////////////
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2015-06-20 01:52:14 +08:00
|
|
|
#ifdef HAVE_OPENCL
|
2013-12-27 17:04:02 +08:00
|
|
|
static void ensureSizeIsEnough(int rows, int cols, int type, UMat &m)
|
|
|
|
{
|
|
|
|
if (m.type() == type && m.rows >= rows && m.cols >= cols)
|
|
|
|
m = m(Rect(0, 0, cols, rows));
|
|
|
|
else
|
|
|
|
m.create(rows, cols, type);
|
|
|
|
}
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
static bool ocl_matchSingle(InputArray query, InputArray train,
|
|
|
|
UMat &trainIdx, UMat &distance, int distType)
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
2014-09-23 19:13:46 +08:00
|
|
|
if (query.empty() || train.empty())
|
2013-12-27 17:04:02 +08:00
|
|
|
return false;
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
const int query_rows = query.rows();
|
|
|
|
const int query_cols = query.cols();
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ensureSizeIsEnough(1, query_rows, CV_32S, trainIdx);
|
|
|
|
ensureSizeIsEnough(1, query_rows, CV_32F, distance);
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ocl::Device devDef = ocl::Device::getDefault();
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
UMat uquery = query.getUMat(), utrain = train.getUMat();
|
|
|
|
int kercn = 1;
|
|
|
|
if (devDef.isIntel() &&
|
|
|
|
(0 == (uquery.step % 4)) && (0 == (uquery.cols % 4)) && (0 == (uquery.offset % 4)) &&
|
|
|
|
(0 == (utrain.step % 4)) && (0 == (utrain.cols % 4)) && (0 == (utrain.offset % 4)))
|
|
|
|
kercn = 4;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
int block_size = 16;
|
|
|
|
int max_desc_len = 0;
|
|
|
|
bool is_cpu = devDef.type() == ocl::Device::TYPE_CPU;
|
2013-12-27 17:04:02 +08:00
|
|
|
if (query_cols <= 64)
|
2014-09-23 19:13:46 +08:00
|
|
|
max_desc_len = 64 / kercn;
|
2013-12-27 17:04:02 +08:00
|
|
|
else if (query_cols <= 128 && !is_cpu)
|
2014-09-23 19:13:46 +08:00
|
|
|
max_desc_len = 128 / kercn;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
int depth = query.depth();
|
|
|
|
cv::String opts;
|
|
|
|
opts = cv::format("-D T=%s -D TN=%s -D kercn=%d %s -D DIST_TYPE=%d -D BLOCK_SIZE=%d -D MAX_DESC_LEN=%d",
|
|
|
|
ocl::typeToStr(depth), ocl::typeToStr(CV_MAKETYPE(depth, kercn)), kercn, depth == CV_32F ? "-D T_FLOAT" : "", distType, block_size, max_desc_len);
|
|
|
|
ocl::Kernel k("BruteForceMatch_Match", ocl::features2d::brute_force_match_oclsrc, opts);
|
|
|
|
if(k.empty())
|
2013-12-27 17:04:02 +08:00
|
|
|
return false;
|
|
|
|
|
2015-10-16 22:10:00 +08:00
|
|
|
size_t globalSize[] = {((size_t)query.size().height + block_size - 1) / block_size * block_size, (size_t)block_size};
|
|
|
|
size_t localSize[] = {(size_t)block_size, (size_t)block_size};
|
2014-09-23 19:13:46 +08:00
|
|
|
|
|
|
|
int idx = 0;
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(uquery));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(utrain));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(trainIdx));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(distance));
|
|
|
|
idx = k.set(idx, uquery.rows);
|
|
|
|
idx = k.set(idx, uquery.cols);
|
|
|
|
idx = k.set(idx, utrain.rows);
|
|
|
|
idx = k.set(idx, utrain.cols);
|
|
|
|
idx = k.set(idx, (int)(uquery.step / sizeof(float)));
|
|
|
|
|
|
|
|
return k.run(2, globalSize, localSize, false);
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_matchConvert(const Mat &trainIdx, const Mat &distance, std::vector< std::vector<DMatch> > &matches)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
if( (trainIdx.type() != CV_32SC1) || (distance.type() != CV_32FC1 || distance.cols != trainIdx.cols) )
|
|
|
|
return false;
|
|
|
|
|
|
|
|
const int nQuery = trainIdx.cols;
|
|
|
|
|
|
|
|
matches.clear();
|
|
|
|
matches.reserve(nQuery);
|
|
|
|
|
|
|
|
const int *trainIdx_ptr = trainIdx.ptr<int>();
|
|
|
|
const float *distance_ptr = distance.ptr<float>();
|
|
|
|
for (int queryIdx = 0; queryIdx < nQuery; ++queryIdx, ++trainIdx_ptr, ++distance_ptr)
|
|
|
|
{
|
|
|
|
int trainIndex = *trainIdx_ptr;
|
|
|
|
|
|
|
|
if (trainIndex == -1)
|
|
|
|
continue;
|
|
|
|
|
|
|
|
float dst = *distance_ptr;
|
|
|
|
|
|
|
|
DMatch m(queryIdx, trainIndex, 0, dst);
|
|
|
|
|
|
|
|
std::vector<DMatch> temp;
|
|
|
|
temp.push_back(m);
|
|
|
|
matches.push_back(temp);
|
|
|
|
}
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_matchDownload(const UMat &trainIdx, const UMat &distance, std::vector< std::vector<DMatch> > &matches)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
Mat trainIdxCPU = trainIdx.getMat(ACCESS_READ);
|
|
|
|
Mat distanceCPU = distance.getMat(ACCESS_READ);
|
|
|
|
|
|
|
|
return ocl_matchConvert(trainIdxCPU, distanceCPU, matches);
|
|
|
|
}
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
static bool ocl_knnMatchSingle(InputArray query, InputArray train, UMat &trainIdx,
|
|
|
|
UMat &distance, int distType)
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
2014-09-23 19:13:46 +08:00
|
|
|
if (query.empty() || train.empty())
|
2013-12-27 17:04:02 +08:00
|
|
|
return false;
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
const int query_rows = query.rows();
|
|
|
|
const int query_cols = query.cols();
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ensureSizeIsEnough(1, query_rows, CV_32SC2, trainIdx);
|
|
|
|
ensureSizeIsEnough(1, query_rows, CV_32FC2, distance);
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
trainIdx.setTo(Scalar::all(-1));
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ocl::Device devDef = ocl::Device::getDefault();
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
UMat uquery = query.getUMat(), utrain = train.getUMat();
|
|
|
|
int kercn = 1;
|
|
|
|
if (devDef.isIntel() &&
|
|
|
|
(0 == (uquery.step % 4)) && (0 == (uquery.cols % 4)) && (0 == (uquery.offset % 4)) &&
|
|
|
|
(0 == (utrain.step % 4)) && (0 == (utrain.cols % 4)) && (0 == (utrain.offset % 4)))
|
|
|
|
kercn = 4;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
int block_size = 16;
|
|
|
|
int max_desc_len = 0;
|
|
|
|
bool is_cpu = devDef.type() == ocl::Device::TYPE_CPU;
|
|
|
|
if (query_cols <= 64)
|
|
|
|
max_desc_len = 64 / kercn;
|
|
|
|
else if (query_cols <= 128 && !is_cpu)
|
|
|
|
max_desc_len = 128 / kercn;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
int depth = query.depth();
|
|
|
|
cv::String opts;
|
|
|
|
opts = cv::format("-D T=%s -D TN=%s -D kercn=%d %s -D DIST_TYPE=%d -D BLOCK_SIZE=%d -D MAX_DESC_LEN=%d",
|
|
|
|
ocl::typeToStr(depth), ocl::typeToStr(CV_MAKETYPE(depth, kercn)), kercn, depth == CV_32F ? "-D T_FLOAT" : "", distType, block_size, max_desc_len);
|
|
|
|
ocl::Kernel k("BruteForceMatch_knnMatch", ocl::features2d::brute_force_match_oclsrc, opts);
|
|
|
|
if(k.empty())
|
2013-12-27 17:04:02 +08:00
|
|
|
return false;
|
|
|
|
|
2015-10-16 22:10:00 +08:00
|
|
|
size_t globalSize[] = {((size_t)query_rows + block_size - 1) / block_size * block_size, (size_t)block_size};
|
|
|
|
size_t localSize[] = {(size_t)block_size, (size_t)block_size};
|
2014-09-23 19:13:46 +08:00
|
|
|
|
|
|
|
int idx = 0;
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(uquery));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(utrain));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(trainIdx));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(distance));
|
|
|
|
idx = k.set(idx, uquery.rows);
|
|
|
|
idx = k.set(idx, uquery.cols);
|
|
|
|
idx = k.set(idx, utrain.rows);
|
|
|
|
idx = k.set(idx, utrain.cols);
|
|
|
|
idx = k.set(idx, (int)(uquery.step / sizeof(float)));
|
|
|
|
|
|
|
|
return k.run(2, globalSize, localSize, false);
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_knnMatchConvert(const Mat &trainIdx, const Mat &distance, std::vector< std::vector<DMatch> > &matches, bool compactResult)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
if(trainIdx.type() != CV_32SC2 && trainIdx.type() != CV_32SC1) return false;
|
|
|
|
if(distance.type() != CV_32FC2 && distance.type() != CV_32FC1)return false;
|
|
|
|
if(distance.size() != trainIdx.size()) return false;
|
|
|
|
if(!trainIdx.isContinuous() || !distance.isContinuous()) return false;
|
|
|
|
|
|
|
|
const int nQuery = trainIdx.type() == CV_32SC2 ? trainIdx.cols : trainIdx.rows;
|
|
|
|
const int k = trainIdx.type() == CV_32SC2 ? 2 : trainIdx.cols;
|
|
|
|
|
|
|
|
matches.clear();
|
|
|
|
matches.reserve(nQuery);
|
|
|
|
|
|
|
|
const int *trainIdx_ptr = trainIdx.ptr<int>();
|
|
|
|
const float *distance_ptr = distance.ptr<float>();
|
|
|
|
|
|
|
|
for (int queryIdx = 0; queryIdx < nQuery; ++queryIdx)
|
|
|
|
{
|
|
|
|
matches.push_back(std::vector<DMatch>());
|
|
|
|
std::vector<DMatch> &curMatches = matches.back();
|
|
|
|
curMatches.reserve(k);
|
|
|
|
|
|
|
|
for (int i = 0; i < k; ++i, ++trainIdx_ptr, ++distance_ptr)
|
|
|
|
{
|
|
|
|
int trainIndex = *trainIdx_ptr;
|
|
|
|
|
|
|
|
if (trainIndex != -1)
|
|
|
|
{
|
|
|
|
float dst = *distance_ptr;
|
|
|
|
|
|
|
|
DMatch m(queryIdx, trainIndex, 0, dst);
|
|
|
|
|
|
|
|
curMatches.push_back(m);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
if (compactResult && curMatches.empty())
|
|
|
|
matches.pop_back();
|
|
|
|
}
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_knnMatchDownload(const UMat &trainIdx, const UMat &distance, std::vector< std::vector<DMatch> > &matches, bool compactResult)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
Mat trainIdxCPU = trainIdx.getMat(ACCESS_READ);
|
|
|
|
Mat distanceCPU = distance.getMat(ACCESS_READ);
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
return ocl_knnMatchConvert(trainIdxCPU, distanceCPU, matches, compactResult);
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_radiusMatchSingle(InputArray query, InputArray train,
|
|
|
|
UMat &trainIdx, UMat &distance, UMat &nMatches, float maxDistance, int distType)
|
|
|
|
{
|
|
|
|
if (query.empty() || train.empty())
|
|
|
|
return false;
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
const int query_rows = query.rows();
|
|
|
|
const int train_rows = train.rows();
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ensureSizeIsEnough(1, query_rows, CV_32SC1, nMatches);
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
if (trainIdx.empty())
|
|
|
|
{
|
2014-09-23 19:13:46 +08:00
|
|
|
ensureSizeIsEnough(query_rows, std::max((train_rows / 100), 10), CV_32SC1, trainIdx);
|
|
|
|
ensureSizeIsEnough(query_rows, std::max((train_rows / 100), 10), CV_32FC1, distance);
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
nMatches.setTo(Scalar::all(0));
|
|
|
|
|
2014-09-23 19:13:46 +08:00
|
|
|
ocl::Device devDef = ocl::Device::getDefault();
|
|
|
|
UMat uquery = query.getUMat(), utrain = train.getUMat();
|
|
|
|
int kercn = 1;
|
|
|
|
if (devDef.isIntel() &&
|
|
|
|
(0 == (uquery.step % 4)) && (0 == (uquery.cols % 4)) && (0 == (uquery.offset % 4)) &&
|
|
|
|
(0 == (utrain.step % 4)) && (0 == (utrain.cols % 4)) && (0 == (utrain.offset % 4)))
|
|
|
|
kercn = 4;
|
|
|
|
|
|
|
|
int block_size = 16;
|
|
|
|
int depth = query.depth();
|
|
|
|
cv::String opts;
|
|
|
|
opts = cv::format("-D T=%s -D TN=%s -D kercn=%d %s -D DIST_TYPE=%d -D BLOCK_SIZE=%d",
|
|
|
|
ocl::typeToStr(depth), ocl::typeToStr(CV_MAKETYPE(depth, kercn)), kercn, depth == CV_32F ? "-D T_FLOAT" : "", distType, block_size);
|
|
|
|
ocl::Kernel k("BruteForceMatch_RadiusMatch", ocl::features2d::brute_force_match_oclsrc, opts);
|
|
|
|
if (k.empty())
|
|
|
|
return false;
|
|
|
|
|
2015-10-16 22:10:00 +08:00
|
|
|
size_t globalSize[] = {((size_t)train_rows + block_size - 1) / block_size * block_size, ((size_t)query_rows + block_size - 1) / block_size * block_size};
|
|
|
|
size_t localSize[] = {(size_t)block_size, (size_t)block_size};
|
2014-09-23 19:13:46 +08:00
|
|
|
|
|
|
|
int idx = 0;
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(uquery));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrReadOnly(utrain));
|
|
|
|
idx = k.set(idx, maxDistance);
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(trainIdx));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(distance));
|
|
|
|
idx = k.set(idx, ocl::KernelArg::PtrWriteOnly(nMatches));
|
|
|
|
idx = k.set(idx, uquery.rows);
|
|
|
|
idx = k.set(idx, uquery.cols);
|
|
|
|
idx = k.set(idx, utrain.rows);
|
|
|
|
idx = k.set(idx, utrain.cols);
|
|
|
|
idx = k.set(idx, trainIdx.cols);
|
|
|
|
idx = k.set(idx, (int)(uquery.step / sizeof(float)));
|
|
|
|
idx = k.set(idx, (int)(trainIdx.step / sizeof(int)));
|
|
|
|
|
|
|
|
return k.run(2, globalSize, localSize, false);
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_radiusMatchConvert(const Mat &trainIdx, const Mat &distance, const Mat &_nMatches,
|
|
|
|
std::vector< std::vector<DMatch> > &matches, bool compactResult)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty() || _nMatches.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
if( (trainIdx.type() != CV_32SC1) ||
|
|
|
|
(distance.type() != CV_32FC1 || distance.size() != trainIdx.size()) ||
|
|
|
|
(_nMatches.type() != CV_32SC1 || _nMatches.cols != trainIdx.rows) )
|
|
|
|
return false;
|
|
|
|
|
|
|
|
const int nQuery = trainIdx.rows;
|
|
|
|
|
|
|
|
matches.clear();
|
|
|
|
matches.reserve(nQuery);
|
|
|
|
|
|
|
|
const int *nMatches_ptr = _nMatches.ptr<int>();
|
|
|
|
|
|
|
|
for (int queryIdx = 0; queryIdx < nQuery; ++queryIdx)
|
|
|
|
{
|
|
|
|
const int *trainIdx_ptr = trainIdx.ptr<int>(queryIdx);
|
|
|
|
const float *distance_ptr = distance.ptr<float>(queryIdx);
|
|
|
|
|
|
|
|
const int nMatches = std::min(nMatches_ptr[queryIdx], trainIdx.cols);
|
|
|
|
|
|
|
|
if (nMatches == 0)
|
|
|
|
{
|
|
|
|
if (!compactResult)
|
|
|
|
matches.push_back(std::vector<DMatch>());
|
|
|
|
continue;
|
|
|
|
}
|
|
|
|
|
|
|
|
matches.push_back(std::vector<DMatch>(nMatches));
|
|
|
|
std::vector<DMatch> &curMatches = matches.back();
|
|
|
|
|
|
|
|
for (int i = 0; i < nMatches; ++i, ++trainIdx_ptr, ++distance_ptr)
|
|
|
|
{
|
|
|
|
int trainIndex = *trainIdx_ptr;
|
|
|
|
|
|
|
|
float dst = *distance_ptr;
|
|
|
|
|
|
|
|
DMatch m(queryIdx, trainIndex, 0, dst);
|
|
|
|
|
|
|
|
curMatches[i] = m;
|
|
|
|
}
|
|
|
|
|
|
|
|
std::sort(curMatches.begin(), curMatches.end());
|
|
|
|
}
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
|
|
|
static bool ocl_radiusMatchDownload(const UMat &trainIdx, const UMat &distance, const UMat &nMatches,
|
|
|
|
std::vector< std::vector<DMatch> > &matches, bool compactResult)
|
|
|
|
{
|
|
|
|
if (trainIdx.empty() || distance.empty() || nMatches.empty())
|
|
|
|
return false;
|
|
|
|
|
|
|
|
Mat trainIdxCPU = trainIdx.getMat(ACCESS_READ);
|
|
|
|
Mat distanceCPU = distance.getMat(ACCESS_READ);
|
|
|
|
Mat nMatchesCPU = nMatches.getMat(ACCESS_READ);
|
|
|
|
|
|
|
|
return ocl_radiusMatchConvert(trainIdxCPU, distanceCPU, nMatchesCPU, matches, compactResult);
|
|
|
|
}
|
2015-06-20 01:52:14 +08:00
|
|
|
#endif
|
2013-12-27 17:04:02 +08:00
|
|
|
|
2010-09-23 21:44:23 +08:00
|
|
|
/****************************************************************************************\
|
|
|
|
* DescriptorMatcher *
|
|
|
|
\****************************************************************************************/
|
2010-11-23 02:27:08 +08:00
|
|
|
DescriptorMatcher::DescriptorCollection::DescriptorCollection()
|
|
|
|
{}
|
|
|
|
|
|
|
|
DescriptorMatcher::DescriptorCollection::DescriptorCollection( const DescriptorCollection& collection )
|
|
|
|
{
|
|
|
|
mergedDescriptors = collection.mergedDescriptors.clone();
|
2012-12-21 23:58:51 +08:00
|
|
|
std::copy( collection.startIdxs.begin(), collection.startIdxs.begin(), startIdxs.begin() );
|
2010-11-23 02:27:08 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
DescriptorMatcher::DescriptorCollection::~DescriptorCollection()
|
|
|
|
{}
|
|
|
|
|
2013-02-25 00:14:01 +08:00
|
|
|
void DescriptorMatcher::DescriptorCollection::set( const std::vector<Mat>& descriptors )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
clear();
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
size_t imageCount = descriptors.size();
|
2010-10-29 16:44:42 +08:00
|
|
|
CV_Assert( imageCount > 0 );
|
|
|
|
|
|
|
|
startIdxs.resize( imageCount );
|
|
|
|
|
|
|
|
int dim = -1;
|
|
|
|
int type = -1;
|
|
|
|
startIdxs[0] = 0;
|
|
|
|
for( size_t i = 1; i < imageCount; i++ )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
int s = 0;
|
2010-11-23 02:27:08 +08:00
|
|
|
if( !descriptors[i-1].empty() )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
dim = descriptors[i-1].cols;
|
|
|
|
type = descriptors[i-1].type();
|
|
|
|
s = descriptors[i-1].rows;
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
startIdxs[i] = startIdxs[i-1] + s;
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
2010-10-29 16:44:42 +08:00
|
|
|
if( imageCount == 1 )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
if( descriptors[0].empty() ) return;
|
2010-10-29 16:44:42 +08:00
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
dim = descriptors[0].cols;
|
|
|
|
type = descriptors[0].type();
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
2013-04-10 19:54:14 +08:00
|
|
|
CV_Assert( dim > 0 );
|
2010-10-29 16:44:42 +08:00
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
int count = startIdxs[imageCount-1] + descriptors[imageCount-1].rows;
|
2010-10-29 16:44:42 +08:00
|
|
|
|
|
|
|
if( count > 0 )
|
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
mergedDescriptors.create( count, dim, type );
|
2010-10-29 16:44:42 +08:00
|
|
|
for( size_t i = 0; i < imageCount; i++ )
|
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
if( !descriptors[i].empty() )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
CV_Assert( descriptors[i].cols == dim && descriptors[i].type() == type );
|
|
|
|
Mat m = mergedDescriptors.rowRange( startIdxs[i], startIdxs[i] + descriptors[i].rows );
|
|
|
|
descriptors[i].copyTo(m);
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
}
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
void DescriptorMatcher::DescriptorCollection::clear()
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
startIdxs.clear();
|
2010-11-23 02:27:08 +08:00
|
|
|
mergedDescriptors.release();
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
const Mat DescriptorMatcher::DescriptorCollection::getDescriptor( int imgIdx, int localDescIdx ) const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
CV_Assert( imgIdx < (int)startIdxs.size() );
|
|
|
|
int globalIdx = startIdxs[imgIdx] + localDescIdx;
|
|
|
|
CV_Assert( globalIdx < (int)size() );
|
|
|
|
|
|
|
|
return getDescriptor( globalIdx );
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
const Mat& DescriptorMatcher::DescriptorCollection::getDescriptors() const
|
|
|
|
{
|
|
|
|
return mergedDescriptors;
|
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
const Mat DescriptorMatcher::DescriptorCollection::getDescriptor( int globalDescIdx ) const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
CV_Assert( globalDescIdx < size() );
|
2010-11-23 02:27:08 +08:00
|
|
|
return mergedDescriptors.row( globalDescIdx );
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
void DescriptorMatcher::DescriptorCollection::getLocalIdx( int globalDescIdx, int& imgIdx, int& localDescIdx ) const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2011-03-13 08:31:17 +08:00
|
|
|
CV_Assert( (globalDescIdx>=0) && (globalDescIdx < size()) );
|
|
|
|
std::vector<int>::const_iterator img_it = std::upper_bound(startIdxs.begin(), startIdxs.end(), globalDescIdx);
|
|
|
|
--img_it;
|
2011-05-17 23:44:01 +08:00
|
|
|
imgIdx = (int)(img_it - startIdxs.begin());
|
2011-03-13 08:31:17 +08:00
|
|
|
localDescIdx = globalDescIdx - (*img_it);
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
int DescriptorMatcher::DescriptorCollection::size() const
|
|
|
|
{
|
|
|
|
return mergedDescriptors.rows;
|
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
/*
|
|
|
|
* DescriptorMatcher
|
|
|
|
*/
|
2013-02-25 00:14:01 +08:00
|
|
|
static void convertMatches( const std::vector<std::vector<DMatch> >& knnMatches, std::vector<DMatch>& matches )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
|
|
|
matches.clear();
|
|
|
|
matches.reserve( knnMatches.size() );
|
|
|
|
for( size_t i = 0; i < knnMatches.size(); i++ )
|
|
|
|
{
|
|
|
|
CV_Assert( knnMatches[i].size() <= 1 );
|
|
|
|
if( !knnMatches[i].empty() )
|
|
|
|
matches.push_back( knnMatches[i][0] );
|
|
|
|
}
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
DescriptorMatcher::~DescriptorMatcher()
|
|
|
|
{}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::add( InputArrayOfArrays _descriptors )
|
2010-11-23 02:27:08 +08:00
|
|
|
{
|
2016-01-23 23:57:46 +08:00
|
|
|
if( _descriptors.isUMatVector() )
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
std::vector<UMat> descriptors;
|
2016-01-23 23:57:46 +08:00
|
|
|
_descriptors.getUMatVector( descriptors );
|
2013-12-27 17:04:02 +08:00
|
|
|
utrainDescCollection.insert( utrainDescCollection.end(), descriptors.begin(), descriptors.end() );
|
|
|
|
}
|
2016-01-23 23:57:46 +08:00
|
|
|
else if( _descriptors.isUMat() )
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
std::vector<UMat> descriptors = std::vector<UMat>(1, _descriptors.getUMat());
|
|
|
|
utrainDescCollection.insert( utrainDescCollection.end(), descriptors.begin(), descriptors.end() );
|
|
|
|
}
|
2016-01-23 23:57:46 +08:00
|
|
|
else if( _descriptors.isMatVector() )
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
std::vector<Mat> descriptors;
|
|
|
|
_descriptors.getMatVector(descriptors);
|
|
|
|
trainDescCollection.insert( trainDescCollection.end(), descriptors.begin(), descriptors.end() );
|
|
|
|
}
|
2016-01-23 23:57:46 +08:00
|
|
|
else if( _descriptors.isMat() )
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
std::vector<Mat> descriptors = std::vector<Mat>(1, _descriptors.getMat());
|
|
|
|
trainDescCollection.insert( trainDescCollection.end(), descriptors.begin(), descriptors.end() );
|
|
|
|
}
|
|
|
|
else
|
2016-01-23 23:57:46 +08:00
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
CV_Assert( _descriptors.isUMat() || _descriptors.isUMatVector() || _descriptors.isMat() || _descriptors.isMatVector() );
|
2016-01-23 23:57:46 +08:00
|
|
|
}
|
2010-11-23 02:27:08 +08:00
|
|
|
}
|
|
|
|
|
2013-02-25 00:14:01 +08:00
|
|
|
const std::vector<Mat>& DescriptorMatcher::getTrainDescriptors() const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
return trainDescCollection;
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
void DescriptorMatcher::clear()
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
utrainDescCollection.clear();
|
2010-10-29 16:44:42 +08:00
|
|
|
trainDescCollection.clear();
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
bool DescriptorMatcher::empty() const
|
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
return trainDescCollection.empty() && utrainDescCollection.empty();
|
2010-11-23 02:27:08 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
void DescriptorMatcher::train()
|
|
|
|
{}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::match( InputArray queryDescriptors, InputArray trainDescriptors,
|
|
|
|
std::vector<DMatch>& matches, InputArray mask ) const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
Ptr<DescriptorMatcher> tempMatcher = clone(true);
|
2013-12-27 17:04:02 +08:00
|
|
|
tempMatcher->add(trainDescriptors);
|
|
|
|
tempMatcher->match( queryDescriptors, matches, std::vector<Mat>(1, mask.getMat()) );
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::knnMatch( InputArray queryDescriptors, InputArray trainDescriptors,
|
|
|
|
std::vector<std::vector<DMatch> >& matches, int knn,
|
|
|
|
InputArray mask, bool compactResult ) const
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
Ptr<DescriptorMatcher> tempMatcher = clone(true);
|
2013-12-27 17:04:02 +08:00
|
|
|
tempMatcher->add(trainDescriptors);
|
|
|
|
tempMatcher->knnMatch( queryDescriptors, matches, knn, std::vector<Mat>(1, mask.getMat()), compactResult );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::radiusMatch( InputArray queryDescriptors, InputArray trainDescriptors,
|
|
|
|
std::vector<std::vector<DMatch> >& matches, float maxDistance, InputArray mask,
|
|
|
|
bool compactResult ) const
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
Ptr<DescriptorMatcher> tempMatcher = clone(true);
|
2013-12-27 17:04:02 +08:00
|
|
|
tempMatcher->add(trainDescriptors);
|
|
|
|
tempMatcher->radiusMatch( queryDescriptors, matches, maxDistance, std::vector<Mat>(1, mask.getMat()), compactResult );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
2014-02-17 07:00:20 +08:00
|
|
|
void DescriptorMatcher::match( InputArray queryDescriptors, std::vector<DMatch>& matches, InputArrayOfArrays masks )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2013-02-25 00:14:01 +08:00
|
|
|
std::vector<std::vector<DMatch> > knnMatches;
|
2010-11-23 02:27:08 +08:00
|
|
|
knnMatch( queryDescriptors, knnMatches, 1, masks, true /*compactResult*/ );
|
2010-10-29 16:44:42 +08:00
|
|
|
convertMatches( knnMatches, matches );
|
|
|
|
}
|
|
|
|
|
2014-02-17 07:00:20 +08:00
|
|
|
void DescriptorMatcher::checkMasks( InputArrayOfArrays _masks, int queryDescriptorsCount ) const
|
2010-11-23 02:27:08 +08:00
|
|
|
{
|
2014-02-17 07:00:20 +08:00
|
|
|
std::vector<Mat> masks;
|
|
|
|
_masks.getMatVector(masks);
|
|
|
|
|
2011-05-02 01:38:52 +08:00
|
|
|
if( isMaskSupported() && !masks.empty() )
|
|
|
|
{
|
|
|
|
// Check masks
|
2013-12-27 17:04:02 +08:00
|
|
|
size_t imageCount = std::max(trainDescCollection.size(), utrainDescCollection.size() );
|
2011-05-02 01:38:52 +08:00
|
|
|
CV_Assert( masks.size() == imageCount );
|
|
|
|
for( size_t i = 0; i < imageCount; i++ )
|
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
if( !masks[i].empty() && (!trainDescCollection[i].empty() || !utrainDescCollection[i].empty() ) )
|
2011-05-02 01:38:52 +08:00
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
int rows = trainDescCollection[i].empty() ? utrainDescCollection[i].rows : trainDescCollection[i].rows;
|
2011-05-02 01:38:52 +08:00
|
|
|
CV_Assert( masks[i].rows == queryDescriptorsCount &&
|
2016-08-18 19:54:43 +08:00
|
|
|
masks[i].cols == rows && masks[i].type() == CV_8UC1);
|
2011-05-02 01:38:52 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2010-11-23 02:27:08 +08:00
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::knnMatch( InputArray queryDescriptors, std::vector<std::vector<DMatch> >& matches, int knn,
|
2014-02-17 07:00:20 +08:00
|
|
|
InputArrayOfArrays masks, bool compactResult )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2011-04-29 22:12:17 +08:00
|
|
|
if( empty() || queryDescriptors.empty() )
|
|
|
|
return;
|
2010-11-23 02:27:08 +08:00
|
|
|
|
2011-04-29 22:12:17 +08:00
|
|
|
CV_Assert( knn > 0 );
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
checkMasks( masks, queryDescriptors.size().height );
|
2010-11-23 02:27:08 +08:00
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
train();
|
2010-11-23 02:27:08 +08:00
|
|
|
knnMatchImpl( queryDescriptors, matches, knn, masks, compactResult );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void DescriptorMatcher::radiusMatch( InputArray queryDescriptors, std::vector<std::vector<DMatch> >& matches, float maxDistance,
|
2014-02-17 07:00:20 +08:00
|
|
|
InputArrayOfArrays masks, bool compactResult )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2011-04-29 22:12:17 +08:00
|
|
|
matches.clear();
|
|
|
|
if( empty() || queryDescriptors.empty() )
|
|
|
|
return;
|
2010-11-23 02:27:08 +08:00
|
|
|
|
2011-04-29 22:12:17 +08:00
|
|
|
CV_Assert( maxDistance > std::numeric_limits<float>::epsilon() );
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
checkMasks( masks, queryDescriptors.size().height );
|
2010-11-23 02:27:08 +08:00
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
train();
|
2010-11-23 02:27:08 +08:00
|
|
|
radiusMatchImpl( queryDescriptors, matches, maxDistance, masks, compactResult );
|
|
|
|
}
|
|
|
|
|
|
|
|
void DescriptorMatcher::read( const FileNode& )
|
|
|
|
{}
|
|
|
|
|
|
|
|
void DescriptorMatcher::write( FileStorage& ) const
|
|
|
|
{}
|
|
|
|
|
2014-02-17 07:00:20 +08:00
|
|
|
bool DescriptorMatcher::isPossibleMatch( InputArray _mask, int queryIdx, int trainIdx )
|
2010-11-23 02:27:08 +08:00
|
|
|
{
|
2014-02-17 07:00:20 +08:00
|
|
|
Mat mask = _mask.getMat();
|
2010-11-23 02:27:08 +08:00
|
|
|
return mask.empty() || mask.at<uchar>(queryIdx, trainIdx);
|
|
|
|
}
|
|
|
|
|
2014-02-17 07:00:20 +08:00
|
|
|
bool DescriptorMatcher::isMaskedOut( InputArrayOfArrays _masks, int queryIdx )
|
2010-11-23 02:27:08 +08:00
|
|
|
{
|
2014-02-17 07:00:20 +08:00
|
|
|
std::vector<Mat> masks;
|
|
|
|
_masks.getMatVector(masks);
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
size_t outCount = 0;
|
|
|
|
for( size_t i = 0; i < masks.size(); i++ )
|
|
|
|
{
|
|
|
|
if( !masks[i].empty() && (countNonZero(masks[i].row(queryIdx)) == 0) )
|
|
|
|
outCount++;
|
|
|
|
}
|
|
|
|
|
|
|
|
return !masks.empty() && outCount == masks.size() ;
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
////////////////////////////////////////////////////// BruteForceMatcher /////////////////////////////////////////////////
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
BFMatcher::BFMatcher( int _normType, bool _crossCheck )
|
2010-11-25 23:59:37 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
normType = _normType;
|
|
|
|
crossCheck = _crossCheck;
|
|
|
|
}
|
|
|
|
|
2016-10-08 06:13:01 +08:00
|
|
|
Ptr<BFMatcher> BFMatcher::create(int _normType, bool _crossCheck )
|
2016-10-08 03:58:57 +08:00
|
|
|
{
|
|
|
|
return makePtr<BFMatcher>(_normType, _crossCheck);
|
|
|
|
}
|
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
Ptr<DescriptorMatcher> BFMatcher::clone( bool emptyTrainData ) const
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
Ptr<BFMatcher> matcher = makePtr<BFMatcher>(normType, crossCheck);
|
2012-03-15 22:36:01 +08:00
|
|
|
if( !emptyTrainData )
|
2010-11-25 23:59:37 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
matcher->trainDescCollection.resize(trainDescCollection.size());
|
|
|
|
std::transform( trainDescCollection.begin(), trainDescCollection.end(),
|
|
|
|
matcher->trainDescCollection.begin(), clone_op );
|
2010-11-25 23:59:37 +08:00
|
|
|
}
|
2012-03-15 22:36:01 +08:00
|
|
|
return matcher;
|
|
|
|
}
|
|
|
|
|
2015-06-20 01:52:14 +08:00
|
|
|
#ifdef HAVE_OPENCL
|
2014-01-24 23:39:05 +08:00
|
|
|
static bool ocl_match(InputArray query, InputArray _train, std::vector< std::vector<DMatch> > &matches, int dstType)
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
UMat trainIdx, distance;
|
2014-01-24 23:39:05 +08:00
|
|
|
if (!ocl_matchSingle(query, _train, trainIdx, distance, dstType))
|
|
|
|
return false;
|
|
|
|
if (!ocl_matchDownload(trainIdx, distance, matches))
|
|
|
|
return false;
|
2013-12-27 17:04:02 +08:00
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
2014-01-24 23:39:05 +08:00
|
|
|
static bool ocl_knnMatch(InputArray query, InputArray _train, std::vector< std::vector<DMatch> > &matches, int k, int dstType, bool compactResult)
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
2014-01-17 14:09:44 +08:00
|
|
|
UMat trainIdx, distance;
|
|
|
|
if (k != 2)
|
|
|
|
return false;
|
2014-01-24 23:39:05 +08:00
|
|
|
if (!ocl_knnMatchSingle(query, _train, trainIdx, distance, dstType))
|
|
|
|
return false;
|
|
|
|
if (!ocl_knnMatchDownload(trainIdx, distance, matches, compactResult) )
|
|
|
|
return false;
|
2013-12-27 17:04:02 +08:00
|
|
|
return true;
|
|
|
|
}
|
2015-06-20 01:52:14 +08:00
|
|
|
#endif
|
2012-03-15 22:36:01 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void BFMatcher::knnMatchImpl( InputArray _queryDescriptors, std::vector<std::vector<DMatch> >& matches, int knn,
|
|
|
|
InputArrayOfArrays _masks, bool compactResult )
|
2012-03-15 22:36:01 +08:00
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
int trainDescType = trainDescCollection.empty() ? utrainDescCollection[0].type() : trainDescCollection[0].type();
|
|
|
|
CV_Assert( _queryDescriptors.type() == trainDescType );
|
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
const int IMGIDX_SHIFT = 18;
|
|
|
|
const int IMGIDX_ONE = (1 << IMGIDX_SHIFT);
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
if( _queryDescriptors.empty() || (trainDescCollection.empty() && utrainDescCollection.empty()))
|
2010-11-25 23:59:37 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
matches.clear();
|
|
|
|
return;
|
2010-11-25 23:59:37 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
std::vector<Mat> masks;
|
|
|
|
_masks.getMatVector(masks);
|
|
|
|
|
|
|
|
if(!trainDescCollection.empty() && !utrainDescCollection.empty())
|
|
|
|
{
|
|
|
|
for(int i = 0; i < (int)utrainDescCollection.size(); i++)
|
|
|
|
{
|
|
|
|
Mat tempMat;
|
|
|
|
utrainDescCollection[i].copyTo(tempMat);
|
|
|
|
trainDescCollection.push_back(tempMat);
|
|
|
|
}
|
|
|
|
utrainDescCollection.clear();
|
|
|
|
}
|
|
|
|
|
2015-06-20 01:52:14 +08:00
|
|
|
#ifdef HAVE_OPENCL
|
2013-12-27 17:04:02 +08:00
|
|
|
int trainDescVectorSize = trainDescCollection.empty() ? (int)utrainDescCollection.size() : (int)trainDescCollection.size();
|
|
|
|
Size trainDescSize = trainDescCollection.empty() ? utrainDescCollection[0].size() : trainDescCollection[0].size();
|
2014-01-15 14:01:40 +08:00
|
|
|
int trainDescOffset = trainDescCollection.empty() ? (int)utrainDescCollection[0].offset : 0;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
if ( ocl::useOpenCL() && _queryDescriptors.isUMat() && _queryDescriptors.dims()<=2 && trainDescVectorSize == 1 &&
|
2014-01-15 14:01:40 +08:00
|
|
|
_queryDescriptors.type() == CV_32FC1 && _queryDescriptors.offset() == 0 && trainDescOffset == 0 &&
|
2013-12-27 17:04:02 +08:00
|
|
|
trainDescSize.width == _queryDescriptors.size().width && masks.size() == 1 && masks[0].total() == 0 )
|
|
|
|
{
|
|
|
|
if(knn == 1)
|
|
|
|
{
|
|
|
|
if(trainDescCollection.empty())
|
|
|
|
{
|
|
|
|
if(ocl_match(_queryDescriptors, utrainDescCollection[0], matches, normType))
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
if(ocl_match(_queryDescriptors, trainDescCollection[0], matches, normType))
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
if(trainDescCollection.empty())
|
|
|
|
{
|
|
|
|
if(ocl_knnMatch(_queryDescriptors, utrainDescCollection[0], matches, knn, normType, compactResult) )
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
if(ocl_knnMatch(_queryDescriptors, trainDescCollection[0], matches, knn, normType, compactResult) )
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2015-06-20 01:52:14 +08:00
|
|
|
#endif
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
Mat queryDescriptors = _queryDescriptors.getMat();
|
|
|
|
if(trainDescCollection.empty() && !utrainDescCollection.empty())
|
|
|
|
{
|
|
|
|
for(int i = 0; i < (int)utrainDescCollection.size(); i++)
|
|
|
|
{
|
|
|
|
Mat tempMat;
|
|
|
|
utrainDescCollection[i].copyTo(tempMat);
|
|
|
|
trainDescCollection.push_back(tempMat);
|
|
|
|
}
|
|
|
|
utrainDescCollection.clear();
|
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
matches.reserve(queryDescriptors.rows);
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
Mat dist, nidx;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
int iIdx, imgCount = (int)trainDescCollection.size(), update = 0;
|
|
|
|
int dtype = normType == NORM_HAMMING || normType == NORM_HAMMING2 ||
|
|
|
|
(normType == NORM_L1 && queryDescriptors.type() == CV_8U) ? CV_32S : CV_32F;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
CV_Assert( (int64)imgCount*IMGIDX_ONE < INT_MAX );
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
for( iIdx = 0; iIdx < imgCount; iIdx++ )
|
2010-11-25 23:59:37 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
CV_Assert( trainDescCollection[iIdx].rows < IMGIDX_ONE );
|
|
|
|
batchDistance(queryDescriptors, trainDescCollection[iIdx], dist, dtype, nidx,
|
|
|
|
normType, knn, masks.empty() ? Mat() : masks[iIdx], update, crossCheck);
|
|
|
|
update += IMGIDX_ONE;
|
2010-11-25 23:59:37 +08:00
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
if( dtype == CV_32S )
|
2011-11-08 20:01:49 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
Mat temp;
|
|
|
|
dist.convertTo(temp, CV_32F);
|
|
|
|
dist = temp;
|
2011-11-08 20:01:49 +08:00
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
for( int qIdx = 0; qIdx < queryDescriptors.rows; qIdx++ )
|
2011-11-08 20:01:49 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
const float* distptr = dist.ptr<float>(qIdx);
|
|
|
|
const int* nidxptr = nidx.ptr<int>(qIdx);
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-02-25 00:14:01 +08:00
|
|
|
matches.push_back( std::vector<DMatch>() );
|
|
|
|
std::vector<DMatch>& mq = matches.back();
|
2012-03-15 22:36:01 +08:00
|
|
|
mq.reserve(knn);
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-04-14 05:50:59 +08:00
|
|
|
for( int k = 0; k < nidx.cols; k++ )
|
2012-03-15 22:36:01 +08:00
|
|
|
{
|
|
|
|
if( nidxptr[k] < 0 )
|
|
|
|
break;
|
|
|
|
mq.push_back( DMatch(qIdx, nidxptr[k] & (IMGIDX_ONE - 1),
|
|
|
|
nidxptr[k] >> IMGIDX_SHIFT, distptr[k]) );
|
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
if( mq.empty() && compactResult )
|
|
|
|
matches.pop_back();
|
2011-11-08 20:01:49 +08:00
|
|
|
}
|
2010-11-25 23:59:37 +08:00
|
|
|
}
|
2010-11-23 02:27:08 +08:00
|
|
|
|
2015-06-20 01:52:14 +08:00
|
|
|
#ifdef HAVE_OPENCL
|
2014-01-24 23:39:05 +08:00
|
|
|
static bool ocl_radiusMatch(InputArray query, InputArray _train, std::vector< std::vector<DMatch> > &matches,
|
2013-12-27 17:04:02 +08:00
|
|
|
float maxDistance, int dstType, bool compactResult)
|
|
|
|
{
|
|
|
|
UMat trainIdx, distance, nMatches;
|
2014-01-24 23:39:05 +08:00
|
|
|
if (!ocl_radiusMatchSingle(query, _train, trainIdx, distance, nMatches, maxDistance, dstType))
|
|
|
|
return false;
|
|
|
|
if (!ocl_radiusMatchDownload(trainIdx, distance, nMatches, matches, compactResult))
|
|
|
|
return false;
|
2013-12-27 17:04:02 +08:00
|
|
|
return true;
|
|
|
|
}
|
2015-06-20 01:52:14 +08:00
|
|
|
#endif
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void BFMatcher::radiusMatchImpl( InputArray _queryDescriptors, std::vector<std::vector<DMatch> >& matches,
|
|
|
|
float maxDistance, InputArrayOfArrays _masks, bool compactResult )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2013-12-27 17:04:02 +08:00
|
|
|
int trainDescType = trainDescCollection.empty() ? utrainDescCollection[0].type() : trainDescCollection[0].type();
|
|
|
|
CV_Assert( _queryDescriptors.type() == trainDescType );
|
|
|
|
|
|
|
|
if( _queryDescriptors.empty() || (trainDescCollection.empty() && utrainDescCollection.empty()))
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
matches.clear();
|
|
|
|
return;
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
std::vector<Mat> masks;
|
|
|
|
_masks.getMatVector(masks);
|
|
|
|
|
|
|
|
if(!trainDescCollection.empty() && !utrainDescCollection.empty())
|
|
|
|
{
|
|
|
|
for(int i = 0; i < (int)utrainDescCollection.size(); i++)
|
|
|
|
{
|
|
|
|
Mat tempMat;
|
|
|
|
utrainDescCollection[i].copyTo(tempMat);
|
|
|
|
trainDescCollection.push_back(tempMat);
|
|
|
|
}
|
|
|
|
utrainDescCollection.clear();
|
|
|
|
}
|
|
|
|
|
2015-06-20 01:52:14 +08:00
|
|
|
#ifdef HAVE_OPENCL
|
2013-12-27 17:04:02 +08:00
|
|
|
int trainDescVectorSize = trainDescCollection.empty() ? (int)utrainDescCollection.size() : (int)trainDescCollection.size();
|
|
|
|
Size trainDescSize = trainDescCollection.empty() ? utrainDescCollection[0].size() : trainDescCollection[0].size();
|
2014-01-15 14:01:40 +08:00
|
|
|
int trainDescOffset = trainDescCollection.empty() ? (int)utrainDescCollection[0].offset : 0;
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
if ( ocl::useOpenCL() && _queryDescriptors.isUMat() && _queryDescriptors.dims()<=2 && trainDescVectorSize == 1 &&
|
2014-01-15 14:01:40 +08:00
|
|
|
_queryDescriptors.type() == CV_32FC1 && _queryDescriptors.offset() == 0 && trainDescOffset == 0 &&
|
2013-12-27 17:04:02 +08:00
|
|
|
trainDescSize.width == _queryDescriptors.size().width && masks.size() == 1 && masks[0].total() == 0 )
|
|
|
|
{
|
2014-01-24 23:39:05 +08:00
|
|
|
if (trainDescCollection.empty())
|
2013-12-27 17:04:02 +08:00
|
|
|
{
|
|
|
|
if(ocl_radiusMatch(_queryDescriptors, utrainDescCollection[0], matches, maxDistance, normType, compactResult) )
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
2014-01-24 23:39:05 +08:00
|
|
|
if (ocl_radiusMatch(_queryDescriptors, trainDescCollection[0], matches, maxDistance, normType, compactResult) )
|
2014-10-03 19:17:28 +08:00
|
|
|
{
|
|
|
|
CV_IMPL_ADD(CV_IMPL_OCL);
|
2013-12-27 17:04:02 +08:00
|
|
|
return;
|
2014-10-03 19:17:28 +08:00
|
|
|
}
|
2013-12-27 17:04:02 +08:00
|
|
|
}
|
|
|
|
}
|
2015-06-20 01:52:14 +08:00
|
|
|
#endif
|
2013-12-27 17:04:02 +08:00
|
|
|
|
|
|
|
Mat queryDescriptors = _queryDescriptors.getMat();
|
|
|
|
if(trainDescCollection.empty() && !utrainDescCollection.empty())
|
|
|
|
{
|
|
|
|
for(int i = 0; i < (int)utrainDescCollection.size(); i++)
|
|
|
|
{
|
|
|
|
Mat tempMat;
|
|
|
|
utrainDescCollection[i].copyTo(tempMat);
|
|
|
|
trainDescCollection.push_back(tempMat);
|
|
|
|
}
|
|
|
|
utrainDescCollection.clear();
|
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
matches.resize(queryDescriptors.rows);
|
|
|
|
Mat dist, distf;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
int iIdx, imgCount = (int)trainDescCollection.size();
|
2017-08-31 23:24:02 +08:00
|
|
|
int dtype = normType == NORM_HAMMING || normType == NORM_HAMMING2 ||
|
2012-03-15 22:36:01 +08:00
|
|
|
(normType == NORM_L1 && queryDescriptors.type() == CV_8U) ? CV_32S : CV_32F;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
for( iIdx = 0; iIdx < imgCount; iIdx++ )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
batchDistance(queryDescriptors, trainDescCollection[iIdx], dist, dtype, noArray(),
|
|
|
|
normType, 0, masks.empty() ? Mat() : masks[iIdx], 0, false);
|
|
|
|
if( dtype == CV_32S )
|
|
|
|
dist.convertTo(distf, CV_32F);
|
2010-09-23 21:44:23 +08:00
|
|
|
else
|
2012-03-15 22:36:01 +08:00
|
|
|
distf = dist;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
for( int qIdx = 0; qIdx < queryDescriptors.rows; qIdx++ )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2012-08-15 19:02:20 +08:00
|
|
|
const float* distptr = distf.ptr<float>(qIdx);
|
|
|
|
|
2013-02-25 00:14:01 +08:00
|
|
|
std::vector<DMatch>& mq = matches[qIdx];
|
2012-08-15 19:02:20 +08:00
|
|
|
for( int k = 0; k < distf.cols; k++ )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2012-03-15 22:36:01 +08:00
|
|
|
if( distptr[k] <= maxDistance )
|
|
|
|
mq.push_back( DMatch(qIdx, k, iIdx, distptr[k]) );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
int qIdx0 = 0;
|
|
|
|
for( int qIdx = 0; qIdx < queryDescriptors.rows; qIdx++ )
|
|
|
|
{
|
|
|
|
if( matches[qIdx].empty() && compactResult )
|
|
|
|
continue;
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
if( qIdx0 < qIdx )
|
|
|
|
std::swap(matches[qIdx], matches[qIdx0]);
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
std::sort( matches[qIdx0].begin(), matches[qIdx0].end() );
|
|
|
|
qIdx0++;
|
|
|
|
}
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
///////////////////////////////////////////////////////////////////////////////////////////////////////
|
2012-08-15 19:02:20 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
/*
|
|
|
|
* Factory function for DescriptorMatcher creating
|
|
|
|
*/
|
2013-03-23 00:37:49 +08:00
|
|
|
Ptr<DescriptorMatcher> DescriptorMatcher::create( const String& descriptorMatcherType )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
Ptr<DescriptorMatcher> dm;
|
2012-03-15 22:36:01 +08:00
|
|
|
if( !descriptorMatcherType.compare( "FlannBased" ) )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<FlannBasedMatcher>();
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
2012-03-15 22:36:01 +08:00
|
|
|
else if( !descriptorMatcherType.compare( "BruteForce" ) ) // L2
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<BFMatcher>(int(NORM_L2)); // anonymous enums can't be template parameters
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
2012-03-15 22:36:01 +08:00
|
|
|
else if( !descriptorMatcherType.compare( "BruteForce-SL2" ) ) // Squared L2
|
2011-08-08 21:18:12 +08:00
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<BFMatcher>(int(NORM_L2SQR));
|
2011-08-08 21:18:12 +08:00
|
|
|
}
|
2012-03-15 22:36:01 +08:00
|
|
|
else if( !descriptorMatcherType.compare( "BruteForce-L1" ) )
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<BFMatcher>(int(NORM_L1));
|
2012-03-15 22:36:01 +08:00
|
|
|
}
|
|
|
|
else if( !descriptorMatcherType.compare("BruteForce-Hamming") ||
|
|
|
|
!descriptorMatcherType.compare("BruteForce-HammingLUT") )
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<BFMatcher>(int(NORM_HAMMING));
|
2012-03-15 22:36:01 +08:00
|
|
|
}
|
|
|
|
else if( !descriptorMatcherType.compare("BruteForce-Hamming(2)") )
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
dm = makePtr<BFMatcher>(int(NORM_HAMMING2));
|
2012-03-15 22:36:01 +08:00
|
|
|
}
|
|
|
|
else
|
2013-04-10 19:54:14 +08:00
|
|
|
CV_Error( Error::StsBadArg, "Unknown matcher name" );
|
2011-08-08 21:18:12 +08:00
|
|
|
|
2012-03-15 22:36:01 +08:00
|
|
|
return dm;
|
2011-08-08 21:18:12 +08:00
|
|
|
}
|
|
|
|
|
2016-10-08 03:58:57 +08:00
|
|
|
Ptr<DescriptorMatcher> DescriptorMatcher::create(int matcherType)
|
|
|
|
{
|
2016-10-08 06:13:01 +08:00
|
|
|
|
|
|
|
|
2016-10-08 03:58:57 +08:00
|
|
|
String name;
|
|
|
|
|
|
|
|
switch(matcherType)
|
|
|
|
{
|
|
|
|
case FLANNBASED:
|
|
|
|
name = "FlannBased";
|
|
|
|
break;
|
|
|
|
case BRUTEFORCE:
|
|
|
|
name = "BruteForce";
|
|
|
|
break;
|
|
|
|
case BRUTEFORCE_L1:
|
|
|
|
name = "BruteForce-L1";
|
|
|
|
break;
|
|
|
|
case BRUTEFORCE_HAMMING:
|
|
|
|
name = "BruteForce-Hamming";
|
|
|
|
break;
|
|
|
|
case BRUTEFORCE_HAMMINGLUT:
|
|
|
|
name = "BruteForce-HammingLUT";
|
|
|
|
break;
|
|
|
|
case BRUTEFORCE_SL2:
|
|
|
|
name = "BruteForce-SL2";
|
|
|
|
break;
|
|
|
|
default:
|
|
|
|
CV_Error( Error::StsBadArg, "Specified descriptor matcher type is not supported." );
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
|
|
|
|
return DescriptorMatcher::create(name);
|
2016-10-08 06:13:01 +08:00
|
|
|
|
2016-10-08 03:58:57 +08:00
|
|
|
}
|
|
|
|
|
2011-08-08 21:18:12 +08:00
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
/*
|
|
|
|
* Flann based matcher
|
|
|
|
*/
|
|
|
|
FlannBasedMatcher::FlannBasedMatcher( const Ptr<flann::IndexParams>& _indexParams, const Ptr<flann::SearchParams>& _searchParams )
|
|
|
|
: indexParams(_indexParams), searchParams(_searchParams), addedDescCount(0)
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
CV_Assert( _indexParams );
|
|
|
|
CV_Assert( _searchParams );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
2016-10-08 06:13:01 +08:00
|
|
|
Ptr<FlannBasedMatcher> FlannBasedMatcher::create()
|
2016-10-08 03:58:57 +08:00
|
|
|
{
|
|
|
|
return makePtr<FlannBasedMatcher>();
|
|
|
|
}
|
|
|
|
|
2014-02-17 05:37:32 +08:00
|
|
|
void FlannBasedMatcher::add( InputArrayOfArrays _descriptors )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2014-02-17 05:37:32 +08:00
|
|
|
DescriptorMatcher::add( _descriptors );
|
|
|
|
|
2016-01-23 23:58:16 +08:00
|
|
|
if( _descriptors.isUMatVector() )
|
|
|
|
{
|
|
|
|
std::vector<UMat> descriptors;
|
|
|
|
_descriptors.getUMatVector( descriptors );
|
|
|
|
|
|
|
|
for( size_t i = 0; i < descriptors.size(); i++ )
|
|
|
|
{
|
|
|
|
addedDescCount += descriptors[i].rows;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
else if( _descriptors.isUMat() )
|
|
|
|
{
|
|
|
|
addedDescCount += _descriptors.getUMat().rows;
|
|
|
|
}
|
|
|
|
else if( _descriptors.isMatVector() )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-01-23 23:58:16 +08:00
|
|
|
std::vector<Mat> descriptors;
|
|
|
|
_descriptors.getMatVector(descriptors);
|
|
|
|
for( size_t i = 0; i < descriptors.size(); i++ )
|
|
|
|
{
|
|
|
|
addedDescCount += descriptors[i].rows;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
else if( _descriptors.isMat() )
|
|
|
|
{
|
|
|
|
addedDescCount += _descriptors.getMat().rows;
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
CV_Assert( _descriptors.isUMat() || _descriptors.isUMatVector() || _descriptors.isMat() || _descriptors.isMatVector() );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
void FlannBasedMatcher::clear()
|
|
|
|
{
|
|
|
|
DescriptorMatcher::clear();
|
|
|
|
|
|
|
|
mergedDescriptors.clear();
|
|
|
|
flannIndex.release();
|
|
|
|
|
|
|
|
addedDescCount = 0;
|
|
|
|
}
|
|
|
|
|
|
|
|
void FlannBasedMatcher::train()
|
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2013-08-13 21:30:14 +08:00
|
|
|
if( !flannIndex || mergedDescriptors.size() < addedDescCount )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2014-03-12 19:54:55 +08:00
|
|
|
// FIXIT: Workaround for 'utrainDescCollection' issue (PR #2142)
|
|
|
|
if (!utrainDescCollection.empty())
|
|
|
|
{
|
|
|
|
CV_Assert(trainDescCollection.size() == 0);
|
|
|
|
for (size_t i = 0; i < utrainDescCollection.size(); ++i)
|
|
|
|
trainDescCollection.push_back(utrainDescCollection[i].getMat(ACCESS_READ));
|
|
|
|
}
|
2010-10-29 16:44:42 +08:00
|
|
|
mergedDescriptors.set( trainDescCollection );
|
2013-08-13 21:30:14 +08:00
|
|
|
flannIndex = makePtr<flann::Index>( mergedDescriptors.getDescriptors(), *indexParams );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
}
|
2010-09-23 21:44:23 +08:00
|
|
|
|
2011-08-04 22:13:07 +08:00
|
|
|
void FlannBasedMatcher::read( const FileNode& fn)
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
if (!indexParams)
|
|
|
|
indexParams = makePtr<flann::IndexParams>();
|
2011-08-04 22:13:07 +08:00
|
|
|
|
|
|
|
FileNode ip = fn["indexParams"];
|
|
|
|
CV_Assert(ip.type() == FileNode::SEQ);
|
|
|
|
|
2012-03-17 05:21:04 +08:00
|
|
|
for(int i = 0; i < (int)ip.size(); ++i)
|
2011-08-04 22:13:07 +08:00
|
|
|
{
|
|
|
|
CV_Assert(ip[i].type() == FileNode::MAP);
|
2013-03-23 00:37:49 +08:00
|
|
|
String _name = (String)ip[i]["name"];
|
2011-08-04 22:13:07 +08:00
|
|
|
int type = (int)ip[i]["type"];
|
|
|
|
|
|
|
|
switch(type)
|
|
|
|
{
|
|
|
|
case CV_8U:
|
|
|
|
case CV_8S:
|
|
|
|
case CV_16U:
|
|
|
|
case CV_16S:
|
|
|
|
case CV_32S:
|
2012-06-09 23:00:04 +08:00
|
|
|
indexParams->setInt(_name, (int) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_32F:
|
2012-06-09 23:00:04 +08:00
|
|
|
indexParams->setFloat(_name, (float) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_64F:
|
2012-06-09 23:00:04 +08:00
|
|
|
indexParams->setDouble(_name, (double) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_USRTYPE1:
|
2013-03-23 00:37:49 +08:00
|
|
|
indexParams->setString(_name, (String) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,2):
|
2012-06-09 23:00:04 +08:00
|
|
|
indexParams->setBool(_name, (int) ip[i]["value"] != 0);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,3):
|
2011-08-15 03:46:39 +08:00
|
|
|
indexParams->setAlgorithm((int) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
};
|
|
|
|
}
|
|
|
|
|
2013-08-13 21:30:14 +08:00
|
|
|
if (!searchParams)
|
|
|
|
searchParams = makePtr<flann::SearchParams>();
|
2011-08-04 22:13:07 +08:00
|
|
|
|
|
|
|
FileNode sp = fn["searchParams"];
|
|
|
|
CV_Assert(sp.type() == FileNode::SEQ);
|
|
|
|
|
2012-03-17 05:21:04 +08:00
|
|
|
for(int i = 0; i < (int)sp.size(); ++i)
|
2011-08-04 22:13:07 +08:00
|
|
|
{
|
|
|
|
CV_Assert(sp[i].type() == FileNode::MAP);
|
2013-03-23 00:37:49 +08:00
|
|
|
String _name = (String)sp[i]["name"];
|
2011-08-04 22:13:07 +08:00
|
|
|
int type = (int)sp[i]["type"];
|
|
|
|
|
|
|
|
switch(type)
|
|
|
|
{
|
|
|
|
case CV_8U:
|
|
|
|
case CV_8S:
|
|
|
|
case CV_16U:
|
|
|
|
case CV_16S:
|
|
|
|
case CV_32S:
|
2012-06-09 23:00:04 +08:00
|
|
|
searchParams->setInt(_name, (int) sp[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_32F:
|
2012-06-09 23:00:04 +08:00
|
|
|
searchParams->setFloat(_name, (float) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_64F:
|
2012-06-09 23:00:04 +08:00
|
|
|
searchParams->setDouble(_name, (double) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_USRTYPE1:
|
2013-03-23 00:37:49 +08:00
|
|
|
searchParams->setString(_name, (String) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,2):
|
2012-06-09 23:00:04 +08:00
|
|
|
searchParams->setBool(_name, (int) ip[i]["value"] != 0);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,3):
|
2011-08-15 03:46:39 +08:00
|
|
|
searchParams->setAlgorithm((int) ip[i]["value"]);
|
2011-08-04 22:13:07 +08:00
|
|
|
break;
|
|
|
|
};
|
|
|
|
}
|
|
|
|
|
|
|
|
flannIndex.release();
|
|
|
|
}
|
|
|
|
|
|
|
|
void FlannBasedMatcher::write( FileStorage& fs) const
|
|
|
|
{
|
2016-03-23 06:19:42 +08:00
|
|
|
writeFormat(fs);
|
2011-08-04 22:13:07 +08:00
|
|
|
fs << "indexParams" << "[";
|
|
|
|
|
2012-12-22 03:27:29 +08:00
|
|
|
if (indexParams)
|
2011-08-04 22:13:07 +08:00
|
|
|
{
|
2013-03-23 00:37:49 +08:00
|
|
|
std::vector<String> names;
|
2011-08-04 22:13:07 +08:00
|
|
|
std::vector<int> types;
|
2013-03-23 00:37:49 +08:00
|
|
|
std::vector<String> strValues;
|
2011-08-04 22:13:07 +08:00
|
|
|
std::vector<double> numValues;
|
|
|
|
|
|
|
|
indexParams->getAll(names, types, strValues, numValues);
|
|
|
|
|
|
|
|
for(size_t i = 0; i < names.size(); ++i)
|
|
|
|
{
|
|
|
|
fs << "{" << "name" << names[i] << "type" << types[i] << "value";
|
|
|
|
switch(types[i])
|
|
|
|
{
|
|
|
|
case CV_8U:
|
|
|
|
fs << (uchar)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_8S:
|
|
|
|
fs << (char)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_16U:
|
|
|
|
fs << (ushort)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_16S:
|
|
|
|
fs << (short)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_32S:
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,2):
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,3):
|
|
|
|
fs << (int)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_32F:
|
|
|
|
fs << (float)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_64F:
|
|
|
|
fs << (double)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_USRTYPE1:
|
|
|
|
fs << strValues[i];
|
|
|
|
break;
|
|
|
|
default:
|
|
|
|
fs << (double)numValues[i];
|
|
|
|
fs << "typename" << strValues[i];
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
fs << "}";
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
fs << "]" << "searchParams" << "[";
|
|
|
|
|
2012-12-22 03:27:29 +08:00
|
|
|
if (searchParams)
|
2011-08-04 22:13:07 +08:00
|
|
|
{
|
2013-03-23 00:37:49 +08:00
|
|
|
std::vector<String> names;
|
2011-08-04 22:13:07 +08:00
|
|
|
std::vector<int> types;
|
2013-03-23 00:37:49 +08:00
|
|
|
std::vector<String> strValues;
|
2011-08-04 22:13:07 +08:00
|
|
|
std::vector<double> numValues;
|
|
|
|
|
|
|
|
searchParams->getAll(names, types, strValues, numValues);
|
|
|
|
|
|
|
|
for(size_t i = 0; i < names.size(); ++i)
|
|
|
|
{
|
|
|
|
fs << "{" << "name" << names[i] << "type" << types[i] << "value";
|
|
|
|
switch(types[i])
|
|
|
|
{
|
|
|
|
case CV_8U:
|
|
|
|
fs << (uchar)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_8S:
|
|
|
|
fs << (char)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_16U:
|
|
|
|
fs << (ushort)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_16S:
|
|
|
|
fs << (short)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_32S:
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,2):
|
|
|
|
case CV_MAKETYPE(CV_USRTYPE1,3):
|
|
|
|
fs << (int)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_32F:
|
|
|
|
fs << (float)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_64F:
|
|
|
|
fs << (double)numValues[i];
|
|
|
|
break;
|
|
|
|
case CV_USRTYPE1:
|
|
|
|
fs << strValues[i];
|
|
|
|
break;
|
|
|
|
default:
|
|
|
|
fs << (double)numValues[i];
|
|
|
|
fs << "typename" << strValues[i];
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
fs << "}";
|
|
|
|
}
|
|
|
|
}
|
|
|
|
fs << "]";
|
|
|
|
}
|
|
|
|
|
2010-11-23 02:27:08 +08:00
|
|
|
bool FlannBasedMatcher::isMaskSupported() const
|
|
|
|
{
|
|
|
|
return false;
|
|
|
|
}
|
|
|
|
|
|
|
|
Ptr<DescriptorMatcher> FlannBasedMatcher::clone( bool emptyTrainData ) const
|
|
|
|
{
|
2013-08-13 21:30:14 +08:00
|
|
|
Ptr<FlannBasedMatcher> matcher = makePtr<FlannBasedMatcher>(indexParams, searchParams);
|
2010-11-23 02:27:08 +08:00
|
|
|
if( !emptyTrainData )
|
|
|
|
{
|
2013-04-10 19:54:14 +08:00
|
|
|
CV_Error( Error::StsNotImplemented, "deep clone functionality is not implemented, because "
|
2010-11-23 02:27:08 +08:00
|
|
|
"Flann::Index has not copy constructor or clone method ");
|
|
|
|
//matcher->flannIndex;
|
|
|
|
matcher->addedDescCount = addedDescCount;
|
|
|
|
matcher->mergedDescriptors = DescriptorCollection( mergedDescriptors );
|
2011-01-12 14:39:08 +08:00
|
|
|
std::transform( trainDescCollection.begin(), trainDescCollection.end(),
|
|
|
|
matcher->trainDescCollection.begin(), clone_op );
|
2010-11-23 02:27:08 +08:00
|
|
|
}
|
|
|
|
return matcher;
|
|
|
|
}
|
|
|
|
|
2010-10-29 16:44:42 +08:00
|
|
|
void FlannBasedMatcher::convertToDMatches( const DescriptorCollection& collection, const Mat& indices, const Mat& dists,
|
2013-02-25 00:14:01 +08:00
|
|
|
std::vector<std::vector<DMatch> >& matches )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
|
|
|
matches.resize( indices.rows );
|
|
|
|
for( int i = 0; i < indices.rows; i++ )
|
|
|
|
{
|
|
|
|
for( int j = 0; j < indices.cols; j++ )
|
|
|
|
{
|
|
|
|
int idx = indices.at<int>(i, j);
|
|
|
|
if( idx >= 0 )
|
2010-09-23 21:44:23 +08:00
|
|
|
{
|
2010-10-29 16:44:42 +08:00
|
|
|
int imgIdx, trainIdx;
|
|
|
|
collection.getLocalIdx( idx, imgIdx, trainIdx );
|
2011-09-22 16:52:40 +08:00
|
|
|
float dist = 0;
|
|
|
|
if (dists.type() == CV_32S)
|
|
|
|
dist = static_cast<float>( dists.at<int>(i,j) );
|
|
|
|
else
|
|
|
|
dist = std::sqrt(dists.at<float>(i,j));
|
|
|
|
matches[i].push_back( DMatch( i, trainIdx, imgIdx, dist ) );
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void FlannBasedMatcher::knnMatchImpl( InputArray _queryDescriptors, std::vector<std::vector<DMatch> >& matches, int knn,
|
|
|
|
InputArrayOfArrays /*masks*/, bool /*compactResult*/ )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
Mat queryDescriptors = _queryDescriptors.getMat();
|
2010-11-23 02:27:08 +08:00
|
|
|
Mat indices( queryDescriptors.rows, knn, CV_32SC1 );
|
|
|
|
Mat dists( queryDescriptors.rows, knn, CV_32FC1);
|
|
|
|
flannIndex->knnSearch( queryDescriptors, indices, dists, knn, *searchParams );
|
2010-10-29 16:44:42 +08:00
|
|
|
|
|
|
|
convertToDMatches( mergedDescriptors, indices, dists, matches );
|
|
|
|
}
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
void FlannBasedMatcher::radiusMatchImpl( InputArray _queryDescriptors, std::vector<std::vector<DMatch> >& matches, float maxDistance,
|
|
|
|
InputArrayOfArrays /*masks*/, bool /*compactResult*/ )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2016-08-18 14:53:00 +08:00
|
|
|
CV_INSTRUMENT_REGION()
|
|
|
|
|
2013-12-27 17:04:02 +08:00
|
|
|
Mat queryDescriptors = _queryDescriptors.getMat();
|
2010-10-29 16:44:42 +08:00
|
|
|
const int count = mergedDescriptors.size(); // TODO do count as param?
|
2010-11-23 02:27:08 +08:00
|
|
|
Mat indices( queryDescriptors.rows, count, CV_32SC1, Scalar::all(-1) );
|
|
|
|
Mat dists( queryDescriptors.rows, count, CV_32FC1, Scalar::all(-1) );
|
|
|
|
for( int qIdx = 0; qIdx < queryDescriptors.rows; qIdx++ )
|
2010-10-29 16:44:42 +08:00
|
|
|
{
|
2010-11-23 02:27:08 +08:00
|
|
|
Mat queryDescriptorsRow = queryDescriptors.row(qIdx);
|
2010-10-29 16:44:42 +08:00
|
|
|
Mat indicesRow = indices.row(qIdx);
|
|
|
|
Mat distsRow = dists.row(qIdx);
|
2011-07-14 07:04:39 +08:00
|
|
|
flannIndex->radiusSearch( queryDescriptorsRow, indicesRow, distsRow, maxDistance*maxDistance, count, *searchParams );
|
2010-10-29 16:44:42 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
convertToDMatches( mergedDescriptors, indices, dists, matches );
|
|
|
|
}
|
2010-09-23 21:44:23 +08:00
|
|
|
}
|