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")
# Objektide piirjooned

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

## 1. piirkonna piir
pk <- pk_piir(obj = objektid[i])
sf::st_geometry(pk) %>% plot()
## 2. piirkonna epk10t ruudud
epk10t <- pk_piir_bb(obj = objektid[i])
sf::st_geometry(epk10t) %>% plot(border = 3, lwd = 0.3, col = "#d3fffb")
sf::st_geometry(pk) %>% plot(add = T)
# 2.1 kaardiruutude nimekiri
epk10t_nr <- epk10t$nr

## 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")



## 2. boundarybox
# ----------------- Piirkonna BB leidmine -----------------
## Piirkonna 'pk' Boundary Box.
result <- ruut::pk_bbox(epk10t, type = 0)
pk_bbox <- sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT"))
sf::st_geometry(pk_bbox) %>% plot(border = 2)
sf::st_geometry(epk10t) %>% plot(add = T, border = "darkblue", lwd = 0.5, lty = 1.8)

# -------------- export bbox to postgis database ---------------
conf$table <- "bbox"
ruut::copy_polygon_to_db(x = result, conf = conf, id = "fid")