使用tmap绘制sf数据时出现“二进制运算符的非数字参数”错误

时间:2019-03-27 12:59:26

标签: r sf tmap

当我尝试使用tm_symbols显示点数据时,出现错误消息: non-numeric argument to binary operator

我已经剥离了代码以试图找出问题所在,当然我也搜索了tmap和其他文档。

与其他人的某些链接正在做我想要做的事情:

另请参阅:

这是我的代表:

library(sf)
library(tmap)
library(leaflet)

item_data <- data.frame(
    name=c("Epping Forest District Citizens Advice (Epping)","Epping Forest District Citizens Advice (Loughton)","Epping Forest District Citizens Advice (Waltham Abbey)"),
    latitude=c("51.696921", "51.649158", "51.687181"),
    longitude=c("0.110474", "0.05899", "-0.004736"),
    stringsAsFactors = FALSE
)
items_sf <- st_as_sf(item_data, coords=c("longitude", "latitude"), crs=3857)

tmap_mode("view")
epmap <- tm_basemap(leaflet::providers$Stamen.TonerBackground) +
  tm_shape(items_sf, name="CA Locations") +
  tm_symbols(shape=21)
epmap

这给了我

## Error in b[3:4] - b[1:2] : non-numeric argument to binary operator

我正在尝试使用tmap,正如我所建议的那样,但是我想我也可以尝试使用其他方法来生成地图……如果我这样做:

plot(items_sf)

...它给出错误:

## Error in r[i1] - r[-length(r):-(length(r) - lag + 1L)] : non-numeric argument to binary operator

如果我这样做:

library(mapview)
mapview(items_sf)

...我得到了一个绘制了三个点的地图视图,但是总范围小于一米,因此由于某种原因,坐标不被处理为纬度和经度。

我很乐于解决问题,但是我真的很困在这里,因为我不知道如何处理这些错误消息。

我希望将三个位置的tmap图绘制为底图上的点(点/符号)。实际结果:错误消息,未呈现任何地图。

**编辑: 好吧,引用数字错误对我来说是很愚蠢的,这是受访者的好地方。这是由于我输入数据框,而不是简单地从我实际使用的数据框中复制而造成的。 解决该问题后,我的脚本仍然存在其他错误,但最终得到了解决。

projection / ESPG很有帮助,因为我还不太了解这些内容,并且基本上在猜测该怎么做。所以我在那里也学到了一些东西。 **

1 个答案:

答案 0 :(得分:1)

您的代码似乎有两个问题:

  • 您的坐标存储为文本(即不是数字)
  • 您使用的度量标准CRS(3857)的坐标在以十进制度(角度单位)查看时更有意义

请仔细考虑以下代码(删除引号并将CRS从3857更改为4326 +更改颜色)

library(sf)
library(tmap)
library(leaflet)

item_data <- data.frame(
  name=c("Epping Forest District Citizens Advice (Epping)","Epping Forest District Citizens Advice (Loughton)","Epping Forest District Citizens Advice (Waltham Abbey)"),
  latitude=c(51.696921, 51.649158, 51.687181),
  longitude=c(0.110474, 0.05899, -0.004736),
  stringsAsFactors = FALSE
)

items_sf <- st_as_sf(item_data, coords = c("longitude", "latitude"), crs = 4326)

tmap_mode("view")
epmap <- tm_shape(items_sf, name="CA Locations") + tm_symbols(shape = 21, col = "red") +
  tm_basemap(leaflet::providers$Stamen.TonerBackground)

epmap

enter image description here