高斯克吕格投影

在网上搜集整理有关高斯克吕格投影的代码并用 python 改写。

from math import *
p = 180.0 / pi * 3600def deg2rad(deg):dms = deg*3600/preturn dms
def gauss_zhengsuan(B,L,L0): # 纬度、经度、中央经线经度l = L-L0            l = deg2rad(l) #经度B = deg2rad(B) #纬度a = 6378245.0 # Beijing1954 Krasovsky_1940',6378245.0,298.3f = 1.0/298.3b = a - a*fc = a**2/be = sqrt(a**2 - b**2)/ae1 = sqrt(a**2 - b**2)/bp = 3600*180/piBeta0 = 1 - (3.0/4)*e1**2 + (45.0/64)*e1**4 - (175.0/256)*e1**6 + (11025.0/16384)*e1**8Beta2 = Beta0-1Beta4 = (15.0/32)*e1**4 - (175.0/384)*e1**6 + (3675.0/8192)*e1**8Beta6 = -(35.0/96)*e1**6 + (735.0/2048)*e1**8Beta8 = (315.0/1024)*e1**8Xb0 = c*(Beta0*(B) + sin(B)*(Beta2*cos(B) + Beta4*cos(B)**3 + Beta6*cos(B)**5 + Beta8*cos(B)**7))Np = a/sqrt(1 - (e*sin(B))**2)m0 = l*cos(B)t = tan(B)in2 = (e1)**2*cos(B)**2X = Xb0 + (1.0/2)*Np*t*m0**2 + (1.0/24)*(5-t**2 + 9.0*in2 + 4.0*in2**2)*Np*t*m0**4 + (1.0/720)*(61 - 58.0*t**2)*Np*t*m0**6.0Y = Np*m0 + (1.0/6)*(1 - t**2 + in2)*Np*m0**3 + (1.0/120)*(5 - 18.0*t**2 + t**4 + 14.0*in2 - 58.0*t**2*in2)*Np*m0**5 + 500000return Y,X # X表示纵坐标，Y表示横坐标result = gauss_zhengsuan(32.34, 120.123, 120)

高斯克吕格投影 python2.x 版本相关推荐

【测绘程序设计】高斯克吕格投影：带号及中央经度计算神器V1.0（附源程序）
[问题描述]:很多情况下,我们知道某一地点的坐标(经纬度),需要计算其在高斯克吕格投影中的带号及中央经度.关于该问题,有具体的公式可言,只是计算过程稍微繁琐一些,当然啦,我们可以写程序来解决,谁叫我们 ...
Win10下python3和python2多版本同时安装并解决pip共存问题
Win10下python3和python2多版本同时安装并解决pip共存问题参考文章: (1)Win10下python3和python2多版本同时安装并解决pip共存问题 (2)https://ww ...
高斯-克吕格投影与UTM投影
高斯-克吕格投影与UTM投影高斯-克吕格(Gauss-Kruger)投影与UTM投影(Universal Transverse Mercator,通用横轴墨卡托投影)都是横轴墨卡托投影的变种,目前一 ...
MATLAB实现高斯-克吕格投影正算
MATLAB实现高斯-克吕格投影正算-即经纬度转为x和y 高斯-克吕格投影简介更新2020-06,重新整理一下脚本函数高斯-克吕格投影,是由德国数学家.物理学家.天文学家高斯于1822年代首次提出 ...
ArcGIS 坐标系统文件 ---beijing 1954 高斯克吕格投影在arcgis中的说明
坐标是GIS数据的骨骼框架,能够将我们的数据定位到相应的位置,为地图中的每一点提供准确的坐标. ArcGIS自带了多种坐标系统,在${ArcGISHome}/Coordinate Systems/目录 ...
高斯-克吕格投影与地形图分带
1.高斯-克吕格投影的概念.为了将地球椭球面上的各种量,如方向.长度归算到地图平面上相应的量,就要采用地图投影的数学方法.一般在大于或等于1:50万比例尺的地形图中我国使用高斯-克吕格投影(或简称高斯 ...
北京54和西安80投影坐标系，高斯-克吕格投影
1.首先理解地理坐标系(Geographic coordinate system),Geographic coordinate system直译为地理坐标系统,是以经纬度为地图的存储单位的.很明显, ...
高斯——克吕格投影正算
高斯--克吕格投影正算
高斯——克吕格投影反算
高斯--克吕格投影反算
北京1954-3度分带-高斯克吕格投影
北京1954年3度分带高斯克吕格投影是一种投影方式,用于将地球表面投影到平面图上.这种投影方式经常用于地图制作,在保留经纬度信息的同时减小了对形状和大小的失真.这种投影方式是在北京基准椭球上进行的,因 ...

高斯克吕格投影 python2.x 版本

高斯克吕格投影

高斯克吕格投影 python2.x 版本相关推荐

最新文章

热门文章