'piirkonnale_ruudustike_lisamine' töötav versioon projektsiooniga 3301.

gpkg/01_piirkonnale_ruudustike_lisamine.R

@@ -61,11 +61,11 @@ piirkonnale_ruudustike_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "
     # .quiet = TRUE
   )
   result
-  pk_grid_bb <- sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT"))
-  sf::st_geometry(pk_grid_bb) %>% plot()
-  sf::st_geometry(pk) %>% plot(add = T)
+  # pk_grid_bb <- sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT"))
+  # sf::st_geometry(pk_grid_bb) %>% plot()
+  # sf::st_geometry(pk) %>% plot(add = T)
 
-  system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s", dsn, tmp_gpkg_file, output_layer_name))
+  system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s -t_srs \"EPSG:3301\"", dsn, tmp_gpkg_file, output_layer_name))
 
   ## ------------- 3. piirkonna epk10t (5x5 km) ruudud --------------------
   # 3.1 kogu ruutvõrgustik
@@ -77,7 +77,7 @@ piirkonnale_ruudustike_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "
     conf <- ruut::get_config()
     pg <- ruut::construct_ogr2ogr_PG_connect_str(conf = conf)
     input <- sprintf(
-      'postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=3301 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"%s\".\"%s\" (geometry)',
+      'postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"%s\".\"%s\" (geometry)',
       conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password,
       "maaamet", ruut
     )
@@ -93,10 +93,10 @@ piirkonnale_ruudustike_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "
       # .quiet = TRUE
     )
     result
-    assign(ruut, sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT")))
-    system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s", dsn, tmp_gpkg_file, output_layer_name))
+    # assign(ruut, sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT")))
+    system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s -t_srs \"EPSG:3301\"", dsn, tmp_gpkg_file, output_layer_name))
   }
-  sf::st_geometry(epk10t_grid) %>% plot(add = T)
+  # sf::st_geometry(epk10t_grid) %>% plot(add = T)
 
   # 3.2 ainult piirkonnaga seotud ning informatsiooni sisaldav ruutvõrgustik
   ruudud <- c("epk10t", "epk2t")
@@ -105,9 +105,9 @@ piirkonnale_ruudustike_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "
     ruut <- ruudud[j]
     output_layer_name <- ruut
     conf <- ruut::get_config()
-    pg <- ruut::construct_ogr2ogr_PG_connect_str(conf = conf)
+    # pg <- ruut::construct_ogr2ogr_PG_connect_str(conf = conf)
     input <- sprintf(
-      'postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=3301 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"%s\".\"%s\" (geometry)',
+      'postgres://dbname=\'%s\' host=%s port=%s user=\'%s\' sslmode=%s password=\'%s\' key=\'fid\' srid=4326 type=Polygon checkPrimaryKeyUnicity=\'1\' table=\"%s\".\"%s\" (geometry)',
       conf$dbname, conf$host, conf$port, conf$user, conf$sslmode, conf$password,
       "maaamet", ruut
     )
@@ -123,8 +123,8 @@ piirkonnale_ruudustike_lisamine <- function(obj = NULL, pk = NULL, gpkg_home = "
       # .quiet = TRUE
     )
     result
-    assign(ruut, sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT")))
-    system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s", dsn, tmp_gpkg_file, output_layer_name))
+    # assign(ruut, sf::read_sf(qgisprocess::qgis_output(result, "OUTPUT")))
+    system(sprintf("ogr2ogr -f GPKG -overwrite  %s %s -nln %s -s_srs \"EPSG:4326\" -t_srs \"EPSG:3301\"", dsn, tmp_gpkg_file, output_layer_name))
   }
   # sf::st_geometry(epk10t) %>% plot(add = T)
   ## ---------------------- vaata layer'id ----------------------