问题:Postgis纬度/经度读反了

我有一个类,Geocoordinateable,它为模型提供了一些处理地理坐标的有用方法。其中包括接受数组的“geocoordinatesu003d”。

https://gist.github.com/mboyle/58dd3add830bbdeef316

您可以调用“地理坐标”来获取数组中的坐标:

irb(main):056:0> b.geocoordinates
=> [-118.25, 34.197]

当我像这样分配坐标时:

irb(main):058:0> b.geocoordinates = [34.197, -118.25]
   (1.8ms)  
 UPDATE users
 SET geocoordinates = 'SRID=4326;POINT(' || 34.197 || ' ' || -118.25 || ')'
 WHERE id = 347708

   (0.7ms)  
 UPDATE activities
 SET geocoordinates = users.geocoordinates
 FROM users
 WHERE activities.user_id = users.id
 AND users.id = 347708
 AND NOT postgis.ST_Equals(activities.geocoordinates, users.geocoordinates)

=> [34.197, -118.25]

事情似乎可以保存并正确显示。但是,如果我再次查询模型,然后读取地理坐标,则纬度/经度会反转:

irb(main):059:0> b = User.find(b.id)
irb(main):060:0> b.geocoordinates
   (0.7ms)  
 SELECT postgis.ST_Y(geocoordinates) AS latitude, postgis.ST_X(geocoordinates) AS longitude
 FROM users
 WHERE id = 347708

=> [-118.25, 34.197]

我一生都无法理解为什么会这样。有人有线索吗?

解答

在 PostGIS 中,坐标按 (X, Y) 顺序设置,或 (longitude, latitude) 用于地理坐标对。所以你应该这样做:

UPDATE users
SET geocoordinates = 'SRID=4326;POINT(' || -118.25 || ' ' || 34.197 || ')'
WHERE id = 347708

然后可能也是:

irb(main):058:0> b.geocoordinates = [-118.25, 34.197]

请注意,PostGIS 以OGC指定的众所周知的文本 (WKT) 格式提取文本数据。该标准规定了 (X,Y) 或 (lon,lat) 顺序。所以这不是一个选择,它需要使 PostGIS 能够摄取由其他 OGC 兼容软件生成的 WKT 文件。摄取 WKT 后,PostGIS 将数据以其内部格式(实际上是稍微修改的 OGC Well-Known Binary (WKB) 格式)存储在 PostGIS 定义的geometry或geography数据类型的列中,然后可以将数据存储在以最终用户所需的任何形式进行分析和呈现。