投影坐标转换

第二节 平面坐标基准转换

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

§3·2·1 欧拉角

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

首先，绕O Z'轴，将O X'轴旋转到OX

0轴，所转的角为z；

其次，绕OY0轴，将O Z'轴旋转到OZ0轴，所转的角为y； 最后，绕OX轴，将O Z

0轴旋转到OZ轴，所转的角为x；

Z

Z0 Z'

X' O X0 X Y0 Y

Y' 图A

因此有

X X'

Y = R1（x）R2（y）R3（z） Y' Z Z'

式中 R1（x）、R2（y）、R3（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年国家大地坐标系和19年北京坐标系属参心坐标系，他们所对应得空间直角坐标系是不同的，这里将讨论不同空间直角坐标系的变换模型。

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

Z O' X' Y'

O Y

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

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

X0

+ Y0 公式一 Z0

式有七个变换参数X0、Y0、Z0、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 - dB = -

sinLcosL'' '' 0 X0

(NH)cosB(NH)cosBsinBcosL''sinBsinL''cosB - '' Y0 +

MHNHMHdH cosBcosL cosBsinL sinB Z0

N(1e2)HN(1e2)HtgBcosL tgBsinL -1 x NHNH(NH)Ne2sin2B(NH)Ne2sin2BsinL cosL 0 y + -MHMH-

Ne2sinBcosBsinL\"

Ne2sinBcosBcosL\" 0 z

0 0 -NNe2sinBcosB\" e2sinBcosB\" m+ (MH)aMHN+H-Ne2sin2B -N(1e2sin2B) a0 da M(2e2sin2B)sinBcosB\" 公式三

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

'.

.

球微分公式。

利用一些公共点上两套大地坐标系坐标值，采用最小二乘法可解得变换参数。

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

在三维不同大地坐标系的变换模型中，当进行WGS—84和我国参心大地坐标系的

变换时，由于后者大地高的精度不高（一般在3m左右的误差），加之又难以确定其方差和协方差，因此，也可以考虑选择二维大地坐标系的变换模型。

所谓二维大地坐标系，即当大地高H为零时的椭球面上的大地坐标系。其变换模型，只要在公式三中，将H=0代入即可得到。

将此公式用于GPS网和地面网联合平差时，如果顾及地面网的系统性观测误差对网的定向的影响时，则可在椭球面上引入一个附加旋转参数dA,以使两网更好的配合。由于dA产生的对dL、dB的影响加于公式右端。

为了便于比较GPS网和地面网的大地坐标，若在将GPS网的X、Y、Z反算L、B、H时，采用了地面网的椭球参数，即两网相应的椭球参数已化为一致，则公式中不再有 da、df项。

§3·2·5 不同二维高斯投影平面坐标系的坐标转换

由高斯投影正算公式可得：

xxdx dL LByydy = dB 公式四

LB式中等号右端偏导数由高斯投影正算公式得：

x =NsinBcosBl

Lx1 =M[1+(1-2sin2 B)l2]

B2y1 =NcosB[1+(1-2sin2 B) l2] 公式五

L2y =-MsinBl

B上式中，l=L-L0, L0为子午线得大地经度，公式五中dB、dL见公式三。对于不同二维高斯投影平面坐标系坐标差的模型，可以由公式三和公式四给出。

§3·2·6 同一参考系统下的高斯直角坐标、 大地坐标、空间直角坐标之间的相互转换

'.

.

当运用了§3·2·1～§3·2·5后，我们就可以将一个系统的坐标，转化到另一个系统的对应结果。可以完成对应之间的坐标转化，但是，如果高斯直角坐标、 大地坐标、空间直角坐标之间的相互转换就需要用本节的内容。 3·2·6·1 高斯直角坐标同大地坐标之间的转化

完成两者之间的转化，要用高斯正、反算方法。这里并不详细介绍，只给出其数学模型：

11正算公式：x=X+N*t*cos2B*l2+N*t(5- t2+92+44)cos4B*l4

2241 +N*t(61-58 t2+ t4+270 2-330 2t2)cos6* l6

72011 y=N*cosB*l+N(1-t2+2)cos3B* l3+N(5-18 t2

6120+ t4+14反算公式：B=Bf-2-58 -y

22t2)cos5B*l5

tf2MfNftf+24

tf24MfN3f422(5+3t2f+2f-9tff)y

-

720MfN5f(61+90t2f+45t4f)y6

l=

11223y-(1+2t+)y ff3NfcosBf6NfcosBf +

1242225(5+28t+24t+6+8t) y fffff5120NfcosBf通过高斯正、反算，可以将高斯直角坐标化算为大地坐标，或将大地坐标化算为高斯直角坐标。

3·2·6·2 大地坐标同空间直角坐标的化算

空间一点的空间直角坐标与大地坐标关系式是： 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则可以采用直接解或迭代解法，解算时对公式做些变化即可。

3·2·6·3高斯直角坐标与空间直角坐标的相互转换

'.

.

关于这两者的转换关系并没有直接给出，但是，我们可以利用两者同大地坐标的转换关系，先把其中的一种化算成大地坐标，然后再由大地坐标转换成另一种坐标。

§3·2·7 不同投影之间的转换

陆地地形图采用高斯投影系统，而海图采用的是墨卡托投影系统，这使得两种地形图不能拼接。在此以由高斯投影系统转化到墨卡托投影系统，给出转化的思路和过程：

首先利用布尔莎七参数变换公式（公式一）将一参考系内的高斯坐标转换为另一参考系内的高斯坐标；再利用高斯反算公式，将转换后的高斯直角坐标换算为大地坐标；然后在同一参考系内进行墨卡托投影将大地坐标化算为该平面内的平面坐标。这里考虑的是高斯坐标同墨卡托平面坐标不在同一参考椭球上，如果在同一参考椭球内，则第一步化算可以省略。

'.