从栅格图层中的数据框中提取的绘图点答案

【问题标题】：Plotting points extracted from a dataframe in a Raster layer从栅格图层中的数据框中提取的绘图点
【发布时间】：2022-01-12 02:55:41
【问题描述】：

我想将数据框中的点绘制到我拥有的栅格图层中。对于每个点，我希望像元的值为 1（初始栅格图层上的所有其他像元的值为零）。

我的数据框（数据）看起来像这样（仅前三行）

Year<-c(2020, 2019, 2018)
Lat<-c(48.3,48.79,48.4)
Long<-c(-123.62, -123.36, -123.29)

我设法用以下代码绘制了这些点

points = st_as_sf(data, coords = c("Long", "Lat"), crs = 4326)
plot(st_geometry(points), pch=16, col="navy")

得到了这张图： Graph points plotted

我现在想将这些点绘制到该区域的栅格图层中。我的栅格图层的参数如下：

class      : RasterLayer 
dimensions : 44, 41, 1804  (nrow, ncol, ncell)
resolution : 0.2916667, 0.2916667  (x, y)
extent     : -133.2625, -121.3042, 41.3875, 54.22083  (xmin, xmax, ymin, ymax)
crs        : NA 
names      : Blank_Map

这个栅格的每个像元的值为 0，这正是我所需要的。现在，我想将我的数据帧中的点添加到这个栅格层，对这些数据点所在的每个单元格使用值 1。我还想将整个内容保存为新的栅格图层（然后将具有 0 和 1 值）。

有人可以帮我实现吗？我已经尝试了好几天，但似乎没有任何效果任何帮助表示赞赏！

【问题讨论】：

嗨@FanCou，我在下面的回答中建议了两种解决方案（一种使用raster 库，另一种使用terra 库）。如果他们满足您的需求，请考虑将此答案标记为“已接受”和/或“赞成”，以便其他 SO 用户更容易找到正确答案。如果没有，请告诉我出了什么问题。干杯

标签： r dataframe raster

【解决方案1】：

1.请在下面找到使用 raster 和 sf 库提供可能解决方案的代表。

library(raster)
library(sf)


# Create from scratch the raster with 0 values
raster0 <- raster(nrows = 44, ncols = 41, 
                  xmn = -133.2625, xmx = -121.3042, 
                  ymn = 41.3875, ymx = 54.22083, 
                  vals=0)


# Convert points 'sf' object into 'sp' object
points_sp <- as(points, "Spatial")

# Extract the cells of raster0 matching the points geometry of the 'sp' object
cells_ID <- extract(raster0, points_sp, cellnumbers=TRUE)[,"cells"]

# Copy raster0 (as you want the final result in another raster)
raster01 <- raster0

# Replace 0 values to 1 for cells matching points geometry in the 'raster01'
raster01[cells_ID] <- 1

# Visualization of the final result
plot(raster01)

^{由reprex package 创建于 2021-12-07 (v2.0.1)}

2。请在下面找到使用 terra 和 sf 库提供另一种解决方案的代表

library(terra)
library(sf)

# Create from scratch the 'SpatRaster' with 0 values
raster0 <- rast(nrows=44, ncols=41, 
                nlyrs=1, 
                xmin=-133.2625, xmax=-121.3042, 
                ymin=41.3875, ymax=54.22083, 
                vals = 0)


# Convert points 'sf' object into 'SpatVector' object
points <- vect(points)

# Extract the cells of raster0 matching the points geometry of the 'sp' object
cells_ID <- terra::extract(raster0, points, cells = TRUE)

# Copy raster0 (as you want the final result in another raster)
raster01 <- raster0

# Replace 0 values to 1 for cells matching points geometry in the 'raster01'
raster01[cells_ID[,"cell"]] <- 1

# Visualization of the final result
plot(raster01)

^{由reprex package (v2.0.1) 于 2021-12-07 创建}

【讨论】：