R语言【raster】——projectRaster()：映射一个Raster对象

Package raster version 3.6-27

Description

用另一个投影(坐标参考系统，(CRS))将一个Raster*对象的值投影到一个新的Raster*对象。

您可以通过将新投影作为单个参数提供来实现这一点，在这种情况下，函数将设置新对象的范围和分辨率。为了对转换有更多的控制，例如，为了确保新对象与其他数据集对齐，您可以提供一个Raster*对象，该对象具有输入数据应该投射到的属性。

projectExtent返回一个带投影范围的RasterLayer，但是没有任何值。然后可以调整这个光栅层(例如，通过设置其分辨率)，并在projectaster中用作模板。

Usage

R 复制代码

projectRaster(from, to, res, crs, method="bilinear", 
             alignOnly=FALSE, over=FALSE, filename="", ...) 

projectExtent(object, crs)

Arguments

参数【from】：Raster*对象。

参数【to】：Raster*对象。

参数【res】：单个或两个数字的向量。如果缺少' To '，可选地设置输出分辨率。

参数【crs】：类'CRS'的字符或对象。project .4坐标参考系的描述。在projectRaster中，这是用来设置输出的CRS，如果'to'缺失，或者如果'to'没有有效的CRS。

参数【method】：方法用于计算新的RasterLayer的值。要么'ngb'(最近邻)，这对分类变量很有用，要么'bilinear'(双线性插值;默认值)，这适用于连续变量。

参数【alignOnly】：逻辑值。使用to或其他参数仅用于对齐输出(即相同的原点和分辨率)，但使用来自from的投影范围。

参数【over】：逻辑值。如果为TRUE，则关闭日期行环绕。这对于全局数据是可取的(以避免两次映射相同的区域)，但在其他情况下不可取

参数【filename】：字符值。输出文件名。

参数【...】：应用于 writeRaster 的其他参数。

参数【object】：Raster* 对象。

Details

有两种方法可以用来投影Raster对象的值。

1）提供crs参数，以及可选的res参数，但不提供to参数。

2）用你想要投影到的CRS创建一个模板栅格。你可以使用一个已经存在的对象，或者使用projectExtent或者一个已经存在的Raster*对象。还可以设置行数和列数(或分辨率)，可能还可以调整范围。输出光栅的分辨率通常应该与输入光栅的分辨率相似。然后使用该对象作为参数，将输入光栅投影到。这是首选的方法，因为你有最大的控制权。例如，您可以确保生成的光栅对象与其他光栅对象对齐。

投影是使用PROJ库执行的。

也请参见projInfo('proj')， projInfo('ellps')和projInfo('datum')获取有效的project.4值。

Value

RasterLayer 或 RasterBrick 对象。

Note

如果输出的分辨率比输入的分辨率大得多，您应该首先聚合输入，使输入的分辨率变得与输出更相似(可能更小)。

用户注意。遗憾的是，这个函数使用的GDAL/PROJ库已经部分弃用了PROJ.4表示法。您仍然可以使用它，但只能使用WGS84基准。其他数据将被静默忽略。

当打印一个Spat *对象时，会显示PROJ.4符号，因为它是最简洁和清晰的格式。但是，在内部使用WKT表示(参见crs)。

向量(点，线，多边形)可以用spTransform变换。

当将投影的环极数据转换为(例如)经度/纬度时，projectExtent不能很好地工作。有了这些数据，您可能需要调整返回的对象。例如:do ymax(object) <- 90。