我正在使用光栅包来显示我已转换为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）。简而言之，我想知道代理在地图上的移动位置，以及最常使用的区域，因此网格单元格的值越高，代理程序使用的越多。
我有许多地图，每个地图都有完全相同的坐标系，但每个栅格地图中的代理移动（以及网格单元格值）分布不同，我想比较每个地图以显示代理移动地图之间存在统计差异。
这可能吗？
由于