Ruby Rails Postgis - 查找多边形中的所有点

Question

我想要一些关于构建 sql 查询的帮助，以便在带有 activerecord-postgis-adapter 的 Rails 中使用。我已经做了很多阅读，但现在有点卡住了，任何帮助将不胜感激。

我有两个模型事件和区域：

事件有一个“几何”列，其类型为点

class Event < ActiveRecord::Base
  self.rgeo_factory_generator = RGeo::Geos.factory_generator    
end

t.spatial  "geometry", :limit => {:srid=>4326, :type=>"polygon", :geographic=>true}

区域有一个“几何”列，其类型为多边形

class Area < ActiveRecord::Base 
  self.rgeo_factory_generator = RGeo::Geos.factory_generator
end

t.spatial  "geometry", :limit => {:srid=>4326, :type=>"point", :geographic=>true}

我可以在谷歌地图上创建和绘制事件和区域，并创建区域通过单击地图并将其保存到数据库。

我希望能够执行以下 2 个查询：

1. @area.events - 显示某个区域中的所有事件
2. @event.areas - 显示单个事件所在的所有区域

我知道我在这里可能问得有点多，但是任何帮助将不胜感激 非常

感谢

Answer 1

这是执行此操作的快速方法。这些将简单地返回 ActiveRecord 对象的数组。

class Area
  def events
    Event.joins("INNER JOIN areas ON areas.id=#{id} AND st_contains(areas.geometry, events.geometry)").all
  end
end

class Event
  def areas
    Area.joins("INNER JOIN events ON events.id=#{id} AND st_contains(areas.geometry, events.geometry)").all
  end
end

您可能应该记住（缓存结果），这样您就不会在每次调用该方法时都查询数据库。这应该很简单；我把它作为练习留给读者。

可能变得复杂，并将其包装在真正的 Rails 关联代理中（这样您就可以获得所有 Rails 关联的好处）。不过我还没有研究过这个。无论如何，它都不是标准的 Rails 关联，因为您不存储 ID。

第十二点是正确的：您应该为两个表创建空间索引。 Activerecord-postgis-adapter 应该可以让您在迁移中轻松完成这些任务。

change_table :events do |t|
  t.index :geometry, :spatial => true
end

change_table :areas do |t|
  t.index :geometry, :spatial => true
end

如果您在安装 postgis 时遇到问题，我最近写了很多关于这方面的博客文章。查看 http://www.daniel-azuma.com/blog/archives /category/tech/georarails。我也是 rgeo 和 activerecord-postgis-adapter 的作者，所以如果您遇到困难，我很乐意为您提供帮助。

Answer 2

这个答案对您来说将是一项正在进行的工作。我对 ruby​​ on Rails 不太了解，但我应该能够帮助您完成数据库部分。

您有两个表，包含多边形的区域和将事件保存为单个点的事件（如果事件也是一个区域并且您试图挑选出重叠区域...如果事件是单个点，则情况会更复杂点这有效）。

Select *
from area a inner join event e on 1=1

这将创建加入每个事件的每个区域的列表...如果您有 500 个事件和 20 个区域，则此查询将返回 10'000 行。现在您想要对此进行过滤，以便仅过滤其所加入的区域内的事件。我们可以将 st_contains 用作 st_contains(polygon,point)：

where st_contains(a.polygon,e.point) = 't'

如果运行此命令，它应该为您提供区域内所有事件的 a.,e.。现在只要数你想数的就可以了。

select a.id, count(1)
from area a inner join event e on 1=1
where st_contains(a.polygon,e.point) = 't'
group by 1

这将为您提供所有区域的列表（按 ID）以及其中的事件计数。用 e.id 替换 a.id 将给出事件 id 的列表以及它们所在的数字区域。

不幸的是，我不知道如何在 Ruby 中表达这些查询，但是您需要的数据库概念都在这里。 出于速度考虑，

您应该查看 Postgres 的 GIStree 索引...索引多边形的性能呈指数级提高。

编辑：

PostGIS 是 Postgres 附带的一个 contrib 文件，但标准安装中不存在...您需要找到此 contrib 文件。这些将在您的数据库中安装一系列 GIS 功能，包括 ST_Contains。 （函数驻留在数据库中，因此请确保在您正在使用的数据库中安装函数）

PostGIS contrib 文件安装的第二件事是 template_postGIS 数据库，这是几何数据类型所需的（geom 作为数据类型将不存在）直到安装完毕）。