C++怎么实现一个KD树_C++高维空间近邻搜索数据结构-C++-PHP中文网

C++怎么实现一个KD树_C++高维空间近邻搜索数据结构

尼克

发布： 2025-11-26 15:50:02

原创

795人浏览过

实现KD树需递归划分高维空间，C++中用模板类定义节点结构，包含坐标、分割维度和子树指针；建树时按轮转维度选中位数分割，确保平衡，利用std::nth_element优化至平均O(n)；搜索时递归下降并回溯剪枝，通过比较查询点与分割面距离判断是否遍历兄弟子树，使用欧氏距离平方避免开方，适用于低维场景，高维可改用Ball Tree等近似方法。

c++怎么实现一个kd树_c++高维空间近邻搜索数据结构

实现KD树的关键在于递归划分高维空间，每次选择一个维度进行分割，使得数据在该维度上左右分布。C++中通过结构体或类来组织节点信息，结合递归建树和剪枝搜索策略，可以高效完成近邻查找。

定义KD树节点结构

每个节点需要存储当前点的坐标、分割维度、以及左右子树指针。坐标的维度可以在编译时用模板确定，也可以运行时动态处理。

使用数组或vector保存多维坐标值
记录当前节点用于划分的维度 axis
左右子树指针 left 和 right

示例代码：

template <size_t K>
struct KDNode {
    std::array<float, K> point;
    int axis;
    KDNode* left;
    KDNode* right;
<pre class="brush:php;toolbar:false;"><pre class="brush:php;toolbar:false;">KDNode(const std::array<float, K>& p) : point(p), axis(0), left(nullptr), right(nullptr) {}

登录后复制

};

构建KD树

建树过程是递归的。每层选择一个维度，按该维度对数据排序后取中位数作为分割点，确保树尽量平衡。

立即学习“C++免费学习笔记（深入）”；

选择划分维度：轮转法（如第d层用d%K维）或方差最大维
找到当前数据在选定维度上的中位数元素
以中位数为根节点，左右部分递归建左子树和右子树

关键操作是快速找到中位数——可用std::nth_element优化到平均O(n)。

INFINITE ALBUM

面向游戏玩家的生成式AI音乐

144

查看详情

template <size_t K>
KDNode<K>* buildTree(std::vector<std::array<float, K>>& points, int depth = 0) {
    if (points.empty()) return nullptr;
<pre class="brush:php;toolbar:false;"><pre class="brush:php;toolbar:false;">int axis = depth % K;
auto mid = points.begin() + points.size() / 2;
std::nth_element(points.begin(), mid, points.end(),
    [axis](const auto& a, const auto& b) { return a[axis] < b[axis]; });

KDNode<K>* node = new KDNode<K>(*mid);
node->axis = axis;

std::vector<std::array<float, K>> leftPoints(points.begin(), mid);
std::vector<std::array<float, K>> rightPoints(mid + 1, points.end());

node->left = buildTree<K>(leftPoints, depth + 1);
node->right = buildTree<K>(rightPoints, depth + 1);

return node;

登录后复制

}

最近邻搜索

从根节点开始，根据查询点与分割面的关系决定优先走哪边，再判断另一边是否有更近的可能。

递归下降到叶子节点，记录当前最短距离
回溯过程中检查兄弟子树是否可能包含更近点（通过距离分割面的距离判断）
维护一个最小距离变量，用于剪枝

距离计算通常用欧氏距离平方避免开方开销。

float distance(const std::array<float, K>& a, const std::array<float, K>& b) {
    float dist = 0;
    for (int i = 0; i < K; ++i)
        dist += (a[i] - b[i]) * (a[i] - b[i]);
    return dist;
}
<p>void nearestNeighbor(KDNode<K><em> node, const std::array<float, K>& query,
KDNode<K></em>& best, float& bestDist, int depth = 0) {
if (!node) return;</p><pre class="brush:php;toolbar:false;"><pre class="brush:php;toolbar:false;">float dist = distance(query, node->point);
if (!best || dist < bestDist) {
    best = node;
    bestDist = dist;
}

int axis = depth % K;
KDNode<K>* nearSide = query[axis] < node->point[axis] ? node->left : node->right;
KDNode<K>* farSide = (nearSide == node->left) ? node->right : node->left;

nearestNeighbor(nearSide, query, best, bestDist, depth + 1);

float planeDist = (query[axis] - node->point[axis]) * (query[axis] - node->point[axis]);
if (planeDist < bestDist) {
    nearestNeighbor(farSide, query, best, bestDist, depth + 1);
}

登录后复制

}