不支持多边形 'contains' 和其他几何操作

Question

我对 postgis 和 rgeo 很陌生。我怀疑我可能以错误的方式处理问题，但我有点惊讶地发现一些操作，特别是包含和内部，在基于球形的对象上是不可能的。

我有一堆地理分布的对象，我想根据邮政编码之类的东西将它们组合在一起。对于这些分组中的每一个，我都有一个边界，我想检查一个对象是否在该边界内，并检查一个边界是否在另一个边界内。我正在使用 rails，这是用于设置我的收集模型的迁移

class CreateGeoCollectionDefinition < ActiveRecord::Migration[5.0]
  def change
     create_table :geo_collection_definitions do |t|
       t.string :name
       t.string :geo_place_id
       t.string :geo_place_types, array: true, default: []
       t.st_polygon :boundary, geographic: true
       t.jsonb :boundary_json
       t.st_point :latlng, geographic: true
     end
   end
end

当前边界来自 google 反向地理编码查找。 NorthEast 和 SouthWest 边界框坐标在对象创建时传递给此方法

GEO_FACTORY = RGeo::Geographic.spherical_factory(srid: 4326)

def self.createBoundary(pointOne, pointTwo)
  point1 = GEO_FACTORY.point(pointOne['lat'], pointOne['lng'])
  point2 = GEO_FACTORY.point(pointTwo['lat'], pointTwo['lng'])
  boundingBox = RGeo::Cartesian::BoundingBox.create_from_points(point1,    point2).to_geometry
  boundingBox
end

我已经编写了一些规范来检查是否一切都按照我期望的方式运行。基于简单距离的测试均按预期通过，但用于测试边界功能的测试存在问题。我运行变成下面

 # ------------------
 # --- Caused by: ---
 # PG::UndefinedFunction:
 #   ERROR:  function st_contains(geography, geography) does not exist
 #   LINE 1: ...COUNT(*) FROM "geo_collection_definitions" WHERE (ST_Contain...
 #                                                                ^
 #   HINT:  No function matches the given name and argument types. You might need to add explicit type casts.

当我尝试运行查询时（我知道它有点傻）

 GeoCollectionDefinition.where("ST_Contains(boundary, boundary)")

或者如果我尝试直接使用 rgeo 对象

it "should be possible to test is points belong in GeoCollection.boundary" do
  factory = RGeo::Geographic.spherical_factory(srid: 4326)
  externalPoint = factory.point(EmptyGeocodeLatLag['lng'], EmptyGeocodeLatLag['lat'])
  expect(someplace_def.boundary.contains?(someplace_def.latlng)).to be_truthy
  expect(someplace_def.boundary.contains?(externalPoint)).to be_falsy
end

我明白了

  RGeo::Error::UnsupportedOperation:
   Method Geometry#contains? not defined.

四处挖掘我发现了这个 rgeo issue 和其他证据表明这些操作不支持基于球形工厂的对象

我只是想知道；

一定是这样吗
按照我在迁移中描述的方式对对象的边界和位置进行建模对我来说似乎很有意义，但我猜我错了？
我应该如何使用 postgis、rgeo 和活动记录适配器找出一堆点是否在多边形内？
是否可以检查一个多边形是否在另一个多边形内部？

编辑：

我跟进了 tilt 的建议，运行直接在我的数据库上进行了一些查询。首先，只是为了验证我的设置

SELECT PostGIS_full_version();

 NOTICE:  Function postgis_topology_scripts_installed() not found. Is topology support enabled and topology.sql installed?                                  postgis_full_version

 POSTGIS="2.1.7 r13414" GEOS="3.5.0-CAPI-1.9.0 r4084" PROJ="Rel.   4.9.2, 08 September 2015" GDAL="GDAL 1.11.5, released 2016/07/01" LIBXML="2.9.2" LIBJSON="UNKNOWN" RASTER

我也运行一些疑问；

  SELECT name  FROM geo_collection_definitions WHERE st_contains(latlng, boundary);
 ERROR:  function st_contains(geography, geography) does not exist
 LINE 1: ...ELECT name  FROM geo_collection_definitions WHERE   st_contain...

所以猜测包含在地理位置上不起作用。我四处寻找，发现 st_covers

  SELECT name  FROM geo_collection_definitions WHERE st_covers(boundary, latlng);
  name
  ------
  (0 rows)

 SELECT name  FROM geo_collection_definitions WHERE st_covers(boundary, ST_GeomFromText('POINT(12.9549709 55.5563043)', 4326));
  name
  ------
  (0 rows)

这真的很令人惊讶，因为 latlng 是边界的中心点。我很困惑，并且确定我在做一些非常愚蠢的事情。任何帮助将不胜感激

Answer 1

我推荐使用 ST_DWithin，它很好地支持 PostGIS 的 geography 类型。对于半径参数，您可以使用 0 或 10（即如果您的数据精度为 10 米）。

没有任何计划使 ST_Contains or ST_Within 可用于 geography 类型。

Answer 2

Mike 和 Tilt 的回答帮助我弄清了真相，这有点令人尴尬，但以防万一这些东西可能对其他人有所帮助...

ST_DWithin 绝对是解决我的特定问题的方法，很高兴能对此大声疾呼。我将一些规范查询更改为这样的

  collectionDefs = GeoCollectionDefinition.where("ST_DWithin(boundary, '#{someplace.latlng}', 10)")
  expect(collectionDefs.count).to eq(2)
  collectionDefs = GeoCollectionDefinition.where("ST_DWithin(boundary, 'Point(0.0 0.0)', 10)")
  expect(collectionDefs.count).to eq(0)

查询运行很好，很好，但我仍然没有得到预期的结果。这就是尴尬的地方

根据 Tilt 的建议，我打印出了多边形

POLYGON((55.4965351 12.8894595,55.6445967 12.8894595,55.6445967 13.151087,55.4965351 13.151087,55.4965351 12.8894595)) and the centre point is POINT(13.0108705 55.5790534)

并很快发现坐标错了。损坏是用我上面概述的 createBoundary 方法完成的。我正在创建边界定义点并传递 lat 值代替 lng，反之亦然。我将代码更新为以下

 def self.createBoundary(pointOne, pointTwo)
    point1 = GEO_FACTORY.point(pointOne['lng'], pointOne['lat'])
    point2 = GEO_FACTORY.point(pointTwo['lng'], pointTwo['lat'])
    boundingBox =    RGeo::Cartesian::BoundingBox.create_from_points(point1, point2).to_geometry
   boundingBox
 end

我会羞愧地低下头:)

非常感谢 Tilt 和 Mike。我不确定我是否会找到 ST_DWithin 并且专业的 postgist 调试技巧确实让我省去了很多时间的挫折

不支持多边形 'contains' 和其他几何操作

polygon 'contains' and other operations on geometries not supported

ruby-on-rails

activerecord

postgis

rgeo