R/GIS: Как создать подмножество шейп-файла по ограничивающей рамке по широте и долготе?

Я хочу разделить шейп-файл (.shp и связанные с ним файлы здесь) в другой, ограниченный набором координат, скажем, между longs [80,90] и lats [20,30], а затем запишите это как еще один шейп-файл. Если я использую пакет maptools:

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

, а затем посмотрю структуру файла с as.data.frame (df), я не могу найти никакого очевидного способа подмножества по координатам. Я могуиспользовать пакет PBSmappingдля подмножества:

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

, но тогда я не могу преобразовать это во фрейм данных SpatialLines, который может экспортировать через writeSpatialShape()в maptools. Я продолжаю получать эту ошибку: Ошибка в PolySet2SpatialLines(df_sub): неизвестная система координат . Конечно, я упускаю что-то очень простое, и должен быть простой способ подмножества геоданных по геокоординатам?