万里写入胸怀间-CSDN博客

原创利用丁香园数据生成疫情分布地图（R语言）

疫情牵动大家，不过就个人而言，目前也没法做得更多。欢迎使用Markdown编辑器你好！这是你第一次使用 Markdown编辑器所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。新的改变我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用...

2020-02-13 12:49:25 6805 1

虎扑上有位网友发起了一个很有意思的话题，《湖北是唯一一个到各省区最多只要跨过两个省的省份吗？》这在地理里是个多阶空间邻接矩阵的问题。我觉得这个问题可以作为学习R语言中spdep包的入门吧。顺手解决了一下。中国地图的矢量数据是我从高德上抓取的。先画一个中国地图library(pacman)p_load(tidyverse,sf,spdep,tmap)tm_shape(china_provinc...

2019-03-16 07:56:25 1039

原创 R语言中的空间最邻近点问题

在空间分析中，Nearest Point Problem是一个十分基本而重要的问题。对于小样本量的点集而言，通过计算距离矩阵并进行排序即可暴力解决。而对于上万甚至几十万以上的点集而言，通过计算距离矩阵显然需要耗费太多的资源，一般采用KD树等方法进行搜寻。背后的算法我也不太懂，这里仅把R语言中面向st_point类型的点集时运用两种方法的代码记录一下。point_set为st_point类型的点集。...

2019-03-04 22:24:05 3568

原创 R语言实现最小凸包算法

#R语言实现最小凸包算法p=co_db[which.min(co_db$Y),] #initial p as the start point#1 order point set by its angular (from min to max)co_db$cos.theta=(co_db$X-p$X)/sqrt((co_db$X-p$X)^2+(co_db$Y-p$Y)^2)co_db_o...

2018-10-31 17:12:45 1288

原创 R语言中的空间插值

最邻近插值实际上是利用点生产沃诺尼多边形，每个点所在的沃诺尼多边形的值就等于点值，这样就实现了由点到面的变化，完成了插值。library(dismo)v &amp;amp;amp;lt;- voronoi(dta)plot(v)ca &amp;amp;amp;lt;- aggregate(cata)vca &amp;amp;amp;lt;- intersect(v, ca)spplot(vca, 'prec', col.regions=r

2018-08-14 20:27:07 6488

原创利用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 10036 7

原创 KML格式转换为SF格式

谷歌地球的数据导入R语言可以直接用st_read命令生成sf类型的对象。通常情况下使用是没有问题的。但是实际上，与在R语言内部直接通过坐标点生成的sf类型对象或者读入shp等矢量数据得到的sf对象不同，kml格式的数据空间维度是XYZ 而并非XY ，也即多了一个Z维度。在与XY维度的数据混在一个数据框中构建一个sf对象时，会发生一些莫名其妙的错误。所以，还是需要对数据进行清洗和重组，以下是我的解决...

2018-08-12 10:16:51 489

原创 R语言空间数据处理笔记—矢量-栅格数据交互

矢量-栅格数据交互根据矢量数据范围裁剪栅格数据统一投影vector = st_transform(vector,projection(raster))设定矢量数据范围外栅格数据的值为NAraster_masked= mask(raster,as(vector,'Spatial')本Markdown编辑器使用StackEdit修改而来，用它写博客，将会带来全新的体...

2018-06-27 22:12:22 15717 3

原创 R导出shp文件不乱码的方法

今天下午试了很多包，sf,sp,rgdal 等等纷纷沦陷。还好有maptools包。函数是writeSpatialShapefile.sf %>% as('Spatial') %>% writeSpatialShape('file.shp')这个函数只认sp包里的那些SpatialDataFrame格式的文件，所以要先把sf格式转换一下。这都不是问题，反正从shp文件可以...

2018-06-04 17:54:36 3179

原创使用arcgisbinding包导入Raster文件

使用arcgisbinding包导入Raster文件由于ARCGIS在国内的领先地位，在用R语言处理空间数据时和ESRI公司独有的以数据集的形式来存储的矢量和栅格数据是不可避免的问题。而R语言现有的raster包和rgdal 包好像都无法直接导入raster文件。好在ESRI公司发布了R-ArcGIS桥接工具并开发了arcgisbinding 包。1 安装安装部分直接参考虾神的这篇...

2018-05-31 17:23:33 839 1

万里写入胸怀间