删除空栅格 - 星星包 R

【问题标题】：Remove emtpy raster - stars package R删除空栅格 - 星星包 R
【发布时间】：2021-01-30 20:51:38
【问题描述】：

我已经使用stars 包加载了一个栅格，然后在栅格上创建了切片。现在我根据这些图块对栅格进行子集化。

tiles[[i]] <- st_bbox(c(xmin=x0,ymin=y0,xmax=x1,ymax=y1),crs=st_crs(r))
crop_tiles[[i]] <- r[tiles[[i]]]

这里r 是使用read_stars("filename.tif") 加载的光栅。

现在我想在 crop_tiles 中删除空的。即栅格值为NA的图块

【问题讨论】：

标签： r-stars

【解决方案1】：

您可以检查名为r 的stars 对象的所有值是否都是NA：

all(is.na(r[[1]]))

然后删除上面等于TRUE的图块。

对于更具体的代码，请在您的问题中提供可重现的示例数据，谢谢。

【讨论】：

反之亦然吗？我的数据中缺少时间步长，需要填写空白图层。
@JoanneDemmler 有趣的问题！首先想到的是遍历矢量日期，对于每个日期i 检查它是否存在于您的栅格中：如果存在则保留它，如果不存在则创建一个空栅格。最后，合并所有生成的波段。如果您可以在单独的问题中发布示例，我将尝试准备示例代码。