R/GIS：シェープファイルを緯線のバウンディングボックスでサブセット化する方法は？

Question

シェイプファイル（.shpとそれに関連するファイルはhereです）を、ロング[80,90]とラット[20,30]の間の座標セットで区切られた別のものにサブセット化して、これを書きたいとします別のシェイプファイルとして。私はmaptoolsパッケージを使用する場合：

df = readShapeLines("/path/asia_rivers.shp")

をし、その後as.data.frame(df)とファイルの構造を見て、私は座標によってサブセットの任意の明白な方法を見つけることができません。

df = importShapefile("/path/asia_rivers.shp") 
df_sub = subset(df, X>=80 & X<=90 & Y >=20 & Y <=30) 

をその後私はmaptoolsにwriteSpatialShape()経由でエクスポートすることができデータフレームにこれを強制することができるように見えることはできません：私はサブセット化PBSmappingパッケージを使用することができます。私はこのエラーが発生し続ける：Error in PolySet2SpatialLines(df_sub) : unknown coordinate reference system。確かに私は何か非常に基本的なものを見逃していますし、地理座標によって地理データをサブセット化する簡単な方法があるはずですか？

Answer 1

次の操作を試みることができる：

library(rgeos) 
rivers <- readWKT("MULTILINESTRING((15 5, 1 20, 200 25), (-5 -8,-10 -8,-15 -4), (0 10,100 5,20 230))") 
bbx <- readWKT("POLYGON((0 40, 20 40, 20 0, 0 0, 0 40))") 

rivers.cut <- gIntersection(rivers, bbx) 

plot(rivers, col="grey") 
plot(bbx, add=T, lty=2) 
plot(rivers.cut, add=T, col="blue") 

Answer 2

私はこれが回答されているけど、私はあなたがPBSmappingを使用したい正確に何ができると思います。 PBSmappingは（ポリゴンやラインデータのための）ポリセットをクリップする機能を有しているので、あなたが試すことができます：

df <- importShapefile("/path/asia_rivers.shp") 
df_sub <- clipLines(df, xlim = c(80 , 90) , ylim = c(20 , 30), keepExtra = TRUE) 
dfSL <- PolySet2SpatialLines(df_sub) 

キープの余分はあなたが（私は属性データのために仮定）クリッピングを行うときに、非標準の列を維持することができます。

Answer 3

別の方法：

library(raster) 

s <- shapefile("/path/asia_rivers.shp") 

sub <- crop(s, extent(80, 90, 20, 30)) 

shapefile(sub, 'cropped_rivers.shp')