r - 如何创建 SpatRaster (terra) 的可重现示例?
问题描述
对于特定于我的特定数据集的问题,如果我将其存储为 R 中的 SpatRaster,如何制作该数据集的可重现示例?
数据结构足够复杂,以至于我不知道如何徒手发明一个更简单的版本并将其读取为 SpatRast(即 x <- rast(????????)
我也无法弄清楚如何使用包或命令来提取足够的信息来提供功能上可重现的示例
有关示例,请参见我之前的问题:如何将类名添加到 terra SpatRaster 中的数字栅格值?
解决方案
您可以像这样从头开始创建对象:
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
s有关系吗?文件在磁盘上?文件格式?一种棘手的情况是,如果实际文件大得多,则会出现差异。terraOptions(todisk=TRUE)
您可以通过在函数中设置和使用参数来模拟大文件steps
,例如
b <- clamp(x, steps=5)
如果这些都不允许您复制错误,那么您最后的手段是提供该文件的链接,以便其他人可以下载它。如果你不能这样做,那么至少show
SpatRaster 的内容x
和show(x)
提供代码来创建一个类似的对象as.character(x)
推荐阅读
- javascript - 我无法在 excel web add in (javascript) 中保存对代码的更改
- r - 如何计算 Quanteda 中多词表达的频率?
- linux - 源子中的父脚本名称
- wordpress - 在 Woocommerce 产品属性中保存 URL
- flutter - 将 DraggableScrollableSheet 添加到 Sliver 页面的底部
- excel - Excel VBA多次重复3行相同的值
- c++ - 在 C++ 中检查枚举变量的无效值
- lumen - 找不到路由时的 Lumen API 自定义错误消息
- python - django python - 如何列出列表的渐进数值
- c# - MediaPlayer 和 RenderTargetBitmap 从视频中获取帧图像的问题