当前位置:文档之家› 不同大地坐标系间坐标转换模型研究

不同大地坐标系间坐标转换模型研究

不同大地坐标系间坐标转换模型研究
不同大地坐标系间坐标转换模型研究

空间大地坐标系与平面直角坐标系转换公式

§2.3.1 坐标系的分类 正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即采用什么方法来表示空间位置。人们为了描述空间位置,采用了多种方法,从而也产生了不同的坐标系,如直角坐标系、极坐标系等。 在测量中常用的坐标系有以下几种: 一、空间直角坐标系 空间直角坐标系的坐标系原点位于参考椭球的中心,Z 轴指向参考椭球的北极,X 轴指向起始子午面与赤道的交点,Y 轴位于赤道面上且按右手系与X 轴呈90°夹角。某点在空间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。空间直角坐标系可用图2-3来表示: 图2-3 空间直角坐标系 二、空间大地坐标系 空间大地坐标系是采用大地经、纬度和大地高来描述空间位置的。纬度是空间的点与参考椭球面的法线与赤道面的夹角;经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角;大地高是空间点沿参考椭球的法线方向到参考椭球面的距离。空间大地坐标系可用图2-4来表示:

图2-4空间大地坐标系 三、平面直角坐标系 平面直角坐标系是利用投影变换,将空间坐标空间直角坐标或空间大地坐标通过某种数学变换映射到平面上,这种变换又称为投影变换。投影变换的方法有很多,如横轴墨卡托投影、UTM 投影、兰勃特投影等。在我国采用的是高斯-克吕格投影也称为高斯投影。UTM 投影和高斯投影都是横轴墨卡托投影的特例,只是投影的个别参数不同而已。 高斯投影是一种横轴、椭圆柱面、等角投影。从几何意义上讲,是一种横轴椭圆柱正切投影。如图左侧所示,设想有一个椭圆柱面横套在椭球外面,并与某一子午线相切(此子午线称为中央子午线或轴子午线),椭球轴的中心轴CC ’通过椭球中心而与地轴垂直。 高斯投影满足以下两个条件: 1、 它是正形投影; 2、 中央子午线投影后应为x 轴,且长度保持不变。 将中央子午线东西各一定经差(一般为6度或3度)范围内的地区投影到椭圆柱面上,再将此柱面沿某一棱线展开,便构成了高斯平面直角坐标系,如下图2-5右侧所示。 图2-5 高斯投影 x 方向指北,y 方向指东。 可见,高斯投影存在长度变形,为使其在测图和用图时影响很小,应相隔一定的地区,另立中央子午线,采取分带投影的办法。我国国家测量规定采用六度带和三度带两种分带方法。六度带和三度带与中央子午线存在如下关系: 366 N L =中; n L 33=中 其中,N 、n 分别为6度带和3度带的带号。

坐标转换工具说明书-1208

§10.2坐标转换工具 HGO 数据处理软件包提供了坐标转换程序,可以进行地方坐标与WGS-84坐标的相互转换,同时具备参数求解功能。 下面对这个工具进行介绍: 10.2.1概述 首先,介绍一下常见的三种坐标表示方法:经纬度和椭球高(BLH),空间直角坐标(XYZ),平面坐标和水准高程(xyh/NEU)。注意:椭球高是一个几何量,而水准高是一个物理量。 我们通常说的WGS-84坐标是经纬度和椭球这一种,北京54坐标是平面坐标和水准高程这一种,实质是有平面基准和高程基准组成的。 此外,再注意一下坐标转换的严密性问题,在同一个椭球里的纯几何转换都是严密的(BLH<->XYZ),而在不同的基准之间的转换是不严密的。举个例子,在WGS-84坐标和北京54坐标之间是不存在一套转换参数可以全国通用的,因为前者是一个地心坐标系,后者是一个参心坐标系。高程转换是由几何高向物理高转换。因此在每个地方必须用椭球进行局部拟合,通常用7参数模型来拟合。 那么,两个椭球间的坐标转换应该是怎样的呢?一般而言比较严密的是用七参数法(或称布尔莎模型),即X平移,Y平移,Z平移,X旋转,Y旋转,Z旋转,尺度变化K。要求得七参数就需要在一个地区需要3个以上的已知点(7个参数至少7个方程可解,所以需要三个点列出9个方程),如果区域范围不大、最远点间的距离不大于30Km(经验值)的情况可以用三参数,即X平移,Y平移,Z平移,而将X旋转,Y旋转,Z旋转,尺度变化K视为0,所以三参数只是七参数的一种特例。 七参数模型的实质是用一个局部椭球去拟合地方坐标系的形态;所以转换后获得的地方椭球高就是水准高。当然我们也可以把平面和高程两个方向分别进行拟合。例如平面用四参数模型拟合,高程方向则用二次曲面等模型来拟合。这样分开处理的模式相对七参数模型自由度更高。但是由于四参数模型参数较少,表达能力较弱,通常只用于小区域坐标转换。 综上所述,从实用的角度出发,坐标转换程序提供了两种转换策略供给客户选择使用: 1.七参数模型,一步得到地方平面和水准数据。 2.四参数加高程拟合模型,分两步得到地方平面和水准数据。 由于各厂家的模型和流程定义可能是不一样的,这里就我们公司的转换流程描述如下:七参数的转换过程是这样的:

