R栅格绘制图像，绘制圆圈并屏蔽圆形外的像素

Question

下面的代码绘制图像，然后在该图像上绘制圆圈。我想让所有落在该圈外的像素变黑。我怎么能这样做？

data-url

Answer 1

使用str（）查看＆＃39; x＆＃39; -object，你会看到：

..@ data    :Formal class '.MultipleRasterData' [package "raster"] with 14 slots
  .. .. ..@ values    : num [1:2500, 1:3] 255 248 221 199 198 210 221 190 104 79 ...
  .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. .. ..$ : NULL
  .. .. .. .. ..$ : chr [1:3] "red" "green" "blue"

....所以1:50乘1:50的值被映射到三列。 X值可能为0:2500 %% 50，y值可能为0:2500 %/% 50

所以记住＆＃34;起源＆＃34;如果栅格对象的左上角，但是绘图函数的左下角，以及20的y值变为50-20或30，这将使您接近所要求的内容（抱歉将y序列放在第一位） ）：

x@data@values[( ((1:2500 %/% 50 )- 30)^2 + ((1:2500 %% 50) - 20)^2 ) >=100, 1] <- 0
 x@data@values[( ((1:2500 %/% 50 )- 30)^2 + ((1:2500 %% 50) - 20)^2 ) >=100, 2] <- 0
 x@data@values[( ((1:2500 %/% 50 )- 30)^2 + ((1:2500 %% 50) - 20)^2 ) >=100, 3] <- 0
 plotRGB(x)
 draw.circle(20,20,10,border="blue")

逻辑是标准的形式为（x-dx）^ 2 +（y-dy）^ 2> r ^ 2其中dx和dy是圆的中心坐标，r是半径== 10。

问题改变后的编辑：

对于每个颜色图层，具有命名参数的代码将类似于使得最暗的＆＃34;红色＆＃34;的代码。这给出了一个大致圆形的蒙版，虽然没有正确处理排列中心：

x@data@values[( ((1:(height*width) %/% (height*5/4) )- (height-circley*5/4) )^2 + 
            ((1:(height*width) %% width) -       circlex )^2 ) >= radius^2, 1] <- 0

进一步的实验表明它非常接近：

x@data@values[( ((1:(height*width) %/% (height) )- (height-circley) *5/4)^2 + 
                 ((1:(height*width) %% width) -       circlex )^2 ) >= radius^2, 1] <- 0
plotRGB(x, asp=5/4, addfun=function() draw.circle(circlex,circley,radius,border="blue") )

显然，你可以用width/height比例因子替换5/4出现的新宽高比。

Answer 2

这是一个更通用的解决方案，也适用于具有不同坐标系的栅格。

函数rasterToPoints()将栅格转换为点。按照您的示例，head(rasterToPoints(x))返回以下内容：

> head(rasterToPoints(x))
       x    y red green blue
[1,] 0.5 39.5 255   255  251
[2,] 1.5 39.5 204   205  199
[3,] 2.5 39.5 171   172  164
[4,] 3.5 39.5 157   159  148
[5,] 4.5 39.5 162   164  151
[6,] 5.5 39.5 187   189  176

然后我们需要找到哪些点落在圆圈之外并将它们的值设置为零：

is_outside_circ = (rasterToPoints(x)[,1] - circlex)^2 + (rasterToPoints(x)[,2] - circley)^2 >= radius^2
x@data@values[ is_outside_circ,] <- 0
plotRGB(x)
draw.circle(circlex,circley,radius,border="blue")

Result: black points outside circle