首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

一个R语言中操纵矢量空间数据的标准化工具—sf

目前还不知道,那些R语言中数百个使用了sp包提供的类方法的包,是否会、以及何时会将修改为依赖sf包的类方法。...该软件包的一个新功能是能够检索空间度量,并使用显式度量单元设置距离参数等 (Pebesma et al., 2016): > st_area(st_transform(nc[1, ], 2264)) #...这种能力主要由可见文本(WKT)二进制文本(WKB)序列化提供,它们是sf标准的一部分,并且受到sf的支持。与GDAL、GEOSliblwgeom库的通信都使用WKB方法。...一篇关于在sf中设置空间索引的博文 描述了如何使用索引操纵大内存的空间数据集。对于球面数据,还需要研究liblwgeom或s2提供的索引。...为了进一步了解sf的全部功能及其原理,读者可以参考软件包附带的六篇小品文。 致谢 如果没有Roger Bivand的前期工作和付出,sf包的面世几乎是不可能的。

4.1K51

R语言可视化——关于ggplot所支持的数据地图素材类型

library("plyr") library("dplyr") library("rgdal") library("sf") library(maptools) library("ggplot2")...而sf对象将这种控件数据格式件进行了更加整齐的布局,使用st_read()导入的空间数据对象完全是一个整齐的数据框,拥有整齐的行列,这些行列中包含着数据描述几何多边形的边界点信息。...然而在sf对象中我们仅需指定一次合并即可,即描述层业务指标数据的合并。...,而且建议使用rgdal::readOGRsf::st_read 来导入。...如果我们需要制作填色地图,那么我们需要先将自己的业务指标描述层数据进行整理并,并将合并后的数据与几何映射层的数据框进行合并。 这里我虚拟一个指标。

2.3K41
您找到你想要的搜索结果了吗?
是的
没有找到

rgdal包readOGR使用

1R语言地图数据分类 R语言使用ggplot2作图,所支持的地图数据对象主要包括两类 sp: SpatialPolygonDataFrame sf: Simple feature list column...SF数据特点 最大特点hi是,他将每一个行政区划所对应的几何边界点封装成一个list对象,这条记录就像其他普通的文本记录一样,被排列在对应行政区划描述单元中 使用sf包的st_read()函数导入的空间数据对象完全是一个整齐的数据结构...,这些行列中包括了描述层几何多边形的边界点信息。...SF对象我们只需要指定一次合并即可,即将描述层和你的分析数据合并,使用sf::st_read()函数读取数据即可得到SF数据对象,其为data.frame对象类型。...image.png 2使用readOGR读取shp绘制地图 library(rgdal) #提供readOGR()函数 library(ggplot2) library(dplyr) p<-choose.files

5.6K20

左手用R右手Python系列12——空间数据可视化与数据地图

今天要讲解的主角是R语言中的sfPython中的geopandas库。...也许以上描述过于抽象,因为涉及到到的内容比较深入,我实在是不知道该如何把这些内容将的通俗易懂,接下来会使用图片辅助演示。...我能告诉你的是,geom_ploygon制作地图的时候,剥离了地理信息边界点数据多边形属性信息,所以你需要同时兼顾、处理两个包含空间信息的数据框,如果是对不同区域进行等值线映射,你还需要对这两个数据框进行合并操作...ggplot() + geom_sf(data=china_map,aes(fill=value)) ? 所以使用sf提供的新技术,制作数据地图通常仅需以上几步。任务量大大缩减。...geopandas包同时支持导入shp素材json素材,导入之后得数据结构与R语言中得sf导入之后得结构是一致得,地理多边形边界点信息都被压缩成了一个非常整齐的列表存储,列表内每一个单独的子项目都代表着一个多边形

2K40

如何处理地图投影转换

但是使用geom_polygon()制作地图成本非常高,因为geom_polygon不直接支持GIS的数据模型(如sp、sf等)。...R语言中支持GIS数据模型的包一共有两个:sp包sf包,在旧版的ggplot2中,geom_polygon高度依赖从sp导入的数据对象(虽然也可以从sf中获取)。...但是这种情况马上会随着sf包的逐步完善以及ggplot2sf包的进一步融合而大有改观。 最新版的ggplot2(github上面的开发版)已经内置了geom_sf()图层。...之前使用simple模型练习的图表 ? 常见的多圆锥视角中国地图投影 ? 投影问题涉及到两个关键环节:地理坐标投影坐标的转换。...由于投影后的投影坐标系已经被投影算法转换,所以在使用geom_text等图层函数时,务必要使用与几何对象投影一致的经纬度点,这里使用sf中的点中心计算函数最为快捷。

