根据输入的地点坐标计算中心点

       在开发中，需要根据输入的地点坐标计算中心点，但是在百度，Google上搜索“根据输入的地点坐标计算中心点”或者“计算地图中心点”等等一系列关键字，折腾了很久，也没有找到任何解决方法。不过还好，最后在Google搜索“Latitude and longitude of the center”得到了解决方案，因为解决方案是在英文网站上找到的，所以将解决方案整理出来，供大家参考（呵呵，看来有些东西还是需要到国外去取取经）。

先给出地点坐标类的定义。

public class GeoCoordinate
{
    private readonly double latitude;
    private readonly double longitude;

    public double Latitude 
    { 
        get 
        { 
            return latitude; 
        } 
    }
    public double Longitude 
    { 
        get 
        { 
            return longitude;
        } 
    }

    public GeoCoordinate(double latitude, double longitude)
    {
        this.latitude = latitude;
        this.longitude = longitude;
    }

    public override string ToString()
    {
        return string.Format("{0},{1}", Latitude, Longitude);
    }
}

下面给出完整的计算方法。

/// <summary>
/// 根据输入的地点坐标计算中心点
/// </summary>
/// <param name="geoCoordinateList"></param>
/// <returns></returns>
public GeoCoordinate GetCenterPointFromListOfCoordinates(List<GeoCoordinate> geoCoordinateList)
{
    int total = geoCoordinateList.Count;
    double X = 0, Y = 0, Z = 0;
    foreach (GeoCoordinate g in geoCoordinateList)
    {
        double lat, lon, x, y, z;
        lat = g.Latitude * Math.PI / 180;
        lon = g.Longitude * Math.PI / 180;
        x = Math.Cos(lat) * Math.Cos(lon);
        y = Math.Cos(lat) * Math.Sin(lon);
        z = Math.Sin(lat);
        X += x;
        Y += y;
        Z += z;
    }
    X = X / total;
    Y = Y / total;
    Z = Z / total;
    double Lon = Math.Atan2(Y, X);
    double Hyp = Math.Sqrt(X * X + Y * Y);
    double Lat = Math.Atan2(Z, Hyp);
    return new GeoCoordinate(Lat * 180 / Math.PI, Lon * 180 / Math.PI);
}

对于400km以下时，可以采用下面的简化计算方法。

/// <summary>
/// 根据输入的地点坐标计算中心点（适用于400km以下的场合）
/// </summary>
/// <param name="geoCoordinateList"></param>
/// <returns></returns>
public GeoCoordinate GetCenterPointFromListOfCoordinates(List<GeoCoordinate> geoCoordinateList)
{
    //以下为简化方法（400km以内）
    int total = geoCoordinateList.Count;
    double lat = 0, lon = 0;
    foreach (GeoCoordinate g in geoCoordinateList)
    {
        lat += g.Latitude * Math.PI / 180;
        lon += g.Longitude * Math.PI / 180;
    }
    lat /= total;
    lon /= total;
    return new GeoCoordinate(lat * 180 / Math.PI, lon * 180 / Math.PI);
}

在以下页面包含有多种实现，大家可以参考。

http://stackoverflow.com/questions/6671183/calculate-the-center-point-of-multiple-latitude-longitude-coordinate-pairs

 详细的算法说明，可以参考。

http://www.geomidpoint.com/calculation.html

扩展阅读。

根据两点经纬度计算距离