2

当我尝试在 Matlab 与 R 中计算多边形面积时，我得到了两组不同的值。在 Matlab 中，我使用了 areaint（参考椭圆体设置为 WGS84）。在 R 中，我在 lat-lon 上尝试了 geosphere::areaPolygon，在网格转换的欧几里得坐标上尝试了 splancs::areapl。R 结果彼此一致，但与 Matlab 计算不一致。两组平均有 25% 的差异（Matlab 较低）。我错过了什么？

4

1 回答 1

1

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

另见： http ://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=reply&node=7581419

于 2012-10-26T00:12:57.220 回答