R에서지도의 면적을 계산하는 방법은 무엇입니까?

Question

패키지 맵을 사용하여 유역 면적을 계산하려하지만 성공하지 못했습니다. 코드에 의해 주어진 https://drive.google.com/open?id=0BwVpSqmgvCe-a01ra2RtcTVZTGs

install.packages("maps"); 
library(maps) 
tmarias <- readShapePoly("TrêsMarias_Poligon.shp") 
m <- map(tmarias, fill = TRUE, plot = TRUE) 
area.map(m) 

대답은 4.37이지만, 무언가에 대한 50000 평방 킬로미터해야한다 : 다각형을 다운로드 할 수

링크입니다.

Answer 1

지리 정보 패키지의 다각형 영역을 계산하는 함수는 areaPolygon입니다. 면적은 평방 미터이지만 sq.km로 변환 할 수 있습니다.

library(geosphere) 
areaPolygon(tmarias) 
[1] 50766929868 
areaPolygon(tmarias)/1000000 
[1] 50766.93  ## sq km 

당신은 CRAN Task View: Analysis of Spatial Data

에 많은 도움이되는 도구를 찾을 수 있습니다