setwd(paste0("/home/ardo/apps/R/packages/", new_pkg_name))
new_pkg_name <- "estmap"
setwd(paste0("/home/ardo/apps/R/packages/", new_pkg_name))
document()
styler:::style_active_file()
library(estmap)
estmap::maaamet_ortofotod(epk10t = 54711)
library(estmap)
estmap::maaamet_ortofotod(epk10t = 54711)
estmap::maaamet_ortofotod(epk10t = 54711)
estmap::maaamet_ortofotod(epk10t = 54711)
estmap::maaamet_ortofotod(epk10t = 54711)
library(estmap)
estmap::maaamet_ortofotod(epk10t = 54711)
ls_long[6]
conf = NULL
# Temp directory
tmp_dir <- "/tmp/maaamet_aadressandmed"
if (!dir.exists(tmp_dir)) {
dir.create(tmp_dir)
}
# Download link
url <- "https://xgis.maaamet.ee/adsavalik/valjav6te/"
# Estonia shapefile (ZIP archive)
map_shapefile <- "aadressandmed.zip"
# Download and save
# Aadressandmete õige lingi tuvastamine
adsavalik <- jsonlite::fromJSON(url)
url_download <- sprintf("%s%s", url, adsavalik$fail[which(adsavalik$vvnr == 1 & is.na(adsavalik$kov))])
saveTo <- sprintf("%s/%s", tmp_dir, map_shapefile)
if (!file.exists(saveTo)) {
utils::download.file(
url = url_download,
destfile = saveTo, method = "curl", extra = "-L"
)
}
# Unzip
utils::unzip(saveTo, overwrite = T, exdir = tmp_dir)
# List of files
ls <- list.files(path = tmp_dir, pattern = ".dbf")
ls_long <- list.files(path = tmp_dir, pattern = ".dbf", full.names = T)
tbl_names <- unlist(strsplit(x = ls, split = ".dbf"))
# Export to postgis
if (is.null(conf)) {
conf <- ruut::get_config()
conf$schema <- "maaamet"
conf$table <- "aadressandmed"
}
# New schema
ruut::db_create_new_schema(conf = conf)
# List of files
ls <- list.files(path = tmp_dir, pattern = ".csv")
ls_long <- list.files(path = tmp_dir, pattern = ".csv", full.names = T)
tbl_names <- unlist(strsplit(x = ls, split = ".csv"))
# Multi layer
for (i in 1:length(tbl_names)) {
print(tbl_names[i])
conf$table <- tbl_names[i]
source <- sprintf('"%s" "%s"', tmp_dir, tbl_names[i])
## Export to postgis database.
ruut::copy_shp_to_db(
dir = tmp_dir, layer = tbl_names[i], conf = conf,
id = "fid", crs_source = "EPSG:3301", crs_target = "EPSG:4326",
geometry_type = "POINT"
)
cat(sprintf(
"\nShp fail %s kopeeriti POSTGIS andmebaasi %s.%s\n\n",
ls[i], conf$schema, conf$table
))
}
tbl_names[i]
library(estmap)