3

我正在使用 Tigris 包下载带有代码的形状文件,

options(tigris_class = "SF") 

虽然我可以轻松地映射多边形

ggplot(data=zctas) +
  geom_sf(aes())

我正在努力创建一个标记/点而不是一个填充的多边形(例如,一个形状 - 就像 zctas 内的一个圆圈),因为我从底格里斯岛拉下来的形状文件没有 x 的纬度/经度和 y AES 映射。形状文件有一个“几何”列,我想知道我是否可以使用它?

这里的文档,https://ggplot2.tidyverse.org/reference/ggsf.html

似乎表明geom_sf可以用来创建点?(我猜想使用 zctas 多边形的质心?)但我找不到一个例子?感谢任何资源来识别从这个 Tigris 生成的形状文件中映射点的代码,和/或提示和/或替代方法。

4

2 回答 2

10

您正在寻找stat_sf_coordinates(),在此处描述

library(ggplot2)

nc <- sf::st_read(system.file("shape/nc.shp", package="sf"))
#> Reading layer `nc' from data source `/Library/Frameworks/R.framework/Versions/3.6/Resources/library/sf/shape/nc.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 100 features and 14 fields
#> geometry type:  MULTIPOLYGON
#> dimension:      XY
#> bbox:           xmin: -84.32385 ymin: 33.88199 xmax: -75.45698 ymax: 36.58965
#> epsg (SRID):    4267
#> proj4string:    +proj=longlat +datum=NAD27 +no_defs

ggplot(nc) +
  geom_sf()

ggplot(nc) +
  geom_sf() +
  stat_sf_coordinates()
#> Warning in st_point_on_surface.sfc(sf::st_zm(x)): st_point_on_surface may
#> not give correct results for longitude/latitude data

ggplot(nc) +
  geom_sf() +
  geom_point(
    aes(color = SID74, size = AREA, geometry = geometry),
    stat = "sf_coordinates"
  ) +
  scale_color_viridis_c(option = "C") +
  theme(legend.position = "bottom")
#> Warning in st_point_on_surface.sfc(sf::st_zm(x)): st_point_on_surface may
#> not give correct results for longitude/latitude data

reprex 包(v0.3.0)于 2019-11-03 创建

关于最后一个示例调用中的geometry = geometry语句的一条评论:如果在数据集中找到几何列,两者都会自动映射几何列。但是,不知道几何列,因此不自动映射,因此我们必须手动映射。aes()geom_sf()stat_sf_coordinates()geom_point()

于 2019-11-03T17:48:26.397 回答
4

如果我正确理解您要实现的目标,您可能只需要使用st_centroid()函数将数据转换为点...最基本,如下所示:

ggplot(data=st_centroid(zctas)) +
  geom_sf()

也许:

ggplot(data=st_centroid(st_geometry(zctas))) +
  geom_sf()

更多细节和示例,请参见: https ://r-spatial.github.io/sf/articles/sf5.html

于 2019-11-03T16:55:50.687 回答