210127_artikkel/00_piirkonna_uute_geomeetriate_loomine.R

#' Uute geomeetriliste joonte genereerimine ja salvestamine andmebaasi.
#'

source("01_funktsioonid.R")

## 1. Piirkonda katvate epk10t ruudud
for (i in 1:length(objektid)) {
  cat(sprintf("\n-----------------\n%s\n\n", objektid[i]))
  result <- create_pk_epk10t(obj = objektid[i])
  epk10t <- sf::read_sf(qgis_output(result, "OUTPUT"))
  sf::st_geometry(epk10t) %>% plot()
  conf$table <- sprintf("%s_epk10t", objektid[i])
  ruut::copy_qgis_object_to_db(x = result, conf = conf, geometry_type = "POLYGON")
}
## 1. Piirkonda katvate epk2t ruudud
for (i in 1:length(objektid)) {
  cat(sprintf("\n-----------------\n%s\n\n", objektid[i]))
  result <- create_pk_epk2t(obj = objektid[i])
  epk2t <- sf::read_sf(qgis_output(result, "OUTPUT"))
  # sf::st_geometry(epk2t) %>% plot()
  conf$table <- sprintf("%s_epk2t", objektid[i])
  ruut::copy_qgis_object_to_db(x = result, conf = conf, geometry_type = "POLYGON")
}

## 2. Ortofotode allalaadimine
# kaardiruutude nimekiri
epk10t <- pk_epk10t(obj = objektid[i])
epk10t_nr <- epk10t$nr
j <- 2
pattern <- sprintf("%s.tif", epk10t_nr[j])
# pattern <- "62084.tif"
ls <- list.files(path = "~/ortofotod/", pattern = pattern)
if (length(ls) == 0) {
  estmap::maaamet_ortofotod(epk10t = epk10t_nr[j])
}

## 3. Kaardiruutude epk200t piirikasti loomine
create_epk200t_bb()

## 4. Kogu Eestit katva võrgustike (grid) loomine
create_epk200t_grid() # 100x100 km
create_epk10t_grid() # 5x5 km
create_epk2t_grid() # 1x1 km
# create_epk02t_grid() # 100x100 m

## 5. Piirkonna piirikastide + vastavate ruudustike loomine
for (i in 1:length(objektid)) {
  cat(sprintf("\n-----------------\n%s\n\n", objektid[i]))
  create_boundarybox_3301_to_pk(obj = objektid[i])
  create_epk10t_grid_to_pk_bb(obj = objektid[i])
  create_epk2t_grid_to_pk_bb(obj = objektid[i])
}

## 6. Piirkonna piirikastide maatriksite loomine