2015-04-02 109 views
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")) 

它没有工作。我只是得到我的纬度/经度变量的副本。

我没有找到一个关于如何做到这一点的好文档。也许别人可以帮助我?

+1

您的CRS字符串不正确:EPSG:4326是WGS84 longlat(请参阅http://spatialreference.org/ref/epsg/wgs-84/);试试EPSG:32632 – rcs 2015-04-02 18:32:33

+0

是的,就是这样!谢谢! – SEMson 2015-04-03 06:52:14

回答

0

您的CRS字符串不正确:EPSG:4326是WGS84 longlat(请参阅here);试试EPSG:32632。

相关问题