0
我用utm32-shapefile制作了一张等值线地图,现在我想用另一个数据框的数据添加第二层。这第二个数据帧只有Lon/Lat信息。 我想将Lon/Lat转换为utm32坐标,所以我可以将它们与ggplot2一起用作另一个 - 会聚层。如何使用ggplot2绘制utm32中的经度/纬度坐标进行绘图?
数据帧是这样的:
GPS_Lat GPS_Lon Index
51,133 14,683 12.75
首先,我用gsub
重新编写了 “” 到 “”然后我将这些变量转换为数字。
#Then I tried to define the variables as Lon/Lat-Coordinates
cord.dec<-SpatialPointsDataFrame(cbind(plot.eeg$GPS_Lon,plot.eeg$GPS_Lat),data=plot.eeg,proj4string=CRS("+proj=longlat"))
#Then i tried to convert them to utm32
cord.UTM<-spTransform(cord.dec,CRS("+init=epsg:4326"))
它没有工作。我只是得到我的纬度/经度变量的副本。
我没有找到一个关于如何做到这一点的好文档。也许别人可以帮助我?
您的CRS字符串不正确:EPSG:4326是WGS84 longlat(请参阅http://spatialreference.org/ref/epsg/wgs-84/);试试EPSG:32632 – rcs 2015-04-02 18:32:33
是的,就是这样!谢谢! – SEMson 2015-04-03 06:52:14