不同空间直角坐标系的转换

不同空间直角坐标系的转换 欧勒角 不同空间直角坐标系的转换,包括三个坐标轴的平移和坐标轴的旋转,以及两个坐标系的尺度比参数,坐标轴之间的三个旋转角叫欧勒角。 三参数法 三参数坐标转换公式是在假设两坐标系间各坐标轴相互平行,轴系间不存在欧勒角的条件下得出的。实际应用中,因为欧勒角不大,可以用三参数公式近似地进行空间直角坐标系统的转换。公共点只有一个时,采用三参数公式进行转换。

七参数法 用七参数进行空间直角坐标转换有布尔莎公式,莫洛琴斯基公式和范氏公式等。下面给出布尔莎七参数公式: 坐标转换多项式回归模型 坐标转换七参数公式属于相似变换模型。大地控制网中的系统误差一般呈区域性,当区域较小时,区域性的系统误差被相似变换参数拟合,故局部区域的坐标转换采用七参数公式模型是比较适宜的。但对全国或一个省区范围内的坐标转换,可以采用多项式回归模型,将各区域的系统偏差拟合到回归参数中,从而提高坐标转换精度。 两种不同空间直角坐标系转换时,坐标转换的精度取决于坐标转换的数学模型和求解转换系数的公共点坐标精度,此外,还与公共点的分布有关。鉴于地面控制网系统误差在???? ??????+??????????=??????????000111222Z Y X Z Y X Z Y X ???? ??????+????????????????????---+??????????+=??????????000111111222000)1(Z Y X Z Y X Z Y X m Z Y X X Y X Z Y Z εεεεεε

不同区域并非是一个常数,所以采用分区进行坐标转换能更好地反映实际情况,提高坐标转换的精度。

坐标转换之计算公式

坐标转换之计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释: A :参心空间直角坐标系: a) 以参心0为坐标原点; b) Z 轴与参考椭球的短轴(旋转轴)相重合; c) X 轴与起始子午面和赤道的交线重合; d) Y 轴在赤道面上与X 轴垂直,构成右手直角坐标系0-XYZ ; e) 地面点P 的点位用(X ,Y ,Z )表示; B :参心大地坐标系: a) 以参考椭球的中心为坐标原点,椭球的短轴与参考椭球旋转轴重合; b) 大地纬度B :以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ; c) 大地经度L :以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ; d) 大地高H :地面点沿椭球法线至椭球面的距离为大地高H ; e) 地面点的点位用(B ,L ,H )表示。 2 参心大地坐标转换为参心空间直角坐标: ?? ???+-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中,N 为椭球面卯酉圈的曲率半径,e 为椭球的第一偏心率,a 、b 椭球的长短半 径,f 椭球扁率,W 为第一辅助系数 a b a e 2 2-= 或 f f e 1*2-= W a N B W e =-=22sin *1( 3 参心空间直角坐标转换参心大地坐标

[]N B Y X H H e N Y X H N Z B X Y L -+=+-++==cos ))1(**)()(*arctan( )arctan(2 2222 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形 高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关; 3. 离中央子午线越远变形越大 为控制投影后的长度变形,采用分带投影的方法。常用3度带或6度带分带,城市或工 程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式: 5 2224253 2236 4254 42232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24 cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++=) 3、高斯投影反算公式:

GPS坐标和国家大地坐标之间的转换

GPS坐标和国家大地坐标之间的转换 一、前言 WGS-84坐标系是目前GPS所采用的坐标系统,GPS所发布的星历参数就是基于此坐标系统的。WGS-84坐标系统的全称是World Geodical System-84(世界大地坐标系-84),它是一个地心地固坐标系统。WGS-84坐标系统由美国国防部制图局建立,于1987年取代了当时GPS所采用的坐标系统-WGS-72坐标系统而成为GPS的所使用的坐标系统。WGS-84坐标系的坐标原点位于地球的质心,轴指向BIH1984.0定义的协议地球极方向,轴指向BIH1984.0的启始子午面和赤道的交点。 采用椭球参数为:a=6 378 137m,f= 1/298.257 223 563。 北京54 坐标系、西安80 坐标系—属于参心坐标系, 北京54 坐标系采用克拉索夫斯基椭球参数,长轴a= 6 3 78 2 4 5 米, 扁率f=l : 2 98.3 ;西安80 大地系坐标系椭球参数采用国际大=地测量和地球物理联合19 7 5 后推荐的地球椭球参数, 长轴a= 6 3 7 8 140 米, 扁率f1 : 298.257,大地原点在我西安市径阳县永乐镇。西安80 坐标系的建立是在54 年北京坐标系的基础上完成的。 在实际的工作中,对于GPS的测量数据。我们需要将其转换成所需要的54或80坐标系,才能够使用。或是将其转换成相应的地方坐标系。在转换的过程中需要进行一系列的变换。本文将对其过程做详细的说明。 二、转换过程 (1)数据测量:在实际操作中,首先进行的是数据的观测。根据实际工作需要,采用相应的观测方法进行观测,得到合格的测量成果。本文主要是针对GPS控制网的转换来说明的。 (2)平差:在GPS控制网的测量工程中,在进行完基线测量(地面坐标和高程)后,需要对测量结果进行平差,得到相应的平差结果。下面对相应的条件平差①做具体说明: AV-W=0 [1] L#=L+V [2] 基础方程和它的解: 设有r个平差线性条件方程: [3]

大地坐标转换成施工坐标公式

大地(高斯平面)坐标系工程坐标系转换大地坐标系--->工程坐标系 ======================== 待转换点为P,大地坐标为:Xp、Yp 工程坐标系原点o: 大地坐标:Xo、Yo 工程坐标:xo、yo 工程坐标系x轴之大地方位角:a dX=Xp-Xo dY=Yp-Yo P点转换后之工程坐标为xp、yp: xp=dX*COS(a)+dY*SIN(a)+xo yp=-dX*SIN(a)+dY*COS(a)+yo 工程坐标系--->大地坐标系 ======================== 待转换点为P,工程坐标为:xp、yp 工程坐标系原点o: 大地坐标:Xo、Yo 工程坐标:xo、yo 工程坐标系x轴之大地方位角:a dx=xp-xo dy=yp-yo P点转换后之工程坐标为xp、yp: xp=Xo+dx*COS(a)-dy*SIN(a)

yp=Yo+dx*SIN(a)+dy*COS(a) 坐标方位角计算程序 置镜点坐标:ZX ZY 后视点坐标:HX HY 方位角:W 两点间距离: S Lb1 0← {A, B, C, D}← A〝ZX=〞:B〝ZY=〞:C〝HX=〞:D 〝HY=〞:W=tg1((D-B)÷(C-A)):(D-B)>0=>(C-A)>0=>W=W:∟∟(D-B)>0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)>0=>W=360+W∟∟W=W◢ S=√((D-B)2+(C-A)2) ◢ Goto 0← CASIO fx-4500p坐标计算程序 根据坐标计算方位角 W=W+360△W:“ALF(1~2)=”L1 A“X1=”:B“Y1=”:Pol(C“X2”-A,D“Y2”-B:“S=”▲W<0 直线段坐标计算 L1 X“X(0)”:Y“Y(0)”:S“S(0)”:A“ALF” L2 Lb1 2 L3 {L}:L“LX”

ARCGIS中坐标转换

ArcGIS 坐标转换 1.坐标分析 问题:对于某地A中心点坐标为455299.845,3223622.525的CAD矩形,CAD施工图。将其转换为WGS-84坐标,如何转换? 分析:分析455299.845为6位,则为东向Y坐标,省去了带号,加上了5000000加常数,其最大为为4,说名在中央子午线的左侧(左侧为负值,加上500万后肯定小于500万,首位为4。若在中央子午线右侧,则最大位数为5);3223622.525为7位,为北向X坐标。 查看“某地A”的经度为92.5度,因为为CAD施工图,比例尺肯定大于1:5万,所以为3度带,所以此点的中央子午线为93E,带号为Beijing_54_Zone_31。 2.CAD转为shp格式并设定坐标系: ArcTool box-Convesion Tools->To Geodatabse->CAD to Geodatabase: 其中空间参考坐标系选择Beijing_1954_3_Degree_GK_CM_93E。 具体原因:选择投影坐标系-Gauss Kruger-Bei Jing54,此时3度带有两种:Beijing_1954_3_Degree_GK_CM_93E和Beijing_54_Zone_31,前者表示中央子午线为93E的3度带,后者表示北京54 31度带,二者意义一样,但选择哪种呢?因为点坐标东向为455299.845为6位,不带带号,因此选择Beijing_1954_3_Degree_GK_CM_93E(若东向坐标

为31455299.845,则选择Beijing_54_Zone_31), 3.北京54到WGS84坐标的转换 1.1加载图层: 打开ArcTool box-Data Management Tools->Project and transformation->feature->Project,加载shp图层,弹出下列窗口: 出现红色“X”号,说明原始图层坐标系没有识别出,则需要首先设定其坐标系后再转换。具体设坐标系参考“9 设置或改变Shp文件坐标系” 1.2选择输出图层地址和名称: 在Out Put Dataset or Feature处输入输出图层名:

直角坐标系下的画图及其转换公式

直角坐标系下的画图及其转换公式 在直角坐标系下我们的圆方程是: 222()()x a y b R -+-= 其中,a 和b 是圆心,R 是半径。但在画圆的时候,你就会发现如果按该公式画圆,多半是不成功的,或者画了一半,所以在matlab 中画圆,一半采用极坐标形式 圆对应的极坐标转换公式为: cos sin x R y R θ θ =?? =?(公式1) 这个很容易理解,你画个单位圆来看看就知道了。 那么上面那个黑色的点的x 坐标和y 坐标用半径和连线与坐标轴x 的夹角来表示,就得到了公式1。 观察这个公式,我们发现,在极坐标系下,圆的半径没变,夹角是在不断变化的,所以,在matlab 中极坐标系下画单位圆的问题可以这样来考虑: 首先将夹角360等分,也就是每一个步长为360度/360; 但需要指出的是,matlab 中正弦预先函数的变量其实是弧度,并不是度。这个你在matlab 命令窗里就可以试: 比如你要得到30度的正弦值,一般是sin (pi/6),而不是sin(30)。这里的pi 是3.1415926的在matlab 中的表示。 所以我们的步长应该是弧度制的,我们知道,1度对应的弧度为360/(2*pi)。也即180/pi; 所以我们的夹角应该是: Theta=0:180/pi:2*pi-180/pi; 注意,由于是从零开始画图的,所以最后一个应该是2*pi-180/pi;而不是2*pi ; 这个时候我们可以开始画图了 X=R*cos(Theta); Y=R*sin(Theta); Plot(x,y,’r.’) axis square %保证画出来的圆是圆的。

大地坐标与直角空间坐标转换计算公式

大地坐标与直角空间坐标转换计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释: A :参心空间直角坐标系: a) 以参心0为坐标原点; b) Z 轴与参考椭球的短轴(旋转轴)相重合; c) X 轴与起始子午面和赤道的交线重合; d) Y 轴在赤道面上与X 轴垂直,构成右手直角坐标系0-XYZ ; e) 地面点P 的点位用(X ,Y ,Z )表示; B :参心大地坐标系: a) 以参考椭球的中心为坐标原点,椭球的短轴与参考椭球旋转轴重合; b) 大地纬度B :以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ; c) 大地经度L :以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ; d) 大地高H :地面点沿椭球法线至椭球面的距离为大地高H ; e) 地面点的点位用(B ,L ,H )表示。 2 参心大地坐标转换为参心空间直角坐标: ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中,N 为椭球面卯酉圈的曲率半径,e 为椭球的第一偏心率,a 、b 椭球的长短半径,f 椭球扁率,W 为第一辅助系数 a b a e 2 2-= 或 f f e 1 *2-= W a N B W e = -=22 sin *1( 西安80椭球参数: 长半轴a=6378140±5(m )

短半轴b=6356755.2882m 扁 率α=1/298.257 3 参心空间直角坐标转换参心大地坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形 高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关; 3. 离中央子午线越远变形越大 为控制投影后的长度变形,采用分带投影的方法。常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式: 52224253 2236 425442232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++ =) 3、高斯投影反算公式:

