Question

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

我正在尝试关注Daniel Azuma的RubyConf2012地理空间分析截屏视频。以下代码始终返回nil：

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

这是我的Region类：

 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时，通常意味着其他东西。

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

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

之内

属于Polygon的点

1 个答案: