关于GPS 坐标系的那些事

一、投影

将GPS坐标投影的方式有很多种，等角度的有两种：高斯-克吕格投影、UTM投影

关于高斯-克吕格投影的资料，这两份相对比较全一些，资料1、资料2

了解的内容主要是几点：

1、X Y轴：中央经线为X轴，赤道为Y轴

2、角度：从北顺时针到直线的夹角为方位角

3、坐标的表示：为了方便区分带号，在每个点横坐标Y值前加所在的带号

二、经纬度的基本知识

1、度分秒坐标与经纬度的转换

经纬度1度=60分=3600秒

例子：109度10分52秒转为10进制109.181111

109

10/60=0.1666666

52/3600=0.014444

相加得：109.181111

109.181111转度分秒 109.181111-109=0.181111 然后再用0.181111*60=10.8 取整10，（0.1811111-10/60）*3600=52

2、坐标格式

A、经纬度。这个是球面坐标，对于北京来说，就是(116.388171,39.935961)这样的坐标。比如腾讯、高德、百度都是这样的经纬度坐标。谷歌是经纬度顺序写反的经纬度坐标。

墨卡托坐标。是一种投影坐标, 又称为高斯-克吕格投影（Gauss-KrugerProjection），它不是用经纬度来表示的，数字一般都比较大，像这样的：(215362.00021333335, 99526.00034912192)，多用于航海图航空图。墨卡托坐标主要用于程序的后台计算。直线距离，加加减减计算方便。

3、大地坐标系与地理坐标系有何不同？

大地坐标系和地理坐标系都是经纬度表示的坐标系，本身并不包含投影信息，很多时候，这两种说法都是相同的。

4、世界大地坐标系统

WGS-84是世界大地坐标系

5、其他

https://www.sohu.com/a/121612743_488161

6、经纬度坐标格式：

经纬度格式分为三种：

度：（ ddd.ddddd °）十进制小数部分（5位）百度地图，谷歌地图，手机自带gps经纬度输出

度 . 分：（ddd°mm.mmm’ ）十进制小数部分（3位） GPGGA

度 . 分 . 秒：（ddd°mm’ss’’）谷歌地图

关系：

一度（°）等于60分钟（’）等于3600秒（“）：

整度（d）等于十进制的度的整数部分（dd）：

分钟（m）等于十进制的度的整数部分（dd）减去整数度（d）的60倍：

秒（s）等于十进制的度（dd）减去整数度（d）减分（M）除以60乘以3600：

三、高斯投影中角度都是弧度

代码：度转弧度

'''
度变成弧度
'''
def deg_to_rad(x):deg,min,sec = x.split(' ')latDeg = float(deg) + float(min)/60 + float(sec)/3600print latDegreturn (float(deg) + float(min)/60 + float(sec)/3600)/180*pi#return float(x)*pi/180def central_meridian_diff(lon):deg,min,sec = lon.split(' ')lonDeg = float(deg) + float(min) / 60 + float(sec) / 3600print lonDeg#lonDeg = float(lon)n = trunc(lonDeg/6)+1lon0 = n*6 - 3l = (lonDeg-lon0)/180*pireturn l,  n