Question

对于特定于我的特定数据集的问题，如果我将它作为 SpatRaster 存储在 R 中，我如何制作该数据集的可重现示例？

数据结构足够复杂，我不知道如何徒手发明一个更简单的版本并将其作为 SpatRast 读取（即 x <- rast(????????)

我也无法弄清楚如何使用包或命令来提取足够的信息来提供功能上可重现的示例

以我之前的问题为例：How can I add a class name to numeric raster values in a terra SpatRaster?

Answer 1

您可以像这样从头开始创建对象：

library(terra)
r <- rast()
s <- rast(ncols=22, nrows=25, nlyrs=5, xmin=0)

有关您可以使用的其他参数和替代方法，请参阅 ?terra::rast。

您也可以使用 R 附带的文件。例如：

f <- system.file("ex/elev.tif", package="terra")
r <- rast(f)

您还可以从头开始创建一个新的 SpatRaster，其（大部分）与返回的属性相同

as.character(r)

然后用类似的东西重新创建它

r <- rast(ncols=95, nrows=90, nlyrs=1, xmin=5.74166666666667, xmax=6.53333333333333, ymin=49.4416666666667, ymax=50.1916666666667, names=c('elevation'), crs='GEOGCRS[\"WGS 84\",DATUM[\"World Geodetic System 1984\",ELLIPSOID[\"WGS 84\",6378137,298.257223563,LENGTHUNIT[\"metre\",1]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433]],CS[ellipsoidal,2],AXIS[\"geodetic latitude (Lat)\",north,ORDER[1],ANGLEUNIT[\"degree\",0.0174532925199433]],AXIS[\"geodetic longitude (Lon)\",east,ORDER[2],ANGLEUNIT[\"degree\",0.0174532925199433]],ID[\"EPSG\",4326]]')
r <- init(r, "cell")

如果您无法使用示例数据复制错误，这可能会给您提供有关问题的提示。它与 NA 有关系吗？文件在磁盘上？文件格式？一种棘手的情况是，如果实际文件大得多，就会存在差异。您可以通过设置 terraOptions(todisk=TRUE) 并在函数中使用 steps 参数来模拟大文件，例如

b <- clamp(x, steps=5)

如果这些都不允许您复制错误，那么最后的方法是提供文件的链接，以便其他人可以下载它。如果您不能这样做，那么至少 show 带有 x 的 SpatRaster show(x) 的内容并提供代码以使用 as.character(x)

创建类似的对象

如何创建可重现的 SpatRaster (terra) 示例？

1 个答案: