6

我目前正在尝试将北卡罗来纳州的所有邮政编码映射到传单地图上。最终，我想用另一个数据集的值对它们进行颜色编码。

虽然我已经设法使用这些值将标记放在地图上，但以彩色方式使用邮政编码这样做确实会为我的项目增加价值。

我在 R 中下载了 USAboundaries 库，所以我认为我有足够的数据在传单中创建多边形，但数据不是 SpatialPointsDataFrame 格式。

USAboundaries 数据集的示例输出如下所示：

USAboundaries::us_zipcodes()

Simple feature collection with 33144 features and 6 fields
geometry type:  POINT
dimension: XY
bbox:  
 xmin: -176.6316 
 ymin: -14.22075 
 xmax: 145.7536 
 ymax: 71.2738

epsg (SRID):    4326

proj4string:    +proj=longlat +datum=WGS84 +no_defs

First 10 features:
  zipcode zcta5ce10     affgeoid10 geoid10   aland10 awater10                   geometry
    21914     21914 8600000US21914   21914   1924479   477175 POINT (-75.98187 39.57303)
    01001     01001 8600000US01001   01001  29731610  2118827 POINT (-72.62585 42.06258)
    34736     34736 8600000US34736   34736 322808220 78588518 POINT (-81.89624 28.55458)
    46151     46151 8600000US46151   46151 530632048  9804480  POINT (-86.4305 39.44806)
    48039     48039 8600000US48039   48039  59592687  4845242 POINT (-82.55054 42.69666)

在我看来，这应该足以转换为 SpatialPointsDataFrame，但我不知道如何在 R 中做到这一点！任何提示将非常感谢！

4

1 回答 1

4

Leaflet 支持sf对象，因此您不应该尝试任何问题。不过，您可以将sf对象转换为sp（例如空间点数据框）as(YourObject, "Spatial")。对于从 sp 到 sf 的更改，您可以使用as_as_sf

library(sf)
library(sp)
data(meuse)
coordinates(meuse) = ~x+y
m.sf = st_as_sf(meuse) 
x = st_sfc(st_point(c(5,5)), st_point(c(6,9)), crs = 4326)
as(x, "Spatial")

看到这个

于 2019-09-10T03:14:37.087 回答