ruut/R/gpkg_piirkonnale_punktide_lisamine.R

#' Piirkonnale andmebaasides olevate punktide lisamine
#'
#' Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud polügoonid. Andmed salvestatakse GPKG faili kihtidena.
#'
#' @param obj str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja lisada.
#' @param gpkg_home path Salvestatavate GPKG faili asukoht.
#' @return Uute GPKG andmebaasi kihtide 'piir_...' ja 'bb_...' loomine.
#' @seealso [sf::st_read()], [sf::write_sf()],[sf::st_transform()],[ruut::gpkg_piirkonnale_ruudustike_lisamine()] ,[ruut::gpkg_piirkonnale_polygoonide_lisamine()],[ruut::gpkg_sellest_alustame_gpkg_loomist()],[ruut::gpkg_piirkonnale_joonte_lisamine()],[ruut::gpkg_piirkonnale_punktide_lisamine()]
#' @keywords GPKG, boundary box, EPSG:3301
#' @export
#' @examples
#' \dontrun{
#'
#' gpkg_home <- "/tmp"
#' obj <- "marja"
#' gpkg_piirkonnale_punktide_lisamine(obj = obj, gpkg_home = gpkg_home)
#'
#' # Layers list.
#' dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
#' sf::st_layers(dsn = dsn)
#' }
gpkg_piirkonnale_punktide_lisamine <- function(obj = NULL, gpkg_home = "/tmp") {
  dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
  if (!file.exists(dsn)) {
    cat(sprintf("\nSellist faili \"%s\" ei leitud.\n", dsn))
    return(NULL)
  }
  ## Konfiguratsiooni muutujale väärtuste omistamine
  conf <- ruut::get_config()
  tmp_gpkg_file <- tempfile(fileext = ".gpkg")
  conf$gpkg_home <- dirname(tmp_gpkg_file)
  conf$gpkg_file <- gsub(paste0(conf$gpkg_home, "/"), "", tmp_gpkg_file)
  conf$gpkg_file <- gsub(".gpkg", "", conf$gpkg_file)

  conf$gpkg_home <- gpkg_home
  conf$gpkg_file <- obj
  tmp_gpkg_file_output <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = F)
  tmp_gpkg_file_input <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = T)

  ## Algorithm
  # ruut::qgis_algorithm_search_by_word(str = "extract")
  # algorithm <- "native:extractbylocation"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))

  ## -------------------- Loop -----------------------
  intersect_layers <- c("piir", "bb")
  intersect_layers <- c("bb")
  andmed <- data.frame("schema" = character(0), "table" = character(0))
  andmed <- rbind(andmed, data.frame("schema" = "maaamet", "table" = "aadressandmed")) # POINT
  andmed <- rbind(andmed, data.frame("schema" = "osm_shp", "table" = "pofw"))
  andmed <- rbind(andmed, data.frame("schema" = "osm_shp", "table" = "pois"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_bussipeatus"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "kilomeetripostid"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_jaotus"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_kandur"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_mahasoit"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_onnetus"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_rdtyl"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_ristmik"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_ristumispunkt"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_sild"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_teeosa_points"))
  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_ylek"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_101_kivi_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_103_pinnavorm_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_202_seisuveekogu_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_301_muu_kolvik_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_305_puittaimestik_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_402_korgrajatis_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_403_muu_rajatis_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_602_tehnopaigaldis_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "kirik_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "kohanimi"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "korgus_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "roostik_p"))
  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "sygavus_p"))
  for (intersect in intersect_layers) {
    conf$gpkg_table <- intersect
    intersect_layer <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = T)
    for (i in 1:nrow(andmed)) {
      ## teeregister_wfs andmebaasis on importimisel geomeetria välja tähistus 'geometry'.
      if (andmed$schema[i] %in% c("teeregister_wfs")) geom <- "geometry" else geom <- "geom"
      conf$gpkg_table <- sprintf("%s_%s", intersect, andmed$table[i])
      conf$table <- andmed$table[i]
      conf$schema <- andmed$schema[i]
      output <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = F)
      input <- ruut::construct_to_gpkg_output_postgres_str(
        conf = conf, geometry_type = "Point", srid = 3301,
        checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = geom
      )
      ## ------------------- QGIS: qgisprocess ----------------------
      str <- paste0("{ 'INPUT' : '", input, "', 'INTERSECT' : '", intersect_layer, "', 'OUTPUT' : '", output, "', 'PREDICATE' : [0] }")
      algorithm <- "native:extractbylocation"
      cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
      cat(sprintf("\n%s\n\n", cmd))
      system(cmd)
      ## Filtreerime aadressandmete kihi eraldi alamkihtideks
      if (andmed$table[i] == "aadressandmed") {
        aadressandmed <- unique(as.data.frame(sf::read_sf(dsn = dsn, layer = sprintf("%s_aadressandmed", intersect), as_tibble = T))[, c("adob_liik"), drop = FALSE])
        parent_table <- conf$gpkg_table
        for (k in 1:nrow(aadressandmed)) {
          table_suffix <- aadressandmed$adob_liik[k]
          conf$gpkg_table <- parent_table
          input <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = TRUE)
          conf$gpkg_table <- sprintf("%s_%s", conf$gpkg_table, table_suffix)
          output <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = FALSE)
          str <- paste0("{ 'INPUT' : '", input, "', 'INTERSECT' : '", intersect_layer, "', 'OUTPUT' : '", output, "', FIELD : 'adob_liik', OPERATOR : 0, VALUE : ", aadressandmed$adob_liik[k], " }")
          algorithm <- "native:extractbyattribute"
          cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
          cat(sprintf("\n%s\n", cmd))
          system(cmd)
        }
      }
    }
  }
  ## Layers list
  sf::st_layers(dsn = dsn)
}