0
我一直在使用以下代碼來讀取人口普查道數據併爲每個人口普查道產生質心數據。我想通過需要經度和緯度的Foursqaure API來抽取一些樣本的緯度和經度。我注意到這段代碼讓我回到了與經度不同的座標系。任何想法如何我可以提取緯度經度是否轉換座標或通過不同的功能?從座標數據轉換爲R中的經度和緯度?
#load libraries-----------------------------
library(rgdal)
library(sp)
library(GISTools)
library(ggplot2)
#load census tact file------------------------
tracts = readOGR("nycb2010.shp", layer="nycb2010")
#extract centroids
centroids = as.data.frame(getSpPPolygonsLabptSlots(tracts))
#map new data
ggplot() + geom_polygon(data=tracts, aes(x=long, y=lat, group=group), fill="black", colour="grey90", alpha = 1)+
geom_point(data=centroids, aes(x=V1, y=V2, group = 1, color = 'red'), size = .1)
沒關係我知道我沒怎麼傻使用新的變量。謝謝! – LoF10