ruut/R/gpkg_piirkonnale_ruudustike_lisamine.R

#' Piirkonna ruudustike leidmine
#'
#' 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 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()]
#' @keywords GPKG, boundary box, EPSG:3301
#' @export
#' @examples
#' \dontrun{
#'
#' gpkg_home <- "/tmp"
#' obj <- "marja"
#' gpkg_piirkonnale_ruudustike_lisamine(obj = obj, gpkg_home = gpkg_home)
#'
#' # Layers list.
#' dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
#' sf::st_layers(dsn = dsn)
#' }
gpkg_piirkonnale_ruudustike_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)
  }
  ## Loeme andmebaasist piiri ja piirikasti.
  pk <- sf::read_sf(dsn = dsn, layer = "piir")
  bb <- sf::read_sf(dsn = dsn, layer = "bb")
  cc <- sf::read_sf(sprintf("%s|layername=%s", dsn, "bb"))
  ## Konfiguratsiooni muutujale väärtuste omistamine
  conf <- ruut::get_config()
  conf$gpkg_home <- gpkg_home
  conf$gpkg_file <- obj
  ## epk02t_grid (100x100 meetrit ruudustik)
  conf$schema <- "maaamet"
  conf$table <- "epk2t_grid"
  input <- ruut::construct_to_gpkg_output_postgres_str(conf = conf)
  # ruut::qgis_algorithm_search_by_word(str = "extract")
  postgres <- sprintf(
    "postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'id\' srid=3301 type=Polygon checkPrimaryKeyUnicity=\'1\' ",
    conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password
  )
  ## Algorithm
  algorithm <- "native:extractbylocation"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))

  ## Loop
  intersect_layers <- c("piir", "bb")
  ruudud <- c("epk10t_grid", "epk2t_grid", "epk02t_grid", "epk10t", "epk2t")
  ruudud <- c("epk02t_grid")
  for (intersect in intersect_layers) {
    for (ruut in ruudud) {
      conf$gpkg_table <- sprintf("%s_%s", intersect, ruut)
      output <- ruut::construct_to_gpkg_output_file_str(conf = conf)
      result <- qgisprocess::qgis_run_algorithm(
        algorithm = algorithm,
        INPUT = sprintf(
          '%s table=\"%s\".\"%s\" (geom)',
          postgres, "maaamet", ruut
        ),
        INTERSECT = sprintf("/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/marja.gpkg|layername=%s", intersect),
        OUTPUT = output,
        PREDICATE = c(0)
      )
    }
  }
  ## Layers list
  sf::st_layers(dsn = dsn)
}