210127_artikkel/geomeetria_teisendused/create_pk_epk10t.R

#' Funktsioon: piirkonna piiri katvate epk10t (5x5km) ruutude loomine
#' Valitud objekti järele leitakse milliste ruutude 'maaamet.epk10t' sees
#' asub valitud piirkond.

create_pk_epk10t <- function(obj) {
  # obj - objekti nimetus (näiteks: valga)
  # ruut::qgis_algorithm_search_by_word(str = "extract")
  algorithm <- "native:extractbylocation"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
  result <- qgisprocess::qgis_run_algorithm(
    algorithm = algorithm,
    INPUT = 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable password=\'osm\' key=\'fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"maaamet\".\"epk10t\" (geometry)',
    INTERSECT = sprintf('postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable password=\'osm\' key=\'ogc_fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"xxx_artikkel_210127\".\"%s_piir\" (geometry)', obj),
    OUTPUT = qgisprocess::qgis_tmp_vector(),
    PREDICATE = c(0, 1),
    .quiet = TRUE
  )
  result
}
# create_pk_epk10t(obj = objektid[i])