maaamet_ehak.R 2.9 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677
  1. #' Maa-ameti haldus- ja asustusjaotus
  2. #'
  3. #' Source: \url{https://geoportaal.maaamet.ee/est/Ruumiandmed/Haldus-ja-asustusjaotus-p119.html} . Andmed salvestatakse postgisi andmebaasi. Schema = 'maaamet'. Koniguratsiooni muutmiseks muuda konfiguratsiooni. Muutujate vaikeväärtused on sellised, et ei ole vaja midagi muuta.
  4. #' @param conf A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.
  5. #' @return No output.
  6. #' @seealso [ruut::get_config()], [ruut::copy_shp_to_db()]
  7. #' @keywords postgis, maps, ESRI Shpfile, OSM
  8. #' @export
  9. #' @examples
  10. #' ## Not run:
  11. #' ##
  12. #' ## maaamet_ehak()
  13. #' ##
  14. #' ## End(**Not run**)
  15. maaamet_ehak <- function(conf = NULL) {
  16. ans <- utils::askYesNo("Do you want to import maps into database?")
  17. if (!ans) {
  18. cat("\n------------------------\n")
  19. cat("Kaardikihte ei lisatud.")
  20. cat("\n------------------------\n")
  21. }
  22. if (ans) {
  23. # Download links
  24. urls <- c(
  25. "https://geoportaal.maaamet.ee/docs/haldus_asustus/maakond_shp.zip",
  26. "https://geoportaal.maaamet.ee/docs/haldus_asustus/omavalitsus_shp.zip",
  27. "https://geoportaal.maaamet.ee/docs/haldus_asustus/asustusyksus_shp.zip",
  28. "https://geoportaal.maaamet.ee/docs/haldus_asustus/omavalitsus_valispiirini_shp.zip"
  29. )
  30. # Estonia shapefile (ZIP archive)
  31. map_shapefile <- "maaamet_shapefile.zip"
  32. # Export to postgis
  33. if (is.null(conf)) {
  34. conf <- ruut::get_config()
  35. conf$schema <- "maaamet"
  36. }
  37. # New schema
  38. ruut::db_create_new_schema(conf = conf)
  39. for (j in 1:length(urls)) {
  40. # Temp directory
  41. tmp_dir <- sprintf("/tmp/maaamet_ehak/%s", j)
  42. if (!dir.exists(tmp_dir)) {
  43. dir.create(tmp_dir, recursive = T)
  44. }
  45. # Download and save
  46. saveTo <- sprintf("%s/%s", tmp_dir, map_shapefile)
  47. utils::download.file(
  48. url = urls[j],
  49. destfile = saveTo, method = "curl", extra = "-L"
  50. )
  51. # Unzip
  52. utils::unzip(saveTo, overwrite = T, exdir = tmp_dir)
  53. # List of files
  54. ls <- list.files(path = tmp_dir, pattern = ".dbf")
  55. ls_long <- list.files(path = tmp_dir, pattern = ".dbf", full.names = T)
  56. tbl_names <- unlist(strsplit(x = ls, split = ".dbf"))
  57. # Multi layer
  58. for (i in 1:length(tbl_names)) {
  59. print(tbl_names[i])
  60. conf$table <- gsub("[^[:alnum:]]", "", gsub("[[:digit:]]", "", tbl_names[i]))
  61. source <- sprintf('"%s" "%s"', tmp_dir, tbl_names[i])
  62. ## Export to postgis database.
  63. ruut::copy_shp_to_db(
  64. dir = tmp_dir, layer = tbl_names[i], conf = conf,
  65. id = "fid", crs_source = "EPSG:3301", crs_target = "EPSG:4326",
  66. geometry_type = "POLYGON"
  67. )
  68. cat(sprintf(
  69. "\nShp fail %s kopeeriti POSTGIS andmebaasi %s.%s\n\n",
  70. ls[i], conf$schema, conf$table
  71. ))
  72. }
  73. # Delete temp directory
  74. # system(sprintf("rm -rf %s/*", tmp_dir))
  75. }
  76. }
  77. }