我设法使VlFeat的SIFT实现得以实现,并且我想尝试匹配两组图像描述符。
SIFT的特征向量是128个元素的浮点数组,我将描述符列表存储在std::vector
中,如下面的代码片段所示:
std::vector<std::vector<float> > ldescriptors = leftImage->descriptors;
std::vector<std::vector<float> > rdescriptors = rightImage->descriptors;
/* KDTree, L1 comparison metric, dimension 128, 1 tree, L1 metric */
VlKDForest* forest = vl_kdforest_new(VL_TYPE_FLOAT, 128, 1, VlDistanceL1);
/* Build the tree from the left descriptors */
vl_kdforest_build(forest, ldescriptors.size(), ldescriptors.data());
/* Searcher object */
VlKDForestSearcher* searcher = vl_kdforest_new_searcher(forest);
VlKDForestNeighbor neighbours[2];
/* Query the first ten points for now */
for(int i=0; i < 10; i++){
int nvisited = vl_kdforestsearcher_query(searcher, &neighbours, 2, rdescriptors[i].data());
cout << nvisited << neighbours[0].distance << neighbours[1].distance;
}
据我所知,这应该行得通,但我得到的只是距离
nan
。描述符数组检出的长度,因此似乎确实有数据进入树中。我已经绘制了关键点,它们看起来也很合理,因此数据相当合理。我想念什么?
此处的文档很少(链接到API):http://www.vlfeat.org/api/kdtree.html
最佳答案
我想念什么?vl_kdforestsearcher_query
的第二个参数采用指向VlKDForestNeighbor
的指针:
vl_size
vl_kdforestsearcher_query(
VlKDForestSearcher *self,
VlKDForestNeighbor *neighbors,
vl_size numNeighbors,
void const *query
);
但是在这里,您声明了
VlKDForestNeighbor neighbours[2];
,然后将&neighbours
作为第二个参数传递了-这是不正确的-您的编译器可能发出了incompatible pointer types
警告。由于声明了数组,因此必须做的就是显式传递指向第一个邻居的指针:
int nvisited = vl_kdforestsearcher_query(searcher, &neighbours[0], 2, qrys[i]);
或者让编译器为您完成:
int nvisited = vl_kdforestsearcher_query(searcher, neighbours, 2, qrys[i]);
编辑
实际上,还有第二个(主要)问题与使用
ldescriptors.data()
构建kd树的方式有关。当VLFeat期望
std::vector<float>*
连续数组包含按行主要顺序排列的所有数据点时,在此处传递float *
指针。因此,您可以按照以下格式复制数据:float *data = new float[128*ldescriptors.size()];
for (unsigned int i = 0; i < ldescriptors.size(); i++)
std::copy(ldescriptors[i].begin(), ldescriptors[i].end(), data + 128*i);
vl_kdforest_build(forest, ldescriptors.size(), data);
// ...
// then, right after `vl_kdforest_delete(forest);`
// do a `delete[] data;`
关于c++ - VLeat Kdtree设置和查询,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/28606011/