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
|
#' @importFrom utils head object.size str tail
#' @importFrom stats aggregate na.omit rbinom runif setNames
#' @importFrom tools file_ext file_path_sans_ext
#' @importFrom methods as new slot slotNames
#' @importFrom grid convertHeight convertUnit convertWidth current.viewport linesGrob nullGrob pathGrob pointsGrob polylineGrob unit viewport
#' @import graphics
#' @importFrom grDevices dev.size rgb
#' @importFrom Rcpp evalCpp
#' @importFrom DBI dbConnect dbDisconnect dbExecute dbGetQuery dbReadTable dbSendQuery dbWriteTable
#' @importFrom units as_units drop_units make_unit_label set_units
#' @importFrom classInt classIntervals
#' @useDynLib sf
NULL
#' @importFrom magrittr %>%
#' @export
magrittr::`%>%`
setOldClass("sf")
setOldClass(c("sfc_POINT", "sfc"))
setOldClass(c("sfc_MULTIPOINT", "sfc"))
setOldClass(c("sfc_LINESTRING", "sfc"))
setOldClass(c("sfc_MULTILINESTRING", "sfc"))
setOldClass(c("sfc_POLYGON", "sfc"))
setOldClass(c("sfc_MULTIPOLYGON", "sfc"))
setOldClass(c("sfc_GEOMETRY", "sfc"))
setOldClass("sfg")
setOldClass("crs")
.sf_cache <- new.env(FALSE, parent=globalenv())
pathGrob <- NULL
.onLoad = function(libname, pkgname) {
if (getRversion() < as.numeric_version("3.6")) { # nocov start
pathGrob <<- function(..., pathId.lengths) {
grid::pathGrob(...)
}
} # nocov end
load_gdal()
assign(".sf.use_s2", Sys.getenv("_SF_USE_S2") == "true", envir=.sf_cache)
}
.onUnload = function(libname, pkgname) {
unload_gdal() # nocov
}
.onAttach = function(libname, pkgname) {
m = paste0("Linking to GEOS ", strsplit(CPL_geos_version(TRUE), "-")[[1]][1],
", GDAL ", CPL_gdal_version(), ", PROJ ", CPL_proj_version())
packageStartupMessage(m)
if (length(grep(CPL_geos_version(FALSE, TRUE), CPL_geos_version(TRUE))) != 1) { # nocov start
packageStartupMessage("WARNING: different compile-time and runtime versions for GEOS found:")
packageStartupMessage(paste(
"Linked against:", CPL_geos_version(TRUE, TRUE),
"compiled against:", CPL_geos_version(FALSE, TRUE)))
packageStartupMessage("It is probably a good idea to reinstall sf, and maybe rgeos and rgdal too")
} # nocov end
}
#' Provide the external dependencies versions of the libraries linked to sf
#'
#' Provide the external dependencies versions of the libraries linked to sf
#' @export
sf_extSoftVersion = function() {
structure(c(CPL_geos_version(), CPL_gdal_version(), CPL_proj_version(),
ifelse(CPL_gdal_with_geos(), "true", "false"),
ifelse(CPL_proj_h(), "true", "false")),
names = c("GEOS", "GDAL", "proj.4", "GDAL_with_GEOS", "USE_PROJ_H"))
}
load_gdal <- function() {
if (file.exists(system.file("proj/nad.lst", package = "sf")[1])) {
# nocov start
prj = system.file("proj", package = "sf")[1]
if (! CPL_set_data_dir(prj)) { # if TRUE, uses C API to set path, leaving PROJ_LIB alone
assign(".sf.PROJ_LIB", Sys.getenv("PROJ_LIB"), envir=.sf_cache)
Sys.setenv("PROJ_LIB" = prj)
}
CPL_use_proj4_init_rules(1L)
assign(".sf.GDAL_DATA", Sys.getenv("GDAL_DATA"), envir=.sf_cache)
gdl = system.file("gdal", package = "sf")[1]
Sys.setenv("GDAL_DATA" = gdl)
# nocov end
}
CPL_gdal_init()
register_all_s3_methods() # dynamically registers non-imported pkgs (tidyverse)
}
unload_gdal <- function() {
CPL_gdal_cleanup_all()
if (file.exists(system.file("proj/nad.lst", package = "sf")[1])) {
# nocov start
if (! CPL_set_data_dir(system.file("proj", package = "sf")[1])) # set back:
Sys.setenv("PROJ_LIB"=get(".sf.PROJ_LIB", envir=.sf_cache))
Sys.setenv("GDAL_DATA"=get(".sf.GDAL_DATA", envir=.sf_cache))
# nocov end
}
#units::remove_symbolic_unit("link")
#units::remove_symbolic_unit("us_in")
#units::remove_symbolic_unit("ind_yd")
#units::remove_symbolic_unit("ind_ft")
#units::remove_symbolic_unit("ind_ch")
}
#' @export
#' @name sf_project
#' @details \code{sf_add_proj_units} loads the PROJ units `link`, `us_in`, `ind_yd`, `ind_ft`, and `ind_ch` into the udunits database, and returns \code{TRUE} invisibly on success.
sf_add_proj_units = function() {
#nocov start
units::install_conversion_constant("m", "link", 0.201168)
units::install_conversion_constant("m", "us_in", 1./39.37)
units::install_conversion_constant("m", "ind_yd", 0.91439523)
units::install_conversion_constant("m", "ind_ft", 0.30479841)
units::install_conversion_constant("m", "ind_ch", 20.11669506)
invisible(TRUE)
#nocov end
}
|
