我正在使用 raster 包来显示已转换为 ASCII 文件的纺织品。
使用 ggplot2 显示地图我没有问题。我拥有的地图包含一个 506 * 608 的网格。每个网格单元都有一个与之关联的值。这是地图中包含的唯一信息，以及 x 和 y 坐标。
我在 R 中使用以下代码显示了地图：
Map <- raster("filename.txt")
map.spdf <- as(Map, "SpatialPixelsDataFrame")
map.df <- as.data.frame(map.spdf)
head(map.df)
gMap <- ggplot(map.df, aes(x=x, y=y)) + geom_tile(aes(fill=layer)) + coord_equal()
gMap <- gMap + theme(panel.background = element_rect(fill='white'))
gMap <- gMap + scale_fill_gradient2(limits = c(0, 60000), low=muted("green"), mid="red", high="white", midpoint=30000)
gMap
我的地图中的每个网格单元值对应于我在基于代理的模拟中模拟的多个代理在该位置的移动量（范围从 0 到 20,000）。总之，我想知道agent在地图上移动了哪些地方，哪些区域使用频率最高，所以grid cell的值越高，agent使用的越多。
我有许多地图，每个地图都具有完全相同的坐标系，但在每个栅格地图中具有不同的代理移动分布（以及网格单元值），我想比较每个地图以显示代理移动在统计上是否不同地图。
这可能吗？
谢谢