题
如何将像素值添加到绘图中?我可以通过使用 click()
但我希望它出现在情节中。
library(raster)
r <- raster(nrow=3, ncol=3)
r[] <- 1:ncell(r)
plot(r)
click(r)
解决方案
试试下面的方法,它是基于由返回的函数拼凑而成的
getMethod("click", signature="Raster")
.
myClick <- function(x, n = Inf, id = FALSE, xy = FALSE, cell = FALSE,
type = "n", show = TRUE, ...) {
i <- 0
n <- max(n, 1)
while (i < n) {
i <- i + 1
loc <- locator(1, type, ...)
xyCoords <- cbind(x = loc$x, y = loc$y)
cells <- na.omit(cellFromXY(x, xyCoords))
if (length(cells) == 0)
break
value <- extract(x, cells)
text(xyCoords, labels = value)
}
}
## Try it out
myClick(r, n=4)
其他提示
如果你想显示所有的值,你可以使用 text
方法:
library(raster)
r <- raster(nrow=3, ncol=3, vals=1:9)
plot(r)
text(r)
对于一个子集,你可以做类似的事情:
z <- rasterToPoints(r, function(x) x > 6 )
plot(r)
text(z[,1], z[,2], z[,3])
我知道这个问题已经被标记为answered,但是基于Josh的解决方案和Eddie的后续问题,这里有一个小小的for-loop,它做了Eddie所要求的(绘制没有十进制数的栅格值并且不 click
):
r <- raster(nrow=3, ncol=3)
r[] <- runif(ncell(r))
plot(r)
for(i in 1:ncell(r)){
xycoords <- xyFromCell(r, cell = i)
value <- extract(r, xycoords)
text(xycoords, labels = round(value))
}
不隶属于 StackOverflow