我从这里下载了绘制爱荷华州的shapefile - https://geodata.iowa.gov/dataset/county-boundaries-iowa/resource/183e782f-2d43-4073-8524-fe8e634cf17a
但是,当我尝试打印地图时,我只能获得底图,而不能获得县级别的数据。
mymap<-st_read("county.shp")
str(mymap)
library(leaflet)
# map the polygons in shp
mymap %>%
leaflet() %>%
addTiles() %>%
addPolygons()
此外,我还收到以下警告-
Warning messages:
1: sf layer is not long-lat data
2: sf layer has inconsistent datum (+proj=utm +zone=15 +datum=NAD83 +units=m +no_defs).
Need '+proj=longlat +datum=WGS84'
我是一个在这个领域的初学者,将非常感谢任何在这个问题上的帮助。
发布于 2020-07-25 03:33:14
使用@jazzurro的评论作为未来搜索者的答案-
mymap<-st_read("county.shp")%>% st_transform(crs = 4326)
str(mymap)
library(leaflet)
# map the polygons in shp
mymap %>%
leaflet() %>%
addTiles() %>%
addPolygons()
显然,我忘了指定导致问题的坐标参考系统。
https://stackoverflow.com/questions/63082820
复制相似问题