Algandmete koopia.

00_algandmed.R

@@ -2,11 +2,7 @@
 #'
 system("ogr2ogr --long-usage")
 
-# Processing algorithm…
-# Algorithm 'Count points in polygon' starting…
-# Input parameters:
-#   { 'CLASSFIELD' : '', 'FIELD' : 'NUMPOINTS', 'OUTPUT' : 'ogr:dbname=\'/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/valga_results.gpkg\' table=\"count_aadressandmed_ee_epk02t\" (geom)', 'POINTS' : '/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/valga_p.gpkg|layername=aadressandmed_ee', 'POLYGONS' : 'memory://Polygon?crs=EPSG:3301&field=fid:long(0,0)&field=id:integer(0,0)&field=cat:double(0,0)&field=row:double(0,0)&field=col:integer(0,0)&field=NUMPOINTS:double(0,0)&uid={b5c46c2b-5663-417c-83d4-6811b2ca538b}', 'WEIGHT' : '' }
-str <- "{ 'INPUT' : '/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/valga.gpkg|layername=piir', 'OUTPUT' : 'ogr:dbname=\'/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/valga_l.gpkg\' table=\"piir_l\" (geom)' }"
+str <- "{ 'INPUT' : 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable authcfg=0i30k14 key=\'id\' srid=3301 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"maaamet\".\"epk2t_grid\" (geom)', 'INTERSECT' : '/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa/marja.gpkg|layername=piir', 'OUTPUT' : 'TEMPORARY_OUTPUT', 'PREDICATE' : [5] }"
 ruut::construct_qgis_output_result_to_beter_format(str = str)
 
 source("functions/delete_existing_variables.R")
@@ -23,6 +19,65 @@ conf <- ruut::get_config()
 conf$schema <- "xxx_artikkel_210127"
 ruut::db_create_new_schema(conf = conf)
 
+## --------------- testi algus 'ruudu' funktsioonidele ------------------
+geojson_txt <- paste('{"type":"MultiPoint","coordinates":[[658300,6474800],[658300,6475000],[658700,6476000],[659600,6475000],[659600,6474800],[658300,6474800]]}')
+# read geojson from string:
+geojson_txt <- paste('{"type":"MultiPoint","coordinates":[[658300,6474800],
+[658300,6475000],[658400,6475100],[658600,6475000],[658600,6474800],
+[658300,6474800]]}')
+x <- sf::read_sf(geojson_txt) %>% sf::st_cast("POLYGON")
+sf::st_crs(x) <- 3301
+x
+sf::st_geometry(x) %>% graphics::plot()
+
+gpkg_home <- "/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa"
+obj <- "mikihiir"
+dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj)
+## GPKG loomine
+ruut::gpkg_sellest_alustame_gpkg_loomist(pk = x, obj = obj, gpkg_home = gpkg_home)
+# Layers list
+sf::st_layers(dsn = dsn)
+## GPKG ruudustike lisamine
+ruut::gpkg_piirkonnale_ruudustike_lisamine(obj = obj, gpkg_home = gpkg_home)
+## GPKG polügoonide lisamine
+ruut::gpkg_piirkonnale_polygoonide_lisamine(obj = obj, gpkg_home = gpkg_home)
+
+
+## Loeme andmebaasist piiri ja piirikasti.
+# Layers list
+gpkg_info <- sf::st_layers(dsn = dsn)
+layer_names <- gpkg_info$name
+for (layer_name in layer_names) {
+  cat(sprintf("\nFailist %s laeti alla: %s", dsn, layer_name))
+  assign(layer_name, sf::read_sf(dsn = dsn, layer = layer_name))
+}
+## Test plot
+sf::st_geometry(bb_epk02t_grid) %>% graphics::plot(border = 0, lwd = 0.3, col = "#ffffff") # valge
+sf::st_geometry(bb) %>% graphics::plot(add = T, border = 3, lwd = 0.3, col = "#d3fffb")
+sf::st_geometry(piir) %>% graphics::plot(add = T, border = 3, lwd = 0.3, col = "#a3fffb")
+
+sf::st_geometry(bb_epk02t_grid) %>% graphics::plot(border = 0, lwd = 0.3, col = "#ffffff") # valge
+sf::st_geometry(bb) %>% graphics::plot(add = T, border = 3, lwd = 0.3, col = "#d3fffb")
+sf::st_geometry(bb_epk02t_grid) %>% graphics::plot(add = T, border = 3, lwd = 0.3)
+sf::st_geometry(bb_shp_katastriyksus) %>% graphics::plot(add = T, border = 2, lwd = 0.3)
+
+sf::st_geometry(bb_epk02t_grid) %>% graphics::plot(border = 0, lwd = 0.3, col = "#ffffff") # valge
+sf::st_geometry(piir) %>% graphics::plot(add = T, border = 3, lwd = 0.3, col = "#d3fffb")
+sf::st_geometry(bb_epk02t_grid) %>% graphics::plot(add = T, border = 3, lwd = 0.3)
+sf::st_geometry(piir_shp_katastriyksus) %>% graphics::plot(add = T, border = 2, lwd = 0.3)
+
+
+
+
+
+
+
+##------------------ testi lõpp -------------------
+
+
+
+
+
 
 # Objektide nimekiri
 objektid <- c("valga", "matsalu", "lahemaa")
@@ -33,28 +88,38 @@ conn <- ruut::db_connect(conf = conf)
 ## Valitud objekti indeks
 i <- 1
 source("geomeetria_teisendused/pk_piir.R")
+source("gpkg/61_kas_on_piirkond.R")
 source("gpkg/01_piirkonnale_ruudustike_lisamine.R")
 source("gpkg/02_piirkonnale_polygoonide_lisamine.R")
 source("gpkg/03_piirkonnale_punktide_lisamine.R")
 source("gpkg/04_piirkonnale_joonte_lisamine.R")
 
-for (i in 1:length(objektid)) {
-  ## ---------------- 1. piirkonna piir ------------------
+
+
+# for (i in 1:length(objektid)) {
+#   ## ---------------- 1. piirkonna piir ------------------
   ## Muutujad: pk - piirkond
   obj <- objektid[i]
   pk <- pk_piir(obj = obj)
-  pk <- sf::st_transform(pk, sf::st_crs("EPSG:3301"))
+  # pk <- sf::st_transform(pk, 3301)
   # sf::st_crs(pk)
 
   sf::st_geometry(pk) %>% plot()
   gpkg_home <- "/data/gpkg/artiklid/artikkel_210127_valga_matsalu_lahemaa"
-  piirkonnale_ruudustike_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
-  piirkonnale_polygoonide_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
-  piirkonnale_punktide_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
-  piirkonnale_joonte_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
-}
-
-
+#
+#
+#   piirkonnale_ruudustike_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
+#   piirkonnale_polygoonide_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
+#   piirkonnale_punktide_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
+#   piirkonnale_joonte_lisamine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
+# }
+
+
+m <- piirkonnale_maatriksi_loomine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)
+m_filename <- sprintf("minu_%s_maatriks.csv", obj)
+write.table(m,file=m_filename, sep = ";", row.names = T)
+y <- read.table(m_filename,header=TRUE, sep = ";")
+unlink(m_filename) # Kustuta fail