Rの座標データから緯度および経度に変換する。

Question

私は以下のコードを使用して国勢調査データを読み込み、各国勢調査の重心データを生成しています。私は緯度と経度を必要とするFoursqaure APIを使っていくつかのサンプリングを行っているので、これから緯度と経度を引き出したいと思います。私は、このコードが私に緯度経度とは異なる座標系を返すことに気付きました。座標を変換しているのか、別の関数を使っているのか、緯度経度をどのように取り出すことができるのでしょうか？

#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) 

Answer 1

あなたはこのように、空間データの座標系を変更する機能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"