2016-11-07 129 views
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) 

回答

0

您可以使用函數spTransform改變你的空間數據的座標系統,就像這樣:

library(sp) 
library(rgdal) 

lat_long = CRS("+init=epsg:4326") 
tracts_lat_long <- spTransform(tracts, lat_long) 

proj4string(tracts_lat_long) 
[1] "+init=epsg:3035 +proj=laea +lat_0=52 +lon_0=10 +x_0=00 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" 
+1

沒關係我知道我沒怎麼傻使用新的變量。謝謝! – LoF10

相關問題