根据标签的网格的颜色正方形

Question

我已经使用以下代码创建了100x100网格：

nb=100
nb1=nb
nb2=nb
minlat=45.18306
maxlat=45.6325
binlat=seq(minlat-0.05, maxlat+0.05, length.out = nb1)
minlon=8.724444
maxlon=9.519722
binlon=seq(minlon-0.05, maxlon+0.05, length.out = nb2)

我的数据如下：

n = 100
coord = matrix(0, n, 2)
coord[,1]= runif(100, min=45.18306, max=45.6325)
coord[,2]= runif(100, min=8.724444, max=9.519722)

lab=matrix(0, n, 1)
m=1
for(i in 1:nb1){
  for(j in 1:nb2){
    lab[coord[,1]>binlat[i] & coord[,1]<=binlat[i+1] & coord[,2]>binlon[j] & coord[,2]<=binlon[j+1]]=m
    m=m+1
  }
}

假设我有一个如下：

a <- c(1000,2536,3045,6654)
b <- c(1,1,2,3)
matrix <- cbind(a,b)

where a是向量识别网格平方ID和b是广场所属的组的向量。

是否可以根据a根据其属于的组彩色的平方可视化网格？ a中未包含的所有正方形必须是非彩色/白色的。

Answer 1

您可以使用explion.grid获取所有构成网格的唯一x和y值的数据框架。

df <- expand.grid(x = binlon, y = binlat)

这将为您提供每个y值的所有x值，以下是：

head(df)
#>          x        y
#> 1 8.674444 45.13306
#> 2 8.683487 45.13306
#> 3 8.692530 45.13306 
#> 4 8.701574 45.13306
#> 5 8.710617 45.13306
#> 6 8.719660 45.13306

...等于10,000行。

组标记正确的单元格，请访问：

df$group <- NA_real_
df$group[matrix[,'a']] <- matrix[,'b']
df$group <- factor(df$group)

如果您的a列用此订购代表单元格的矩阵位置，则可以通过执行以下的 explion.grid中的参数

绘制结果可能最容易使用ggplot

library(ggplot2)

ggplot(df, aes(x, y, fill = group)) + 
  geom_tile(colour = 'gray90') +
  scale_fill_brewer(palette = 'Set1', na.value = 'white')