ruut/R/pk_sellest_alustame_db_loomist.R

#' Piirkonna postgis andmebaasi loomine
#'
#' Funktsioon loob piirkonna geojoonest ja nimest postgis andmebaasi. Samas lisatakse faili piir ja piiri piirkast. Ülejäänud funktsioonid kasutavad selle funktsiooni poolt loodud objekti nime ja piiri ega piirikasti pole vaja enam teiste funktsioonide rakendamisel lisada.
#'
#' @importFrom magrittr %>%
#' @importFrom rlang .data
#' @param pk Piirkonna geomeetriline joon.
#' @param obj str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja edaspidi lisada.
#' @param conf A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.
#' @return Andmebaasi loomine koos kihtidega 'piir' ja 'bb' (boundary box, piirikast).
#' @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{
#'
#' # read geojson from string:
#' geojson_txt <- paste('{"type":"MultiPoint","coordinates":[[658300,6474800],
#' [658300,6475000],[658700,6476000],[659600,6475000],[659600,6474800],
#' [658300,6474800]]}')
#' x <- sf::read_sf(geojson_txt) %>% sf::st_cast("POLYGON")
#' sf::st_crs(x) <- 3301
#' x
#' sf::st_geometry(x) %>% graphics::plot()
#'
#' obj <- "marja"
#' pk_sellest_alustame_db_loomist(pk = x, obj = obj)
#'
#' # Layers list.
#' ruut::db_schema_tablenames(conf = conf)
#' }
pk_sellest_alustame_db_loomist <- function(pk = NULL, obj = NULL, conf = NULL) {
  ## ----------- argumentide vastavuse kontroll -----------
  if (is.null(pk) || !sf::st_is_valid(pk)) {
    cat("\nPalun kontrolli geomeetrilise kujundi \u00F%igsust.\n")
    return(NULL)
  }
  attributes(pk)$names <- attributes(pk)$sf_column <- "geom"
  if (is.null(obj)) {
    cat("\nPalun sisesta objekti nimi.\n")
    return(NULL)
  }
  obj <- gsub(" ", "_", tolower(obj))
  ## Konfiguratsiooni muutujale väärtuste omistamine, kui seda pole antud.
  if (is.null(conf)) {
    conf <- ruut::get_config()
    conf$table <- obj
  }
  ## ------------- muutujad ja teisendused ---------------
  tmp_dir <- tempdir()
  tmp_gpkg_file <- sprintf("%s/%s.gpkg", tmp_dir, obj)
  system(sprintf("touch %s", tmp_gpkg_file))
  conf$gpkg_home <- tmp_dir
  conf$gpkg_file <- obj
  conf$gpkg_table <- "test_layer"
  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)
  # write to gpkg
  sf::write_sf(pk, dsn = tmp_gpkg_file, layer = "piir", driver = "gpkg", fid_column_name = "id", delete_dsn = T)
  # unlink(tmp_gpkg_file)

  ## ------------------- programmi osa -----------------

  ## 1. Piirjoone salvestamine
  pk <- sf::st_transform(pk, 3301)
  # write to postgis
  conf$schema <- conf$gpkg_file
  ruut::db_create_new_schema(conf = conf)
  conn <- ruut::db_connect(conf = conf)
  sf::write_sf(pk, dsn = conn, layer = c(conf$schema, "piir"), driver = "PostgreSQL", fid_column_name = "id", delete_dsn = F)


  ## 2. Piirikasti salvestamine
  ## konstrueerime EXTENT, HSPACING ja VSPACING väärtused piirikasti konstrueerimiseks.
  pk_attributes <- attributes(pk$geometry)
  extent <- sprintf("%s,%s,%s,%s [EPSG:3301]", round((pk_attributes$bbox["xmin"] / 100), digits = 0) * 100, ceiling((pk_attributes$bbox["xmax"] / 100)) * 100, round((pk_attributes$bbox["ymin"] / 100), digits = 0) * 100, ceiling((pk_attributes$bbox["ymax"] / 100)) * 100)
  hspacing <- (ceiling((pk_attributes$bbox["xmax"] / 100)) * 100) - (round((pk_attributes$bbox["xmin"] / 100), digits = 0) * 100)
  vspacing <- (ceiling((pk_attributes$bbox["ymax"] / 100)) * 100) - (round((pk_attributes$bbox["ymin"] / 100), digits = 0) * 100)
  conf$table <- "bb"
  output <- ruut::construct_to_gpkg_output_postgres_str(conf = conf, geometry_field = "geom", geometry_type = NULL, srid = 3301, checkPrimaryKeyUnicity = FALSE, key = "id")
  # ruut::qgis_algorithm_search_by_word(str = "grid")
  algorithm <- "native:creategrid"
  # cat(ruut::qgis_show_help(algorithm = algorithm))
  cmd <- sprintf("qgis_process run %s --CRS=\"EPSG:3301\" --EXTENT=\"%s\" --HOVERLAY=0 --HSPACING=%s --TYPE=2 --VOVERLAY=0 --VSPACING=%s --OUTPUT=\"%s\"", algorithm, extent, hspacing, vspacing, output)
  cat(sprintf("\n%s\n", cmd))
  system(cmd)


  ## 3. epk2t piirikasti konstrueerimine ja salvestamine s.t kuhu mahuvad sisse
  ## meie piirkonna epk2t (1x1km) ruudud.
  ## 3.1 Leiame piirikastiga seotud epk2t ruudud
  # str <- "{ 'INPUT' : 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable  password=\'osm\' key=\'id\' srid=3301 type=MultiPolygon checkPrimaryKeyUnicity=\'1\' table=\"maaamet\".\"epk2t_grid\" (geom)', 'INTERSECT' : '/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/marja.gpkg|layername=bb', 'OUTPUT' : 'TEMPORARY_OUTPUT', 'PREDICATE' : [5,6] }"
  # ruut::qgis_algorithm_search_by_word(str = "extract")
  algorithm <- "native:extractbylocation"
  # cat(ruut::qgis_show_help(algorithm = algorithm))
  conf$table <- "epk2t_grid"
  conf$schema <- "maaamet"
  input <- ruut::construct_to_gpkg_output_postgres_str(
    conf = conf, geometry_type = "Polygon", srid = 3301,
    checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = "geom"
  )
  conf$table <- "bb"
  conf$schema <- obj
  intersect <- ruut::construct_to_gpkg_output_postgres_str(conf = conf, geometry_type = "Polygon", srid = 3301, checkPrimaryKeyUnicity = TRUE, key = "id")
  output <- tmp_gpkg_file_output # ajutine fail
  str <- sprintf("{ 'INPUT' : '%s', 'INTERSECT' : '%s', 'OUTPUT' : '%s', 'PREDICATE' : [5,6] }", input, intersect, output)
  cmd <- ruut::construct_qgis_output_result_to_better_format(str, algorithm)
  system(cmd)

  ## Ajutise faili abil piirikasti konstrueerimine
  test_layer_obj <- sf::read_sf(dsn = tmp_gpkg_file, layer = conf$gpkg_table)
  pk_attributes <- attributes(test_layer_obj$geom)
  extent <- sprintf("%s,%s,%s,%s [EPSG:3301]", round((pk_attributes$bbox["xmin"] / 100), digits = 0) * 100, ceiling((pk_attributes$bbox["xmax"] / 100)) * 100, round((pk_attributes$bbox["ymin"] / 100), digits = 0) * 100, ceiling((pk_attributes$bbox["ymax"] / 100)) * 100)
  hspacing <- (ceiling((pk_attributes$bbox["xmax"] / 100)) * 100) - (round((pk_attributes$bbox["xmin"] / 100), digits = 0) * 100)
  vspacing <- (ceiling((pk_attributes$bbox["ymax"] / 100)) * 100) - (round((pk_attributes$bbox["ymin"] / 100), digits = 0) * 100)
  conf$table <- "bb2"
  output <- ruut::construct_to_gpkg_output_postgres_str(conf = conf, geometry_field = "geom", geometry_type = NULL, srid = 3301, checkPrimaryKeyUnicity = FALSE, key = "id")
  # ruut::qgis_algorithm_search_by_word(str = "grid")
  algorithm <- "native:creategrid"
  # cat(ruut::qgis_show_help(algorithm = algorithm))
  cmd <- sprintf("qgis_process run %s --CRS=\"EPSG:3301\" --EXTENT=\"%s\" --HOVERLAY=0 --HSPACING=%s --TYPE=2 --VOVERLAY=0 --VSPACING=%s --OUTPUT=\"%s\"", algorithm, extent, hspacing, vspacing, output)
  cat(sprintf("\n%s\n", cmd))
  system(cmd)
  # Layers list
  conf$schema <- obj
  ruut::db_schema_tablenames(conf = conf)
  sf::st_layers(dsn = tmp_gpkg_file)
}