如何解决计算R中的栅格/ tiff标准差
我想知道如何计算tiff / raster文件的标准偏差值,以及如何从中找到值的范围?
背景信息: 我有一个喜马拉雅地区总初级生产力的tiff文件,我想知道该tiff文件的标准差值。假设该tiff文件的平均值为3.23,如果标准偏差为0.11,范围将为3.12-3.34吗?
我要附上一个栅格示例:
ras1 <- raster(matrix(c(1,1,2,2)))
解决方法
要从栅格计算任何统计信息,您需要像这样使用[]
:
> library(raster)
> ras1 <- raster(matrix(c(1,1,2,2)))
> xmean <- mean(ras1[])
> xsd <- sd(ras1[])
> xmean
[1] 1.5
> xsd
[1] 0.5477226
要获取范围,您可能可以做到
> c(xmean - xsd,xmean + xsd)
[1] 0.9522774 2.0477226
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。