SpatialPolygons - Erstellung einer Menge von Polygonen in R aus Koordinaten

Ich versuche mich zu nehmen, erstellen Sie eine Reihe von Polygonen, die aus vertex-Standorte, gespeichert in X,Y format.

Hier ist ein Beispiel meiner Daten - jede Zeile stellt die Eckpunkte für ein polygon. die Polygone sind Vierecke

square <- rbind(c(255842.4, 4111578, 255862.4, 4111578, 255862.4, 4111558, 
                  255842.4, 4111558, 255842.4, 4111578, 255842.4, 4111578),
                c(257397.0, 4111309, 257417.0, 4111309, 257417.0, 4111289, 
                  257397.0, 4111289, 257397.0, 4111309, 257397.0, 4111309))

ID <- c("SJER1", "SJER2")'

Ich bin mit SpatialPolygonsdamit meine Daten müssen in einer Liste. so habe ich eine Schleife, um zu versuchen, um meine Daten in Form einer Liste, die aus einer matrix.

Ich eine Schleife erstellen, die folgenden code, fand ich in einigen anderen Fragen auf dieser Website. Ich brach jeden Schritt, um zu versuchen zu verstehen, warum bin ich nur immer ein polygon als Ausgabe auch durch ich habe 2 sets von Punkten.

for (i in 1:2) {  
  pts <- rbind(c(square[i,1], square[i,2]), c(square[i,3], square[i,4]), 
               c(square[i,5],square[i,6]), c(square[i,7],square[i,8]), 
               c(square[i,9],square[i,10]))
  sp1 <- list(Polygon(pts))
  sp2 <- list(Polygons(sp1,i))
  sp = SpatialPolygons(sp2)  
}
plot(sp)

Können Sie mir bitte helfen zu verstehen, wie ich das einstellen den code zu schreiben, aus zwei Polygone anstelle von nur einem? Und auch, wie ich weisen Sie die ID für jedes polygon gegeben, ich bin mit einer matrix (quadratisch) als meine ab-dataset, und wenn ich das zuweisen einer Charakter-id, verwandelt es alle meine Daten in ein Zeichen.

Mein Ziel ist es zwei Polygone in der SpatialPolygons Objekt, die erste mit der ID SJER1 und die zweite mit der ID SJER2 gespeichert in der SpatialPolygons Objekt.

Dann Schreibe ich dies aus einem shapefile.

InformationsquelleAutor der Frage L Wasser | 2014-10-28

r rgdal sp spatial

25

Gibt es einige Informationen zu ?'SpatialPolygons-class'aber Sie mehr oder weniger möchten, tun Sie Folgendes:
```
polys <- SpatialPolygons(list(
  Polygons(list(Polygon(matrix(square[1, ], ncol=2, byrow=TRUE))), ID[1]),
  Polygons(list(Polygon(matrix(square[2, ], ncol=2, byrow=TRUE))), ID[2])
))

plot(polys)
```
Der Grundtenor ist, dass Sie brauchen, um zu erstellen Polygon Objekte (z.B. vom 2-Spalten-Matrizen mit x - Koordinaten in die erste Spalte und die y Koordinaten in der zweiten). Diese sind kombiniert Listen zu erstellen Polygons Objekte (von denen jeder sollte haben eine eindeutige ID). Diese Polygons Objekte werden zusammengefasst in einer Liste, um eine SpatialPolygons Objekt. Sie können eine CRS wenn Sie möchten, sich mit der proj4string argument SpatialPolygons (siehe ?SpatialPolygons).

Schreiben es sich um ein ESRI-shape-Datei benötigen Sie zu konvertieren, um ein SpatialPolygonsDataFrame Objekt durch die Kombination der polys Objekt, das wir erstellt und einige Daten. Wir werden nur hinzufügen, dass die IDs-Daten für den Mangel an etwas interessanter.
```
polys.df <- SpatialPolygonsDataFrame(polys, data.frame(id=ID, row.names=ID))
```
und dann schreiben Sie es auf...
```
writeOGR(polys.df, '.', 'fancysquares', 'ESRI Shapefile')
```
Das zweite argument ('.') sagt, es zu schreiben, aus dem aktuellen Arbeitsverzeichnis.

BEARBEITEN

Schnell eine SpatialPolygonsDataFrame wenn Sie viele Zeilen beschreiben Polygone, könnten Sie Folgendes verwenden:
```
# Example data
square <- t(replicate(50, {
  o <- runif(2)
  c(o, o + c(0, 0.1), o + 0.1, o + c(0.1, 0), o)
}))
ID <- paste0('sq', seq_len(nrow(square)))

# Create SP
polys <- SpatialPolygons(mapply(function(poly, id) {
  xy <- matrix(poly, ncol=2, byrow=TRUE)
  Polygons(list(Polygon(xy)), ID=id)
}, split(square, row(square)), ID))

# Create SPDF
polys.df <- SpatialPolygonsDataFrame(polys, data.frame(id=ID, row.names=ID))

plot(polys.df, col=rainbow(50, alpha=0.5))
```
InformationsquelleAutor der Antwort jbaums

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.