1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
Sys.setenv(TZ="UTC")
suppressPackageStartupMessages(library(sf))
if ("GPKG" %in% st_drivers()$name) {
tst = st_read(system.file("gpkg/nc.gpkg", package="sf"), "nc.gpkg", crs = 4267, quiet = TRUE)
tst = st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE) # default layer name
}
tst = st_read(system.file("shape/nc.shp", package="sf"), quiet = TRUE) # no layer specified
# data from https://github.com/edzer/sfr/issues/6
tst = st_read(system.file("shape/storms_xyz.shp", package="sf"), quiet = TRUE)
class(st_geometry(tst))
class(st_geometry(tst)[[1]])
tst = st_read(system.file("shape/storms_xyzm.shp", package="sf"), quiet = TRUE)
class(st_geometry(tst))
class(st_geometry(tst)[[1]])
tst = st_read(system.file("shape/storms_xyz_feature.shp", package="sf"), quiet = TRUE)
class(st_geometry(tst))
class(st_geometry(tst)[[1]])
tst = st_read(system.file("shape/storms_xyzm_feature.shp", package="sf"), quiet = TRUE)
class(st_geometry(tst))
class(st_geometry(tst)[[1]])
if ("GPKG" %in% st_drivers()$name) { # read Int64
print(st_read(system.file("gpkg/tl.gpkg", package="sf"), quiet = TRUE)$AWATER)
print(st_read(system.file("gpkg/tl.gpkg", package="sf"), quiet = TRUE, int64_as_string = TRUE)$AWATER)
}
# see https://github.com/edzer/sfr/issues/45 :
if ("OSM" %in% st_drivers()$name && Sys.info()['sysname'] != "Darwin") {
osm = system.file("osm/overpass.osm", package="sf")
print(st_layers(osm))
suppressWarnings(print(st_layers(osm, do_count = TRUE)))
suppressWarnings(print(st_read(osm, "multipolygons", quiet = TRUE)))
}
# layer opening option:
st_read(system.file("shape/nc.shp", package="sf"), quiet = TRUE,
options = c("ADJUST_TYPE=YES", "ENCODING=CPL_ENC_UTF8"))
x <- st_sf(a = 1:2, geom = st_sfc(st_point(0:1), st_multipoint(matrix(1:4,2,2))))
try(st_layers("foo")) # cannot open datasource
try(st_read("foo")) # cannot open datasource
try(st_read("x.gpkg", "xyz")) # cannot open layer
try(st_write(c("foo", "bar")))
try(st_write(x, c("foo", "bar")))
try(st_write(x, "foo", driver = "foo"))
if (Sys.getenv("USER") == "travis") {
try(st_write(x, "/x", driver = "ESRI Shapefile"))
st_write(x, "xxx.gpkg")
st_write(x, "xxx.gpkg", append = TRUE, quiet = FALSE) # appends to layer
y <- st_sf(b = 1:2, geom = st_sfc(st_point(0:1), st_multipoint(matrix(1:4,2,2))))
try(st_write(y, "xxx.gpkg"))
}
geom = st_sfc(st_point(0:1), st_multipoint(matrix(1:4,2,2)))
st_write(geom, "geom.gpkg")
st_write(geom, "geom1.gpkg", layer = "foo")
df <- data.frame(
a = c(0, 1, NA, -Inf, Inf),
b = c("a", "b", NA, "c", ""),
c = c(as.Date("2001-01-01"),
as.Date(c(NA, -99, 0, 1), origin = "1970-01-01")),
d = c(as.POSIXct("2001-01-01"),
as.POSIXct(c(NA, -99, 0, 1), origin = "1970-01-01")),
x = 1:5,
y = 1:5)
x = st_as_sf(df, coords = c("x", "y"))
if ("GPKG" %in% st_drivers()$name) {
st_write(x, "x2.gpkg", quiet = TRUE)
y = st_read("x2.gpkg", quiet = TRUE)
print(y)
}
if ("SQLite" %in% st_drivers()$name && require(RSQLite)) {
db = system.file("sqlite/meuse.sqlite", package = "sf")
dbcon <- dbConnect(dbDriver("SQLite"), db)
m = dbReadTable(dbcon, "meuse.sqlite")
m$GEOMETRY = st_as_sfc(m$GEOMETRY, spatialite = FALSE) # ISO wkb
print(st_sf(m), n = 3)
# or:
(s = st_read(dbcon, "meuse.sqlite"))[1:3,]
dbDisconnect(dbcon)
db = system.file("sqlite/nc.sqlite", package = "sf")
dbcon <- dbConnect(dbDriver("SQLite"), db)
m = dbReadTable(dbcon, "nc.sqlite")
m$GEOMETRY = st_as_sfc(m$GEOMETRY, spatialite = FALSE) # ISO wkb
print(st_sf(m), n = 3)
dbDisconnect(dbcon)
db = system.file("sqlite/b.sqlite", package = "sf") # has an INT8 field
b = st_read(db, quiet = TRUE)
print(b)
b = st_read(db, int64_as_string = TRUE, quiet = TRUE)
print(b)
if (.Platform$endian == "little") {
db = system.file("sqlite/test3.sqlite", package = "sf")
dbcon <- dbConnect(dbDriver("SQLite"), db)
m = dbReadTable(dbcon, "HighWays")
m$Geometry = st_as_sfc(m$Geometry, spatialite = TRUE) # spatialite wkb
print(st_sf(m), n = 3)
m = dbReadTable(dbcon, "Towns")
m$Geometry = st_as_sfc(m$Geometry, spatialite = TRUE) # spatialite wkb
print(st_sf(m), n = 3)
m = dbReadTable(dbcon, "Regions")
m$Geometry = st_as_sfc(m$Geometry, spatialite = TRUE) # spatialite wkb
print(st_sf(m), n = 3)
}
}
csv = system.file("csv/pt.csv", package = "sf")
identical(st_read(quiet = TRUE, csv, options = "AUTODETECT_TYPE=Yes")$Int64[3], NA_real_)
identical(st_read(quiet = TRUE, csv, int64_as_string = TRUE, stringsAsFactors = FALSE,
options = "AUTODETECT_TYPE=Yes")$Int64[3], NA_character_)
identical(st_read(quiet = TRUE, csv, options = "AUTODETECT_TYPE=Yes")$Int32[3], NA_integer_)
if ("GML" %in% st_drivers()$name) {
gml = system.file("gml/fmi_test.gml", package = "sf")
print(dim(st_read(gml, quiet = TRUE)))
gml = system.file("gml/20170930_OB_530964_UKSH.xml.gz", package = "sf")
print(dim(st_read(gml, layer = "Parcely", quiet = TRUE)))
print(dim(st_read(gml, layer = "Parcely", int64_as_string=TRUE, quiet = TRUE)))
}
# logical:
if ("GPKG" %in% st_drivers()$name) {
tst = st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE) # default layer name
tst$bool = tst$NWBIR79 > 800 # logical
tst$bool[1:3] = NA
st_write(tst, "tst__.gpkg")
tst2 = st_read("tst__.gpkg")
print(identical(tst$bool, tst2$bool))
}
# spatial filter:
nc = read_sf(system.file("gpkg/nc.gpkg", package="sf"))
wkt = st_as_text(st_geometry(nc[1,]))
wkt
nc_filtered = read_sf(system.file("gpkg/nc.gpkg", package="sf"), wkt_filter = wkt)
try(read_sf(system.file("gpkg/nc.gpkg", package="sf"), wkt_filter = "wrong"))
nc_filtered$NAME
|
