将坐标多边形内的区域转换为米

Question

给出一组坐标

 lat <- c(47.2325618, 47.2328269, 47.2330041, 47.2330481, 47.2330914, 
          47.2331172, 47.2331291, 47.2331499)
 lon <- c(11.3707441, 11.3707791, 11.3708087, 11.3708031, 11.3707818, 
          11.3707337, 11.3706588, 11.370284)
coords <- cbind(lon,lat)

我想计算多边形的面积。我使用包 splancs ：

中的函数areapl()

library(splancs)
areapl(coords)
# [1] 1.4768e-07

这给我留下了平方度（？）的维数。所以我的问题是：我如何将其转换为米/公里？

非常感谢。

Answer 1

使用spTransform包中的sp转换为笛卡尔网格系统（例如UTM区域），或在areaPolygon包中尝试geosphere。

> areaPolygon(coords)
[1] 7688.568