投影坐标转换

第二节 平面坐标基准转换

由于海上和陆地上在测量时，使用不同的坐标系和不同参考椭球，而且采用的投影也不同，使得我们获得的数据不统一，必须进行坐标转换。

§3·2·1 欧拉角

设有两个空间直角坐标系，分别为O-XYZ 和O-X 'Y 'Z '，为了便于讨论其相应坐标轴间的变换，设其原点相同如图所示，选择εx 、y ε、z ε为欧拉角，又称旋转参数，经过三次旋转，使两个坐标系重合，既：（图见下页A ）

首先，绕O Z '轴，将O X '轴旋转到OX 0轴，所转的角为z ε；

其次，绕OY 0轴，将O Z '轴旋转到OZ 0轴，所转的角为y ε；

最后，绕OX 轴，将O Z 0轴旋转到OZ 轴，所转的角为εx ；

Z

Z 0 Z '

X ' O

X 0

X Y 0 Y

Y '

图A

因此有

X X '

Y = R 1（εx ）R 2（y ε）R 3（z ε） Y '

Z Z '

式中 R 1（εx ）、R 2（y ε）、R 3（z ε）为旋转矩阵，其表达式在ε、y ε、z ε很小时可以最终表示为：

X 1 z ε y ε X '

Y = -z ε 1 εx Y ' 公式1

Z y ε - εx 1 Z '

§3·2·2 不同三维空间直角坐标系的变换模型

GPS 测量的WGS —84属地心坐标系，而1980年国家大地坐标系和1954年北京坐标系属参心坐标系，他们所对应得空间直角坐标系是不同的，这里将讨论不同空间直角坐标系的变换模型。

如图B 两个空间直角坐标系分别为O-XYZ 和O '-X 'Y 'Z '，其坐标系原点不同则存在三个平移参数∆X 0、∆Y 0、∆Z 0，他们表示O '- X 'Y 'Z '坐标系原点O '相对于O-XYZ 坐标系原点O 在三个坐标轴上的分量；又当各坐标轴相互不平行时，既存在三个旋转参数εx 、

Y '

考虑到两个坐标系的平移和旋转以及尺度参数可得公式如下：

X X ' 1 z ε y ε X '

Y =（1+m ） Y ' -z ε 1 εx Y ' Z Z ' y ε - εx 1 Z '

∆X 0

+ ∆Y 0 公式一

∆Z 0

式中共有七个变换参数∆X 0、∆Y 0、∆Z 0、εx 、y ε、z ε、m,简称此公式为布

尔莎七参数变换公式，是坐标变换中一个非常重要的公式。七参数变换公式，除了布尔莎公式外，还有莫洛琴斯基公式和范氏公式。这三种公式，它们之间的七个参数相差很大，但各自构成完整的数学模型，参数间存在着明确的解析关系，可以相互间转换。分别用它们来换算点的坐标时，其结果是完全相同的。因此，这三个公式是等价的。我国的地心坐标变换参数地心二号是七个变换参数，既采用布尔莎公式。

当公式一中ε

x =

y

ε=

z

ε=m=0,既称之为三参数公式。三参数公式表明两个空间

直角坐标系尺度一致，且两个坐标轴相互平行。我国地心坐标变换参数地心一号系三个变换参数。同理在公式一中，略去某些参数，可分别得到四参数、五参数、六参数等坐标变换参数。公式一中的变换参数，一般利用公共点上的两套空间直角坐

标系坐标值（X，Y，Z）

i 和(X',Y', Z')

i

即可采用最小二乘法解得。

应该指出，当进行两种不同空间直角坐标系变换时，坐标变换的精度除取决于坐标变换的数学模型和求解变换参数的公共点坐标精度外，还和公共点的多少、几何形状结构有关。鉴于地面网可能存在一定的系统误差，且在不同区域并非完全一样，所以采用分区变换参数，分区进行坐标转换，可以提高坐标变换精度。无论是从我国的多普勒网还是GPS网，利用布尔莎公式求解和地面大地网间得变换参数，分区变换均较明显地提高了坐标变换的精度。

§3·2·3 不同三维大地坐标系的变换模型

对于不同的三维大地坐标系的变换模型，除了上节的七个变换参数外，还应增加两个变换参数，，这就是两个大地坐标系所对应的地球椭球参数的不同。不同大地坐标的变换公式，又称大地坐标微分公式或变换椭球微分公式。当包括旋转参数和尺度参数时，称为广义大地坐标微分公式或广义变换椭球微分公式。

空间一点的空间直角坐标与大地坐标关系式是：

X (N+H)cosBcosL

Y = (N+H)cosBsinL 公式二

Z [N (1-e2)+H]sinB

式中N为卯酉圈曲率半径。在这个公式中当已知L，B，H时，求X，Y，Z是非常简单的，只要代入公式即可。当已知X，Y，Z时反求L，B，H则可以采用直接解或迭代解法，解算时对公式做些变化即可。

由公式二最终我们可以得到不同三维大地坐标系的变换公式；

dL -B H N L cos )(sin +''ρ B H N L cos )(cos +''ρ 0 ∆X 0 dB = -

H M L B +cos sin ''ρ -H N L B +sin sin ''ρ H

M B +cos ''ρ ∆Y 0 + dH cosBcosL cosBsinL sinB ∆Z 0

L tgB H N H e N cos )1(2++- L tgB H

N H e N sin )1(2++- -1 εx -L H M B Ne H N sin sin )(22+-+ L H

M B Ne H N cos sin )(22+-+ 0 y ε + -"2sin cos sin ρL

B B Ne "2cos cos sin ρL B B Ne 0 z ε

0 0

-"2cos sin ρB B e H

M N + m+ "2cos sin )(ρB B e a H M N + N+H-Ne 2sin B 2 -)sin 1(22B e a

N - 0 da

"22cos sin )

1)(()sin 2(ρB B f H M B e M -+- 公式三 B B e f

M 222sin )sin 1(1-- df 式中dL 、dB 以弧度秒为单位，等式右端L 、B 、H 均以换算前坐标值代入。公式三也就是顾及七个参数和椭球大小变化的广义大地坐标微分公式或广义变换椭球微分公式。略去旋转参数和尺度变化参数的影响，即为一般的大地坐标微分公式或椭