我是 R 新手并使用人口普查数据,我正在尝试构建一个 csv 文件,该文件可以传递给另一个团队,以显示人口普查块的纬度/经度轮廓。
有了这个,我可以到达佛罗里达州的 3 个特定人口普查区。
FL_blocks <- blocks("FL", year = 2010, )
FL_blocks_Alachua <- filter(FL_blocks, COUNTYFP == "001")
NTIA_FL_CB <- FL_blocks_Alachua %>%
filter(FL_blocks_Alachua$GEOID10 %in% c ("120010019071008", "120010019071007", "120010019071009"))
现在我想制作一个表格,只显示纬度/经度和纬度/经度关联的相应人口普查块。这将为我提供每个人口普查块的纬度/经度列表,我可以用人口普查块的多边形绘制它们,但是为了将其传递给不熟悉 R 的团队,我需要 csv 数据输出。
NTIA_FL_CB_Shape_xy <- as.data.frame(st_coordinates(NTIA_FL_CB$geometry))
NTIA_FL_CB_Shape_xy <- NTIA_FL_CB_Shape_xy %>%
rename( Longitude = X, Latitude = Y) %>%
select(Latitude,Longitude )
# save lat/long as csv
st_write(NTIA_FL_CB_Shape_xy, "fl_3cb_ntia_latlong.csv", coords = TRUE)
# plot the 3 census blocks with the outline of the shapefiles marked as green circles
leaflet(NTIA_FL_CB) %>%
addTiles() %>%
addPolygons(popup = ~GEOID10) %>%
addCircleMarkers(data = NTIA_FL_CB_Shape_xy, color = "green")
这张图片显示了我的想法,如果我可以通过概述每个人口普查块的纬度/经度,那么使用我的 csv 文件输出的团队应该能够将其加载到他们的 GIS 软件中进行叠加。