210127_artikkel/01_ajutine.R

source("functions/delete_existing_variables.R")
source("functions/omavalitsused_nimekiri.R")
source("functions/spatial_indeksite_lisamine.R")

library(dplyr)
# library(qgisprocess)
library(sf)
library(ruut)

start.time <- Sys.time()

# Plot'i marginite määramine.
par(oma = c(0, 0, 0, 0)) # outer margin
par(mar = c(0, 0, 0, 0) + 0.0)

# Baaskonfiguratsiooni lugemine.
conf <- ruut::get_config()

#' Omavalitsused
onimi <- omavalitsused_nimekiri()
#' Valime omavalitsuse
i <- 1 # Ruhnu vald
# i <- 5 # Hiiumaa vald
pk_nimi <- onimi[i]
#' Omavalitsuse piir.
conn <- ruut::db_connect()
q <- sprintf("SELECT geom
FROM maaamet.omavalitsus
WHERE onimi ='%s';",
             pk_nimi)
cat(sprintf("\n-----------------\n%s\n\n", q))
x <- sf::st_read(conn, query = q)
sf::st_crs(x) <- 3301
x
# sf::st_geometry(x) %>% graphics::plot()

obj <- paste0("pk_", tolower(strsplit(pk_nimi, split = " ")[[1]][1]))
conf$schema <- obj
# ruut::db_create_new_schema(conf = conf)
## --------------------------------------

## postgisi loomine
ruut::pk_sellest_alustame_db_loomist(pk = x, obj = obj)
## spatial indeksite lisamine
lisa_spatial_indeksid(schema = obj)
# Layers list
ruut::db_schema_tablenames(conf = conf)
## ruudustike lisamine
ruut::pk_lisame_ruudustikud(obj = obj)
## spatial indeksite lisamine
lisa_spatial_indeksid(schema = obj)

## polügoonide lisamine
ruut::pk_lisame_polygoonid(obj = obj)
# joonte lisamine
ruut::pk_lisame_jooned(obj = obj)
# punktide lisamine
ruut::pk_lisame_punktid(obj = obj)
# spatial indeksite lisamine
lisa_spatial_indeksid(schema = obj)

## polügoonide teisendamine
ruut::pk_teisendame_polygoone(obj = obj)
## joonte teisendamine
ruut::pk_teisendame_jooni(obj = obj)
# punktide teisendamine
ruut::pk_teisendame_punkte(obj = obj)
# spatial indeksite lisamine
lisa_spatial_indeksid(schema = obj)

## liiklussageduse kihtide eraldamine
ruut::pk_teisendame_liiklussagedusi(obj = obj)

##
# tabels_list <- ruut::db_schema_tablenames(conf = conf)
# tabels_list[grepl("^data_p_", tabels_list)]
# tabels_list[grepl("^grid_p_", tabels_list)]

# rahvastikuarvude lisamine
ruut::pk_lisame_rahvaarvud(obj = obj)

# Raster tabelid raster T/F tabeliteks.
ruut::pk_tif_to_true(obj = obj)

## GPKG failide nimekiri
# gpkg_home <- "/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa"
# dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
# # Layers list
# sf::st_layers(dsn = dsn)

end.time <- Sys.time()
time.taken <- end.time - start.time
time.taken