210127_artikkel/gpkg/10_piirkonnale_maatriksi_loomine.R

#' Maatriksi loomine
#'
#' Etteantud piirkonna geomeetrilise piirjoone ('piir') järele leitakse selle piirkonna EPSG:3301 projektsioonile vastav piirikast ('boundarybox_3301'). Piirikasti järele leitakse selle kastiga kaetud 5x5 km ('epk10t_grid'), 1x1 km ('epk2t_grid') ja 100x100 meetrit ('epk02t_grid') ruudustikud. Samuti leitakse piirkonna piiri sees olevad ja statistilist ning muud infot sisaldavad 5x5 km ('epk10t') ja 1x1 km ('epk2t') sisaldavad ruudu. Andmed salvestatakse GPKG faili kihtidena. Faili nimeks on objekti nimi.
#'
#' @param obj str Objekti nimi.
#' @param pk Piirkonna geomeetriline joon.
#' @param gpkg_home path Salvestatavate GPKG faili asukoht.
#' @param obj_only TRUE/FALSE Kas geomeetriad leitakse ainult objekti või kogu piirikastiga määratud ala jaoks.
#'
#'
piirkonnale_maatriksi_loomine <- function(obj = NULL, pk = NULL, gpkg_home = "/tmp", obj_only = TRUE) {
  if (!kas_on_piirkond(obj = obj, pk = pk)) {
    return()
  }


  ## Algandmed
  dsn <- sprintf("%s/%s.gpkg", gpkg_home, obj) # Polygon
  layer <- "bb_epk02t_grid"
  dsn_info <- sf::st_layers(dsn = dsn)

  ## 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("\n%s", layer_name))
    assign(layer_name, sf::read_sf(dsn = dsn, layer = layer_name))
  }

  # Valge piirikast
  sf::st_geometry(bb) %>% plot(border = 0, lwd = 0.3, col = "#ffffff")
  # Must piir
  sf::st_geometry(piir) %>% plot(add = T, border = 3, lwd = 0.3, col = "#63fffb")



  sf::st_geometry(bb) %>% plot(border = 3, lwd = 0.3, col = "#d3fffb")
  sf::st_geometry(bb_epk02t_grid) %>% plot(border = 3, lwd = 0.3, col = "#d3fffb")
  sf::st_geometry(piir_epk02t_grid) %>% plot(add = T, border = 3, lwd = 0.3, col = "#a3fffb")
  sf::st_geometry(piir) %>% plot(add = T, border = 3, lwd = 0.3, col = "#63fffb")

  sp <- sf::read_sf(dsn = dsn, layer)
  sp$value <- 0
  unique(sp$left)
  unique(sp$bottom)
  length(sp$value)
#
#   pk <- sf::st_transform(pk, sf::st_crs(sp))

  sf::st_geometry(sp) %>% plot()
  sf::st_geometry(pk) %>% plot(add = T)

  # 0 - intersect
  # 1 - contain
  # 2 - disjoint
  # 3 - equal
  # 4 - touch
  # 5 - overlap
  # 6 - are within
  # 7 - cross

  x <- sf::st_contains(pk, sp, sparse = TRUE)
  sp$value[x[[1]]] <- 1
  x[[1]]
  max(sp$value)
  sp
  if (length(x[[1]]) == 0) {
    cat("\n--------------\nAndmebaasis olev ruudustik ei kattu piirkonnaga.\nPalu genereeri eelnevalt piirikast\n funktsiooniga 'piirkonnale_ruudustike_lisamine()'.\n\n")
    return()
  }
  # plot(sp)
  # sf::st_geometry(sp) %>% plot()
  rows <- sort(unique(sp$bottom), decreasing = f)
  cols <- sort(unique(sp$left))
  length(rows)*length(cols)

  (m0 <- matrix(0L, nrow = length(rows), ncol = length(cols), byrow = F)) # 0-maatriks
  (m <- matrix(sp$value, nrow = length(rows), ncol = length(cols), byrow = F))

  > mat1.data <- c(1,2,3,4,5,6,7,8,9)
  > mat1 <- matrix(mat1.data,nrow=3,ncol=3,byrow=TRUE)
  > mat1
  # colnames(m) <- paste0("X",cols)
  # rownames(m) <- paste0("Y",rows)
  # m["802","3480"]<- 100
  cat("\n---------------\nTagastati maatriks.\n")
  m
}
# piirkonnale_maatriksi_loomine(obj = obj, pk = pk, gpkg_home = gpkg_home, obj_only = TRUE)