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

Question

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

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

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

有关示例，请参见我之前的问题：如何将类名添加到 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")

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

b <- clamp(x, steps=5) 

如果这些都不允许您复制错误，那么您最后的手段是提供该文件的链接，以便其他人可以下载它。如果你不能这样做，那么至少showSpatRaster 的内容x和show(x)提供代码来创建一个类似的对象as.character(x)