编辑:我已经添加了一个演示SQL Server 2017的登录详细信息。
我正在使用FreeTDS驱动程序(它支持Linux中的Windows身份验证,因此比Microsoft odbc驱动程序更受欢迎)从R查询SQL Server。
我无法使用st_read
解释几何列。
连接是这样初始化的:
library(odbc)
library(DBI)
library(sf)
username <- 'SO-user'
sql_server_ip <- '35.214.169.110'
password <- 'SQLaskingfortrouble?!'
con <- dbConnect(odbc(),
Driver = "FreeTDS",
Server = sql_server_ip,
Database = "stackoverflow-example",
UID = username,
PWD = password,
Port = 1433)
# I have also tried Driver = "ODBC Driver 17 for SQL Server"
以下是查询结果:
data <- dbGetQuery(con, 'SELECT TOP 2
objectid, geom,
geom.STGeometryType() geom_type,
geom.STSrid STSrid
FROM spatialdata')
objectid geom geom_type stsRID
864 blob[416 B] LineString 25832
865 blob[416 B] LineString 25832
class(data$geom)
"blob" "vctrs_list_of" "vctrs_vctr"
原始/二进制几何图形看起来像十六进制:
rawToHex(data$geom[1])
e86400000104180000004003780b38331f416b4df3ea9ecf5741a04c150c38331f4159ca32f09ecf574160ba...
我正在尝试使用 st_read
,以便从中获取有用的 sf 数据框架,但是这取决于我传递给函数的选项,会给我带来一系列错误消息:
st_read(con, geometry_column = "geom",
query = "SELECT TOP 10 objectid, geom FROM \"spatialdata\"")
Error: nanodbc/nanodbc.cpp:1617: 00000: [FreeTDS][SQL Server]Invalid object name 'spatial_ref_sys'. [FreeTDS][SQL Server]Statement(s) could not be prepared.
<SQL> 'select srtext from spatial_ref_sys where srid = 70647808'
In addition: Warning message:
In CPL_crs_from_input(x) :
GDAL Error 6: EPSG PCS/GCS code 70647808 not found in EPSG support files. Is this a valid EPSG coordinate system?
看起来EPSG(25832)没有被正确识别。
在我的情况下,st_read的正确构造是什么?
额外问题:如何高效地将边界框附加到st_read中?