之前做过两个空间数据的叠加:
微笑牛油果:R语言空间数据处理:intersectionzhuanlan.zhihu.com从那之后就想做一个系列,主要想研究sf的工具能不能替代我在ArcGIS里头常用的工具。本次先尝试做一些简单工作,看看可行性。
这里先声明一下,本文中出现的warnings并没有被解决,一些输出数据可能不精确。
- 读入shp文件并加入FID
代码:
grid<-read_sf('BJ_1km_132x138.shp', fid_column_name = 'FID')
grid_my<-read_sf('GRID_WGS84_ID.shp')
说明:
fid_column_name为字符串则会在数据框里头新增一列以该字符串命名的FID,shp文件的FID会从0开始编号。如果要把处理数据跟原始文件匹配的话,建议加上FID。
2. 实现ArcGIS中的join功能
代码:
raw<-fread(filepath,
colClasses = c('character',rep('numeric',2)),
col.names = c('grid','hour','value'