在已知坐标和半径的栅格上绘制圆

Question

我正在尝试将具有定义半径的圆绘制到光栅图像上。我已经成功地绘制了我的栅格并将点放在图像上：

# open image to save    
png("RasterImage.png",
        width=10, height=10, units="in", res=144)

# plot raster
    plot(VI.SOS.mean, main="VI Mean SOS",
         col=terrain.colors(length(seq(100,220,20))-1), axes=F, breaks=seq(100,200,20))

# add points
    points(sensors$X, sensors$Y)

# close png file
    dev.off()

生成的图像如下所示：

我的栅格中的像素是 30 m：

class       : RasterLayer 
dimensions  : 2871, 3205, 9201555  (nrow, ncol, ncell)
resolution  : 30, 30  (x, y)
extent      : 254265, 350415, 4731885, 4818015  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs 
data source : in memory
names       : layer 
values      : 38.6, 309  (min, max)

我想做的是在每个点上放置一个半径为 x 的圆（例如 300 m）。我可以通过弄乱点的大小来手动执行此操作（例如，使用cex），但是如果我更改图像的尺寸或裁剪我的光栅，那么圆圈不再是正确的大小。有没有办法使用栅格中定义的单位在地图上放置圆圈？

谢谢！

Answer 1

你可以试试这个：

library(dismo)
cs <- circles(sensors[, c('X','Y')], d=150)
plot(VI.SOS.mean)
plot(polygons(cs), add=TRUE)

重叠的圆圈被合并。如果你不想这样，你可以使用内部函数 dismo:::.generateCircles

pls <- dismo:::.generateCircles(sensors[, c('X','Y')], d=150, lonlat=FALSE, crs=NA)
plot(pls, add=TRUE)