投影坐标转换.doc

摘自武测毕业论文（赖增先）——我也是从网上弄来，转载时请务必保留此此信息。
第二节平面坐标基准转换
由于海上和陆地上在测量时，使用不同的坐标系和不同参考椭球，而且采用的投影也不同，使得我们获得的数据不统一，必须进行坐标转换。
§3·2·1欧拉角
设有两个空间直角坐标系，分别为O-XYZ和O-XYZ，为了便于讨论其相应坐标轴间的变换，设其原点相同如图所示，选择、、为欧拉角，又称旋转参数，经过三次旋转，使两个坐标系重合，既：（图见下页A）
首先，绕OZ轴，将OX轴旋转到OX轴，所转的角为；
其次，绕OY轴，将OZ轴旋转到OZ轴，所转的角为；
最后，绕OX轴，将OZ轴旋转到OZ轴，所转的角为；
Z

ZZ

XO
X
XYY

Y
图A

因此有
XX
Y=R（）R（）R（）Y
ZZ
式中R（）、R（）、R（）为旋转矩阵，其表达式在、、很小时可以最终表示为：
X1X
Y=-1Y公式1
Z-1Z
§3·2·2不同三维空间直角坐标系的变换模型
GPS测量的WGS—84属地心坐标系，而1980年国家大地坐标系和1954年北京坐标系属参心坐标系，他们所对应得空间直角坐标系是不同的，这里将讨论不同空间直角坐标系的变换模型。
如图B两个空间直角坐标系分别为O-XYZ和O-XYZ，其坐标系原点不同则存在三个平移参数X、Y、Z，他们表示O-XYZ坐标系原点O相对于O-XYZ坐标系原点O在三个坐标轴上的分量；又当各坐标轴相互不平行时，既存在三个旋转参数、、。Z
Z
O
X

Y
OY

X
考虑到两个坐标系的平移和旋转以及尺度参数可得公式如下：
XX1X
Y=（1+m）Y-1Y
ZZ-1Z
X
+Y公式一
Z
式中共有七个变换参数X、Y、Z、、、、m,简称此公式为布尔莎七参数变换公式，是坐标变换中一个非常重要的公式。七参数变换公式，除了布尔莎公式外，还有莫洛琴斯基公式和范氏公式。这三种公式，它们之间的七个参数相差很大，但各自构成完整的数学模型，参数间存在着明确的解析关系，可以相互间转换。分别用它们来换算点的坐标时，其结果是完全相同的。因此，这三个公式是等价的。我国的地心坐标变换参数地心二号是七个变换参数，既采用布尔莎公式。
当公式一中===m=0,既称之为三参数公式。三参数公式表明两个空间直角坐标系尺度一致，且两个坐标轴相互平行。我国地心坐标变换参数地心一号系三个变换参数。同理在公式一中，略去某些参数，可分别得到四参数、五参数、六参数等坐标变换参数。公式一中的变换参数，一般利用公共点上的两套空间直角坐标系坐标值（X，Y，Z）和(X,Y,Z)即可采用最小二乘法解得。
应该指出，当进行两种不同空间直角坐标系变换时，坐标变换的精度除取决于坐标变换的数学模型和求解变换参数的公共点坐标精度外，还和公共点的多少、几何形状结构有关。鉴于地面网可能存在一定的系统误差，且在不同区域并非完全一样，所以采用分区变换参数，分区进行坐标转换，可以提高坐标变换精度。无论是从我国的多普勒网还是GPS网，利用布尔莎公式求解和地面大地网间得变换参数，分区变换均较明显地提高了坐标变换的精度。

§3·2·3不同三维大地坐标系的变换模型
对于不同的三维大地坐标系的变换模型，除了上节的七个变换参数外，还应增加两个变换参数，，这就是两个大地坐标系所对应的地球椭球参数的不同。不同大地坐标的变换公式，又称大地坐标微分公式或变换椭球微分公式。当包括旋转参数和尺度参数时，称为广义大地坐标微分公式或广义变换椭球微分公式。
空间一点的空间直角坐标与大地坐标关系式是：
X(N+H)cosBcosL

Y=(N+H)cosBsinL公式二
Z[N(1-e)+H]sinB
式中N为卯酉圈曲率半径。在这个公式中当已知L，B，H时，求X，Y，Z是非常简单的，只要代入公式即可。当已知X，Y，Z时反求L，B，H则可以采用直接解或迭代解法，解算时对公式做些变化即可。
由公式二最终我们可以得到不同三维大地坐标系的变换公式；

dL-0X
dB=--Y+

dHcosBcosLcosBsinLsinBZ

-1
-0+
-0
00
-m+
N+H-Nesin-
0da
公式三
df
式中dL、dB以弧度秒为单位，等式右端L、B、H均以换算前坐标值代入。公式三也就是顾及七个参数和椭球大小变化的广义大地坐标微分公式或广义变换椭球微分公式。略去旋转参数和尺度变化参数的影响，即为一般的大地坐标微分公式或椭球微分公式。
利用一些公共点上两套大地坐标系坐标值，采用最小二乘法可解得变换参数。

§3·2·4不同两维大地坐标系的变换模型
在三维不同大地坐标系的变换模型中，当进行WGS—84和我国参心大地坐标系的变换时，由于后者大地高的精度不高（一般在3m左右的误差），加之又难以确定其方差和协方差，因此，也可以考虑选择二维大地坐