光栅投影到 UTM 到 R 中的 Lat-Lon

Question

我是 R 的初学者。我有一个 rasterlayer，我想更改它的坐标参考系统（从 UTM 到 lan-lot）。我使用了以下代码：

r <-projectRaster(from = tiff,
   crs=crs("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))

（这里的 tiff 是要重新投影的 rasterlayer）

但我遇到了一个错误

“错误：无法分配大小为 439.0 Mb 的向量”

有人，请提供解决方案

Answer 1

这很奇怪。你把tiff中的所有值都读入内存了吗？

问题可能会消失

• 在 projectraster 中使用文件名参数。例如。 projectRaster(tiff, filename="out.tif", crs= ...)

• 将一些 rasterOptions 设置为较低的值，例如，rasterOptions(chunksize=1e+06, maxmemory=1e+08)

此外，将projectRaster 与crs= 参数一起使用通常不是一个好主意。最好提供一个 Raster* 对象作为要投影到的模板。这样，您可以设置范围和分辨率，以便输出栅格与您正在使用的其他栅格数据对齐。