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

Question

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

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

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

示例见我之前的问题：

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)

如果其中的 none 允许您复制错误，您最后的选择是为该文件提供 link 以便其他人可以下载它。如果你不能这样做，那么至少 show SpatRaster x 的内容 show(x) 并提供代码来创建一个类似的对象 as.character(x)

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

How can I create a reproducible example of a SpatRaster (terra)?

gis

r

raster

terra