Muudetud joontega seotud asjad.

NAMESPACE

@@ -26,9 +26,11 @@ export(gpkg_sellest_alustame_gpkg_loomist)
 export(gpkg_teisendame_jooni)
 export(gpkg_teisendame_polygoone)
 export(gpkg_teisendame_punkte)
+export(pk_lisame_jooned)
 export(pk_lisame_polygoonid)
 export(pk_lisame_ruudustikud)
 export(pk_sellest_alustame_db_loomist)
+export(pk_teisendame_jooni)
 export(pk_teisendame_polygoone)
 export(qgis_algorithm_search_by_word)
 export(qgis_show_help)

R/pk_lisame_jooned.R

@@ -0,0 +1,103 @@
+#' Piirkonnale andmebaasides olevate joonte lisamine
+#'
+#' Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud polügoonid. 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_lisame_piirkonnale_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_jooned(obj = obj, conf = NULL)
+#'
+#' # Layers list.
+#' ruut::db_schema_tablenames(conf = conf)
+#' }
+pk_lisame_jooned <- function(obj = NULL, conf = NULL) {
+  obj <- gsub(" ", "_", tolower(obj))
+  ## Konfiguratsiooni muutujale väärtuste omistamine, kui seda pole antud.
+  if (is.null(conf)) {
+    conf <- ruut::get_config()
+    conf$schema <- obj
+  }
+  if (!any("piir" %in% ruut::db_schema_tablenames(conf = conf))) {
+    cat("\nAndmebaas loomata. Palun funktsiooniga ruut::pk_sellest_alustame_db_loomist() andmebaasi loomist.\n")
+    return(NULL)
+  }
+
+  ## ------------- muutujad ja teisendused ---------------
+  conn <- ruut::db_connect(conf = conf)
+  try(piir <- sf::st_read(dsn = conn, layer = c(conf$schema, "piir")))
+  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(piir, dsn = tmp_gpkg_file, layer = "piir", driver = "gpkg", fid_column_name = "id", delete_dsn = T)
+  # unlink(tmp_gpkg_file)
+
+  ## Algorithm
+  # ruut::qgis_algorithm_search_by_word(str = "extract")
+  # algorithm <- "native:extractbylocation"
+  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
+
+  ## -------------------- Loop -----------------------
+  intersect_layers <- c("bb2") # c("piir", "bb", "bb2")
+  andmed <- data.frame("schema" = character(0), "table" = character(0))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_alusdokumendid_ja_lepingud"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "pohimaantee"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "korvalmaantee"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "tugimaantee"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "muutee"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "ramp"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "teeosa"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_kergliiklustee"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_liiklussagedus"))
+  andmed <- rbind(andmed, data.frame("schema" = "teeregister_wfs", "table" = "n_omand"))
+  andmed <- rbind(andmed, data.frame("schema" = "gtfs", "table" = "shapes"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_203_vooluveekogu_j"))
+  # andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_305_puittaimestik_j"))
+  # andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_405_piire_j"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_501_tee_j"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_502_roobastee_j"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_505_liikluskorralduslik_rajatis_j"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "e_601_elektriliin_j"))
+  # andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "korgus_j"))
+  andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "sild_j"))
+  # andmed <- rbind(andmed, data.frame("schema" = "eesti", "table" = "sygavus_j"))
+  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 (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$table <- andmed$table[i]
+      conf$schema <- andmed$schema[i]
+      input <- ruut::construct_to_gpkg_output_postgres_str(
+        conf = conf, geometry_type = "Linestring", srid = 3301,
+        checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = geom
+      )
+      conf$table <- sprintf("%s", andmed$table[i])
+      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")
+      str <- paste0("{ 'INPUT' : '", input, "', 'INPUT_FIELDS' : [], 'OUTPUT' : '", output, "', 'OVERLAY' : '", intersect_layer, "', 'OVERLAY_FIELDS' : [], 'OVERLAY_FIELDS_PREFIX' : '' }")
+      algorithm <- "native:intersection"
+      cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
+      cat(sprintf("\n%s\n\n", cmd))
+      system(cmd)
+    }
+  }
+  ## Layers list
+  conf$schema <- obj
+  ruut::db_schema_tablenames(conf = conf)
+  sf::st_layers(dsn = tmp_gpkg_file)
+}

R/pk_teisendame_jooni.R

