如何在R中插入网格空间对象

install.packages("sp")
install.packages("raster")
library(sp)
library(raster)

# 创建一个空的网格空间对象
grid <- GridTopology(cellcentre.offset = c(0, 0), cellsize = c(1, 1), cells.dim = c(10, 10))

# 将网格空间对象转换为SpatialGrid对象
grid_sp <- SpatialGrid(grid)

# 创建一个空的SpatialPointsDataFrame对象
points <- SpatialPointsDataFrame(coords = data.frame(x = numeric(), y = numeric()), data = data.frame(id = character()))

# 生成一些随机点坐标
set.seed(123)
random_points <- data.frame(x = runif(10, 0, 10), y = runif(10, 0, 10))

# 将随机点坐标插入到网格空间对象中
for (i in 1:nrow(random_points)) {
  point <- SpatialPoints(coords = random_points[i,])
  points <- rbind(points, point)
}

# 将插入的点对象与网格空间对象合并
grid_with_points <- SpatialGridDataFrame(grid_sp, data = data.frame(id = character()), point = points)

# 可视化网格空间对象
plot(grid_with_points)