csv到kml文件新西兰

时间:2013-10-07 23:57:25

标签: r csv gis kml google-earth

我正在尝试从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

0 个答案:

没有答案