2

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

我正在尝试关注 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
4

1 回答 1

2

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

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

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

于 2013-08-31T22:29:23.360 回答