我正在尝试从csv文件生成kml文件,以便可以在Google地球中绘制。我正在使用我在另一个帮助部分找到的代码,它似乎适用于kml的生成,但绘制的点是错误的。
library("sp")
library("rgdal")
data(meuse)
coordinates(meuse) <- c("x", "y")
proj4string(meuse) <- CRS("+init=epsg:28992")
meuse_ll <- spTransform(meuse, CRS("+proj=longlat +datum=WGS84"))
writeOGR(meuse_ll["zinc"], "meuse.kml", layer="zinc", driver="KML")
第一个解决方法是摆脱小数位,这使得我从法国到新西兰,但积分仍然大幅下降。他们不是新西兰奥克兰的北部,而是查塔姆群岛的南部。我怀疑我正在处理投射失误匹配或软糖因素问题。我也试过从这里复制zip文件:http://trac.osgeo.org/proj/ 要使用NZGD49,希望这样可行,但它会引发这个错误:
meuse_ll2 <- spTransform(test, CRS("+proj=longlat +datum=NZGD49"))
#Error in spTransform(test, CRS("+proj=longlat +datum=NZGD49")) :
# error in evaluating the argument 'CRSobj' in selecting a method for function 'spTransform': Error in CRS("+proj=longlat +datum=NZGD49") :
#unknown elliptical parameter name
所以我有什么想法可以解决?
Date Time Latitude1 Longitude1
27/09/2013 12:18:02 -36.67593 174.429688
27/09/2013 12:18:03 -36.675892 174.42981