参考系坐标系及转换汇总

1 天球坐标系、地球坐标系和卫星测量中常用的坐标系的建立方法。天球直角坐标系 天球坐标系 天球球面坐标系 坐标系 地球直角坐标系 地球坐标系 地球大地坐标系 常用的天球坐标系:天球赤道坐标系、天球地平坐标系和天文坐标系。在天球坐标系中,天体的空间位置可用天球空间直角坐标系或天球球面坐标系两种方式来描述。 1 天球空间直角坐标系的定义 地球质心O为坐标原点,Z轴指向天球北极,X轴指向春分点,Y轴垂直于XOZ平面,与X轴和Z轴构成右手坐标系。则在此坐标系下,空间点的位置由坐标(X,Y,Z)来描述。 春分点:当太阳在地球的黄道上由天球南半球进入北半球,黄道与赤道的交 点).

2 天球球面坐标系的定义 地球质心O为坐标原点,春分点轴与天轴(天轴:地球自转的轴)所在平面为天球经度(赤经)测量基准——基准子午面,赤道为天球纬度测量基准而建立球面坐标。空间点的位置在天球坐标系下的表述为(r,α,δ)。

表示:2-1天球空间直角坐标系与天球球面坐标系的关系可用图

岁差和章动的影响 岁差:地球实际上不是一个理想的球体,地球自转轴方向不再保持不变,这 使春分点在黄道上产生缓慢的西移,这种现象在天文学中称为岁差。章动:在日月引力等因素的影响下,瞬时北天极将绕瞬时平北天极旋转,大致呈椭圆,这种现象称为章动。 极移:地球自转轴相对地球体的位置并不是固定的,因而,地极点在地球表面上的位置,是随时间而变化的,这种现象称为极移。地球的自转轴不仅受日、月引力作用而使其在空间变化,而且还受地球内部质量不均匀影响在地球内部运动。前者导致岁差和章动,后者导致极移。 协议天球坐标系:为了建立一个与惯性坐标系统相接近的坐标系,人们通常选择某一时刻,作为标准历元,并将此刻地球的瞬时自转轴(指向北极)和地心至瞬时春分点的方向,经过瞬时的岁差和章动改正后,分别作为X轴和Z轴的指向,。协议天球坐标系由此建立的坐标系称为 3 地球坐标系

