210127_artikkel/geomeetria_teisendused/create_epk02t_grid.R

#  Võrgustike genereerimine ja salvestamine andmebaasi
#' Funktsioon gnereerib kogu Eestit katva ruudustike piirikasti
#' 'maaamet.epk200t_bb' sisse ruutvõrgustiku ruudu suurusega 100x100 m.
#' Võrgustik salvestatakse andmebaasi nime 'epk02t_grid' alla.
#' !!!!!!! Võrgustiku genereerimine ja salvestamine andmebaasi
#' võtab aega üle 12 tunni. Salvestatavate andmete maht ca 5 GB.

create_epk02t_grid <- function() {
  # ruut::qgis_algorithm_search_by_word(str = "v.mkgrid")
  algorithm <- "grass7:v.mkgrid"
  # cat(qgisprocess::qgis_show_help(algorithm = algorithm))
  if (!dir.exists("tmp/epk02t_grid")) dir.create("tmp/epk02t_grid", recursive = T)
  conf <- ruut::get_config()
  conf$table <- "epk02t_grid"
  conf$schema <- "maaamet"
  ## Loe instaleeritud tabelid andmebaasis
  installed_tables <- ruut::db_schema_tablenames(conf = conf)
  if (!conf$table %in% installed_tables) {
    # 1:40
    for (i in 1:40) {
      coordinates <- sprintf("300000,6%s0000 [EPSG:3301]", (29 + i))
      map <- sprintf("tmp/epk02t_grid/%s.gpkg", i)
      cat(sprintf("\n%s/40 1. coordinates: %s, 2. map: %s", i, coordinates, map))
      # ------------------
      result <- qgisprocess::qgis_run_algorithm(
        algorithm = algorithm,
        "-h" = 0,
        "-p" = 0,
        GRASS_OUTPUT_TYPE_PARAMETER = 3,
        GRASS_REGION_PARAMETER = "300000,800000,6320000,6700000  [EPSG:3301]",
        GRASS_VECTOR_DSCO = "",
        GRASS_VECTOR_EXPORT_NOCAT = 0,
        GRASS_VECTOR_LCO = "",
        angle = 0,
        box = "100,100",
        breaks = 0,
        coordinates = coordinates,
        grid = "100,5000", # "4000,5000"
        map = map,
        position = 0,
        .quiet = TRUE
      )
      # result
      ## Muudame rea numbreid, kuna tsükliga lähevad need kaduma.
      layers <- rgdal::ogrListLayers(map)
      x <- rgdal::readOGR(map, layers[1])
      x$row <- x$row + (i - 1) * 100
      x$cat <- rev(x$cat + (i - 1) * 100)
      rgdal::writeOGR(obj = x, dsn = map, layer = layers[1], driver = "GPKG", overwrite_layer = T)
      ## Lisame kihi andmebaasi
      cmd <- sprintf(
        paste0(
          "ogr2ogr -append --config PG_USE_COPY YES -f PostgreSQL ",
          "PG:\" dbname='%s' host=%s port=%d user='%s' password='%s' ",
          "sslmode=%s active_schema=%s \" -lco DIM=2 \"%s\" -nlt GEOMETRY ",
          "-lco GEOMETRY_NAME=geometry -lco FID=%s -nln %s.%s ",
          "-s_srs %s -t_srs %s -nlt %s"
        ),
        conf$dbname, conf$host, conf$port, conf$user, conf$password, conf$sslmode,
        conf$schema, map, "id", conf$schema, conf$table, "EPSG:3301", "EPSG:4326",
        "POLYGON"
      )
      cat(sprintf("\n-----------------\n%s\n", cmd))
      system(cmd)
    }
    cat(sprintf("\n-----------------\nAlgoritm: %s valmis.\n\n", algorithm))
    # result
  } else {
    cat(sprintf(
      "\n-----------------\nUue andmebaasi '%s' installimiseks pead vana enne kustutama.\n",
      conf$table
    ))
  }
}