Implement spatial joins in R
发布时间
阅读量：
阅读量
今日完成了一项具体任务：
所涉及的代码逻辑并不繁琐，但实际操作过程中包含多个较为复杂的流程环节。
1）导入以CSV格式存储的points数据集，并对其进行重采样处理，从而生成grid数据集；
2）加载具备zone属性信息的多边形矢量数据；
3）将grid数据与多边形数据进行空间叠加分析，进而为每个grid单元赋予对应的多边形zone属性特征。
library(rgdal)
library(gstat) # Use gstat's idw routine
library(sp) # Used for the spsample function
library(raster) # Used to clip out thiessen polygons
imos=read.csv('D://research//GHM//imos.csv') ## 导入points的CSV文件
coordinates(imos) <- c("lon","lat")#定义坐标
vector <- readOGR('D:/research/3_GHM/Result/Zone/shp2.shp') ## 导入泰森三角形矢量文件
# Add P's projection 全部评论 (0)
还没有任何评论哟~