1

我正在阅读 Daniel Azuma 的 Geo-Rails 博客系列，我遇到了一个我试图复制的部分，但使用了不同的 SRID。具体来说，我引用了http://blog.daniel-azuma.com/archives/69和标题为使用位置数据的部分。

Daniel Azuma 有一个locations带有latlon属性的表：

class CreateLocations < ActiveRecord::Migration
  def change
    create_table :locations do |t|
      t.string :name
      t.point :latlon, :geographic => true
      t.timestamps
    end
  end
end

相应的Location类使用spherical_factorywith SRID 4326：

class Location < ActiveRecord::Base
  set_rgeo_factory_for_column(:latlon,
    RGeo::Geographic.spherical_factory(:srid => 4326))
end

创建两个示例点，并确定它们之间的距离：

ruby-1.9.3-p0 :001 > loc = Location.create
ruby-1.9.3-p0 :003 > loc.latlon = "POINT(-122.193963 47.675086)"
ruby-1.9.3-p0 :005 > loc2 = Location.create(:name => 'Space Needle',
                        :latlon => 'POINT(-122.349341 47.620471)')
ruby-1.9.3-p0 :006 > puts "Distance is %.02f meters" %
                      loc.latlon.distance(loc2.latlon)
Distance is 13143.18 meters

我正在做类似的事情，但我得到了截然不同的结果（超过6公里的差异）。

我正在使用SRID 3785- 但我的印象是，如果我从 long/lat ( ) 到进行适当的转换，结果应该相当接近：43263785

class AddGeopointToLocations < ActiveRecord::Migration
  def change
    add_column :locations, :geopoint, :point, srid: 3785
    add_index :locations, :geopoint, spatial: true
  end
end

我的Location模型：

class Location < ActiveRecord::Base

  RGEO_FACTORY = RGeo::Geographic.simple_mercator_factory

  set_rgeo_factory_for_column(:geopoint, RGEO_FACTORY.projection_factory)

  def geopoint_geographic
    RGEO_FACTORY.unproject(self.geopoint)
  end
  def geopoint_geographic=(value)
    self.geopoint = RGEO_FACTORY.project(value)
  end

end

通过设置点parse_wkt并计算距离：

2.0.0p247 :001 > loc.geopoint_geographic = 
            ::Location::RGEO_FACTORY.parse_wkt("POINT (-122.193963 47.675086)")
2.0.0p247 :002 > loc2.geopoint_geographic = 
            ::Location::RGEO_FACTORY.parse_wkt("POINT (-122.349341 47.620471)")
2.0.0p247 :003 > loc.geopoint.distance(loc2.geopoint)
 => 19509.352351913036

通过设置点point并计算距离：

2.0.0p247 :004 > loc.geopoint_geographic = 
            ::Location::RGEO_FACTORY.point(-122.193963, 47.675086)
2.0.0p247 :005 > loc2.geopoint_geographic = 
            ::Location::RGEO_FACTORY.point(-122.349341, 47.620471)
2.0.0p247 :006 > loc.geopoint.distance(loc2.geopoint)
 => 19509.352351913036

我的结果是19.5 kilometersDaniel Azuma 的结果是13.1. 为什么差距这么大？我转换不正确吗？

数据库版本：PostgreSQL 9.3.1和POSTGIS="2.1.0 r11822" GEOS="3.3.8-CAPI-1.7.8" PROJ="Rel. 4.8.0, 6 March 2012"

4

1 回答 1

0

流行可视化 CRS 墨卡托 ( EPSG:3785 ) 的单位是米 (m)。WGS84 ( EPSG:4326 ) 的单位是度。

于 2015-01-05T01:10:36.250 回答