代码之家 › 专栏 › 技术社区 › BetterCallMe

使用R中光栅包中的分位数函数计算光栅堆栈的第一个和第三个四边形

r-raster r

BetterCallMe · 技术社区 · 7 年前

我使用以下代码在R中使用光栅包中的分位数函数进行四分位数计算。我只想要输出,不考虑任何数据值。我尝试了以下代码:

    library(raster)    
    raster_data<-list.files(path= getwd() , pattern= "\\.tif$", all.files=FALSE,
    full.names=TRUE,recursive=TRUE) 
    s <- stack(raster_data)
    quantile(s, na.rm = TRUE)

输出:

分位数

                   0% 25% 50% 75% 100%

    X2008001h25v06 NA  NA  NA  NA   NA

    X2008002h25v06 NA  NA  NA  NA   NA

    X2008003h25v06 NA  NA  NA  NA   NA

    X2008004h25v06 NA  NA  NA  NA   NA

    X2008005h25v06 NA  NA  NA  NA   NA

    X2008006h25v06 NA  NA  NA  NA   NA
      --------------   
     -----------------  so on

1 回复 | 直到 7 年前

Robert Hijmans 7 年前

下面是一个有效的示例:

library(raster)
fn <- system.file("external/test.grd", package="raster")
s <- stack(fn, fn)
quantile(s)

#            0%      25%     50%      75%    100%
#test.1 128.434 293.2325 371.412 499.8195 1805.78
#test.2 128.434 293.2325 371.412 499.8195 1805.78

可以通过以下方式获得相同的结果(使用较小的数据集):

t(apply(values(s), 2, quantile, na.rm=TRUE))

因此,您需要找出文件的不同之处。

我能想到的唯一返回分位数为NA的方法是当所有值都为NA时:

quantile(rep(NA, 10), na.rm=TRUE)
# 0%  25%  50%  75% 100% 
# NA   NA   NA   NA   NA

因此,文件中的所有值似乎都是 NA . 情况可能并非如此,但如果没有其中一个文件,则很难对其进行故障排除。

推荐文章

Amp · 使用R ggplot2删除geom_radial中axis.line和panel.border之间的空格

6 月前

Hard_Course · 用另一列中的值替换行的最后一个非NA条目

6 月前

Mark R · 使用geom_sf()删除地球仪上不需要的网格线

6 月前

Joe · 根据对工作日和本周早些时候的日期的了解,找到一个日期

6 月前

Ben · 统计向量中的单词在字符串中出现的频率

6 月前

TheCodeNovice · R中符号格式的尾随零和其他问题[重复]

6 月前

katefull06 · 在R中使用terra修改范围时,会为单独的SpatRaster重写范围

6 月前

dez93_2000 · 在R管道子功能中引用管道对象的当前状态

6 月前

accibio · 在ggplot2中为同一变量创建两个连续的颜色渐变比例

6 月前

Mankka · 如何在Ggplot2中绘制均匀的径向图

6 月前