R 包 `terra` 如何确定在没有模板的情况下投影的栅格的维数？

Question

我想知道

terra

如何确定在没有模板栅格的情况下投影的栅格中的行数和列数。例如：

library(terra)

from <- rast()
mollweide <- 'PROJCRS[\"Sphere_Mollweide\",\n    BASEGEOGCRS[\"Unknown datum based upon the Authalic Sphere\",\n        DATUM[\"Not specified (based on Authalic Sphere)\",\n            ELLIPSOID[\"Sphere\",6371000,0,\n                LENGTHUNIT[\"metre\",1]]],\n        PRIMEM[\"Greenwich\",0,\n            ANGLEUNIT[\"Degree\",0.0174532925199433]]],\n    CONVERSION[\"Sphere_Mollweide\",\n        METHOD[\"Mollweide\"],\n        PARAMETER[\"Longitude of natural origin\",0,\n            ANGLEUNIT[\"Degree\",0.0174532925199433],\n            ID[\"EPSG\",8802]],\n        PARAMETER[\"False easting\",0,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8806]],\n        PARAMETER[\"False northing\",0,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8807]]],\n    CS[Cartesian,2],\n        AXIS[\"(E)\",east,\n            ORDER[1],\n            LENGTHUNIT[\"metre\",1]],\n        AXIS[\"(N)\",north,\n            ORDER[2],\n            LENGTHUNIT[\"metre\",1]],\n    USAGE[\n        SCOPE[\"Not known.\"],\n        AREA[\"World.\"],\n        BBOX[-90,-180,90,180]],\n    ID[\"ESRI\",53009]]'

to <- project(from, mollweide)

dim(from)
dim(to)

从

project()

的帮助推断，似乎

terra

正在使用

gdalwarp

，除非用户禁用它，但gdalwarp的帮助

page

没有阐明。

Answer 1

选择输出栅格几何形状，使其范围和分辨率与输入栅格紧密匹配。