estmap/R/minu_teed_riigiteed_ls.R

#' Riigiteedevõrgu liiklussagedused
#'
#' Riigiteede liiklussagedused. Maanteeameti liiklussageduse kaart \url{https://maanteeamet.maps.arcgis.com/apps/View/index.html?appid=293d200a16454c1c84f2cfe35720149f}, liiklussageduse mõiste \url{https://www.mnt.ee/et/tee/liiklussagedus}. Andmestik teederegistri WFS (Web Feature Service). WFS andmestiku kirjeldus: \url{https://teeregister-api.mnt.ee/teenus/wfs?request=GetCapabilities&service=WFS}.
#'
#' @param conf A list() of configuration variables. Default values \code{\link[ruut]{get_config}}.
#' @return No output.
#' @seealso [ruut::get_config()], [ruut::copy_shp_to_db()], [estmap::minu_teed_riigiteed()], [estmap::minu_teed_muud_teed()], [estmap::minu_teed_eesti()], [estmap::minu_teed_riigiteed_valjapaasud()], [estmap::minu_teed_riigiteed_ls()]
#' @keywords postgis
#' @export
#' @examples
#' \dontrun{
#'
#' conf <- ruut::get_config()
#' minu_teed_riigiteed_ls(conf = conf)
#' }
minu_teed_riigiteed_ls <- function(conf = NULL) {
  ## --------------------- muutujad ja teisendused -----------------------------
  vars <- ajutised_muutujad(conf = conf)
  conf <- vars$conf

  # conf$schema <- "minu_teed"
  # conf$table <- "zzz"
  # output <- ruut::construct_to_gpkg_output_postgres_str(conf = conf, geometry_field = "geom")
  # cat(paste(paste0('"', ruut::db_table_colnames(conf = conf)$column_name, '"'), collapse = ","))

  ## -------- 1.1 Kontrollime andmebaasis vajalike tabelite olemasolu ----------
  conf$schema <- "teeregister_wfs"
  tbl_names <- ruut::db_schema_tablenames(conf = conf)
  if (!all(c("n_liiklussagedus") %in% tbl_names)) {
    cat("\n------------------\nPuuduvad vajalikud andmebaasitabelid.\n")
    return(NULL)
  }

  conf$schema <- "minu_teed"
  ## -------------- 1.1 Liilussageduste summeerimine (AKÖL) --------------------
  # ruut::qgis_algorithm_search_by_word("calculator")
  algorithm <- "native:fieldcalculator"
  # ruut::qgis_show_help(algorithm = algorithm)
  conf$schema <- "teeregister_wfs"
  conf$table <- "n_liiklussagedus"
  input <- ruut::construct_to_gpkg_output_postgres_str(
    conf = conf, geometry_type = "MultiLineString", srid = 3301,
    checkPrimaryKeyUnicity = TRUE, key = "id", geometry_field = "geometry"
  )
  output <- vars$tmp_gpkg_file_output_1
  cmd <- sprintf("qgis_process run %s --FIELD_NAME='akol' --FIELD_TYPE=1 --FIELD_LENGTH=6 --FIELD_PRECISION=0 --GLOBAL=\"\" --FORMULA=' \"sapaarv\" + \"vaabarv\" + \"ararv\" ' --INPUT='%s' --OUTPUT='%s' ", algorithm, input, output)
  system(cmd)
  cat(sprintf("\n-------------\nLoodi andmebaas '%s'.\n", output))

  ## ------------------- 1.2 Jätame alles vajalikud veerud ---------------------
  # ruut::qgis_algorithm_search_by_word("retain")
  algorithm <- "native:retainfields"
  # ruut::qgis_show_help(algorithm = algorithm)
  input <- vars$tmp_gpkg_file_input_1
  output <- vars$tmp_gpkg_file_output_2
  cmd <- sprintf("qgis_process run %s  --FIELDS='sapa' --FIELDS='vaab' --FIELDS='ar' --FIELDS='sapaarv' --FIELDS='vaabarv' --FIELDS='ararv' --FIELDS='akol' --FIELDS='tee_number' --FIELDS='tee_nimi' --INPUT='%s' --OUTPUT='%s' ", algorithm, input, output)
  system(cmd)
  cat(sprintf("\n-------------\nLoodi andmebaas '%s'.\n", output))

  ## ------------------ 1.3 Copy temp gpkg file to postgis ---------------------
  conf$schema <- "minu_teed"
  conf$table <- "teedevork_riigiteed_ls"
  PG <- ruut::construct_ogr2ogr_PG_connect_str(conf = conf)
  cmd <- sprintf("ogr2ogr -progress --config PG_USE_COPY YES -f PostgreSQL %s -lco DIM=2 %s %s -overwrite -nlt %s -lco GEOMETRY_NAME=geom -lco FID=id -nln %s.%s -nlt PROMOTE_TO_MULTI", PG, vars$tmp_gpkg_file, vars$layer_2, "MULTILINESTRING", conf$schema, conf$table)
  system(cmd)
  cat(sprintf("\n-------------\nLoodi andmebaas '%s'.\n", conf$table))
}