|
|
@@ -0,0 +1,99 @@
|
|
|
+#' Piirkonnale andmebaasides olevate polügoonide teisendamine
|
|
|
+#'
|
|
|
+#' Etteantud piirkonna geomeetrilise piirjoone ('piir') ja selle joone piirikasti ('bb') järele leitakse nende aladega kaetud polügoonid. Andmed salvestatakse GPKG faili kihtidena.
|
|
|
+#'
|
|
|
+#' @param obj str Objekti nimi. Edaspidi on oluline ainult see nimi. Piirkonna geomeetrilist joont ei ole vaja lisada.
|
|
|
+#' @param gpkg_home path Salvestatavate GPKG faili asukoht.
|
|
|
+#' @return Uute GPKG andmebaasi kihtide 'piir_...' ja 'bb_...' loomine.
|
|
|
+#' @seealso [sf::st_read()], [sf::write_sf()],[sf::st_transform()],[ruut::gpkg_piirkonnale_ruudustike_lisamine()] ,[ruut::gpkg_piirkonnale_polygoonide_lisamine()],[ruut::gpkg_sellest_alustame_gpkg_loomist()],[ruut::gpkg_piirkonnale_joonte_lisamine()],[ruut::gpkg_piirkonnale_punktide_lisamine()],[ruut::gpkg_teisendame_polygoone()],[ruut::gpkg_teisendame_jooni()],[ruut::gpkg_teisendame_punkte()]
|
|
|
+#' @keywords GPKG, boundary box, EPSG:3301
|
|
|
+#' @export
|
|
|
+#' @examples
|
|
|
+#' \dontrun{
|
|
|
+#'
|
|
|
+#' gpkg_home <- "/tmp"
|
|
|
+#' obj <- "marja"
|
|
|
+#' gpkg_teisendame_polygoone(obj = obj, gpkg_home = gpkg_home)
|
|
|
+#'
|
|
|
+#' # Layers list.
|
|
|
+#' dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
|
|
|
+#' sf::st_layers(dsn = dsn)
|
|
|
+#' }
|
|
|
+gpkg_teisendame_polygoone <- function(obj = NULL, gpkg_home = "/tmp") {
|
|
|
+ dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
|
|
|
+ if (!file.exists(dsn)) {
|
|
|
+ cat(sprintf("\nSellist faili \"%s\" ei leitud.\n", dsn))
|
|
|
+ return(NULL)
|
|
|
+ }
|
|
|
+ ## Konfiguratsiooni muutujale väärtuste omistamine
|
|
|
+ conf <- ruut::get_config()
|
|
|
+ tmp_gpkg_file <- tempfile(fileext = ".gpkg")
|
|
|
+ conf$gpkg_home <- dirname(tmp_gpkg_file)
|
|
|
+ conf$gpkg_file <- gsub(paste0(conf$gpkg_home, "/"), "", tmp_gpkg_file)
|
|
|
+ conf$gpkg_file <- gsub(".gpkg", "", conf$gpkg_file)
|
|
|
+
|
|
|
+ conf$gpkg_home <- gpkg_home
|
|
|
+ conf$gpkg_file <- obj
|
|
|
+ tmp_gpkg_file_output <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = F)
|
|
|
+ tmp_gpkg_file_input <- ruut::construct_to_gpkg_output_file_str(conf = conf, is_input_str = T)
|
|
|
+
|
|
|
+ ## ====================== teisendused =====================
|
|
|
+ # Layers list
|
|
|
+ gpkg_info <- sf::st_layers(dsn = dsn)
|
|
|
+ layer_names <- gpkg_info$name
|
|
|
+ layer_names <- layer_names[!layer_names %in% c("bb", "bb_epk10t_grid", "bb_epk02t_grid", "bb_epk2t_grid", "bb_epk10t", "bb_epk02t")]
|
|
|
+ for (layer_name in layer_names) {
|
|
|
+ ruumiline_obj <- sf::read_sf(dsn = dsn, layer = layer_name)
|
|
|
+ ## Kontlrollime kas geomeetriline objekt on polügoon.
|
|
|
+ ## Arvutame sel juhul pindala.
|
|
|
+ is_polygon <- any(grepl("polygon", tolower(attributes(ruumiline_obj$geom)$class), fixed = TRUE))
|
|
|
+ if (is_polygon) {
|
|
|
+ ## Lõikame objekti ruudistikga tükkideks. Ühte ruutu võib jääda mitu tükki.
|
|
|
+ ## Peame need pärast ühendama
|
|
|
+ # ruut::qgis_algorithm_search_by_word("Intersection")
|
|
|
+ algorithm <- "native:intersection"
|
|
|
+ conf$gpkg_table <- sprintf("%s", layer_name)
|
|
|
+ input <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = T)
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer")
|
|
|
+ output <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = F)
|
|
|
+ conf$gpkg_table <- sprintf("%s", "bb_epk02t_grid")
|
|
|
+ overlay <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = T)
|
|
|
+ str <- sprintf("{ 'INPUT' : '%s', 'INPUT_FIELDS' : ['fid'], 'OUTPUT' : '%s', 'OVERLAY' : '%s', 'OVERLAY_FIELDS' : [], 'OVERLAY_FIELDS_PREFIX' : '' }", input, output, overlay)
|
|
|
+ cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
|
|
|
+ system(cmd)
|
|
|
+ ## Ühendame ruudus olevad pinnad
|
|
|
+ # ruut::qgis_algorithm_search_by_word("Dissolve")
|
|
|
+ algorithm <- "native:dissolve"
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer")
|
|
|
+ input <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = T)
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer_2")
|
|
|
+ output <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = F)
|
|
|
+ str <- sprintf("{ 'FIELD' : ['%s'], 'INPUT' : '%s', 'OUTPUT' : '%s' }", "id", input, output)
|
|
|
+ cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
|
|
|
+ system(cmd)
|
|
|
+ ## Lisame pindade pindalad ja perimeetri
|
|
|
+ # ruut::qgis_algorithm_search_by_word("attributes")
|
|
|
+ algorithm <- "qgis:exportaddgeometrycolumns"
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer_2")
|
|
|
+ input <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = T)
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer")
|
|
|
+ output <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = F)
|
|
|
+ str <- sprintf("{ 'CALC_METHOD' : 0, 'INPUT' : '%s', 'OUTPUT' : '%s' }", input, output)
|
|
|
+ cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
|
|
|
+ system(cmd)
|
|
|
+ ## Add field
|
|
|
+ ## Arvutame pindala protsendi ruudust
|
|
|
+ # ruut::qgis_algorithm_search_by_word("pythonfield")
|
|
|
+ algorithm <- "qgis:advancedpythonfieldcalculator"
|
|
|
+ conf$gpkg_table <- sprintf("%s", "test_layer")
|
|
|
+ input <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = T)
|
|
|
+ conf$gpkg_table <- sprintf("%s_epk02t", layer_name)
|
|
|
+ output <- ruut::construct_to_gpkg_output_file_str(conf = conf, geometry_field = "geom", is_input_str = F)
|
|
|
+ str <- sprintf("{ 'FIELD_NAME' : 'area_pc', 'FIELD_TYPE' : 1, 'FIELD_LENGTH' : 4, 'FIELD_PRECISION' : 2, 'GLOBAL' : '', 'FORMULA' : 'value = round($geom.area()/10000,4)', 'INPUT' : '%s', 'OUTPUT' : '%s' }", input, output)
|
|
|
+ cmd <- ruut::construct_qgis_output_result_to_better_format(str = str, algorithm = algorithm)
|
|
|
+ system(cmd)
|
|
|
+ }
|
|
|
+ }
|
|
|
+ ## Layers list
|
|
|
+ sf::st_layers(dsn = dsn)
|
|
|
+}
|