高斯平面直角坐标系与大地坐标系相互转化

高斯平面直角坐标系与大地坐标系相互转化 高斯平面直角坐标系与大地坐标系转换 1. 高斯投影坐标正算公式(1) 高斯投影正算:已知椭球面上某点的大地坐标(L,B),求该点在高斯投影平面上的直角坐标(x,y),即(L,B)->(x,y)的坐标变换。(2) 投影变换必须满足的条件中央子午线投影后为直线; 中央子午线投影后长度不变; 投影具有正形性质,即正形投影条件。(3) 投影过程在椭球面上有对称于中央子午线的两点P 1 和P 2 ,它们的大地坐标分别为(L,B)及(l,B),式中l 为椭球面上P 点的经度与中央子午线(L 0 )的经度差:l=L-L 0 ,P 点在中央子午线之东,l 为正,在西则为负,则投影后的平面坐标一定为P 1 ’(x,y)和P 2 ’(x,-y)。(4) 计算公式 4 ' ' 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 9 5 ( cos sin 2 sin 2 l t B B N Bl N X x 5 ' ' 4 2 5 5 ' ' 3 ' ' 2 2 3 ' ' ' ' ' ' ) 18 5 ( cos 120 ) 1 ( 6 cos l t t B N l t B N Bl N y 当要求转换精度精确至0.001m时,用下式计算: 6 ' ' 4 2 5 6 ' ' 4 ' ' 4 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 58 61 ( cos sin 720 ) 4 9 5 ( cos sin 24 sin 2 l t t B B N l t B B N Bl N X x 5 ' ' 2 2 2 4 2 5 5 ' ' 3 ' ' 2 2 3 3 ' ' ' ' ' ' ) 58 14 18 5 ( cos 720 ) 1

坐标转换模型

坐标转换模型 1.空间直角坐标系间的转换模型(七参数模型) ①公式(布尔莎模型): ②分析: (1)将O-XYZ中的长度单位缩放l+m倍,使其与O'-X'Y'Z'的长度单位一致; (2)从X反向看向原点O,以O为旋转点,让O-XYZ绕X轴顺时针旋转Wx角,使经过旋转后的Y轴与O'-X'Y'Z’平面平行; (3)从Y反向看向原点O,以O为旋转点,让O-XYZ绕Y轴顺时针旋转Wy角,使经过旋转后的X轴与O'-X'Y'Z'平面平行。显然,此时Z轴也与Z'轴平行; (4)从Z反向看向原点O,以O点为旋转点,O-XYZ绕Z轴顺时针旋转Wz角,使经过旋转后的X轴与X’轴平行。显然,此时O-XYZ的三个坐标轴己与O'-X'Y'Z’中相应的坐标轴平行; 原坐标为O-XYZ,转换到新坐标O-X’Y’Z’.(两坐标系都为空间直角坐标系)其中(dX dY dZ)为坐标原点的平移参数,即将坐标O-XYZ的原点分别沿三个坐标轴平移-dX,-dY,-dZ,使原坐标轴与O-X’Y’Z’的点重合。m为尺度参数,(w1 w2 w3)分别为坐标轴的旋转参量(角度),构成的旋转矩阵分别为: 分别将R1 R2 R3代入上式,可得:

当旋转角度w1 w2 w3很小时(<=10),cos(w)=1,sin(w)=0;在误差允许范围内可以将模型简化为:(同样七参数模型) 四参数模型是在七参数模型的特例,没有考虑坐标轴的旋转量,只考虑坐标轴的平移。 总结: 类似布尔莎模型(以坐标原点为参考点),还有莫洛金斯基坐标模型(以目标点为变换中心)、武测转换模型和范士转换模型(以控制网参考点的站心地平坐标系的三个坐标轴为旋转轴),这些坐标转换模型很容易实现相关坐标在不同坐标系的转换,但是参考位置的偏移向量的相关参数,在实际运用中这些参量是很难测定的,并且受地球重力等物理因素的影响,两个坐标系统即使经过相似变换,仍可能存在较大的残差,所以这些模型适用于简单且规则模型中。 ④程序: clc clear all dX=input('please input value of dX=');

参考系坐标系及转换

1天球坐标系、地球坐标系和卫星测量中常用的坐标系的建立方法。 L天球直角坐标系 厂天球坐标系 天球球面坐标系 地球直角坐标系地球大地坐标系 常用的天球坐标系:天球赤道坐标系、天球地平坐标系和天文坐标系。 在天球坐标系中,天体的空间位置可用天球空间直角坐标系或天球球面坐标系两种方式来描述。 1天球空间直角坐标系的定义 地球质心0为坐标原点,Z轴指向天球北极,X轴指向春分点,丫轴垂直于XOZ 平面,与X轴和Z轴构成右手坐标系。则在此坐标系下,空间点的位置由坐标(X,丫Z)来描述。 春分点:当太阳在地球的黄道上由天球南半球进入北半球,黄道与赤道的交点)

