У меня есть набор координат точки в зоне 48 UTM.

  x           y       
615028.3  2261614    
615016.3  2261635    
614994.4  2261652    

CSV-файлВот.

Я хотел бы загрузить CSV и создать шейп-файл, используя R. Мой код:

library(maptools)
library(rgdal)
library(sp)

    UTMcoor=read.csv(file="https://dl.dropboxusercontent.com/u/549234/s1.csv")
    coordinates(UTMcoor)=~X+Y
    proj4string(UTMcoor)=CRS("++proj=utm +zone=48") # set it to UTM
    LLcoor<-spTransform(UTMcoor,CRS("+proj=longlat")) #set it to Lat Long
    plot(LLcoor)
    points(LLcoor$X,LLcoor$Y,pch=19,col="blue",cex=0.8) #to test if coordinate can be plot as point map
    writeOGR(UTMcoor, dsn="c:/todel" ,layer="tsb",driver="ESRI Shapefile")
    writeSpatialShape("LLcoor","test")

В последней команде (writeSpatialShape) R выдают следующую ошибку:

Error in writeSpatialShape("LL2", "test") : 
  x is acharacterobject, not a compatible Spatial*DataFrame

Когда я читаю LLcoor из консоли, кажется, что это уже Spatial DataFrame. Запись файла формы с использованием writeOGR (пакет RGdal) также дает похожую ошибку. Любой намек очень ценится.