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

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)

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

Plot Circle on Raster with Known Coordinates and Radius

plot

r

raster

r-raster