我的问题几乎与此类似。但就我而言,多边形不一定相互接触/重叠。它们遍布整个空间。
我有一大堆这样的多边形。同样,我有很多点。我目前正在运行一个 RoR 模块,它一次需要 1 个点,并一次检查相对于 1 个多边形的交集。数据库是 PostGIS。性能相当缓慢。
有没有更快或最佳的方法来做到这一点?
我的问题几乎与此类似。但就我而言,多边形不一定相互接触/重叠。它们遍布整个空间。
我有一大堆这样的多边形。同样,我有很多点。我目前正在运行一个 RoR 模块,它一次需要 1 个点,并一次检查相对于 1 个多边形的交集。数据库是 PostGIS。性能相当缓慢。
有没有更快或最佳的方法来做到这一点?
可以作为一个选择语句完成,但为了性能....查看多边形上的要点索引。为简单起见,假设我有一个包含多边形字段(geom 数据类型)和点字段(geom 数据类型)的表。如果您正在制作多边形列表中的点列表,请进行交叉连接,以便比较每个多边形和每个点。
select *
from t1 inner join t2 on 1=1
where st_contains(t1.poly,t2.point) = 't'
(修改为包含表连接示例。我使用的是交叉连接,这意味着每个多边形都将连接到每个点并进行比较。如果我们谈论的是一个大型记录集,那么让那些 GIS 树索引继续运行)
我目前这样做是为了在几百个多边形内定位几百万个点。如果您有重叠的多边形,这将为位于 2 个或更多多边形中的每个点返回多行。
您的积分存储为的数据类型可能会失败。如果他们在 geom 领域,它会流动的很好。如果您使用文本值,则需要使用 st.geomfromtext 语句将字符转换为点。这看起来更像:
st_contains(poly, st_geomfromtext('POINT('||lon||' ' ||lat ||')')) = 't'
我使用了一个纬度/经度示例...这里唯一要注意的是 geomfromtext 要求您使用 || 创建点 从您的字段创建字符串。如果您需要有关 st_geomfromtext 概念的帮助,请告诉我。