공간 객체에서 위도와 경도 가져 오기 R

Question

셰이프 파일에서 위도와 경도를 얻고 싶습니다. 지금까지는 셰이프 파일을 읽는 방법 만 알고 있습니다.

library(rgdal) 
centroids.mp <- readOGR(".","35DSE250GC_SIR") 

하지만 centroids.mp에서 위도와 경도를 어떻게 추출 할 수 있습니까?

Answer 1

이 질문에 대한 몇 가지 수준이있다.

경도와 위도를 묻지 만이 좌표계는이 개체에서 사용하지 않을 수 있습니다. 당신이이 SpatialPointsDataFrame, 모든 라인의 목록입니다 이것은 SpatialLinesDataFrame 경우 좌표, 그냥 무게 중심 경우 경우 "무게 중심"이 좌표의 모든 것이라고이

coordinates(centroids.mp) 

주와 같은 좌표를 얻을 수 있습니다 이것은 SpatialPolygonsDataFrame입니다.

좌표는 경도와 위도 일 수 있지만 개체는이를 알지 못할 수 있습니다. 사용

proj4string(centroids.mp) 

"NA"인 경우 개체는 (A)를 알지 못합니다. "+ proj = longlat"을 포함하면 객체는 알 수 있으며 경도/위도 (B)입니다. "+ proj ="와 다른 이름 ("longlat"이 아님)이 포함되어 있으면 객체가 알기 때문에 경도/위도 (C)가 아닙니다.

(A) 알아야 할 경우 또는 값에서 분명한 경우 일 수 있습니다.

(B) 완료 한 경우 (먼저 가정을 확인해야하지만,이 메타 데이터가 잘못 될 수 있음).

하는 경우 (C) 할 수 있습니다 (꽤 안정적으로 먼저 가정을 확인해야하지만) 이런 식으로 (데이텀 WGS84에) 경도 위도로 변환 :

그것은 당신의 질문은 전혀 "GDAL"를 포함 분명하지 않다

coordinates(spTransform(centroids.mp, CRS("+proj=longlat +datum=WGS84"))) 

Answer 2

사용과 같은 coordinates() :

library(maptools) 
xx <- readShapePoints(system.file("shapes/baltim.shp", package="maptools")[1]) 
coordinates(xx) 
#  coords.x1 coords.x2 
# 0  907.0  534.0 
# 1  922.0  574.0 
# 2  920.0  581.0 
# 3  923.0  578.0 
# 4  918.0  574.0 
#  [.......]