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

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米：

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米）。我可以通过弄乱点的大小来手动执行此操作（例如，使用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)