C#实现墨卡托投影坐标系经纬度与米制单位之间的互转

using System;
using GeoJSON.Net.Geometry;
 
namespace GISWebService.Common
{
    /// <summary>
    /// 墨卡托坐标转换（经纬度与米之间互转）
    /// </summary>
    public class MercatorCoordinatesConverter
    {
        /// <summary>
        /// 经纬度转Web墨卡托（单位：米）
        /// </summary>
        /// <param name="longitude">经度</param>
        /// <param name="latitude">纬度</param>
        /// <returns>转换后的位置</returns>
        public static Position Degree2WebMercatorMeter(double longitude, double latitude)
        {
            var xValue = longitude * 20037508.34 / ;
            var y = Math.Log(Math.Tan(( + latitude) * Math.PI / )) / (Math.PI / );
            var yValue = y * 20037508.34 / ;
            return new Position(xValue, yValue);
        }
 
        /// <summary>
        /// 经纬度转World墨卡托（单位：米）
        /// </summary>
        /// <param name="longitude">经度</param>
        /// <param name="latitude">纬度</param>
        /// <returns>转换后的位置</returns>
        public static Position Degree2WorldMercatorMeter(double longitude, double latitude)
        {
            const int radius = ;
            const double minorRadius = 6356752.314245179;
 
            const double d = Math.PI / ;
            const double r = radius;
            var y = latitude * d;
            const double tmp = minorRadius / r;
            double e = Math.Sqrt( - tmp * tmp),
                con = e * Math.Sin(y);
 
            var ts = Math.Tan(Math.PI /  - y / ) / Math.Pow(( - con) / ( + con), e / );
            y = -r * Math.Log(Math.Max(ts, 1E-));
 
            var xValue = longitude * d * r;
            var yValue = y;
 
            return new Position(xValue, yValue);
        }
 
        /// <summary>
        /// Web墨卡托转经纬度
        /// </summary>
        /// <param name="x">X坐标值（单位：米）</param>
        /// <param name="y">Y坐标值（单位：米）</param>
        /// <returns>转换后的位置</returns>
        public static Position WebMercatorMeter2Degree(double x, double y)
        {
            var xValue = x / 20037508.34 * ;
            var yValue = y / 20037508.34 * ;
            yValue =  / Math.PI * ( * Math.Atan(Math.Exp(yValue * Math.PI / )) - Math.PI / );
            var longitude = xValue;
            var latitude = yValue;
            return new Position(longitude, latitude);
        }
    }
}

代码原文连接 https://blog.csdn.net/a_dev/article/details/80990492