使用折线提取栅格 - 去除角

如何解决使用折线提取栅格 - 去除角

我有以下栅格

r <- raster(nrows=6,ncols=7,xmn=0,xmx=7,ymn=0,ymx=6,crs="+proj=utm +units=m")

r[] <- c(2,2,1,5,8,7,3,9,9)

我想沿着以下折线提取单元格值

pts <- data.frame(x= c(2.5,4.5),y=c(3.5,5.5))
coordinates(pts) <- c("x","y")
proj4string(pts) <- crs(r)
lns <- spLines(pts,crs=crs(r))

绘制一切

plot(r)
points(pts)
lines(lns)

提取单元格值

extract(r,lns,cellnumbers=T)
[[1]]
     cell layer
[1,]    5     5
[2,]   11     8
[3,]   12     5
[4,]   17     1
[5,]   18     8

问题：尽管该线仅通过 3 个单元格，但仍有 5 个值。我知道这是因为提取函数还提取线穿过的 2 个角处的值。但是，我不想包括角点。我该怎么做？

我失败的尝试描述如下。提前致谢。

尝试 1：使用坐标提取

extract(r,coordinates(lns)[[1]][[1]],cellnumbers=T)
    cells layer
[1,]    17     1
[2,]     5     5

它只提取用于创建 lns 的点，因此它缺少创建线时未包含的单元格，因此它不起作用。

尝试 2：使用掩码

我想我可以创建一个蒙版栅格，然后从该蒙版中提取值。然而，掩码也会提取与角对应的单元格，如下所示

lns.m <- mask(r,lns)
plot(lns.m)