使用Python实现了一下我们同事的C++高斯投影正反算,实际跑通,可用。 #!/ usr/bin/python # -*- coding:utf-8 -*- import math def LatLon2XY(latitude, longitude): a = 6378137.0 # b = 6356752.3142 # c = 6399593.6258 # alpha = 1 / 298.257223563 e2 = 0.0066943799013 # epep = 0.00673949674227 #将经纬度转换为弧度 latitude2Rad = (math.pi / 180.0) * latitude beltNo = int((longitude + 1.5) / 3.0) #计算3度带投影度带号 L = beltNo * 3 #计算中央经线 l0 = longitude - L #经差 tsin = math.sin(latitude2Rad) tcos = math.cos(latitude2Rad) t = math.tan(latitude2Rad) m = (math.pi / 180.0) * l0 * tcos et2 = e2 * pow(tcos, 2) et3 = e2 * pow(tsin, 2) X = 111132.9558 * latitude - 16038.6496 * math.sin(2 * latitude2Rad) + 16.8607 * math.sin( 4 * latitude2Rad) - 0.0220 * math.sin(6 * latitude2Rad) N = a / math.sqrt(1 - et3) x = X + N * t * (0.5 * pow(m, 2) +
python实现高斯投影正反算
最新推荐文章于 2025-06-02 13:03:38 发布