R语言【raster】——rasterize()：栅格化点、线、面

news/发布时间2024/5/15 7:10:39

Package raster version 3.6-27

Description

将与“对象”类型空间数据(点、线、多边形)相关的值转移到栅格单元。

对于多边形，如果多边形覆盖栅格单元的中心，则传输值。对于行，将值传输到与行接触的所有单元格。您可以通过先将多边形光栅化为直线，然后再光栅化为多边形来组合此行为。

如果x代表点，每个点被分配到一个网格单元。落在单元格之间边界上的点被放置在右边的单元格和/或下面的单元格中。网格单元格的值由与点和函数fun相关的值决定。

Usage

## S4 method for signature 'matrix,Raster'
rasterize(x, y, field, fun='last', background=NA,mask=FALSE, update=FALSE, updateValue='all', filename="", na.rm=TRUE, ...)## S4 method for signature 'SpatialPoints,Raster'
rasterize(x, y, field, fun='last', background=NA,mask=FALSE, update=FALSE, updateValue='all', filename="", na.rm=TRUE, ...)## S4 method for signature 'SpatialLines,Raster'
rasterize(x, y, field, fun='last', background=NA,mask=FALSE, update=FALSE, updateValue='all', filename="", ...)## S4 method for signature 'SpatialPolygons,Raster'
rasterize(x, y, field, fun='last', background=NA,mask=FALSE, update=FALSE, updateValue='all', filename="",getCover=FALSE, silent=TRUE, ...)

Arguments

参数【x】：点(一个SpatialPoints*对象，或一个两列矩阵(或data.frame))， SpatialLines*， SpatialPolygons*，或一个Extent对象。

参数【y】：Raster*对象。

参数【field】：数字或字符。要传输的值。这可以是单个数字，也可以是与空间特征(点、线、多边形)数量长度相同的数字向量。如果x是一个Spatial*DataFrame，这可以是要传输的变量的列名。如果缺少，则使用属性索引(即从1到特征数的数字)。您还可以提供与空间特征数量相同长度的向量，或者提供行数与空间特征数量匹配的矩阵。

参数【fun】：功能或特征。确定要为被多个空间特征覆盖的单元格分配哪些值。您可以使用min, max或mean等函数，或以下字符值之一:'first'， 'last'， 'count'。默认值是'last'。在SpatialLines*的情况下，'length'也是允许的(目前仅用于平面坐标系统)。

如果x代表点数，fun必须接受na.rm参数，无论是显式地还是通过省略号('dots')。这意味着fun=length失败了，但是fun=function(x，…)length(x)有效，尽管它忽略了na.rm参数。你可以使用这样的函数:fun=function(x, na.rm){if (na.rm) length(na.omit(x)) else (length(x)}，或者使用一个在所有情况下都删除NA值的函数，就像这个函数来计算每个网格单元“丰富度”的唯一值的数量:fun=function(x，…){length(unique(na.omit(x)))}。如果你想计算每个网格单元的点数，你可以使用fun='count'或fun=function(x，…){length(x)}。你也可以使用fun=function(x，…)c(length(x)，mean(x))这样的语句传递多个函数，在这种情况下返回的对象是一个RasterBrick(多层)。

参数【background】：数字。值放置在未被x的任何特征覆盖的单元格中。默认值为NA。

参数【mask】：逻辑值。如果为TRUE，则输入栅格对象的值被x的空间特征“掩盖”。也就是说，与空间特征在空间上重叠的单元格保留其值，其他单元格成为NA。默认为FALSE。当update=TRUE时，此选项不能使用。

参数【update】：逻辑值。如果为TRUE，则更新与x的空间特征重叠的单元格的Raster*对象的值。默认为FALSE。当mask=TRUE时不能使用。

参数【updateValue】：数字(通常是整数)或字符。仅在update=TRUE时相关。根据单元格的值选择要用空间特征的值进行更新的单元格。有效的字符值是'all'， 'NA'和'!NA'。默认为“all”。

参数【filename】：字符值。输出文件名（可选）。

参数【na.rm】：如果为TRUE，则如果fun使用na.rm参数，则删除NA值。

参数【getCover】：逻辑值。如果为TRUE，则返回被多边形覆盖的每个网格单元的百分比(并且忽略field, fun, mask和update的值)。通过将每个细胞划分为100个子细胞并确定每个子细胞中心多边形的存在/不存在来估计覆盖的分数。

参数【silent】：逻辑值。如果为TRUE，对多边形数量的反馈将被抑制。默认为FALSE。

参数【...】：文件写入的附加参数，如writeRaster。