To plot the bei dataset in Google Earth use:
library(spatstat) library(maptools) library(rgdal) data(bei) lon0lat0 <- data.frame(lat=9.15125, lon=-79.8553) coordinates(lon0lat0) <- ~lon+lat proj4string(lon0lat0) <- CRS("+proj=longlat +datum=WGS84") x0y0 <- data.frame(spTransform(lon0lat0, CRS("+init=epsg:32617"))) bei.pnt <- data.frame(x=x0y0[1,1]+bei$x, y=x0y0[1,2]+bei$y, no=rep(1, length(bei$x))) coordinates(bei.pnt) <- ~x+y proj4string(bei.pnt) <- CRS("+init=epsg:32617") bei.ll <- spTransform(bei.pnt, CRS("+proj=longlat +datum=WGS84")) writeOGR(bei.ll, "bei.kml", "bei", "KML")