0

我已经在 python 中看到了此代码，但无法弄清楚如何在 R 中执行此操作。我有一些代码，但它对我不起作用。

我的数据有一个纬度列和一个经度列，但它们在国家网格系统中。我不知道如何发布数据集以提供帮助

我正在尝试使用 rgdal。我取 lat 和 long 并制作一个 lat, long 变量。我发现应该有帮助的代码是这个

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

我不知道我是否正确设置了威士忌。坐标是基本问题。当我运行第三个命令时，出现此错误

(function (classes, fdef, mtable) 中的错误：无法为签名 '"character", "CRS"' 的函数 'proj4string<-' 找到继承的方法</p>

RowID 酒厂纬度经度 1 Aberfeldy 286580 749680 2 Aberlour 326340 842570 3 AnCnoc 352960 839320

4

1 回答 1

0

我将假设您正在使用点数据。

我认为您还没有将威士忌定义为空间点数据框？

代替您的第二行，我将使用以下内容

coordinates(whiskies)<-~Longitude+Latitude
proj4string(whiskies)<-CRS("+init=epsg:27700")
于 2014-03-05T00:57:03.167 回答