A <空闵直笥坐瑟厂K V : z 丿的楚辽” 2天球球面坐标系的定义 地球质心0为坐标原点,春分点轴与天轴(天轴:地球自转的轴)所在平面为天 球经度(赤经)测量基准一一基准子午面,赤道为天球纬度测量基准而建立球面 坐标。空间点的位置在天球坐标系下的表述为(r ,a,S )。 天欢申诗与地球质?M 重合T 赤礙刊为舍天黏 和感分点的天球子牛面 与过天体$的天球子牛面 之间的夾角,未纬 S 为 原点Mi 天体£的连規与 天球击道面之间的夹角, 旬題丫为展点Mi 天体S 球球】?坐抚1就,S 1 r )的C 义: 天球空间直角坐标系与天球球面坐标系的关系可用图 2-1表示: 感鼻—地I 球质心M 一孑塾一指向天球北奴Pn 、 ¥菇'一垂直于XMZ 平面, 与X 抽和Z 抽枸成右 手坐 标系统。 Pn A Z y X 1 \y X 奋 My\5 Ps / /

对同一空间点,直角坐标糸与其著效的球面坐标糸参教间有如下转换关务: C X - /cos a cos S < Y= / sin cos -Z = ysin 5 Y V a = arctan —— L Xz d -arctail . 岁差和章动的影响 岁差:地球实际上不是一个理想的球体,地球自转轴方向不再保持不变,这使春分点在黄道上产生缓慢的西移,这种现象在天文学中称为岁差。 章动:在日月引力等因素的影响下,瞬时北天极将绕瞬时平北天极旋转,大致呈椭圆,这种现象称为章动。 极移:地球自转轴相对地球体的位置并不是固定的,因而,地极点在地球表面上的位置,是随时间而变化的,这种现象称为极移。地球的自转轴不仅受日、月引力作用而使其在空间变化,而且还受地球内部质量不均匀影响在地球内部运动。 前者导致岁差和章动,后者导致极移。 协议天球坐标系:为了建立一个与惯性坐标系统相接近的坐标系,人们通常选择某一时刻,作为标准历元,并将此刻地球的瞬时自转轴(指向北极)和地心至瞬 时春分点的方向,经过瞬时的岁差和章动改正后,分别作为 X轴和Z轴的指向, 由此建立的坐标系称为协议天球坐标系。天味奋 5 y X X Ps

