2015-04-02 94 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。

相關問題