您最初倾向于围绕每个点绘制一个圆并计算该圆中其他点的数量是一个好的方法,并且如unutbu所述,KDTree将是一种快速方法解决这个问题。
这可以很容易地用PySAL来完成,PySAL使用scipy的kdtree。
import pysal
import numpy
pts = numpy.random.random((100,2)) #generate some random points
radius = 0.2 #pick an arbitrary radius
#Build a Spatial Weights Matrix
W = pysal.threshold_continuousW_from_array(pts, threshold=radius)
# Note: if your points are in Latitude and Longitude you can increase the accuracy by
# passing the radius of earth to this function and it will use arc distances.
# W = pysal.threshold_continuousW_from_array(pts, threshold=radius, radius=pysal.cg.RADIUS_EARTH_KM)
print W.cardinalities
#{0: 10, 1: 15, ..... }
如果你的数据是在shape文件,只需用threshold_continuousW_from_shapefile取代threshold_continuousW_from_array,请参阅该文档的详细信息。