1.7K30

贝叶斯地理统计模型R-INLA-1

贝叶斯地理统计模型INLA 本次博客主要讲述如何使用R-INLA软件进行空间分析,通过随机嵌套偏微分方程方法集成的嵌套Laplace渐进法可为潜在高斯随机场模型中的边际分布提供准确而有效的估计。...由于笔者水平有限,关于理论部分,可前往link,针对数学公式及理论部分,这里不赘述,简化数学公式,强调如何应用,及在R语言里面如何实现。...这里主要介绍如何判断空间位置的点是否存相互独立,更多关于Variogram,请参考此处Variograms & Kriging # set the train and test df_rain = st_as_sf...::filter(ID %in% id ) train=df_rain %>% dplyr::filter(!...(详情请见) inla.mesh.2d() 常用参数: loc 用作初始三角剖分节点的坐标矩阵 loc.domain 定义空间域的多边形的坐标 max.edge 内部(外部)区域的最大允许三角形边缘长度

1.6K20

详谈R语言构建地理投影系统绘制高端地图

❝本节来详细介绍如何使用R语言来构建地理投影系统绘制世界地图,细节挺多的小编做了详细的注释;结果仅供参考❞ 加载R包 library(tidyverse) library(sf) library(camcorder...投影系统用来确定如何在平面上把地球表面的地理空间信息进行投影,以便更好地展示分析数据。...这个投影系统叫做 "Winkel Tripel",它是一种广泛用于地图制作的投影方式,可以保证较为均衡地表示地球上的经纬度信息。...使用 st_graticule 函数来生成经纬线。该函数接受一个参数 lat,用来指定经纬线的纬度间隔。 在此将纬度间隔设为 -89.9,seq(-80, 80, 20) 89.9。...st_sfc函数将这个空间几何对象包装在一个简单空间几何集合sfc对象中,并为其指定坐标参考系统 st_sf函数将这个简单空间几何集合对象转换为一个空间数据框对象,并使用st_transform_proj

1.7K20

空间地理数据可视化之 leaflet 包及其拓展

1.基本画图设置 Leaflet 包是制作交互式地图非常流行的开源 JavaScript 库,可以很容易地在 R 中合成控制地图。...我们可以调用 leaflet() 函数来创建地图,并可以使用 addTiles() (添加背景地图)、 addPolygons() (添加多边形)、 addLegend() (添加图例) 等来添加图层。...在使用 leaflet包前,要求先将地图数据转化为 EPSG4326 下的投影,使用的是 sf 包中的 st_transform() 函数。...leaflet(map) %>% addTiles() %>% addPolygons( color = "black", fillColor = ~ pal(SID74), ##指定多边形边界填充的颜色...本篇是空间地理数据可视化系列的第四期,主要由 林华师 制作。本系列的宗旨是带你系统学习如何使用 R 对空间地理数据进行可视化。下一期将会继续介绍 mapview 包的使用,敬请期待。

2.5K10

空间地理数据可视化之 ggplot2 包及其拓展

基本画图设置 ggplot2[2]是一个基于图形语法来创建图形的包,因此我们可以使用 ggplot() 函数以下元素创建一个图: 想要可视化的数据; 指定数据的几何形状,如点或条。...其中,aes() 用于将数据中的变量映射为对象的视觉属性; 可选的元素,如标尺、标题、标签、图例主题等。 我们可以使用 geom_sf() 函数一个简单特征对象( sf 类)来创建地图。...更多设置 在 ggplot() 中,离散变量的默认色标是 scale_*_hue() ,这里 * 表示颜色(为点线等特征着色)或填充(为多边形或柱状图着色); scale_*_grey() 用来改变灰色颜色的默认比例...内容扩展 gganimate 包[3] plotly 包[4]可以与 ggplot2 包结合使用,分别创建动画交互式绘图。...本篇是空间地理数据可视化系列的第二期,主要由 林华师 制作。本系列的宗旨是带你系统学习如何使用 R 对空间地理数据进行可视化。下一期将会继续介绍其他可视化的 R 包,敬请期待。

3K30

一篇文章教你搞定JSON素材,从此告别SHP时代~

mymapdata<-merge(mapdata,city) 合并边界点数据各区名称与分组依据(主要是ggplot映射时作为分组变量使用) 因为各区的行政中心经纬度未知,这里暂时提取多边形中心作为其参考值...(long,lat))) 以上过程展示了如何从json格式的数据文件中提取我们制作数据地图所需要的指标(核心指标由三个:lon、lat、group),但是以上只够我们画出一幅单色地图,因为没有指定任何指标...id基于多边形编号的group才会一一对应,否则不会出现严格对应关系)。...,安徽省的各市级单位经纬度信息数据看起来在list不是同级的,即有些城市是单独一个list,有些城市是一个list里面嵌套好几个子list(这就解释了上面所讲过的,有些城市辖区不接壤,需要分别对其进行多边形描述定义...下一篇,跟大家细讲关于ggplot在制作数据地图过程中的变量映射规则注意事项。

1.7K60

结构方程模型 SEM 多元回归模型诊断分析学生测试成绩数据与可视化

p=24694 本文首先展示了如何将数据导入 R。然后,生成相关矩阵,然后进行两个预测变量回归分析。最后,展示了如何将矩阵输出为外部文件并将其用于回归。 数据输入清理 首先,我们将加载所需的包。...library(dplyr) #用于清理数据 library(Hmisc) #相关系数的显着性 然后,我们将使用 Fortran 读入数据文件并稍微清理数据文件。...NA 是默认值 # 使用 dplyr 对特定测试进行子集化 select(sub, c(T1, T2, T4)) # 使用 psych 包获取描述 请注意,R 将原始数据中的空白单元格视为缺失,...如果我们能找出一个异常的案例,我们在有没有这个案例的情况下进行分析,以确定其影响。输出的变化将是对杠杆的测试。 现在我们制作测试之间关系的 3d 散点图。...plot(T1,T2, T4, 3d(model) #使用我们先前的模型来绘制一个回归平面 使用相关矩阵的多元回归 现在我们将展示如何使用相关矩阵进行回归。

3K20

如何基于一张图片来创建3D模型?包含哪些步骤?

首先,基于一张图片(一个角度)来制作一个3D模型是不大现实的。因为三维物体是立体的、多维的,单从一个角度去观察很难判断物体其他视角的构造,制作出来的模型也就一个平面化的物体。...一般制作一个三维模型,用来参考的图片最少要包含正、侧、顶三个视图,这样就可以制作出对称性的三维模型了,比如人体建模。而对于不具有对称性的物体的建模,则需要更多角度的参考图了。...我们就以最常用的三维制作软件maya来简单说一下参照图片建模的大致步骤。 我们就以建人头模型为例。首先,我们需要用ps处理一下,保证正面侧面的照片脸部大小基本都对上。...分别将正面侧面导入到Maya的前视图窗口侧视图窗口 创建模型有多重方法,如画多边形、用box来修改等,基于对称性,我们用box来造型,把box左右分两等分,删掉左边的部分(或者右边的部分,根据个人喜好来决定...接着对照参考图,通过加点、线的方式对多边形进行编辑,最终完成头像的建模工作。

2K00

R语言有多强大?十个你不知道的功能

R语言的确提供了很全面的统计分析的软件包,比如CRAN,Bioconductor,Neuroconductor,以及ROpenSci;并且提供了优秀的包管理功能。...1.R的标记语言可以制作可重复生成的WordPowerpoint文档 R语言中的rmarkdown包可以制作可重复生成的Word文档Powerpoint幻灯片,而这只需要改变一行YAML的代码。...4.通过使用R语言的dplyr/dbplyr,几乎各种数据库都可以连接 使用dbplyr包,用R语言连接各种数据库,无论是本地的还是远程的,都非常方便。...5.本地或多个不同的数据存储,在R语言里可以利用相同的dblyr语法来操作 当你学会如何利用dplyr来转换数据,本地远程的数据库、数据存储都可以利用相同的代码来操作。...6.你可以用KerasTensorflow训练深度学习模型 使用keras包或TensorFlow接口,你可以利用R语言来学习预训练或者开发全新的深度学习模型。

1K30

一篇小短文助你打开数据可视化的任督二脉!

本文主要讨论ggplot2是如何通过颜色信号来对多边形进行填充的底层理念,这也是想要进阶R语言数据可视化过程中必须搞明白的关键环节。...边界点不声明顺序,软件如何知道该按照什么样的顺序来进行打印呢?...当这种grouporder顺序定义之后,软件首先将所有的经纬度坐标点按照group顺序打印,即先打印group顺序排在第一的多边形,group内部按照order的顺序,依次打印左边点,单个group但因完毕之后...(至于这两种颜色映射方法的具体规则原理,以后会讲到)。...(这里说明以下,所用地图仅用于案例演示使用,限于经济能力,无法获得高精度地图,如果出现边界线纠纷问题,不代表官方观点)。 当然,你可以将颜色映射在行政区划的轮廓线上。

1.3K40
领券