根据经纬度计算多边形面积

计算方法比较简单，主要是求出多边形外接矩形已米为单位面积和已经纬度为单位面积比值，然后用这个比值乘以多边形经纬度为单位面积，即可得出这个多边形以米为单位面积。

double GetArea(const vector<Coordinate>& ls)
{
    if (ls.size() < 4)
        return 0;

    double sum = 0;

    for (size_t i=0; i<ls.size()-1; ++i)
    {
        const Coordinate& p = ls[i];
        const Coordinate& q = ls[i+1];
        sum += (p.x + q.x) * (q.y - p.y);
    }

    return sum/2;
}

double GetPrjArea(const vector<Coordinate> &ls)
{
    if (ls.size() < 4)
        return 0;

    double dArea = GetArea(ls);
    dArea = abs(dArea);
    if (dArea == 0)
        return 0;

    double xmin, ymin, xmax, ymax;
    xmin = xmax = ls[0].x;
    ymax &