Question

我在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

我不知道我是否正确设置了whiskies.coord是基本问题。当我运行第三个命令时，我收到此错误

（函数（classes，fdef，mtable）中的错误：无法找到函数'proj4string＆lt; - '的签名'“字符”，“CRS”'

的继承方法

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

Answer 1

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

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

我将使用以下

代替你的第二行

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

试图将国家电网系统转换为长期转发

1 个答案: