如何将矩阵转换为光栅文件，以及如何计算R中不同区域的平均值

我有一个数据集，它是

dat<-matrix(seq(1:5290),c(115,46))

如果我想计算每个州的平均值，我应该怎么做？还是特定地区

非常感谢您的反馈

---

您可以尝试以下方法：

library(raster)
library(maps)

dat<-matrix(seq(1:5290),c(115,46))  

取整个光栅的平均值

mean(getValues(r))

区域平均数

r.sub <- crop(r, c(-100, -54, 40, 45))
mean(getValues(r.sub))

并提取平均值：

dat <- extract(r, poly)
sapply(dat, mean)

---

dat您可以尝试以下方法：

library(raster)
library(maps)

dat<-matrix(seq(1:5290),c(115,46))  

取整个光栅的平均值

mean(getValues(r))

区域平均数

r.sub <- crop(r, c(-100, -54, 40, 45))
mean(getValues(r.sub))

并提取平均值：

dat <- extract(r, poly)
sapply(dat, mean)

dat这很有效。有没有一种方法，我们可以在地图上添加代表不同生态区域的多边形，计算不同生态区域的平均值，然后用平均值映射生态区域。这正是我想要的。我将使用中的shapefile，我将测试是否可以完成它。谢谢。ecoregions.shp您的proj4string建议不同的预测。你能在哪里画出它们？绘图是否匹配？我想在绘图（r）的顶部绘制它们。他们彼此不匹配。这就是形状文件的来源，我如何在同一个绘图上绘制它们？这非常有用。有没有一种方法，我们可以在地图上添加代表不同生态区域的多边形，计算不同生态区域的平均值，然后用平均值映射生态区域。这正是我想要的。我将使用中的shapefile，我将测试是否可以完成它。谢谢。ecoregions.shp您的proj4string建议不同的预测。你能在哪里画出它们？绘图是否匹配？我想在绘图（r）的顶部绘制它们。他们彼此不匹配。这就是形状文件的来源，如何在同一个绘图上绘制它们？