函数rather：：getData的子集状态看起来非常奇怪

xfb7svmp 于 2023-06-19 发布在其他

关注(0)|答案(1)|浏览(69)

我使用raster::getData获得Large SpatialPolygonsDataFrame的详细信息

states <- raster::getData("GADM", country = "United States", level = 1

我把IL和WI分了...

statesilwi <- states[states$NAME_1 == "Illinois" | states$NAME_1 == "Wisconsin",]

但是当我绘制它的时候，它伊利诺斯州看起来有点时髦，为什么会发生这种情况，有没有办法解决它？

plot(statesilwi)

来源：https://stackoverflow.com/questions/76450905/subsetted-states-look-very-weird-with-function-rastergetdata

1条答案

按热度按时间

lyfkaqu11#

“看起来有点时髦”不是很清楚。但在我看来，@rawr是正确的。密歇根湖的一部分属于伊利诺伊州。
你可以像这样删除密歇根湖：

library(terra)
library(geodata)
us <- geodata::gadm("United States", level = 2, path=".")
us <- us[us$TYPE_2 != "Water body", ]
s <- us[us$NAME_1 %in% c("Illinois", "Wisconsin"),]
s <- aggregate(s, "NAME_1")
plot(s)

赞(0）回复(0）举报 2023-06-19

我来回答

函数rather：：getData的子集状态看起来非常奇怪

1条答案

相关问题

热门标签

最新问答