空间数据处理
万里写入胸怀间
这个作者很懒,什么都没留下…
展开
-
R语言空间数据处理笔记—矢量-栅格数据交互
矢量-栅格数据交互根据矢量数据范围裁剪栅格数据 统一投影vector = st_transform(vector,projection(raster))设定矢量数据范围外栅格数据的值为NAraster_masked= mask(raster,as(vector,'Spatial')本Markdown编辑器使用StackEdit修改而来,用它写博客,将会带来全新的体...原创 2018-06-27 22:12:22 · 15970 阅读 · 3 评论 -
利用R语言从高德地图生成AOI(Area Of Interests)
p_load(“httr”,’tidyverse’,’jsonlite’,’sf’,’rgdal’,’rlist’,’rgdal’) library(‘Rgctc2’,lib.loc=’~/GitHub/R_coordination_transformation’)search for poi’s ID through amap WEB_apioptions(digits=11) ...原创 2018-08-13 12:02:32 · 10333 阅读 · 7 评论 -
KML格式转换为SF格式
谷歌地球的数据导入R语言可以直接用st_read命令生成sf类型的对象。通常情况下使用是没有问题的。但是实际上,与在R语言内部直接通过坐标点生成的sf类型对象或者读入shp等矢量数据得到的sf对象不同,kml格式的数据空间维度是XYZ 而并非XY ,也即多了一个Z维度。在与XY维度的数据混在一个数据框中构建一个sf对象时,会发生一些莫名其妙的错误。所以,还是需要对数据进行清洗和重组,以下是我的解决...原创 2018-08-12 10:16:51 · 540 阅读 · 0 评论 -
R语言中的空间插值
最邻近插值实际上是利用点生产沃诺尼多边形,每个点所在的沃诺尼多边形的值就等于点值,这样就实现了由点到面的变化,完成了插值。library(dismo)v <- voronoi(dta)plot(v)ca <- aggregate(cata)vca <- intersect(v, ca)spplot(vca, 'prec', col.regions=r原创 2018-08-14 20:27:07 · 6562 阅读 · 0 评论 -
R语言中的空间最邻近点问题
在空间分析中,Nearest Point Problem是一个十分基本而重要的问题。对于小样本量的点集而言,通过计算距离矩阵并进行排序即可暴力解决。而对于上万甚至几十万以上的点集而言,通过计算距离矩阵显然需要耗费太多的资源,一般采用KD树等方法进行搜寻。背后的算法我也不太懂,这里仅把R语言中面向st_point类型的点集时运用两种方法的代码记录一下。point_set为st_point类型的点集。...原创 2019-03-04 22:24:05 · 3673 阅读 · 0 评论 -
哪个省能称为“九省通衢”?
虎扑上有位网友发起了一个很有意思的话题,《湖北是唯一一个到各省区最多只要跨过两个省的省份吗?》这在地理里是个多阶空间邻接矩阵的问题。我觉得这个问题可以作为学习R语言中spdep包的入门吧。顺手解决了一下。中国地图的矢量数据是我从高德上抓取的。先画一个中国地图library(pacman)p_load(tidyverse,sf,spdep,tmap)tm_shape(china_provinc...原创 2019-03-16 07:56:25 · 1111 阅读 · 0 评论