210127_artikkel/00_algandmed.R

#' Andmebaasi loomine ja ühendamine
#'

source("functions/delete_existing_variables.R")

library(dplyr)
library(qgisprocess)
library(sf)

par(oma = c(0, 0, 0, 0)) # outer margin
par(mar = c(0, 0, 0, 0) + 0.0)

source("01_funktsioonid.R")

# Uue schema loomine
conf <- ruut::get_config()
conf$schema <- "xxx_artikkel_210127"
ruut::db_create_new_schema(conf = conf)


# Objektide nimekiri
objektid <- c("valga", "matsalu", "lahemaa")

# ----------------- Loe piirkond (objekt) -----------------
## Piirkonna 'pk'
conn <- ruut::db_connect(conf = conf)
## Valitud objekti indeks
i <- 1

## 1. piirkonna piir
pk <- pk_piir(obj = objektid[i])
sf::st_geometry(pk) %>% plot()
## 2. piirkonna 3301 projektsiooniga piirikast
conn <- ruut::db_connect()
q <- sprintf("SELECT * FROM %s.%s_bb", "xxx_artikkel_210127", objektid[i])
cat(sprintf("\n-----------------\n%s\n\n", q))
pk_bb_3301 <- sf::st_read(conn, query = q)
sf::st_geometry(pk) %>% plot()
sf::st_geometry(pk_bb_3301) %>% plot(add = T)

## 3. piirkonna epk10t ruudud
# 3.1 kogu ruutvõrgustik
epk10t_grid <- pk_epk10t_grid(obj = objektid[i])
# 3.1 ainult piirkonna ruutvõrgustik
epk10t <- pk_epk10t(obj = objektid[i])
sf::st_geometry(epk10t_grid) %>% plot(border = 3, lwd = 0.3, col = "#d3fffb")
sf::st_geometry(epk10t) %>% plot(add = T, border = 3, lwd = 0.3, col = "#a3fffb")
sf::st_geometry(pk) %>% plot(add = T)

# 4 piirkonna epk10t kaardiruutude nimekiri ortofotode allalaadimiseks
epk10t_nr <- pk_epk10t_ruutude_nimekiri(objektid[i])

## 5. piirkonna epk2t ruudud
# 3.1 kogu ruutvõrgustik
epk2t_grid <- pk_epk2t_grid(obj = objektid[i])
# 3.1 ainult piirkonna ruutvõrgustik
epk2t <- pk_epk2t(obj = objektid[i])
sf::st_geometry(epk2t_grid) %>% plot(border = 3, lwd = 0.3, col = "#d3fffb")
sf::st_geometry(epk2t) %>% plot(add = T, border = 3, lwd = 0.3, col = "#a3fffb")
sf::st_geometry(pk) %>% plot(add = T)




## -------------- Muud ruudustikega seotud demo joonised ----------------
## 2. Kaardiruudustiku epk200t (100x100km) piirikast
conn <- ruut::db_connect()
q <- sprintf("SELECT * FROM %s.%s", "maaamet", "epk200t_bb")
cat(sprintf("\n-----------------\n%s\n\n", q))
epk200t_bb <- sf::st_read(conn, query = q)
sf::st_geometry(epk200t_bb) %>% plot()

## 3. Kaardiruudustiku epk200t (100x100km)
conn <- ruut::db_connect()
q <- sprintf("SELECT * FROM %s.%s", "maaamet", "epk200t")
cat(sprintf("\n-----------------\n%s\n\n", q))
epk200t <- sf::st_read(conn, query = q)
sf::st_geometry(epk200t) %>% plot(add = T, border = 3, lwd = 0.3, col = "#d3fffb")