210127_artikkel/geomeetria_teisendused/create_epk10t_grid.R

#  Võrgustike genereerimine ja salvestamine andmebaasi
#' Funktsioon gnereerib kogu Eestit katva ruudustike piirikasti
#' 'maaamet.epk200t_bb' sisse ruutvõrgustiku ruudu suurusega 5x5 km.
#' Võrgustik salvestatakse andmebaasi nime 'epk10t_grid' alla.

create_epk10t_grid <- function() {
  # ruut::qgis_algorithm_search_by_word(str = "v.mkgrid")
  algorithm <- "grass7:v.mkgrid"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
  result <- qgisprocess::qgis_run_algorithm(
    algorithm = algorithm,
    "-h" = 0,
    "-p" = 0,
    GRASS_OUTPUT_TYPE_PARAMETER = 0,
    GRASS_REGION_PARAMETER = "",
    GRASS_VECTOR_DSCO = "",
    GRASS_VECTOR_EXPORT_NOCAT = 0,
    GRASS_VECTOR_LCO = "",
    angle = 0,
    box = "5000,5000", # 5x5 km
    breaks = 0,
    coordinates = "300000,6300000 [EPSG:3301]",
    grid = "80,100", # ridade ja veergude arv
    map = "tmp/epk10t_grid.gpkg",
    position = 0,
    .quiet = TRUE
  )
  # result
  epk10t_grid <- sf::read_sf(qgisprocess::qgis_output(result, "map"))
  sf::st_geometry(epk10t_grid) %>% plot()
  conf <- ruut::get_config()
  conf$table <- "epk10t_grid"
  conf$schema <- "maaamet"
  ruut::copy_qgis_object_to_db(
    x = result, conf = conf, geometry_type = "POLYGON",
    crs_source = "EPSG:3301"
  )
  cat(sprintf("\n-----------------\nAlgoritm: %s valmis.\n\n", algorithm))
}
# create_epk10t_grid()