@@ -0,0 +1,102 @@
+#' Piirkonnale andmebaasides olevate joonte teisendamine
+#'
+#' Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud jooned. 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_lisame_piirkonnale_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_teisendame_jooni(obj = obj, conf = NULL)
+#'
+#' # Layers list.
+#' ruut::db_schema_tablenames(conf = conf)
+#' }
+pk_teisendame_jooni <- function(obj = NULL, conf = NULL) {
+  obj <- gsub(" ", "_", tolower(obj))
+  ## Konfiguratsiooni muutujale väärtuste omistamine, kui seda pole antud.
+  if (is.null(conf)) {
+    conf <- ruut::get_config()
+    conf$schema <- obj
+  }
+  if (!any("piir" %in% ruut::db_schema_tablenames(conf = conf))) {
+    cat("\nAndmebaas loomata. Palun funktsiooniga ruut::pk_sellest_alustame_db_loomist() andmebaasi loomist.\n")
+    return(NULL)
+  }
+
+  ## ------------- muutujad ja teisendused ---------------
+  conn <- ruut::db_connect(conf = conf)
+  try(piir <- sf::st_read(dsn = conn, layer = c(conf$schema, "piir")))
+  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)
+  conf$gpkg_table <- "test_layer_2"
+  tmp_gpkg_file_output_2 <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = F)
+  tmp_gpkg_file_input_2 <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = T)
+  # write to gpkg
+  sf::write_sf(piir, dsn = tmp_gpkg_file, layer = "piir", driver = "gpkg", fid_column_name = "id", delete_dsn = T)
+  # unlink(tmp_gpkg_file)
+
+  ## ====================== teisendused =====================
+  # Layers list
+  conf$schema <- obj
+  layer_names <- ruut::db_schema_tablenames(conf = conf)
+  layer_names <- layer_names[!layer_names %in% c("bb2", "epk10t_grid", "epk02t_grid", "epk2t_grid")] # c("bb", "bb2", "bb_epk10t_grid", "bb_epk02t_grid", "bb_epk2t_grid", "bb_epk10t", "bb_epk02t")]
+  for (layer_name in layer_names) {
+    try(ruumiline_obj <- sf::st_read(dsn = conn, layer = c(conf$schema, layer_name)))
+    ## Kontlrollime kas geomeetriline objekt on polügoon.
+    ## Arvutame sel juhul pindala.
+    is_line <- any(grepl("line", tolower(attributes(ruumiline_obj$geom)$class), fixed = TRUE))
+    if (is_line) {
+      ## Lõikame objekti ruudistikga tükkideks. Ühte ruutu võib jääda mitu tükki.
+      ## Peame need pärast ühendama
+      # ruut::qgis_algorithm_search_by_word("Intersection")
+      algorithm <- "native:intersection"
+      conf$table <- sprintf("%s", layer_name)
+      input <- ruut::construct_to_gpkg_output_postgres_str(
+        conf = conf, geometry_type = "MultiLineString", srid = 3301,
+        checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = "geom"
+      )
+      conf$table <- sprintf("%s", "epk02t_grid")
+      overlay <- ruut::construct_to_gpkg_output_postgres_str(
+        conf = conf, geometry_type = "Polygon", srid = 3301,
+        checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = "geom"
+      )
+      output <- tmp_gpkg_file_output # ajutine fail
+      str <- sprintf("{ 'INPUT' : '%s', 'INPUT_FIELDS' : ['fid'], 'OUTPUT' : '%s', 'OVERLAY' : '%s', 'OVERLAY_FIELDS' : [], 'OVERLAY_FIELDS_PREFIX' : '' }", input, output, overlay)
+      cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
+      system(cmd)
+      ## Ühendame ruudus olevad pinnad
+      # ruut::qgis_algorithm_search_by_word("Dissolve")
+      algorithm <- "native:dissolve"
+      input <- tmp_gpkg_file_input # ajutine fail
+      output <- tmp_gpkg_file_output_2 # ajutine fail
+      str <- sprintf("{ 'FIELD' : ['%s'], 'INPUT' : '%s', 'OUTPUT' : '%s' }", "id", input, output)
+      cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
+      system(cmd)
+      ## Lisame joonte pikkused
+      # ruut::qgis_algorithm_search_by_word("attributes")
+      algorithm <- "qgis:exportaddgeometrycolumns"
+      input <- tmp_gpkg_file_input_2 # ajutine fail
+      conf$table <- sprintf("%s_epk02t", layer_name)
+      output <- ruut::construct_to_gpkg_output_postgres_str(conf = conf, geometry_field = "geom", geometry_type = "Polygon", srid = 3301, checkPrimaryKeyUnicity = FALSE, key = "id")
+      str <- sprintf("{ 'CALC_METHOD' : 0, 'INPUT' : '%s', 'OUTPUT' : '%s' }", input, output)
+      cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
+      system(cmd)
+    }
+  }
+  ## Layers list
+  conf$schema <- obj
+  ruut::db_schema_tablenames(conf = conf)
+  sf::st_layers(dsn = tmp_gpkg_file)
+}

