Matlab中的areaint与R中的areaPolygon或areapl

当我尝试在Matlab中计算多边形面积时，得到了两组不同的值。在Matlab中，我使用了areaint（参考椭球设置为WGS84）。在R中，我尝试了lat lon上的geosphere:：areaPolygon或网格转换的欧几里德坐标上的splancs:：areapl。R结果相互吻合，但与Matlab计算不符。这两组平均相差25%（Matlab更低）。我遗漏了什么？

显然，来自Matlab的结果来自格式不正确的输入矩阵：areaint（）需要（lat，lon），而R输入格式正确为（lon，lat）。这适用于前面提到的任何R函数

另见：

---

试试R-Sig-Geo邮件列表。不确定你会在这里得到答案。如果你尝试一个简单的多边形，你知道该区域会发生什么？谢谢，交叉张贴在R-Sig-Geo上@DennisJaheruddin使用？areaPolygon:p>>中的自包含示例，结果相差约20%（~40e6 km^2）。