Päevakoopia.

00_algandmed.R

@@ -22,28 +22,28 @@ 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]]}')
+# 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"
+obj <- "marja"
 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)
-# ## GPKG joonte lisamine
-# ruut::gpkg_piirkonnale_joonte_lisamine(obj = obj, gpkg_home = gpkg_home)
-## GPKG punktide lisamine
+## 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)
+## GPKG joonte lisamine
+ruut::gpkg_piirkonnale_joonte_lisamine(obj = obj, gpkg_home = gpkg_home)
+# GPKG punktide lisamine
 ruut::gpkg_piirkonnale_punktide_lisamine(obj = obj, gpkg_home = gpkg_home)
 
 
@@ -94,31 +94,32 @@ 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)) {
+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, 3301)
-  # sf::st_crs(pk)
+  x <- pk_piir(obj = obj)
 
-  sf::st_geometry(pk) %>% plot()
+  sf::st_geometry(x) %>% 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)
-# }
+  ## 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)
+  ## GPKG joonte lisamine
+  ruut::gpkg_piirkonnale_joonte_lisamine(obj = obj, gpkg_home = gpkg_home)
+  # GPKG punktide lisamine
+  ruut::gpkg_piirkonnale_punktide_lisamine(obj = obj, gpkg_home = gpkg_home)
+
+
+
+  }
 
 
 m <- piirkonnale_maatriksi_loomine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)

help/join_attributes_by_location_(summary).txt

@@ -0,0 +1,19 @@
+
+
+## Võtsin genereeritud 5x5 km ruudud ja lugesin (count) palju nendes ruutudes on elamuid (osm_id).
+
+Processing algorithm…
+Algorithm 'Join attributes by location (summary)' starting…
+Input parameters:
+{ 'DISCARD_NONMATCHING' : True, 'INPUT' : 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable authcfg=0i30k14 key=\'ogc_fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"maaamet\".\"epk10t_genereeritud_ruudud\" (geometry)', 'JOIN' : 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable authcfg=0i30k14 key=\'fid\' srid=4326 type=MultiSurface checkPrimaryKeyUnicity=\'1\' table=\"osm_shp\".\"buildings_a\" (geometry)', 'JOIN_FIELDS' : ['osm_id'], 'OUTPUT' : 'TEMPORARY_OUTPUT', 'PREDICATE' : [1], 'SUMMARIES' : [0] }
+
+
+## Selle kaardi rasterdamine
+
+Processing algorithm…
+Algorithm 'Rasterize (vector to raster)' starting…
+Input parameters:
+{ 'BURN' : 255, 'DATA_TYPE' : 0, 'EXTENT' : '21.776430335,28.230611099,57.488980215,59.831538657 [EPSG:4326]', 'EXTRA' : '-l wrgn', 'FIELD' : 'osm_id_count', 'HEIGHT' : 5000, 'INIT' : None, 'INPUT' : 'postgres://dbname=\'data\' host=localhost port=6432 user=\'osm\' sslmode=disable authcfg=0i30k14 key=\'ogc_fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"results\".\"epk10t_count_buildings\" (geometry)', 'INVERT' : False, 'NODATA' : 0, 'OPTIONS' : '', 'OUTPUT' : '/home/ardo/Documents/doktoriõpe/210127 - Valga artikkel/R/210127_artikkel/tmp/tiff/majade_arv_5x5_km_ruudus.tif', 'UNITS' : 0, 'WIDTH' : 5000 }
+
+GDAL command:
+gdal_rasterize -l results.epk10t_count_buildings -a osm_id_count -ts 5000.0 5000.0 -a_nodata 0.0 -te 21.776430335 57.488980215 28.230611099 59.831538657 -ot Byte -of GTiff  "PG:dbname='data' host=localhost port=6432 user='osm' sslmode=disable password='osm'" "/home/ardo/Documents/doktoriõpe/210127 - Valga artikkel/R/210127_artikkel/tmp/tiff/majade_arv_5x5_km_ruudus.tif"