用经纬度计算多边形的面积

Question

我有这段代码，是用这个写的：source1 and this: source 2

public static double CalculatePolygonArea(IList<GpsLocation> coordinates)
    {
        double area = 0;

        if (coordinates.Count > 2)
        {
            for (var i = 0; i < coordinates.Count-1; i++)
            {
                GpsLocation p1, p2;
                p1 = coordinates[i];
                p2 = coordinates[i + 1];
                area += ToRad(p2.Longitude - p1.Longitude) * (2 + Math.Sin(ToRad(p1.Latitude))
                   + Math.Sin(ToRad(p2.Latitude)));

                area = area * R * R / 2;
            }
        }

        return Math.Abs(area);
    }

这是我的测试代码：

[Fact]
    public void GpsPolygonAreaTest()
    {
        var poly = new List<GpsLocation>();
        var p1 = new GpsLocation(0, 0);
        poly.Add(p1);
        var p2 = GpsHelper.CreateLocationBasedOnBearingDistance(p1, 5, 100);
        poly.Add(p2);
        var p3 = GpsHelper.CreateLocationBasedOnBearingDistance(p2, 95, 100);
        poly.Add(p3);
        var p4 = GpsHelper.CreateLocationBasedOnBearingDistance(p3, 185, 100);
        poly.Add(p4);
        poly.Add(p1);

        var area = GpsHelper.CalculatePolygonArea(poly);

        area.Should().Be(10000);
    }

我确认我的多边形是 100m x 100m，见图：

我的测试结果是：期望值是10000，结果是1.28153883377486E+48。

知道我的代码有什么问题吗？

Answer 1

我很确定这个说法：

area = area * R * R / 2;

应该放在顶点循环之后，而不是在每次迭代中重复。

用经纬度计算多边形的面积

Calculate the area of a polygon with latitude and longitude

google-maps

polygon

latitude-longitude

area