00_piirkonna_uute_geomeetriate_loomine.R 1.0 KB

123456789101112131415161718192021222324252627282930313233
  1. #' Uute geomeetriliste joonte genereerimine ja salvestamine andmebaasi.
  2. #'
  3. ## 1. Piirkonda katvate epk10t ruudud
  4. for (i in 1:length(objektid)) {
  5. cat(sprintf("\n-----------------\n%s\n\n", objektid[i]))
  6. result <- create_pk_epk10t(obj = objektid[i])
  7. epk10t <- sf::read_sf(qgis_output(result, "OUTPUT"))
  8. sf::st_geometry(epk10t) %>% plot()
  9. conf$table <- sprintf("%s_epk10t", objektid[i])
  10. ruut::copy_qgis_object_to_db(x = result, conf = conf, geometry_type = "POLYGON")
  11. }
  12. ## 2. Ortofotode allalaadimine
  13. # kaardiruutude nimekiri
  14. epk10t <- pk_epk10t(obj = objektid[i])
  15. epk10t_nr <- epk10t$nr
  16. j <- 2
  17. pattern <- sprintf("%s.tif", epk10t_nr[j])
  18. # pattern <- "62084.tif"
  19. ls <- list.files(path = "~/ortofotod/", pattern = pattern)
  20. if (length(ls) == 0) {
  21. estmap::maaamet_ortofotod(epk10t = epk10t_nr[j])
  22. }
  23. ## 3. Kaardiruutude epk200t piirikasti loomine
  24. create_epk200t_bb()
  25. ## 4. Kogu Eestit katva võrgustike (grid) loomine
  26. create_epk200t_grid() # 100x100 km
  27. create_epk10t_grid() # 5x5 km
  28. create_epk2t_grid() # 1x1 km
  29. create_epk02t_grid() # 100x100 m