R：如何读取oracle空间数据库？

Question

我想从r中读取一个oracle spatial数据库。该数据库包含一个geometry : geometry1 : MDSYS列。SDO _纵坐标_数组(655081，618975579，6861651，4067522，0，655083，444176393，6861649，74471328，0，655094，644962858，6861639，54520924，0，6555109，121755

我可以访问数据库，但无法读取此列

library(tidyverse)
library(DBI)
library(odbc)

con <- dbConnect(odbc::odbc(), database ="Database_oracle", UID="id", PWD="pw")

这项工作：

data <- dbGetQuery(con, "SELECT ID
               FROM TABLE_1") %>% 
as_tibble()

但事实并非如此：

data <- dbGetQuery(con, "SELECT ID, geometry1
               FROM TABLE_1") %>% 
as_tibble()

ORA-01024:OCI调用中的数据类型无效

有没有办法用sf：：st_read读取坐标？我是否必须在SQL查询中使用一些“SDO_GEOMETRY（）”函数来提取坐标？

Answer 1

据我所知，sf 包仅支持 Postgres 数据库直接检索几何列。

您可以在查询中附加< code>Get_WKT()方法，将几何图形转换为“已知文本”，如下所示:

data <- dbGetQuery(con, "SELECT ID, geometry1.Get_WKT() FROM TABLE1")

# Then
st_as_sf(data, wkt = 2)  # set `wkt` parameter to column index

请注意，投影信息在这里丢失了，您需要在以后添加它。希望这有所帮助！

编辑：这里是OracleSDO_GEOMETRY方法的链接-https://docs.oracle.com/database/121/SPATL/sdo_geometry-methods.htm#SPATL498