我使用scipy.spatial.cKDTree的query_ball_point方法查找一个点周围给定半径R内的所有对象。 我想要定期增加R以包含越来越多的点并保存输出。目前我这样做:
for i in range(0,lmax+sphererads,sphererads):
gals = kdt.query_ball_point(origin, i)
DF[n]=pd.DataFrame(gals)
n=n+1
然而,这样我的代码会计算每次迭代中的所有距离,包括已经计算过的距离。有没有办法为查询指定最小半径Rmin?
提前致谢!