estmap/R/maaamet_ortofotod.R

#' Maa-ameti ortofotod
#'
#' Maa-ameti ortofotode allalaadimine. Andmestruktuur: raster. Mõõtkava 1:10 000 (vaata ruudustiku infot \url{https://geoportaal.maaamet.ee/est/ruumiandmed/kaardilehtede-susteemid-p224.html}). Loe: \url{https://geoportaal.maaamet.ee/est/andmed-ja-kaardid/ortofotod-p99.html}. Andmed salvestatakse kataloogi (default = '~/ortofotod'). Eesti topokaardistuse tarbeks toodetud ortofotodes on piksli suurusega 20-40 cm ja katavad kogu riigi territooriumi. Andmete kasutamisel palume viidata Maa-ametile. Näiteks: 'Ortofoto, Maa-amet 2021'. Käsitsi allalaadimine /url{https://geoportaal.maaamet.ee/index.php?lang_id=1&page_id=610}.
#' @importFrom magrittr %>%
#' @param epk10t 5-kohaline number, mis vastab kaardi 1:10 000 ruudu numbrile \url{https://geoportaal.maaamet.ee/est/ruumiandmed/kaardilehtede-susteemid-p224.html}.
#' @param dir GEOTIFF failide salvestamise kataloog.
#' @param conf A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.
#' @return GEOTIFF file.
#' @seealso [ruut::get_config()], [ruut::copy_shp_to_db()]
#' @keywords maa-amet, ortofoto
#' @export
#' @examples
#' ## Not run:
#' ##
#' ## maaamet_ortofotod()
#' ##
#' ## End(**Not run**)
maaamet_ortofotod <- function(epk10t = NULL, dir = "~/ortofotod", conf = NULL) {
  # epk10t <- 54711
  # epk10t numbri õigsuse kontroll
  epk10t_all <- utils::read.csv(file = system.file("csv", "epk10t.csv", package = "estmap"), header = T)[, "nr"]
  if (!epk10t %in% epk10t_all) {
    cat(sprintf("\n---------------\n\nValitud epk10t ruudu number %s ei ole \u00F5ige.\n", epk10t))
    return()
  }
  # GEOTIFF fotode salvestamise kataloog.
  if (!dir.exists(dir)) dir.create(dir, recursive = T)

  ## Geotiff failide allalaadimine
  # Näidis lehelt mille kaudu saab geotiff faile alla laadida.
  url_demo <- paste0(
    "https://geoportaal.maaamet.ee/index.php?",
    "lang_id=1&plugin_act=otsing&page_id=610&=&",
    "kaardiruut=62294&andmetyyp=ortofoto_eesti_rgb&_=1612349208090"
  )
  url_source <- urltools::url_decode(url_demo)
  params <- urltools::param_get(urls = url_source) # get params
  parsed_address <- urltools::url_parse(url_source)
  # Kui geotiff fail on olemas, siis me seda enam ei alla ei lae.
  if (length(list.files(path = dir, pattern = as.character(epk10t))) > 0) {
    cat(sprintf("\nFail %s on juba olemas.\n", epk10t))
    return()
  }
  url_new <- urltools::param_set(url_source, key = "kaardiruut", value = epk10t)
  ## Allalaadimine
  html_file <- readLines(url_new)
  page <- xml2::read_html(html_file)
  website <- ""
  website <- page %>%
    rvest::html_nodes("li a") %>%
    rvest::html_attr("href")
  parsed_address$parameter <- stringr::str_split(website[1], "[?]")[[1]][2]
  geotiff_download_url <- urltools::url_compose(parsed_urls = parsed_address)
  save_to_file <- sprintf("%s/%s.zip", dir, epk10t)
  utils::download.file(url = geotiff_download_url, method = "curl", extra = "-L", destfile = save_to_file)
  # zip failide lahtpakkimine
  cmd <- sprintf("unzip %s -d %s/; mv %s /tmp", save_to_file, dir, save_to_file)
  system(cmd)
  cmd <- sprintf("mv %s/*.tab /tmp", dir)
  system(cmd)
  cmd <- sprintf("mv %s/*.tfw /tmp", dir)
  system(cmd)
  ls <- list.files(path = dir, pattern = sprintf("%s.tif", epk10t))
  if (length(ls) > 0) {
    cat(sprintf("\n---------------\n\nValitud epk10t ruudu number %s eba\u00F5nnestus.\n", epk10t))
  } else {
    cat(sprintf("\n---------------\n\nValitud epk10t ruudu number %s \u00F5nnestus.\n", epk10t))
  }
  return(geotiff_download_url)
}