- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在编写一个算法,要求我搜索点的最近邻居。我从这篇文章( Using Google's C KD Tree Library )中找到了 kdtree 库,但它没有从树中删除单个节点的功能。所以我开始实现我自己的使用www(点)geeksforgeeks.org/k-Dimension-tree-set-3-delete/作为模板。这一切都贯穿始终,但不幸的是有时节点会重复。我的测试用例如下:
#include <stdio.h>
#include <assert.h>
#include <stdlib.h>
#include <math.h>
#include <errno.h>
#include <string.h>
#include <stdarg.h>
#include "kdtree.h"
/* (hopefully) platform independent directory creation */
#if defined(_WIN32) || defined(WIN32) /* this should be defined under windows, regardless of 64 or 32 bit*/
#include <direct.h>
#include <sys/stat.h>
#define GetWorkingDir _getcwd
#define MakeDir(str) _mkdir(str)
#else /* unix based system */
#include <unistd.h>
#include <sys/stat.h>
#define GetWorkingDir getcwd
#define MakeDir(str) mkdir(str, 0777)
#endif
#ifndef MAX_PATH
#define MAX_PATH 260
#endif
void GetLogDir(char* strPath, int nBufSize)
{
if(GetWorkingDir(strPath, nBufSize))
{
strncat(strPath, "/log/", 5);
MakeDir(strPath);
}
else
{
fprintf(stderr, "Could not get working directory");
exit(ENOENT);
}
}
FILE* GetOpenFileHandle(const char* strFilenamePlusPath, const char* strOpenMode)
{
if(strOpenMode == NULL) // too bad we dont have default arguments in C :(
{
strOpenMode = "a+";
}
return(fopen(strFilenamePlusPath, strOpenMode));
}
int CloseFile(FILE* pFile)
{
if(pFile != NULL)
{
fprintf(pFile, "\r\n"); // append a new line before closing!
return(fclose(pFile));
}
fprintf(stderr, "Invalid file handle");
exit(EFAULT);
}
void NodeLabelToFile(FILE* pFile, kdnode* node, const char* strName)
{
fprintf(pFile, "%s [label=\"(%.3f, %.3f)\"] \n", strName, node->pos[0], node->pos[1]);
}
char* NodeToString(kdnode* node, int* num)
{
char* strName = (char*) malloc(MAX_PATH);
if(*num == 0)
{
sprintf(strName, "%s","root");
}
else
{
sprintf(strName, "node%d", *num);
}
return strName;
}
void NodesToFile(FILE* pFile, kdnode* node, const char* strParentname, int* num)
{
if(node && pFile)
{
char* strLeft = NULL;
char* strRight = NULL;
if(node->left)
{
(*num)++;
strLeft = NodeToString(node->left, num);
NodeLabelToFile(pFile, node->left, strLeft);
fprintf(pFile, "%s -> %s \n", strParentname, strLeft);
}
if(node->right)
{
(*num)++;
strRight = NodeToString(node->right, num); // name of the current node
NodeLabelToFile(pFile, node->right, strRight);
fprintf(pFile, "%s -> %s \n", strParentname, strRight);
}
if(strLeft)
{
NodesToFile(pFile, node->left, strLeft, num);
free(strLeft);
}
if(strRight)
{
// (*num)++;
NodesToFile(pFile, node->right, strRight, num);
free(strRight);
}
}
}
FILE* MakeOpenLogFile(const char* strFilename, const char* strOpenMode)
{
if(strOpenMode == NULL)
{
strOpenMode = "a+";
}
char* strFilenamePlusPath = (char*) malloc(MAX_PATH);
GetLogDir(strFilenamePlusPath, MAX_PATH);
strncat(strFilenamePlusPath, strFilename, strlen(strFilename));
FILE* pFile = GetOpenFileHandle(strFilenamePlusPath, strOpenMode);
free(strFilenamePlusPath);
return(pFile);
}
void KDTreeToDotFile(kdtree* Tree, const char* strFilename)
{
if(Tree)
{
FILE* pFile = MakeOpenLogFile(strFilename, "w");
fprintf(pFile, "%s", "digraph d { \n"); // print opening statement for the graph in dot language
// traverse the tree and print the nodes
int* num = (int*) malloc(sizeof(int)); // make this a unique location to make sure numbers can't occur twice
*num = 0;
char* strRoot = NodeToString(Tree->root, num);
NodeLabelToFile(pFile, Tree->root, strRoot);
NodesToFile(pFile, Tree->root, "root", num);
if(strRoot)
{
free(strRoot);
}
free(num);
fprintf(pFile,"%s", "}"); // close the digraph environment
CloseFile(pFile);
}
}
int main(int argc, const char * argv[])
{
int numel = 20;
int toRemove = 19;
double dMax = 3000;
int nNumDim = 2;
printf("init rng");
srand(1234); // seed the rng // srand((unsigned) time(&t));
printf("creating kdtree");
kdtree* TreeRoot = kd_create(nNumDim); // construct the kd tree for the nearest neighbor search
kd_data_destructor(TreeRoot, free); // set free as data destructor
double* pos = (double*) malloc(nNumDim * numel * sizeof(double));
int retval;
for (int ii = 0; ii < numel; ii++)
{
pos[nNumDim * ii] = floor((double)rand()/(double)(RAND_MAX/dMax));
pos[nNumDim * ii + 1] = floor((double)rand()/(double)(RAND_MAX/dMax));
int* randint = (int*) malloc(sizeof(int));
*randint = rand();
retval = kd_insert2(TreeRoot,
pos[nNumDim * ii],
pos[nNumDim * ii + 1],
randint, sizeof(int));
assert(retval == 0);
}
KDTreeToDotFile(TreeRoot, "original.dot");
double* dRemovePos = (double*) malloc(sizeof(double)*nNumDim);
for (int ii = 0; ii < toRemove; ii++)
{
dRemovePos[0] = pos[2*ii];
dRemovePos[1] = pos[2*ii + 1];
kd_remove(TreeRoot, dRemovePos);
}
KDTreeToDotFile(TreeRoot, "removed.dot");
kd_free(TreeRoot); // free kdtree
return 0;
}
删除节点的函数是这样实现的:(我不认为代码太多,所以我只会将我的更改发布到 kd 库。如果我应该添加其余的代码(不幸的是超过 1000 行),请在注释中告诉我。 )
int kd_remove(kdtree* tree, const double* pos)
{
printf("removing node %.3f, %.3f \n", pos[0], pos[1]);
if(tree->root != NULL)
{
assert(tree->dim != 0); // prevent division by 0 (error code 136)
assert(pos != NULL); // make sure a valid position is passed
tree->root = remove_rec(tree->root, pos, tree->dim, tree->destr, 0);
}
return(0);
}
kdnode* remove_rec(kdnode* node, const double* pos, int dim, void (*destr)(void*), int depth)
{
if(node == NULL)
{
return(NULL);
}
int curdim = depth % dim;
if(same_pos(node->pos, pos, dim))
{
// we found the droid we're looking for
if(node->right)
{
// find the minimum in the right subtree
kdnode* node_min = find_min(node->right, curdim, dim);
if(node_min)
{
copy_node_data(node_min, node, dim);
node->right = remove_rec(node->right, node_min->pos, dim, destr, depth + 1);
}
}
else if(node->left)
{
// find the minimum in the left subtree
kdnode* node_min = find_min(node->left, curdim, dim);
if(node_min)
{
copy_node_data(node_min, node, dim);
node->left = remove_rec(node->left, node_min->pos, dim, destr, depth + 1);
}
}
else
{
// no subtrees -> delete the found node
clear_rec(node, destr);
return(NULL);
}
return node; // return the newly filled node to the recursion step one "above"
}
else
{
// points are not the same, look further
if(pos[curdim] < node->pos[curdim])
{
// position we're looking for is smaller -> go left
node->left = remove_rec(node->left, pos, dim, destr, depth + 1);
}
else
{
// go right, position we're looking for is greater
node->right = remove_rec(node->right, pos, dim, destr, depth + 1);
}
return node;
}
}
void copy_node_data(const kdnode* src, kdnode* dst, int dim)
{
if(src && dst)
{
int nNumBytes = dim * sizeof(double);
memcpy(dst->pos, src->pos, nNumBytes);
if(dst->data != NULL)
{
free(dst->data);
dst->data = malloc(src->databytes);
}
memcpy(dst->data, src->data, src->databytes);
dst->databytes = src->databytes;
}
}
int same_pos(const double* pos1, const double* pos2, int dim)
{
for (int i = 0; i < dim; ++i)
{
if(pos1[i] != pos2[i])
{
return 0; // false
}
}
return 1; // true
}
kdnode* find_min(kdnode* node, int dir, int numdim)
{
return find_min_rec(node, dir, 0, numdim);
}
kdnode* find_min_rec(kdnode* node, int dir, int depth, int numdim)
{
if(!node)
{
return NULL;
}
if(node->left == NULL && node->right == NULL)
{
return node; // is leaf node
}
int curdim = depth % numdim;
if(curdim == numdim)
{
if(node->left == NULL)
{
// no smaller node in tree
return node;
}
else
{
// left subtree is populated -> we need to go deeper
return find_min_rec(node->left, node->dir, depth + 1, numdim);;
}
}
// we have to search both subtrees and find the smallest value compared to the current node
return min_node(node, find_min_rec(node->left, node->dir, depth + 1, numdim),
find_min_rec(node->right, node->dir, depth + 1, numdim), node->dir);
}
kdnode* min_node(kdnode* a, kdnode* left, kdnode* right, int dir)
{
if(a == NULL)
{
// node a is the only one that can't be NULL!
fprintf(stderr, "Error: invalid node passed! \n");
exit(EFAULT);
}
kdnode* result = a;
if(left != NULL)
{
if(left->pos[dir] < result->pos[dir])
{
result = left;
}
}
if(right != NULL)
{
if(right->pos[dir] < result->pos[dir])
{
result = right;
}
}
return result;
}
original.dot看起来像这样 removed.dot像那样。从昨天开始我就一直在调试这个问题,我感觉我在这里遗漏了一些非常明显的东西......预先感谢任何愿意提供帮助的人:)
最佳答案
您正在创建 40
个元素
int numel = 20;
int nNumDim = 2;
double* pos = (double*) malloc(nNumDim * numel * sizeof(double)); // Don't cast
但仅删除 38
int toRemove = 19;
for (int ii = 0; ii < toRemove; ii++)
{
dRemovePos[0] = pos[nNumDim * ii];
dRemovePos[1] = pos[nNumDim * ii + 1];
kd_remove(TreeRoot, dRemovePos);
}
在最后一次迭代中:
pos[nNumDim * ii];
= pos[2 * 18];
= pos[36];
pos[nNumDim * ii + 1];
= pos[2 * 18 + 1];
= pos[37];
pos[38]
和 pos[39]
仍然存在。
更改为int toRemove = 20;
。
由于平面数组,您的代码被混淆了,为什么不声明一些类型,例如
struct data {
double el1;
double el2;
};
或
typedef double data[2];
然后
data *value = malloc(numel * sizeof(*value));
关于c - 从 kdtree 中删除时出现重复节点,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/38893517/
到 N 点的列表 [(x_1,y_1), (x_2,y_2), ... ] 我试图找到每个点的最近邻居基于距离的点。我的数据集太大,无法使用蛮力方法,因此 KDtree 似乎是最好的。 我看到 skl
我尝试创建一个 WGS84 坐标的 KD 树并在一定半径内找到邻居 from sklearn.neighbors.dist_metrics import DistanceMetric from skl
以下类从 g++ 返回以下错误: tree.h:9: error: expected unqualified-id before 'template' 如果第 5 行没有注释,同样的错误发生,只在第
我有一个经纬度坐标数组,我正在尝试使用 KDTree 和 scipy 的 query_ball_point 返回指定纬度和经度 1 英里半径范围内的所有数据点。 问题是 query_ball_poin
给定一个样本屏蔽数组和一个 KDTree 查询, In [6]: data = np.arange(9).reshape((3,3))*1.0 In [7]: madata = np.ma.array
我正在用 java 实现 KdTree。我已经完成了程序的大部分其余部分,但我似乎无法让我的最近邻搜索算法正常工作。无论如何,它总是返回根节点的值。这是我的代码: public Point2D n
我一直在尝试从头开始实现 KdTree。成功实现添加、查找最近邻居和在范围方法中查找节点后,我现在被困在删除节点上。 维基百科上描述的方法含糊不清,而且毫无用处。相反,我使用 these slides
我知道如何构建 kd 树。但是我面临的问题是如何使用 KD 树查找最近邻居。我在 google 上搜索过,但无法找到查找最近邻居的代码,尽管给出了算法。但由于语言的原因,我在将该算法转换为代码时遇到了
我正在编写一个算法,要求我搜索点的最近邻居。我从这篇文章( Using Google's C KD Tree Library )中找到了 kdtree 库,但它没有从树中删除单个节点的功能。所以我开始
我正在处理一个核心转储问题,但我找不到解决方案。任何形式的帮助将不胜感激,因为我已经绝望了。 我假设在第二次执行构建函数时出现错误,但我不知道有什么大惊小怪的。 http://pastebin.com
例如,看起来我无法在 sklearn KDTree 中使用此相似性度量,但我需要,因为我正在使用测量单词向量相似性。对于这种情况,快速鲁棒定制算法是什么?我了解本地敏感度哈希,但它应该经过大量调整和测
我正在尝试在我的类(class)中使用 nanoflann 的 kdtree 结构。我有以下 kdtree 适配器的 header ,它只是 this 的一个小(可能是坏的)变体。 ,应该允许以 3d
我使用以下示例: from scipy import spatial x, y = np.mgrid[0:5, 2:8] tree = spatial.KDTree(list(zip(x.ravel(
给定点列表,如何在 KDTree 中获取它们的索引? from scipy import spatial import numpy as np #some data x, y = np.mgrid[0
我正在研究 scipy 库中的 KDTree 实现,发现自己对这行 https://github.com/scipy/scipy/blob/master/scipy/spatial/kdtree.py
所以,我正在实现 KD-Tree进行最近邻搜索。我已经开始构建树部分,但我不认为我完全理解搜索部分。 关于遍历树寻找邻居,维基百科文章是这样说的: Starting with the root nod
所以对于一个项目,我已经实现了一个平衡的 KD 树和 k 最近邻搜索算法,作为标准树的“传统”方式——即我有如下的节点结构 struct KDNode{ int level; //used f
我有一个非常大的 ndarray A 和一个排序的点列表 k(一个小列表,大约 30 个点)。 对于 A 的每个元素,我想确定点列表 k 中最近的元素以及索引。所以像这样: >>> A = np.as
我有一个表示道路网络的图形数据结构(节点是道路中的点/交叉点,边缘是道路)。 Node 对象具有与其关联的纬度和经度。 我正在使用 Accord 的 KDTree 类查找给定 GPS 坐标的附近节点。
我目前正在寻找一种方法来构建几个 kd 树以快速查询一些 n 维数据。但是,我在使用 scipy KD 树算法时遇到了一些问题 我的数据由 id -> {data: somedata, coordin
我是一名优秀的程序员,十分优秀!