2014-02-18 24 views
0

我已經在python中看到了這個代碼,但是不知道如何在R中執行它。我有一些代碼,但它不適用於我。試圖將國家網格系統轉換爲lat long

我的數據有一列緯度和一列經度,但他們在國家電網系統。我不知道如何發佈數據集來幫助

我正在嘗試使用rgdal。我取得經濟和拉長,並做出一個拉特變量。代碼我發現應該有幫助的是這個

library(rgdal) 
whiskies$whiskies.coord = paste0(whiskies$Latitude, ", ", whiskies$Longitude) 
proj4string(whiskies$whiskies.coord) = CRS("+init=epsg:27700") # Specify that our coords are in osgb grid coord 
whiskies.coord <- spTransform(whiskies.coord, CRS("+init=epsg:4326")) # spTransform to convert osgb grid to lat/lon 

我不知道我是否正確設置whiskies.coord是基本問題。當我運行的第三個命令我得到這個錯誤

錯誤(函數(類,FDEF,mtable): 無法找到函數繼承的方法「proj4string < - 」簽字「‘性格’,‘CRS’ 「

行ID釀酒緯度經度 1阿伯費爾迪286580 749680 2亞伯勞爾326340 842570 3 AnCnoc 352960 839320

+0

有多大你的數據?如果它不是特別大,你可以編輯你的問題來發布dict(威士忌),否則後dput(威士忌[1:15,])。很高興有一個小數據可以玩。 – user20650

回答

0

我會假設你使用的是點數據。

我認爲你沒有將威士忌定義爲空間點數據框?

代替你的第二個線,我會用以下

coordinates(whiskies)<-~Longitude+Latitude 
proj4string(whiskies)<-CRS("+init=epsg:27700") 
相關問題