坐标系统的转换

《坐标系统的转换》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1.不同空间直角坐标系之间的转换不同空间直角坐标系之间的转换有三种模型，即布尔莎模型、范士模型、莫洛金斯基模型转换参数包括三个平移参数、三个旋转参数与一个尺度参数。对于坐标换算而言等价，推导布尔莎公式如下：如图所示，P在不同坐标系中的坐标XT＝⊿X0＋(1＋dK)R(e)X（1）式中XT——P在坐标系，OT—XTYTZT中的坐标向量X——P在坐标系O—XYZ中的坐标向量⊿X0——原点平移向量，⊿X0＝(⊿X⊿Y⊿Z)TdK——尺度变化系数，R(e)——旋转矩阵当已知转换参数⊿X0、dK、R(e)时，可按上式将Pi点的X坐标系坐标换

2、算为XT坐标系的坐标。按最小二乘原则求解转换参数⊿X0、dK、R(e)如下。因旋转角e很小，有sine=e和cose=1，若忽略e二阶微小量，则旋转阵代入（10-28）式，忽略二阶微小量dKQXi得XTi＝⊿X0＋R(e)dKXi＋R(e)Xi＝⊿X0＋(E＋Q)dKXi＋(E＋Q)Xi＝⊿X0＋dKXi＋Xi＋QXi顾及则（1）式为（此即用于两空间直角坐标系相互变换的布尔莎七参数公式）若上式中eX＝eY＝0，eZ≠0，则上式为五参数转换模型。若再有eZ＝0，则上式为四参数转换模型。若尺度比参数亦为零，则得三参数转换模型三参数转换

3、公式是在假设两坐标系间各坐标轴相互平行，即轴系间不存在欧勒角的条件下导出的，这在实际情况中往往是不可能的。在欧勒角不大，求得欧勒角误差和欧勒角本身数值属同一数量级时，可以近似地这样处置。此种情况在国内外一些坐标换算中屡见不鲜，如北美坐标系相对于地心坐标系的三参数是X0=-22m，Y0＝157m，Z0=176；欧洲坐标系相对于地心坐标系的三参数是X0＝-84m，Y0＝-103m，Z0＝-127m等。我国地心坐标系转换参数（DX-1）也属三个转换参数。当根据多个公共点按最小二乘法求解转换参数时，对每个点有观测方程设则误差方程法方程单位

4、权方差不同空间大地坐标系之间的转换对于不同大地坐标系的换算，除包含三个平移参数、三个旋转参数和一个尺度变化参数外，还包括两个地球椭球元素变化参数，以下推导不同大地坐标系的换算公式由式取全微分得式中（1）上式两端乘以并加以整理得：(2)式中顾及（1）式及（2）式可写为：(3)上式通常称为广义大地坐标微分公式或广义变换椭球微分公式。如略去旋转参数和尺度变化参数的影响，即简化为一般的大地坐标微分公式。根据3个以上公共点的两套大地坐标值，可列出9个以上(10-24)式的方程，可按最小二乘法求得8个转换参数。2.空间直角坐标与空间大地坐标的

5、转换图２－６表示了空间直角坐标系与空间大地坐标系之间的关系。图2-6　地球空间直角坐标系与大地坐标系在相同的基准下空间大地坐标系向空间直角坐标系的转换公式为：(2-1)式中，　，为椭球的长半轴，为椭球的卯酉圈曲率半径=6378.137km，为椭球的第一偏心率，为椭球的短半轴=6356.7523141km在相同的基准下空间直角坐标系向空间大地坐标系的转换公式为式中2.3不同平面坐标系转换众所周知，坐标系之间的差异主要取决于定位与定向、椭球参数以及坐标系的尺度定义。从原理上讲，严密方法是将旧网的全部观测资料重新归算到新坐标系中，重新平

6、差计算出各点的新坐标。对于二维转换模型，参数的选取依赖于工程项目的需要，在大多数的平面坐标转换应用中，常常使用四参数模型、直接参数法、六参数模型等。四参数模型四参数模型是从布尔莎公式演化而来的，其计算公式为:（1）式中，△x、△y、、θ、m分别为平面上的平移、旋转、尺度参数。不难看出，要求出4个参数至少需要2个已知公共点。当有2个以上转换公共点时，将此模型转换为线性模型用最小二乘求解：（2）式中,a=∆x，b=∆y，c=mcosθ,d=msinθ，m=a2+b2对于n个公共点可列如下误差方程:设所有转换点带有误差的观测值等权，则由

7、式(3)的误差方程，通过间接平差法求得转换参数向量x的最小二乘解为：x=BTPB-1BTPL（其中P为单位权），从而求出a、b、c、d。则平移参数为∆x=a，∆y=b，再用以下两式计算旋转参数θ和尺度因子m：θ=a2+b2，m=arctan(-c/d)直接参数法直接参数法就是利用两套坐标系两个已知公共点的坐标（X1,Y1）、（X2,Y2）、（x1,y1）、（x2,y2）求出坐标转换平移参数、尺度因子、旋转参数。数学模型∆X∆Y=X2Y2-X1Y1，∆x∆y=x2y2-x1y1S=∆X2+∆Y2,s=∆x2+∆y2,A=在此处键入公

8、式。,∂=平移参数(Dx)Dy=X1Y1-x1y1,尺度因子m=S-sS，旋转参数θ=A-∂则其他点（Xi,Yi）的坐标转换公式为：∆Xi∆Yi=XiYi-X1Y1∆Xi∆Yi=(1+m)cosθsinθ-sinθcosθ∆Xi∆Yixiyi=X1