我是R語言的新手，也是處理人口普查資料的新手，我試圖建立一個可以傳遞給另一個團隊的csv檔案，顯示人口普查區的長/緯度輪廓。
有了這個檔案，我可以找到佛羅里達州的3個特定的人口普查區。
FL_blocks < -塊("FL"/span>。 年= 2010, ）
FL_blocks_Alachua < -過濾器(FL_blocks, COUNTYFP == "001"）
NTIA_FL_CB <- FL_blocks_Alachua %>%
filter(FL_blocks_Alachua$GEOID10 %in% c ("120010019071008"/span>。 "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 %>%
重命名(經度= X, 緯度= Y） %> %
select(Latitude,Longitude)
# save lat/long as csv
st_write(NTIA_FL_CB_Shape_xy, "fl_3cb_ntia_latlong. csv", coords = TRUE）
# 繪制3個普查區的輪廓，并將shapefiles的輪廓標為綠色圓圈。
leaflet(NTIA_FL_CB) ％>；％
addTiles(/span>) %>%
addPolygons(popup = ~GEOID10) ％> %
addCircleMarkers(data =NTIA_FL_CB_Shape_xy。 color = "green"）
這個圖片顯示了我的想法，如果我可以把勾勒出每個人口普查區的緯度/長度傳給使用我的csv檔案輸出的團隊，他們應該可以把這個檔案加載到他們的GIS軟體中進行覆寫。
uj5u.com熱心網友回復： 訣竅是投到點上，然后輸出。或者，如果你的團隊正在使用GIS軟體，你可以將你的多多邊形寫成一個shapefile？
標籤：library(sf)
nc <- st_read(system.file("shape/nc. shp",包="sf"）。 安靜 = TRUE)
nc_pt <- st_cast(nc, "POINT")
st_write(nc_pt, "nc_pt。 csv", layer_options="GEOMETRY=AS_XY"）
# st_write(nc, "nc.shp")
