我从 随附的栅格文件开始，该文件是 5 个气候区的世界地图。然后我将其转换为一个简单的功能。我的目标是去掉那些包含在较大多边形中的小多边形，这样我就只能看到五个不同区域的较大斑块的轮廓。作为奖励，使这些线条更平滑会很好，因为边缘在视觉上并不令人愉悦，但作为一个分类量，我无法真正插入这些值。
library(stars)
x=read_stars('climi.tif')
#convert into sf
tes=st_as_sf(x,as_points=FALSE,merge=TRUE,use_integer=TRUE)
#load a simple background map
map.background=sf::st_as_sf(rworldmap::countriesCoarseLessIslands)
ggplot()+ geom_sf(data=map.background,size=0.1) +geom_sf(data=tes,aes(geometry=geometry,color=factor(climi.tif)),fill=NA,size=0.7)
我明白了，但是正如您所看到的，在非常大的多边形内有许多几个像素的小多边形，这使得可视化变得混乱。我想将这些小多边形与围绕它们的较大多边形合并。
我试图计算面积并删除较小的多边形，然后使用 sf_remove_holes，这可以改善结果，但我仍然得到太多较小的多边形或 NA 区域，这些区域应该被着色。
tes$area= st_area(tes);tes$area=as.numeric(tes$area)
subtes=subset(tes,tes$area>199426036400)
subtes=sfheaders::sf_remove_holes(subtes)
但是一些小多边形仍然可见，另外一些区域似乎不属于任何区域。我也尝试过 st_contains 和 raster::clump() 但没有成功。例如，Clump 删除了五个类别之一。我想形成易于区分的连续区域，即使我必须绘制双边框。
