已知四点经纬度位置计算球面面积

1 先把四点坐标的经纬度都转换为角度，四点坐标为：
    A(121.7    31.65)
    B(121.783   31.63)
    C(121.817   31.567)
    D(121.783    31.583)
2 要计算四边形的面积，实际就是计算两个三角形的面积和，即：
    S(四边形ABCD)=S(三角形ABC)+S(三角形BCD)
3 然后要计算球面三角形ABC的三条边的长度，使用球面两点的距离公式：
     l=R*(cos(a1-a2)*cos(b1)*cos(b2)+sin(b1)*sin(b2))
  其中R为地球半径，a1 a2为两点的经度，b1 b2为两点的纬度。