Python SciPy KDTree.query_ball_tree用法及代碼示例

本文簡要介紹 python 語言中 scipy.spatial.KDTree.query_ball_tree 的用法。

用法:

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

找到 self 和 other 之間距離最大為 r 的所有點對。

參數 ：：

other： KDTree 實例

包含要搜索的點的樹。

r： 浮點數

最大距離，必須為正。

p： 浮點數，可選

使用哪個 Minkowski 範數。p必須滿足條件1 <= p <= infinity.

eps： 浮點數，可選

近似搜索。如果最近的點比樹的分支更遠，則不會探索樹的分支r/(1+eps), 如果分支的最遠點接近於r * (1+eps).每股收益必須是非負的。

返回 ：：

results： 列表列表

對於這棵樹的每個元素 self.data[i]，results[i] 是其在 other.data 中的鄰居的索引列表。

例子：

您可以搜索距離內兩個kd-trees之間的所有點對：

>>> 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()