我正在使用 spatstat 包在指定区域上创建六边形镶嵌,并计算每个六边形内的点数。之后,我想将细分保存为 shapefile 以在 GIS 中使用。
如何将 tess 对象转换为 sf 或其他兼容格式,以便将其另存为 shapefile?任何指导将不胜感激!
下面是特定数据集的代码,不幸的是,该数据集未公开,但已显示该图,并且它是对象
QC_empresastur_hex
应该转换为 sf
格式,以便随后导出为 shapefile。
library(spatstat)
hexagonos <- hextess (win_natales, 1000)
QC_ empresastur_hex <- quadratcount(empresastur_ppp, tess = hexagonos)
plot (empresastur_ppp, main = "Conteo de empresas por hexágonos", cols = "orange")
plot(QC_empresastur_hex, add = TRUE, cex = 1)
这是使用内置数据集的示例:
library(sf)
library(spatstat)
hex <- hextess(Window(chorley), s = 4)
counts <- quadratcount(chorley, tess = hex)
plot(counts)
hex_sf_list <- lapply(tiles(hex), st_as_sf)
hex_sf <- Reduce(rbind, hex_sf_list)
counts_sf <- cbind(counts = as.integer(counts), hex_sf)
plot(counts_sf)