代码之家 › 专栏 › 技术社区 › geo_dd

在R中使用sf包映射区域

sf area grid dictionary r

geo_dd · 技术社区 · 6 年前

我花了几天时间想找到一种与 sf 用R打包,但没有成功。我想绘制我的区域,类似于中的示例2 this example 但没有任何意义。我将我的区域加载为 SpatialpolygonsDataFrame 然后我用防御来获得lat和long,如下所示:

area<-readOGR(dsn="/home/ubuntu/..",layer="area")
f_area<-fortify(area)
head(f_area) 

long      lat order  hole piece id group
1 116.1045 57.23717     1 FALSE     1  0   0.1
2 116.1551 57.21548     2 FALSE     1  0   0.1
3 116.2420 57.14505     3 FALSE     1  0   0.1
4 116.1706 57.12011     4 FALSE     1  0   0.1
5 116.1222 57.12006     5 FALSE     1  0   0.1
6 116.0756 57.09926     6 FALSE     1  0   0.1

从这一点上说,我很困惑我必须做什么才能得到我想要的结果。

谢谢你的帮助。

1 回复 | 直到 6 年前

sebdalgarno 6 年前

有了ggplot2::geom_sf,就不需要加固。

# install ggplot2 development version
install.packages('devtools')
devtools::install_github('tidyverse/ggplot2')

# load libraries
library(sf)
library(ggplot2)

# read in data (probably you need to change the path here to read your file)
area <- sf::st_read("/home/ubuntu/area.shp")

# ggplot
ggplot() + geom_sf(area)

# or you can simply use plot
plot(st_geometry(area))

推荐文章

Marc B. · 使用ggplot2创建条形图时“缺少值”

1 年前

deschen · tidyverse与外部向量发生突变,该外部向量的元素是数据帧中的列值

1 年前

Laura · 在Shiny中使用可排序的包拖放名称,这些名称将成为图表

1 年前

Mallikarjun M · 如何使用随机森林进行时间序列预测?

1 年前

ly li · 模型摘要:当表格形状改变时,拟合优度消失

1 年前

C.Robin · 将marginaffects::predictions()的结果连接回main df?

1 年前

monotonic · 如何将格式为“col1+col3+col4”的数据帧的行名转换为一列数字向量“c(1,3,4)”?

2 年前

Shawn Hemelstrand · 为什么我的自定义errorbar函数不能在R中工作?

2 年前

RoyBatty · 统计每个字符在整个数据集中出现的次数

2 年前

stats_noob · R: 记录某个“行为”发生的循环的索引?

2 年前