如何找到R中给定区域内的voronoi细胞的面积？

# 安装和加载必要的包
install.packages("sp")
install.packages("rgeos")
install.packages("deldir")
library(sp)
library(rgeos)
library(deldir)

# 创建点数据
points_df <- data.frame(x = c(1, 2, 5, 6), y = c(2, 5, 3, 7))
coordinates(points_df) <- ~ x + y

# 计算Voronoi图
voronoi <- deldir(points_df$x, points_df$y)

# 将Voronoi图转换为空间多边形对象
voronoi_polygons <- as(voronoi, "SpatialPolygons")

# 定义给定区域
region_polygon <- Polygon(cbind(x = c(0, 2, 2, 0), y = c(0, 0, 7, 7)))
region_polygons <- Polygons(list(region_polygon), ID = "region")
region_sp <- SpatialPolygons(list(region_polygons))

# 裁剪Voronoi图到给定区域
voronoi_in_region <- gIntersection(voronoi_polygons, region_sp, byid = TRUE)

# 计算每个Voronoi细胞的面积
voronoi_areas <- sapply(slot(voronoi_in_region, "polygons"), function(poly) {
  area(poly@coords)
})

# 输出结果
print(voronoi_areas)