R/GIS: Wie Teilmenge einer shape-Datei durch eine lat-long-bounding-box?
Ich will Teilmenge aus einer shape-Datei (die .shp und die zugehörigen Dateien sind hier) in eine andere, begrenzt durch eine Reihe von Koordinaten, sagen wir zwischen longs [80,90] und lats [20,30], und dann schreiben Sie dies als ein weiteres shapefile. Wenn ich die maptools
Paket:
df = readShapeLines("/path/asia_rivers.shp")
anschaut und dann die Struktur der Datei mit as.data.frame(df)
ich finde keine offensichtliche Möglichkeit, Untergruppen von Koordinaten. Ich kann Verwendung der PBSmapping
Paket Teilmenge:
df = importShapefile("/path/asia_rivers.shp")
df_sub = subset(df, X>=80 & X<=90 & Y >=20 & Y <=30)
aber dann kann ich nicht scheinen, um in der Lage sein zu zwingen, diese in eine SpatialLines
Daten-Frames, die exportiert werden können, über writeSpatialShape()
im maptools
. Ich bekomme immer diese Fehlermeldung: Error in PolySet2SpatialLines(df_sub) : unknown coordinate reference system
. Ja, ich bin fehlt etwas sehr basic und sollte es eine einfache Möglichkeit, ausgewählte geo-Daten geo-Koordinaten?
Du musst angemeldet sein, um einen Kommentar abzugeben.
Können Sie Folgendes versuchen:
Ich weiß, das beantwortet wurde, aber ich denke, man kann genau das tun, was Sie wollen mit
PBSmapping
.PBSmapping
hat eine Funktion, clip Polysets (für die polygon-und Linien-Daten), so Sie versuchen können:Die halten extra ermöglicht es Ihnen, um nicht-standard-Spalten, wenn Sie die clipping (ich nehme an, für die Attribut-Daten).
Einen anderen Weg: