我已经围绕一系列坐标点创建了矩形缓冲区,我想做的就是使用2种情况在这个蓝色“足迹”周围创建多边形;
1)在所有至少接触一个顶点的矩形周围绘制多边形 2)在所有至少接触两个顶点的矩形周围绘制多边形
例如,在此处的图像中,场景1中的左下角第二个蓝色正方形将包含在方案1中,但方案2中不包含。
这是我用来绘制数据的代码,对不起,它是不可复制的;
dat$ref <- as.factor(dat$ref) # stipulates whether present/absent blue/red
mysf <- st_as_sf(dat, coords = c("long","lat")) # make df into spatial object
st_crs(mysf) = 4326 # set coordinate reference system
present <- mysf %>% filter(ref==1)
absent <- mysf %>% filter(ref==0)
pt_buffer <- st_buffer(present, dist = 0.025,endCapStyle="SQUARE")
pt2_buffer <- st_buffer(absent, dist = 0.025,endCapStyle="SQUARE")
m <- ggplot() + geom_sf(data = pt_buffer,col="blue",fill="blue") +
geom_sf(data = pt2_buffer,col="red",fill="red")
m + coord_sf(xlim = c(-2, 0), ylim = c(60, 62))
#inspect initial polygon output for number of features
pt_buffer #e.g., features: 4198
在以前的问题中,我找不到任何涉及矩形和连接/触摸点数的东西。