如何计算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来计算geosphere包中多边形的面积。它以平方米为单位，但您可以转换为平方公里。

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

你可以找到许多有用的工具 CRAN Task View: Analysis of Spatial Data