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?

InformationsquelleAutor user702432 | 2012-03-22
Schreibe einen Kommentar