R：删除地图地图中的岛屿答案

【问题标题】：R: Remove island in a tmap mapR：删除地图地图中的岛屿
【发布时间】：2020-10-14 04:08:38
【问题描述】：

我正在尝试使用 tmap 移除岛屿。

我经常使用 lims 来删除 ggplot 中的岛屿：

library(tmap)
library(ggplot2)
library(raster)

chile  <- getData("GADM",country="Chile",level=1)

ggplot()+
    geom_polygon(data=chile, aes(x=long, y=lat, group = group),
         color = "#666666") +
         coord_equal() + 
         lims(x = c( -77.2,-60), y = c(-60, -15))

ggplot Chile map

但是 tmap 不允许使用 lims 和 island 从大陆干扰地图：

tm_shape(chile) +
tm_borders()

tmap Chile with islands

我将 facet 与 tmap 一起使用，因为它可以为我提供更好的结果，但孤立它对我来说是个问题：

    chile2= chile[c(2,4:5,7,8,12:16),]
    chile2$group=1
    chile3= chile[c(1,3,6,9:11),]
    chile3$group=2
    map_split=rbind(chile2, chile3)

    tm_shape(map_split) +
    tm_borders() +
    tm_facets(by = "group")

tmap Chile facet

您可以看到岛屿不允许获得好的构面图。

你知道如何移除 tmap 地图中的岛屿吗？

谢谢

【问题讨论】：

标签： r facet tmap

【解决方案1】：

您可以从数据集中删除不需要的区域。

library(raster)
chile  <- getData("GADM",country="Chile",level=1)
xchile <- crop(chile, extent(-85,-59, -61, -14))

运行需要一段时间，因为智利是一个相当复杂的数据集，因为它有很多岛屿

【讨论】：