属于哪个多边形的点

Question

我正试图找出某个点所在的区域（多边形）。

我正在努力关注Daniel Azuma的RubyConf2012地理空间分析屏幕截图。以下代码总是返回nil:

 Region.where{st_contains(poly, Point.last.coords)}.first

这是我的区域类：

 class Region < ActiveRecord::Base
attr_accessible :poly, :name, :multi
has_many :points
GEOFACTORY = RGeo::Geographic.simple_mercator_factory
set_rgeo_factory_for_column(:poly, GEOFACTORY.projection_factory)

def self.load
path = File.join(Rails.root, 'lib/tasks/uk/districts')
factory = GEOFACTORY
RGeo::Shapefile::Reader::open(path, :factory => factory) do |file|
    file.each do |record|
        name = record['NAME']
        region = Region.where(:name => name).first ||
         record.geometry.projection.each do |poly|
         Region.create(:name => name, :poly => poly)
        end
    end
end
end

积分等级：

class Point < ActiveRecord::Base
attr_accessible :coords
belongs_to :region
def self.find_region
  region = Region.where{st_contains(poly, Point.last.coords)}.first
  puts region
end
end

Answer 1

因此，首先，您说您的查询返回nil。那不是真的。查询结果的first方法返回nil。您的查询返回了一个空数组。即它没有抛出错误，但没有产生任何结果。在rgeo中，当一个方法返回nil时，通常意味着其他东西。

您需要确保区域和点数据都来自相同的投影。这是可能的，如果不可能的话，你的点是地理坐标（朗拉特）。PostGIS只能在平面上进行空间连接。如果您将地理坐标格式化为地理类型（使用RGeo::Geographic工厂之类的东西），那么尝试为投影（平面）数据列执行带有地理形状的空间联接应该会引发错误。但是，如果您没有这样做，并且生成了具有笛卡尔或几何类型的Point实例，那么将地理点扔到方法中肯定会返回[]。毕竟，32, -81也是几何投影上的点。

也有可能最后一点根本不在区域形状内