北京54坐标系转换工具

北京54坐标系转换工具 利用ARCGIS进行自定义坐标系和投影转换 ARCGIS种通过三参数和其参数进行精确投影转换 注意:投影转换成54坐标系需要下载无偏移卫星图像进行转换,有偏移的转换将导致转换后的卫星图像扭曲,坐标错误,无法配准。 第一步:选择无偏移地图源,下载你所需要的卫星图像。 第二步:选择BIGEMAP软件右边工具栏,选择【投影转换】,如下图所示: 2.1 选择说明: 1. 源文件:选择下载好的卫星图像文件(下载目录中后缀为tiff的文件) 2. 源坐标系:打开的源文件的投影坐标系(自动读取,不需要手动填写) 3. 输出文件:选择转换后你要保持文件的文件路径和文件名 4. 目标坐标系:选择你要转换成的目标坐标系,如下图:

选择上图的更多,如下图所示: 1:选择 -Beijing 1954 2:选择地区3:选择分度带对应的带号(一般默认,也可以手动修改)

选择对应的分度带或者中央子午线(请参看:如何选择分度带?),点击【确定】 5. 重采样算法:投影转换需要将影像的像素重新排列,一次每种算法的效率不一样,一般选择【立方卷积采样】,以达到最好的效果。如下图: 6. 指定变换参数:在不知道的情况下,可以不用填此处信息,如果√上,则如下图:

此参数为【三参数】或者【七参数】,均为国家保密参数,需要到当地的测绘部门或者国土部门,以单位名义签保密协议进行购买,此参数各地都不一样,是严格保密的,请不要随便流通。 第三步:点击【确定】,开始转换,如下图:

第四步:完成后,打开你刚才选择的输出文件夹,里面就是转换后的卫星图像。 第五步:如果你需要套合你手里已经有的矢量文件,请参看:【BIGEMAP无偏移影像叠加配准】

不同坐标系之间的变换

不同坐标系之间的变换 SANY GROUP system office room 【SANYUA16H-

§10.6不同坐标系之间的变换 10.6.1欧勒角与旋转矩阵 对于二维直角坐标,如图所示,有: ?? ? ?????????-=??????1122cos sin sin cos y x y x θθθθ(10-8) 在三维空间直角坐标系中,具有相同原点的两坐标系间的变换一般需要在三个坐标平面上,通过三次旋转才能完成。如图所示,设旋转次序为: ①绕1OZ 旋转Z ε角,11,OY OX 旋 转至0 0,OY OX ; ②绕0 OY 旋转Y ε角 10 ,OZ OX 旋转至0 2 ,OZ OX ; ③绕2OX 旋转X ε角, 0,OZ OY 旋转至22,OZ OY 。 Z Y X εεε,,为三维空间直角坐标变换的三个旋转角,也称欧勒角,与 它相对应的旋转矩阵分别为: ???? ? ?????-=X X X X X R εεεεεcos sin 0sin cos 00 01 )(1 (10-10)

????? ?????-=Y Y Y Y Y R εεεεεcos 0sin 010sin 0cos )(2 (10-11) ???? ? ?????-=10 0cos sin 0sin cos )(3Z Z Z Z Z R εεεεε (10-12) 令 )()()(3210Z Y X R R R R εεε= (10- 13) 则有: ???? ? ?????=??????????=??????????1110111321222)()()(Z Y X R Z Y X R R R Z Y X Z Y X εεε (10-14) 代入: ???? ??? ??? +-+++--=Y X Z Y X Z X Z Y X Z X Y X Z Y X Z X Z Y X Z X Y Z Y Z Y R εεεεεεεεεεεεεεεεεεεεεεεεεεεεεcos cos sin sin cos cos sin cos sin cos sin sin cos sin sin sin sin cos cos cos sin sin sin cos sin sin cos cos cos 0一般Z Y X εεε,,为微小转角,可取: sin sin sin sin sin sin sin ,sin ,sin 1cos cos cos =========Z Y Z X Y X Z Z Y Y X X Z Y X εεεεεεεεεεεεεεε 于是可化简

坐标转换之计算公式

创作编号: GB8878185555334563BT9125XW 创作者: 凤呜大王* 坐标转换之计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释: A :参心空间直角坐标系: a) 以参心0为坐标原点; b) Z 轴与参考椭球的短轴(旋转轴)相重合; c) X 轴与起始子午面和赤道的交线重合; d) Y 轴在赤道面上与X 轴垂直,构成右手直角坐标系0-XYZ ; e) 地面点P 的点位用(X ,Y ,Z )表示; B :参心大地坐标系: a) 以参考椭球的中心为坐标原点,椭球的短轴与参考椭球旋转轴重合; b) 大地纬度B :以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ; c) 大地经度L :以过地面点的椭球子午面与起始子午面之间的夹角为大地经度 L ; d) 大地高H :地面点沿椭球法线至椭球面的距离为大地高H ; e) 地面点的点位用(B ,L ,H )表示。 2 参心大地坐标转换为参心空间直角坐标: ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中,N 为椭球面卯酉圈的曲率半径,e 为椭球的第一偏心率,a 、b 椭球的长短半径,f 椭球扁率,W 为第一辅助系数

