2

我对 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

目前边界来自谷歌反向地理编码查找。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 问题和其他证据表明这些操作不支持基于球形工厂的对象

我是在想;

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

编辑:

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

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...

所以guess contains在地理上是行不通的。我四处挖掘,发现了 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 是边界的中心点。我很困惑,肯定我在做一些很愚蠢的事情。任何帮助将不胜感激

4

2 回答 2

2

我推荐使用ST_DWithin,它很好地支持 PostGIS 的geography类型。对于半径参数,您可以使用 0 或 10(即,如果您的数据具有 10 m 的精度)。

没有任何计划使ST_ContainsST_Within可用于geography类型。

于 2016-09-25T23:28:19.637 回答
1

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

我会羞愧地低下头:)

感谢一百万倾斜和迈克。我不确定我是否会找到 ST_DWithin,而专业的 postgist 调试技巧无疑为我节省了很多小时的挫败感

于 2016-09-26T09:52:03.853 回答