python2代码如下,
# -*- encoding: utf-8 -*-
from rtree import index
EPS = 5e-5 #1纬度约等于1.11e5米,1经度约等于1.11e5*cos(纬度)米。
# 创建一个R树索引
idx = index.Index()
# 添加一些空间对象到索引
idx.insert(0, (0, 0, 1, 1)) # 空间对象1,换成list也是ok的
idx.insert(1, (1, 1, 2, 2)) # 空间对象2
idx.insert(2, (2, 2, 3, 3)) # 空间对象3
# 定义给定矩形范围
query_rect = (1.5, 1.5, 2.5, 2.5) # 给定的矩形范围
# 查询与给定矩形范围相交的对象
query_result = list(idx.intersection(query_rect))
print(query_result)
返回结果为,
[1, 2]