a b a e 2 2-= 或 f f e 1 *2-= W a N B W e = -=22 sin *1( 3 参心空间直角坐标转换参心大地坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形 高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关; 3. 离中央子午线越远变形越大 为控制投影后的长度变形,采用分带投影的方法。常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式:

坐标转换器使用说明

大地坐标(BLH) 平面直角坐标(XYZ) 四参数:X 平移、Y 平移、旋转角和比例 七参数:X平移,Y平移,Z 平移,X 轴旋转,Y 轴旋转,Z 轴旋转,缩放比例(尺度比) GPS控制网是由相对定位所求的的基线向量而构成的空间基线基线向量网,在GPS控制网的平差中,是以基线向量及协方差为基本观测量。 图3-1表示为HDS2003数据处理软件进行网平差的基本步骤,从图中可以看到,网平差实际上可以分为三个过程: l、前期的准备工作,这部分是用户进行的。即在网平差之前,需要进行坐标系的设置、并输入已知点的经纬度、平面坐标、高程等。 2、网平差的实际进行,这部分是软件自动完成的; 3、对处理结果的质量分析与控制,这部分也是需要用户分析处理的过程。 图3-1 平差过程 坐标系选择 针对不同的平差,要相应选择不同的坐标系,是否输入相应信息。在笔者接触过的项目中,平差时先通过三维无约束平差后,再进行二维约束平差。由于先进行的时三维无约束平差,是在WGS84坐标系统下进行的。 首先更改项目的坐标系统。在菜单“项目”->“坐标系统”或在工具栏“坐标系统”,则弹出“坐标

系统”对话框,选择WGS-84坐标。 图3-2 坐标系统 这里注意的是,在“投影”下见图,中央子午线是114°。很多情况下这里需要进行修改。 图3-3 WGS84投影 软件中自带的“中国-WGS 84”是允许修改的,我们换种方法:就是新建一个坐标文件,其他参数都和“中国-WGS84”一致,仅仅将中央子午线修改下。 在上图中,点击“新建”,得到“COORD GM”对话框,在“文件”->“新建”,如图

图3-4 新建坐标系统 然后在“设置”->“地图投影”,直接修改中央子午线,这里以81°为例,点击确定后,返回“COORD GM”对话框。 图3-5 投影设置 将输入源坐标和输入目标坐标的椭球,均改为WGS84。在“文件”->“保存”,输入名称和国家(中国),退出操作。

CORS坐标转换软件使用说明

坐标转换软件使用说明  1、功能介绍  在南京进行测量的同行一直受到坐标系统和已知控制点的困扰, 所以往往许多测量成果因坐标系统问题得不到承认,浪费了大量的人 力物力。基于此:本公司集全部精干技术力量,研发本款坐标转换软 件,可以说:它是全体测量工作者的福音。  南京CORS因为其免费,应用十分广泛,但是使用南京CORS在 很多情况下,因为已知控制点原因无法实地取得平面坐标而限制了 CORS优势的发挥。本软件可以实现基于南京CORS测量的WGS84 坐标与92南京地方坐标双向自由转换,转换精度与权威部门转换成 果比较(在南京市6800平方公里范围内,包括高淳、溧水、六合、 浦口):平面残差中误差优于±5mm、高程残差中误差均优于±1cm。精度完全具有保障,免去到处寻找控制点带来的人力、财力和时间浪费。按照最新城市规范规定,这种模式可以实现城市E级GPS控制 点的平面测量。  本软件是一款后处理软件,即:内业处理软件,它不能在实地计 算坐标,通过事后(采集)或事前(放样)数据处理,同样可以让你 在野外无忧无障碍开展工作。  适用平台:Windows 32位所有系统平台。  2、外业采集数据转换操作介绍  外业测量数据从RTK手簿中以WGS84坐标格式导出,导出以后 将文件复制到计算机,假设文件名为0513.dat。在电脑中启动软件,

界面如下:  图一:程序启动界面  首先选择转换方向下拉列表框,此时选择“WGS84—>NJ92”,表示将WGS84坐标转向92南京地方坐标,此时软件会出现一个按钮 键读入数据并转换,点击该按钮,在弹出的文件对话框中选择从手簿 导出的外业坐标文件。如:0513.dat,点击打开按钮即可完成转换。如图二:  图二:选择原始数据文件  记得一定要选择你的原始数据文件格式在点击打开按钮。转换完 成以后又会在对话框中再出现一个按钮导出转换成果,点击它即可将

常用坐标系之间的关系与转换

7.5 常用坐标系之间的关系与转换 一、大地坐标系和空间大地直角坐标系及其关系 大地坐标系用大地纬度企丈地经度L 和丈地髙H 来表示点的位置°这种坐标系是经 典大地 测量甬:両用座标紊7屜据地图投影的理论,大地坐标系可以通过一定的投影转 化为投影平面上的直角坐标系,为地形测图和工程测量提供控制基础。同时,这种坐标系 还是研究地球形状和大小的 种有用坐标系°所以大地坐标系在大地测量中始终有着重要 的作用. 空间大地直角坐标系是-种以地球质心为原点购亘墮?坐标系,一般用X 、化Z 表 示点 BSSTSTT 逐碇SS 範菇飞両H 绕禎扭转冻其轨道平面随时通过 地球质心。对它们的跟踪观测也以地球质心为坐标原点,所以空间大地直角坐标系是卫星 大地测量中一种常用的基本坐标系。现今,利用卫星大地测量的手段*可以迅速地测定点 的空间大地直角坐拯,广泛应用于导航定位等空间技术。同时经过数学变换,还可求岀点 的大地坐标I 用以加强和扩展地面大地网,进行岛屿和洲际联测,使传统的大地测量方法 发生了深刻的变化,所以空间大地宜角坐标系对现今大地测量的发展’具有重要的意义。 、大地坐标系和空间大地直角坐标系的转换 如图7- 23所示’尸点的位置用空间 大地 直角坐标〔X, Y, Z)表示,其相应 的大地坐 标为(E, L)a 将该图与图?一5 上式表明了 2种基本坐标系之间的关系。 加以比较可见,图7-5中的子午椭圆平面 相 当于图7-23中的OJVP 平面.其中 PPz=Z.相 当于图7-5中的j7;OP 3相当 丫于图7-5中的 仏两平面的经度乙可视为 相同,等于"叽 于是可以直接写岀 X=jrcQsi f Y=jrsinL, Z=y 将式(7-21).式(7-20)分别代入上式, 井考虑 式(7-26)得 X=Ncos^cosZr ” Y =NcQsBsinL > (7—78) Z=N (1—护〉sin^ ; BB 7-23