man/pk_lisame_jooned.Rd

@@ -0,0 +1,36 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/pk_lisame_jooned.R
+\name{pk_lisame_jooned}
+\alias{pk_lisame_jooned}
+\title{Piirkonnale andmebaasides olevate joonte lisamine}
+\usage{
+pk_lisame_jooned(obj = NULL, conf = NULL)
+}
+\arguments{
+\item{obj}{str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja lisada.}
+
+\item{conf}{A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.}
+}
+\value{
+Uute andmebaasi kihtide 'piir_...' ja 'bb_...' loomine.
+}
+\description{
+Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud polügoonid. Andmed salvestatakse postgis andmebaasi.
+}
+\examples{
+\dontrun{
+
+obj <- "marja"
+pk_lisame_jooned(obj = obj, conf = NULL)
+
+# Layers list.
+ruut::db_schema_tablenames(conf = conf)
+}
+}
+\seealso{
+\code{\link[sf:st_read]{sf::st_read()}}, \code{\link[sf:st_write]{sf::write_sf()}},\code{\link[sf:st_transform]{sf::st_transform()}},\code{\link[=pk_sellest_alustame_db_loomist]{pk_sellest_alustame_db_loomist()}},\code{\link[=pk_lisame_ruudustikud]{pk_lisame_ruudustikud()}} ,\code{\link[=pk_lisame_polygoonid]{pk_lisame_polygoonid()}},\code{\link[=pk_lisame_jooned]{pk_lisame_jooned()}},\code{\link[=pk_lisame_punktid]{pk_lisame_punktid()}},\code{\link[=pk_lisame_piirkonnale_punktid]{pk_lisame_piirkonnale_punktid()}},\code{\link[=pk_teisendame_polygoone]{pk_teisendame_polygoone()}},\code{\link[=pk_teisendame_jooni]{pk_teisendame_jooni()}},\code{\link[=pk_teisendame_punkte]{pk_teisendame_punkte()}}
+}
+\keyword{EPSG:3301}
+\keyword{boundary}
+\keyword{box,}
+\keyword{postgis,}

man/pk_teisendame_jooni.Rd

@@ -0,0 +1,36 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/pk_teisendame_jooni.R
+\name{pk_teisendame_jooni}
+\alias{pk_teisendame_jooni}
+\title{Piirkonnale andmebaasides olevate joonte teisendamine}
+\usage{
+pk_teisendame_jooni(obj = NULL, conf = NULL)
+}
+\arguments{
+\item{obj}{str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja lisada.}
+
+\item{conf}{A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.}
+}
+\value{
+Uute andmebaasi kihtide 'piir_...' ja 'bb_...' loomine.
+}
+\description{
+Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud jooned. Andmed salvestatakse postgis andmebaasi.
+}
+\examples{
+\dontrun{
+
+obj <- "marja"
+pk_teisendame_jooni(obj = obj, conf = NULL)
+
+# Layers list.
+ruut::db_schema_tablenames(conf = conf)
+}
+}
+\seealso{
+\code{\link[sf:st_read]{sf::st_read()}}, \code{\link[sf:st_write]{sf::write_sf()}},\code{\link[sf:st_transform]{sf::st_transform()}},\code{\link[=pk_sellest_alustame_db_loomist]{pk_sellest_alustame_db_loomist()}},\code{\link[=pk_lisame_ruudustikud]{pk_lisame_ruudustikud()}} ,\code{\link[=pk_lisame_polygoonid]{pk_lisame_polygoonid()}},\code{\link[=pk_lisame_jooned]{pk_lisame_jooned()}},\code{\link[=pk_lisame_punktid]{pk_lisame_punktid()}},\code{\link[=pk_lisame_piirkonnale_punktid]{pk_lisame_piirkonnale_punktid()}},\code{\link[=pk_teisendame_polygoone]{pk_teisendame_polygoone()}},\code{\link[=pk_teisendame_jooni]{pk_teisendame_jooni()}},\code{\link[=pk_teisendame_punkte]{pk_teisendame_punkte()}}
+}
+\keyword{EPSG:3301}
+\keyword{boundary}
+\keyword{box,}
+\keyword{postgis,}