r - 如何在纬度和经度上从 R 中的 PostGis 转换 Geom 对象？

Question

我有一个具有“位置”属性的数据库，它被保存为几何对象，如下所示:0101000020E6100000000000603D1D5EC0000000A06D424740

有没有办法在 R 中从中提取坐标？这些值作为字符串加载到 R 中。它还给出了警告:

Warning message:
In postgresqlExecStatement(conn, statement, ...) :
  RS-DBI driver warning: (unrecognized PostgreSQL field type geometry (id:18832) in column 2)

感谢任何帮助。

Answer 1

您是否尝试过使用 rgeos 函数？一种选择是直接从数据库中读取 WKT(使用 ST_AsText )...

SELECT ST_AsText('0101000020E6100000000000603D1D5EC0000000A06D424740');
                st_astext                 
------------------------------------------
 POINT(-120.456871032715 46.518970489502)
(1 Zeile)

在 R 中，您可以使用 readWKT 提取坐标对:

> readWKT("POINT(-120.456871032715 46.518970489502)")
SpatialPoints:
          x        y
1 -120.4569 46.51897
Coordinate Reference System (CRS) arguments: NA 

另一种选择是使用 ST_X 和 ST_Y 直接从数据库中读取 x,y:

 SELECT ST_X('0101000020E6100000000000603D1D5EC0000000A06D424740'),
        ST_Y('0101000020E6100000000000603D1D5EC0000000A06D424740');
       st_x        |      st_y       
-------------------+-----------------
 -120.456871032715 | 46.518970489502
(1 Zeile)