Wie berechnet man den Bereich einer Karte in R?

Question

Ich versuche, die Fläche einer Wasserscheide mit Hilfe der Paketkarten zu berechnen, aber keinen Erfolg.

Der Link das Polygon zum Download ist: 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) 

Die Antwort durch den Code gegeben ist 4,37, aber es sollte etwas über 50000 km² sein.

Answer 1

Es gibt eine Funktion areaPolygon, um die Fläche eines Polygons im Geosphärenpaket zu berechnen. Es gibt die Fläche in Quadratmetern, aber Sie können in Quadratkilometer umwandeln.

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

Sie können viele nützliche Werkzeuge an der finden CRAN Task View: Analysis of Spatial Data