推导坐标旋转公式

推导坐标旋转公式 数学知识2010-09-12 21:03:53 阅读151 评论0 字号:大中小订阅 在《Flash actionScript 3.0 动画教程》一书中有一个旋转公式: x1=cos(angle)*x-sin(angle)*y; y1=cos(angle)*y+sin(angle)*x; 其中x,y表示物体相对于旋转点旋转angle的角度之前的坐标,x1,y1表示物体旋转angle 后相对于旋转点的坐标 从数学上来说,此公式可以用来计算某个点绕另外一点旋转一定角度后的坐标,例如:A(x,y)绕B(a,b)旋转β度后的位置为C(c,d),则x,y,a,b,β,c,d有如下关系式: 1。设A点旋转前的角度为δ,则旋转(逆时针)到C点后角度为δ+β 2。求A,B两点的距离:dist1=|AB|=y/sin(δ)=x/cos(δ) 3。求C,B两点的距离:dist2=|CB|=d/sin(δ+β)=c/cos(δ+β) 4。显然dist1=dist2,设dist1=r所以: r=x/cos(δ)=y/sin(δ)=d/sin(δ+β)=c/cos(δ+β) 5。由三角函数两角和差公式知: sin(δ+β)=sin(δ)cos(β)+cos(δ)sin(β) cos(δ+β)=cos(δ)cos(β)-sin(δ)sin(β) 所以得出:

c=r*cos(δ+β)=r*cos(δ)cos(β)-r*sin(δ)sin(β)=xcos(β)-ysin(β) d=r*sin(δ+β)=r*sin(δ)cos(β)+r*cos(δ)sin(β)=ycos(β)+xsin(β) 即旋转后的坐标c,d只与旋转前的坐标x,y及旋转的角度β有关 从图中可以很容易理解出A点旋转后的C点总是在圆周上运动,圆周的半径为|AB|,利用这点就可以使物体绕圆周运动,即旋转物体。 上面公式是相对于B点坐标来的,也就是假如B点位(0,0)可以这么做。现在给出可以适合任意情况的公式: x0 = dx * cos(a) - dy * sin(a) y0 = dy * cos(a) + dx * sin(a) 参数解释: x0,y0是旋转后相对于中心点的坐标,也就是原点的坐标,但不是之前点旋转后的实际坐标,还要计算一步,a旋转角度,可以是顺时针或者逆时针。 dx是旋转前的x坐标-旋转后的x坐标 dy是旋转前的y坐标-旋转后的y坐标 x1=b+x0; y1=c+y0; 上面才是旋转后的实际坐标,其中b,c是原点坐标 下面是上面图的公式解答: x0=(x-b)*cos(a)-(y-c)*sin(a); y0=(y-c)*cos(a)+(x-b)*sin(a); x1=x0+b; y1=y0+c;

相关主题
相关文档 最新文档