210127_artikkel/gpkg/04_piirkonnale_joonte_lisamine.R

#' Piirkonna joonte lisamine
#'
#' Etteantud piirkonna geomeetrilise piirjoone ('piir') järele leitakse selles piirkonnas ja piirkonna piiriga piirnevad joonobjektid ().
#'
#' @param obj str Objekti nimi.
#' @param pk Piirkonna geomeetriline joon.
#' @param gpkg_home path Salvestatavate GPKG faili asukoht.
#' @param obj_only TRUE/FALSE Kas geomeetriad leitakse ainult objekti või kogu piirikastiga määratud ala jaoks.
#'
#' @examples
#' \dontrun{
#'
#' }
piirkonnale_joonte_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "/tmp", obj_only = TRUE) {
  if (is.null(pk) || !sf::st_is_valid(pk)) {
    cat("\nPalun kontrolli geomeetrilise kujundi õigsust.\n")
    return()
  }
  if (is.null(obj)) {
    cat("\nPuudu on objekti nimi.\n")
    return()
  }
  cat(sprintf("\nAlgparameetrid: objekti nimi %s ja GPKG faili kataloog %s\n", obj, gpkg_home))
  postgres <- sprintf(
    "postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=4326 type=%s checkPrimaryKeyUnicity=\'1\'",
    conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password, "Polygon"
  )
  postgres_point <- sprintf(
    "postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=4326 type=%s checkPrimaryKeyUnicity=\'1\'",
    conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password, "Point"
  )
  postgres_line <- sprintf(
    "postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=4326 type=%s checkPrimaryKeyUnicity=\'1\'",
    conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password, "Linestring"
  )
  dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj) # Polygon
  ogr <- sprintf(
    "ogr:dbname='%s'",
    dsn
  )
  tmp_gpkg_file <- tempfile(fileext = ".gpkg")


  ## ----------------- punktfaili loomine -------------------
  ## -------------- epk2t_grids tsentroidid -----------------
  # ruut::qgis_algorithm_search_by_word(str = "polygons")
  algorithm <- "native:polygonstolines"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
  result <- qgisprocess::qgis_run_algorithm(
    algorithm = algorithm,
    INPUT = sprintf("%s|layername=%s", dsn, "piir"),
    OUTPUT = sprintf('%s table=\"%s\" (geometry)', ogr, "piir_l")
    # .quiet = TRUE
  )

  ## ----------------- piiri sisse jäävad punktid -------------------
  ## --------------------------- algandmed --------------------------
  andmed.df <- data.frame("schema" = character(0), "table" = character(0))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "n_alusdokumendid_ja_lepingud"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "korvalmaantee"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "muutee"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "n_kergliiklustee"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "n_liiklussagedus"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "n_omand"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "pohimaantee"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "ramp"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "teeosa"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "teeregister_wfs", "table" = "tugimaantee"))
  andmed.df <- rbind(andmed.df, data.frame("schema" = "gtfs", "table" = "shapes"))
  j <- 1
  for (j in 1:nrow(andmed.df)) {
    output_layer_name <- as.character(andmed.df$table[j])
    # ruut::qgis_algorithm_search_by_word(str = "intersect")
    algorithm <- "native:intersection"
    # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
    result <- qgisprocess::qgis_run_algorithm(
      algorithm = algorithm,
      INPUT = sprintf('%s table="%s"."%s" (geometry)', postgres_line, andmed.df$schema[j], andmed.df$table[j]),
      INPUT_FIELDS = "",
      # OUTPUT = sprintf('%s table=\"%s\" (geometry)', ogr, andmed.df$table[j]),
      OUTPUT = tmp_gpkg_file,
      OVERLAY = sprintf("%s|layername=%s", dsn, "piir"),
      OVERLAY_FIELDS = "",
      OVERLAY_FIELDS_PREFIX = ""
      # # .quiet = TRUE
    )
    result
    # assign(as.character(andmed.df$table[j]), sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT")))
    system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s -t_srs \"EPSG:3301\"", dsn, tmp_gpkg_file, output_layer_name))
    ## AINULT ELAMUD aadressandmetest!!
    # if (output_layer_name == "aadressandmed") {
    #   system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s_ee -t_srs \"EPSG:3301\" -where \"adob_liik = 'EE'\"  ", dsn, tmp_gpkg_file, output_layer_name))
    # }
  }


  ## ---------------------- vaata layer'id ----------------------
  # Vaata layer'eid
  sf::st_layers(dsn = dsn)
}
# piirkonnale_joonte_lisamine(obj = NULL, pk = NULL, gpkg_home = "/tmp")