微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

删除R中较大多边形内的小多边形

如何解决删除R中较大多边形内的小多边形

我从 attached raster 文件开始,它是 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)

我明白了,但是正如您所看到的,在非常大的多边形内有许多由几个像素组成的小多边形,这使可视化变得混乱。我想将这些小多边形与围绕它们的较大多边形合并。

enter image description here

我尝试计算面积并移除较小的多边形,然后使用 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 会删除五个类别之一。我想形成易于相互区分的连续区域,即使我必须绘制双边框。

enter image description here

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?