Início / coding / Perguntas / 79203135
Accepted
chill
Asked: 2024-11-19 18:59:23 +0800 CST

Como faço para amostrar pontos dentro de vários polígonos carregados de um arquivo kml no R?

  • 772

Tenho um arquivo kml que descreve quatro polígonos (o arquivo kml é chamado de "core.kml", e coloquei uma cópia neste link ). Gostaria de escolher pontos aleatórios dentro da área definida pelos polígonos. Idealmente, gostaria de escolher pontos uniformemente em toda a área, então escolher mais pontos dos polígonos maiores do que dos pequenos, em proporção à sua área. Peguei emprestado o código da resposta 1 aqui , que funciona perfeitamente com o polígono da Carolina do Norte (pelo menos quando o atualizo inserindo nc <- st_geometry(nc)). Mas com meus polígonos, a st_sampleetapa sempre dá este erro:

st_as_s2(): descartando a coordenada Z e/ou M Erro em wk_handle.wk_wkb(wkb, s2_geography_writer(oriented = oriented, :
O loop 0 não é válido: A aresta 15 cruza a aresta 19

Eu tentei inserir sf_use_s2(FALSE), mas ainda recebo o mesmo erro. O plotcomando funciona, e os polígonos parecem como esperado. Meu código:

core<-st_read("core.kml")
core<-st_geometry(core) 
sf_use_s2(FALSE) 
core.points<-st_sample(core,5) 
plot(core)

1 respostas

  1. Best Answer

    Se você quiser garantir que haja pontos em cada polígono, talvez você queira tentar algo como

     library(sf)
 sf_use_s2(FALSE) 
#> Spherical geometry (s2) switched off
 core = st_read("/Users/cara/Desktop/core.kml") |> st_geometry() |> st_make_valid() 

 core 
#> Geometry set for 4 features 
#> Geometry type: GEOMETRY
#> Dimension:     XY
#> Bounding box:  xmin: -79.06133 ymin: 33.69862 xmax: -78.90788 ymax: 33.8351
#> Geodetic CRS:  WGS 84
#> MULTIPOLYGON (((-79.014 33.80998, -79.00692 33....
#> POLYGON ((-79.06133 33.82931, -79.05514 33.8306...
#> POLYGON ((-78.99585 33.78962, -78.99496 33.7901...
#> POLYGON ((-78.91353 33.70372, -78.91197 33.7052...
 
 ppwa = \(sfc, pts) { 
   # points per weighted area 
   a = st_area(sfc)
   p = proportions(a)
   as.numeric(p * pts) 
 }
 
 # ensure sum(n) == pts, see below
 floor_preserve = \(x) {
   fx = floor(x)
   i = tail(order(x - fx), round(sum(x)) - sum(fx))
   fx[i] = fx[i] + 1
   fx
 }
 
 n = ppwa(core, 25) |> floor_preserve()
 n 
#> [1] 14  2  6  3
 sum(n)
#> [1] 25
 
 set.seed(1119)
 p = st_sample(core, n)
#> although coordinates are longitude/latitude, st_intersects assumes that they
#> are planar
#> ... gets repeated ...

 plot(core) 
 points(p, pch = 4, col = palette.colors()[7])
#> Warning in st_centroid.sfc(x, of_largest_polygon = of_largest_polygon):
#> st_centroid does not give correct centroids for longitude/latitude data

    Observe que (1) não corrigimos problemas com core. Dê uma olhada nos avisos. Esta pode ser uma boa discussão para começar sua pesquisa; ilhas são difíceis; (2) o número total de pontos não é aleatório, pois o definimos como 25. O número de pontos por polígono, portanto, não é aleatório randônico: ele é 25multiplicado pela área relativa de cada polígono, ou seja, ponderado.

    • 1

relate perguntas