复杂矩阵操作

Question

我的数据是这样的，X和Y是缺陷的中心。我想在矩阵中指定缺陷。

defect_ID X Y X_range Y_range
1         2 3  5        0
2         7 12 1        2

我创建了一个只有0的矩阵200*200。我想通过以下方式将1放入矩阵中：

每个坐标X Y都是1。例如，我们可以看到ID 1，它将允许1到坐标（2,3）的单元格。ID 2将允许1进入我的手机（7,12）。

我已经用代码完成了这项工作

img<-matrix(0,200,200)
img[cbind(data[,X],data[,Y])]<-1

现在我想做一些棘手的事情。我defect_ID，我想使用我的X_range和Y_range值将值1分配给这个缺陷周围的所有点。对于1个缺陷，我们取中心（X，Y），并使用X和Y_range移动。

伪码

for defect_ID 1
img[X,Y]<-1
img[X+ 1:X_range,Y]<-1
img[X,Y+1:Y_range]<-1

例如，如果我们看ID 1，我想给点（2,3）、（3,3）、（4,3）、（5,3）、（6,3）、（7,3）赋值1

如果我有X_范围和Y_范围

Answer 1

这里有一个选择：

data <- tibble(defect_ID=c(1,2),X=c(2,7),Y=c(3,12),X_range=c(5,1),Y_range=c(0,2))

输入数据：

> data
# A tibble: 2 x 5
  defect_ID     X     Y X_range Y_range
      <dbl> <dbl> <dbl>   <dbl>   <dbl>
1         1     2     3       5       0
2         2     7    12       1       2

创建矩阵：

data_mat <- matrix(0,200,200)

分配一个：

for (i in 1:2) { # Change this to the length of your dataset
  data_mat[data[[i,'X']]:data[[i,'X']]+data[[i,'X_range']],data[[i,'Y']]] <- 1
  data_mat[data[[i,'X']],data[[i,'Y']]:data[[i,'Y']]+data[[i,'Y_range']]] <- 1
}

输出：

> data_mat[1:10,1:14]
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14]
 [1,]    0    0    0    0    0    0    0    0    0     0     0     0     0     0
 [2,]    0    0    1    0    0    0    0    0    0     0     0     0     0     0
 [3,]    0    0    1    0    0    0    0    0    0     0     0     0     0     0
 [4,]    0    0    1    0    0    0    0    0    0     0     0     0     0     0
 [5,]    0    0    1    0    0    0    0    0    0     0     0     0     0     0
 [6,]    0    0    1    0    0    0    0    0    0     0     0     0     0     0
 [7,]    0    0    1    0    0    0    0    0    0     0     0     1     1     1
 [8,]    0    0    0    0    0    0    0    0    0     0     0     1     0     0
 [9,]    0    0    0    0    0    0    0    0    0     0     0     0     0     0
[10,]    0    0    0    0    0    0    0    0    0     0     0     0     0     0

Answer 2

img = matrix(0,20,15)
for (i in 1:nrow(data)){
  sx = data$X[i]:(data$X[i]+data$Y_range[i])
  sy = data$Y[i]:(data$Y[i]+data$X_range[i])

  img[sx,sy] = 1
}

这至少应该适用于第一部分。第二个问题我不太清楚。如果范围是奇数呢？我应该选择哪个方向？