坐标、投影及坐标转换

空间坐标(球面)

地心坐标:WGS84,国家2000坐标。

参心坐标:北京54,西安80,地方坐标。

表示方式: 大地坐标(L,B,H);空间直角坐标(x,y,z)。这两种坐标都是基于球面的,前者单位是度,是常见的表示方式,后者单位是米,一般是为了同其他坐标进行转化时才出现。大地坐标和空间直角坐标之间可以相互转化,其转化公示及参数是公开。



投影(平面)

为了测绘及制图展示等需求,常把球面坐标投影到平面上,就需要用到地图投影,投影后的坐标采用平面直角坐标来表示。

表示方式: 平面直角坐标(x,y,h),单位是米。这里的平面直角坐标和空间直角坐标虽然表示类似,但是他们是不同的坐标。常见投影包括:高斯克里格投影、墨卡托投影、web墨卡托投影。

高斯克里格投影:像拨开的橘子瓣,我国大中比例尺地图均采用这种投影。

墨卡托投影:像棋盘,航海图一般采用这种投影。

web墨卡托投影:与墨卡托投影类似,互联网地图采用的投影。墨卡托投影与web墨卡托投影关系详见博文



常见数据的坐标和投影

GPS数据:一般是WGS84坐标,以大地坐标(经纬度)的形式展现,无投影。

地形图:一般是北京54或西安80坐标,以平面直角坐标(米)展现,高斯克里格投影。

互联网地图:WGS84坐标(google earth、google map)、加密WGS84坐标(对WGS84坐标进行加密,如GCJ-02坐标或百度坐标),web墨卡托投影。

市政管线数据:地方坐标,以平面直角坐标(米)展现,高斯克里格投影。



坐标转换过程

坐标转换涉及到空间坐标、平面坐标及投影。

1)三维转换:不同椭球体间转换,例如:WGS84<</span>——>北京54, WGS84<</span>——>西安54, 西安54<</span>——>北京54。常采用七参数法、三参数法。

2)二维转换:同一椭球体下转换,如:北京54<</span>——>地方坐标(该地方坐标是基于北京54建立的),常采用四参数法。

不同椭球体的坐标转换,需要用空间直角坐标来计算,采用三参数(1个点)或七参数(3个点);同一椭球体的坐标转换,基于高斯平面坐标来计算,采用四参数(2个点),如果公共点多可以采用最小二乘法拟合,求取最佳参数。

以WGS84坐标转地方坐标系为例,其流程为:WGS84——>北京54——>地方坐标,详细步骤如下:

1)(B,L,H)84——(X,Y,Z)84,空间大地坐标到空间直角坐标的转换。
2)(X,Y,Z)84——(X,Y,Z)54,坐标基准的转换,即Datum转换。通常有三种转换方法:七参数、简化三参数、Molodensky。
3)(X,Y,Z)54——(B,L,H)54,空间直角坐标到空间大地坐标的转换。
4)(B,L)54——(x,y)54, 高斯投影正算。
5)(x,y)54——(x0,y0)54,地方坐标转换,常用四参数法。

上述1-5步骤中,所有的转换公式都是公开的,但是转换参数有些是公开的,有些是非公开的。

公开参数步骤:1、3、4步骤需要的参数都是公开的。

待解参数步骤:2、4的参数是非公开的,需要找测绘部门去转换申请,或利用同名点求取转换的参数,步骤2,采用三参数(1个点)或七参数(3个点);步骤4,采用四参数(2个点),如果同名点多可以采用最小二乘法拟合,求取最佳参数。


坐标转换公开的参数:同一空间坐标系统,大地坐标和直角坐标的转换,空间坐标进行变换为投影坐标的参数。

坐标转换非公开的参数:不同空间坐标系统,各坐标之间的转化。


lua坐标转换 坐标值转换_最小二乘法



WGS84坐标转地方坐标流程图(图片来自宏大测绘的博客)

lua坐标转换 坐标值转换_坐标系转换_02



地方坐标转大地坐标流程图(蒋小军,2010)


转换工具

坐标转换过程中,坐标转换模型的参数是很关键的,需要找测绘部门购买或者自己计算。找同名点、计算转换参数、利用公式及参数进行坐标计算。

Excel,可以用了计算参数、坐标转换,可以实现批量转换,缺点有时有些公式编写比较困难,最小二乘法比较难实现。利用解一元多次方程,实现参数计算;利用普通公式计算,实现坐标批量转化。

ArcGIS,有投影工具,七参数的坐标转换工具,但是缺少转换参数,缺少参数计算工具。

Matlab,可实现最小二乘法拟合,求取最佳坐标转换参数。

此外,网上有一些坐标转换相关软件,很多实现了坐标转换的所有功能,但是有些软件要收费。

使用excel+arcGIS实现坐标转换。Excel中计算坐标转换参数。ArcGIS中提供了, Beijing_1954_To_WGS_1984工具,高斯投影工具、墨卡托投影工具、任意坐标转换工具(需要自己计算参数)。