scipy.spatial.KDTree.

query_ball_tree#

KDTree.query_ball_tree(other, r, p=2.0, eps=0)[source]#

selfother 之间查找所有距离小于等于 r 的点对。

参数::
otherKDTree 实例

包含要搜索的点的树。

r浮点数

最大距离,必须为正数。

p浮点数,可选

使用哪个闵可夫斯基范数。 p 必须满足条件 1 <= p <= infinity

eps浮点数,可选

近似搜索。如果树的最近点距离大于 r/(1+eps),则不会探索树的这些分支,如果树的最远点距离小于 r * (1+eps),则会批量添加分支。 eps 必须是非负数。

返回::
results列表的列表

对于此树的每个元素 self.data[i]results[i] 是其在 other.data 中的邻居索引的列表。

示例

您可以在两个 kd 树之间搜索所有距离在特定距离内的点对

>>> import matplotlib.pyplot as plt
>>> import numpy as np
>>> from scipy.spatial import KDTree
>>> rng = np.random.default_rng()
>>> points1 = rng.random((15, 2))
>>> points2 = rng.random((15, 2))
>>> plt.figure(figsize=(6, 6))
>>> plt.plot(points1[:, 0], points1[:, 1], "xk", markersize=14)
>>> plt.plot(points2[:, 0], points2[:, 1], "og", markersize=14)
>>> kd_tree1 = KDTree(points1)
>>> kd_tree2 = KDTree(points2)
>>> indexes = kd_tree1.query_ball_tree(kd_tree2, r=0.2)
>>> for i in range(len(indexes)):
...     for j in indexes[i]:
...         plt.plot([points1[i, 0], points2[j, 0]],
...             [points1[i, 1], points2[j, 1]], "-r")
>>> plt.show()
../../_images/scipy-spatial-KDTree-query_ball_tree-1.png