ruut/R/pk_lisame_ruudustikud.R

#' Piirkonna ruudustike leidmine ja lisamine
#'
#' Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud 5x5 km ('epk10t_grid'), 1x1 km ('epk2t_grid') ja 100x100 meetrit ('epk02t_grid') ruudustikud. Andmed salvestatakse postgis andmebaasi.
#'
#' @param obj str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja lisada.
#' @param conf A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.
#' @return Uute andmebaasi kihtide 'piir_...' ja 'bb_...' loomine.
#' @seealso [sf::st_read()], [sf::write_sf()],[sf::st_transform()],[ruut::pk_sellest_alustame_db_loomist()],[ruut::pk_lisame_ruudustikud()] ,[ruut::pk_lisame_polygoonid()],[ruut::pk_lisame_jooned()],[ruut::pk_lisame_punktid()],[ruut::pk_teisendame_polygoone()],[ruut::pk_teisendame_jooni()],[ruut::pk_teisendame_punkte()]
#' @keywords postgis, boundary box, EPSG:3301
#' @export
#' @examples
#' \dontrun{
#'
#' obj <- "marja"
#' pk_lisame_ruudustikud(obj = obj, conf = NULL)
#'
#' # Layers list.
#' ruut::db_schema_tablenames(conf = conf)
#' }
pk_lisame_ruudustikud <- function(obj = NULL, conf = NULL) {
  ## ------------- muutujad ja teisendused ---------------
  vars <- ajutised_muutujad(pk = NULL, obj, conf)
  obj <- vars$obj
  piir <- vars$pk
  conf <- vars$conf
  ## Konfiguratsiooni muutujale väärtuste omistamine, kui seda pole antud.
  if (!any("a00_piir" %in% ruut::db_schema_tablenames(conf = conf))) {
    cat(
      "\nAndmebaas loomata. Palun funktsiooniga ruut::pk_sellest_alustame_db_loomist() andmebaasi loomist.\n"
    )
    return(NULL)
  }

  ## Algorithm
  # ruut::qgis_algorithm_search_by_word(str = "extract")
  algorithm <- "native:extractbylocation"
  # cat(ruut::qgis_show_help(algorithm = algorithm))
  ## -------------------- Loop -----------------------
  intersect_layers <- c("a00_bb2") # c("piir", "bb", "bb2")
  ruudud <-
    c("epk10t_grid",
      "epk2t_grid",
      "epk02t_grid",
      "epk10t",
      "epk2t") #
  conf$schema <- "maaamet"
  ## Tabelid mis on juba andmebaasis ja mida ümber ei kirjutata.
  conf$schema <- obj
  olemasolevad_tabelid <- ruut::db_schema_tablenames(conf = conf)
  # ruudud <- c("epk02t_grid")
  for (intersect in intersect_layers) {
    conf$table <- intersect
    conf$schema <- obj
    intersect_layer <-
      ruut::construct_to_gpkg_output_postgres_str(
        conf = conf,
        geometry_type = "Polygon",
        srid = 3301,
        checkPrimaryKeyUnicity = TRUE,
        key = "id"
      )
    for (ruut in ruudud) {
      ## Kui olemasolev tabel eksisteerib andmebaasis, siis jätame arvutused
      ## selle tabeliga vahele.
      uus_tabel <- sprintf("bb2_%s", ruut)
      if (uus_tabel %in% olemasolevad_tabelid) {
        cat(sprintf("\nTabel %s on juba andmebaasis olemas.\n", uus_tabel))
      } else {
        conf$table <- ruut
        conf$schema <- "maaamet"
        input <-
          ruut::construct_to_gpkg_output_postgres_str(
            conf = conf,
            geometry_type = "MultiPolygon",
            checkPrimaryKeyUnicity = T,
            key = "id",
            srid = 3301
          )
        conf$table <- uus_tabel
        conf$schema <- obj
        output <-
          ruut::construct_to_gpkg_output_postgres_str(
            conf = conf,
            geometry_field = "geom",
            geometry_type = NULL,
            srid = 3301,
            checkPrimaryKeyUnicity = FALSE,
            key = "id"
          )
        ## ------------------- QGIS: qgisprocess ----------------------
        str <-
          paste0(
            "{ 'INPUT' : '",
            input,
            "', 'INTERSECT' : '",
            intersect_layer,
            "', 'OUTPUT' : '",
            output,
            "', 'PREDICATE' : [5,6] }"
          )
        cmd <-
          ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
        cat(sprintf("\n%s\n", cmd))
        system(cmd)
      }
    }
  }
  ## Layers list
  conf$schema <- obj
  ruut::db_schema_tablenames(conf = conf)
  sf::st_layers(dsn = vars$tmp_gpkg_file)
}