范文一:七参数坐标转换
实验一 七参数坐标的转换
七参ωx ωy ωz Tx Ty Tz m 数
-5.855977597879 -63.421779679192 9.188365321063 0.000004463236 -0.000002473494 -0.000005115077 -0.000011107611 值
54坐标/m 80坐标/m 转换后坐标差?/mm 备注 点名 控制点/转换点 ?X ?Y ?H X Y H X Y H
3516952.028 519078.188 36.393 AD01 3516904.452 519026.990 34.565 0.080 -0.052 0.152 GD01
3516556.925 519081.738 31.720 GD02 3516509.353 519030.538 29.893 0.440 -0.133 -0.835 GD08
3516569.509 518700.605 31.097 AD04 3516521.939 518649.409 29.272 0.248 0.162 -0.158 GD10
3516587.828 518636.550 30.246 AD07 3516540.258 518585.355 28.421 0.370 -0.037 0.092
3516304.881 518841.764 29.427 GD03 3516257.313 518790.565 27.602 0.000 0.000 -0.001
3516580.081 518699.370 31.056 GD07 3516532.511 518648.174 29.231 0.137 0.229 -0.178
3516360.719 518402.782 34.964 GD09 3516313.153 518351.588 33.141 0.100 0.419 0.645
精度评定:a = 3.475389486229349e-004
Matlab程序编码如下:
clear
format long
B=[1,0,0,0,-36.307,519081.665,3516958.654;
0,1,0,36.307,0,-3516958.654,519081.665;
0,0,1,-519081.665,3516958.654,0,36.307;
1,0,0,0,-32.894,518582.319,3516256.471;
0,1,0,32.894,0,-3516256.471,518582.319;
0,0,1,-518582.319,3516256.471,0,32.894;
1,0,0,0,-36.821,518397.422,3516573.822;
0,1,0,36.821,0,-3516573.822,518397.422;
0,0,1,-518397.422,3516573.822,0,36.821] l1=[3516958.654;
519081.655;
36.307;
3516256.471;
518582.319;
32.894;
3516573.822;
518397.422;
36.821];
l2=[3516911.078;
519030.457;
34.479;
3516208.905;
518531.123;
31.070;
3516526.254;
518346.230;
34.997];
l=l2-l1;
NBB=B'*B;
W=B'*l;
x=inv(NBB)*W
V=B*x-l;
r=2;
a=sqrt(V'*V/r)
g=[3516952.028,3516556.925,3516569.509,3516587.828,3516304.881,3516580.081,3516360.719;
519078.188,519081.738,518700.605,518636.550,518841.764,518699.370,518402.782;
36.393,31.720,31.097,30.246,29.427,31.056,34.964];
for i=1:1:7
c=[x(1,1);x(2,1);x(3,1)];
d=[1,x(6,1),-x(5,1);
-x(6,1),1,x(4,1);
x(5,1),-x(4,1),1];
e=[g(1:3,i)];
f=c+(1+x(7,1))*d*e
end
范文二:坐标转换参数
rtkTransformMethod,0
rtkTargetEllipsoid,0
longOfOrigin,117.0
latOfOrigin,0.0
projectionScale,1.0
falseEasting,500000.0
falseNorthing,0.0
projectionHeight,0.0
sevenParametersDx,0.0
sevenParametersDy,0.0
sevenParametersDz,0.0
sevenParametersRx,0.0
sevenParametersRy,0.0
sevenParametersRz,0.0
sevenParametersPpm,0.0
fourParametersDx,0.0
fourParametersDy,0.0
fourParametersRotate,0.0
fourParametersK,1.0
elevationFittingA0,0.0
elevationFittingA1,0.0
elevationFittingA2,0.0
elevationFittingA3,0.0
elevationFittingA4,0.0
elevationFittingA5,0.0
elevationFittingAverageLatitude,45.000000000000
elevationFittingAverageLongitude,90.000000000000
rtkCorrectX,0.0
rtkCorrectY,0.0
rtkCorrectZ,0.0
phoneGpsCorrectDx,0.0
phoneGpsCorrectDy,0.0
phoneGpsCorrectDz,0.0
范文三:各地坐标转换参数
各地坐标转换参数
GPS 参数
DX=11.9
DY=-120.8
DZ=-62.4
DA=-108.0
DF=0.00000050
E=93°00.000 +1.0000000 +5000000.0 0.0
藏东可用99°,其它参数不变,可对照地形图校对。
广东省GPS 参数:这是WGS84转北京54的,适宜河源、惠州、深圳、东莞地区
DX=-19
DY=-112
DZ=-55
DA=-108.0
DF=0.00000050
E=114°00.000 +1.0000000 +5000000.0 0.0
,WGS84转西安80的是
DX=-96
DY=-51
DZ=12
DA=-3
DF=0.00000000
E=114°00.000 +1.0000000 +5000000.0 0.0
适宜整个广东。
广东? 河源GPS 参数转换参数/
DX=12
DY=-121
DZ=-62
DA=-108
DF=0.00000050
E=114°00.000 +1.0000000 +5000000.0 0.0
坐标参数
海南坐标转换参数:
dx=-9.8
dy=-114.6
dz=-62.7
da=-108.0
df=0.0000005
中央子午线:111
DX = -18
DY = -104.5
DZ = -57.5
DA= -108;
DF= 0.0000005
中央子午经度:117或123(东为123,西为117)
新疆乌鲁木齐地区坐标转换参数:
DX = 19
DY = -33
DZ = 5
DA= -108;
DF= 0.0000005
中央子午经度:87
各地WGS84坐标系转换BJ54坐标系参数(不断加入中... )
注:以下参数仅供参考!!
以下为四川盆地坐标系转换参数
Dx=-4
Dy=-104
Dz=-45
Da=-108
Df=+0.0000005
中央子午经度:105
以下为包头地区坐标系转换参数
Dx=-92
Dy=-49
Dz=-4
Da=-108
Df=+0.0000005
中央子午经度:114
安徽省坐标转换区域化参数:
DY = -120
DZ = -48
DA= -108;
DF= 0.0000005
中央子午经度:117
鄂尔多斯市省坐标转换区域化参数:
DX = 16
DY = -147
DZ = -77
DA= -108;
DF= 0.0000005
中央子午经度:111
新疆阿克苏地区坐标转换参数:
DX = 18
DY = -152
DA= -108;
DF= 0.0000005
中央子午经度:81
西藏坐标转换区域化参数:
DX = 11.9
DY = -120.8
DZ = -62.4
DA= -108;
DF= 0.0000005
中央子午经度:93
赤峰地区坐标转换参数:
DX = -18
DY = -104.5
DZ = -57.5
DA= -108;
DF= 0.0000005
中央子午经度:117或123(东为123,西为117)
注:以下参数仅供参考!!
以下为四川盆地坐标系转换参数
Dx=-4
Dy=-104
Dz=-45
Da=-108
Df=+0.0000005
中央子午经度:105
安徽省坐标转换区域化参数:
DX = -15
DY = -120
DZ = -48
DA= -108;
DF= 0.0000005
中央子午经度:117
鄂尔多斯市省坐标转换区域化参数: DX = 16
DY = -147
DZ = -77
DA= -108;
DF= 0.0000005
中央子午经度:111
新疆阿克苏地区坐标转换参数: DX = 18
DY = -152
DZ = -76
DA= -108;
DF= 0.0000005
中央子午经度:81
西藏坐标转换区域化参数:
DX = 11.9
DY = -120.8
DZ = -62.4
DA= -108;
DF= 0.0000005
中央子午经度:93
坐标系转换问题
对于坐标系的转换,给很多GPS 的使用者造成一些迷惑,尤其是对于刚刚接触的人,搞不明白到底是怎么一回事。我对坐标系的转换问题,也是一知半解,对于没学过测量专业的人来说,各种参数的搞来搞去实在让人迷糊。在我有限的理解范围内,我想在这里简单介绍一下,主要是抛砖引玉,希望能引出更多的高手来指点迷津。
我们常见的坐标转换问题,多数为WGS84转换成北京54或西安80坐标系。其中WGS84坐标系属于大地坐标,就是我们常说的经纬度坐标,而北京54或者西安80属于平面直角坐标。对于什么是大地坐标,什么是平面直角坐标,以及他们如何建立,我们可以另外讨论。这里不多罗嗦。 那么,为什么要做这样的坐标转换呢?
因为GPS 卫星星历是以WGS84坐标系为根据而建立的,我国目前应
用的地形图却属于1954年北京坐标系或1980年国家大地坐标系;因为不同坐标系之间存在着平移和旋转关系(WGS84坐标系与我国应用的坐标系之间的误差约为80),所以在我国应用GPS 进行绝对定位必须进行坐标转换,转换后的绝对定位精度可由80提高到5-10米。简单的来说,就一句话,减小误差,提高精度。
下面要说到的,才是我们要讨论的根本问题:如何在WGS84坐标系和北京54坐标系之间进行转换。
说到坐标系转换,还要罗嗦两句,就是上面提到过的椭球模型。我们都知道,地球是一个近似的椭球体。因此为了研究方便,科学家们根据各自的理论建立了不同的椭球模型来模拟地球的形状。而且我们刚才讨论了半天的各种坐标系也是
建立在这些椭球基准之上的。比如北京54坐标系采用的就是克拉索夫斯基椭球模型。而对应于 WGS84坐标系有一个WGS84椭球,其常数采用 IUGG 第17届大会大地测量常数的推荐值。WGS84椭球两个最常用的几何常数:长半轴:6378137±2(m);扁率:1:298.257223563
之所以说到半长轴和扁率倒数是因为要在不同的坐标系之间转换,就需要转换不同的椭球基准。这就需要两个很重要的转换参数dA 、dF 。 dA 的含义是两个椭球基准之间半长轴的差;dF 的含义是两个椭球基
准之间扁率倒数的差。
在进行坐标转换时,这两个转换参数是固定的,这里,我们给出在进行84—〉54,84—〉80坐标转换时候的这两个参数如下:
WGS84>北京54:DA:-108;DF:0.0000005
WGS84>西安80:DA: -3 ;DF: 0
椭球的基准转换过来了,那么由于建立椭球的原点还是不一致的,还需要在dXdYdZ 这三个空间平移参量,来将两个不同的椭球原点重合,这样一来才能使两个坐标系的椭球完全转换过来。而由于各地的地理位置不同,所以在各个地方的这三个坐标轴平移参量也是不同的,因此需要用当地的已知点来计算这三个参数。具体的计算方法是:
第一步:搜集应用区域内GPS “B ”级网三个以上网点WGS84坐标系
B 、L 、H 值及我国坐标系(BJ54或西安80)B 、L 、h 、x 值。(注:B 、L 、H 分别为大地坐标系中的大地纬度、大地经度及大地高,h 、x 分别为大地坐标系中的高程及高程异常。各参数可以通过各省级测绘局或测绘院具有“A ”级、“B ”级网的单位获得。)
第二步:计算不同坐标系三维直角坐标值。计算公式如下: X=(N+H)cosBcosL
Y=(N+H)cosBsinL
Z=[N(1-e2)+H]sinB
不同坐标系对应椭球的有关常数详见下表:
项目 WGS84坐标系 BJ54坐标系 西安80坐标系
A
e2 0.00669437999013 0.006693427 0.006694385
(注:X 、Y 、Z 为大地坐标系中的三维直角坐标;A 为大地坐标系对应椭球之长半轴;e2为大地坐标系对应椭球第一偏心率; N 为该点的卯酉圈曲率半径,N=A/(1-e2sin2B)1/2;H=h+x,该处H 为BJ54或西安80坐标系中的大地高)
第三步:求出DX ,DY ,DZ 。即利用WGS84坐标系的X 、Y 、Z 值,减去我国坐标系的对应值,得出实现坐标系统转换的三个参数。(应算出WGS84与北京和西安坐标系两套参数。)
第四步:参数验证。参数计算之后必须对其进行验证。验证的方法是在应用区域内选择5个以上水准点进行实测,实测值与测绘部门提供的理论值对比,如果最大误差不大于15米,平均误差不大于10米,则计算出的参数可以使用,否则要重新计算或查找出现问题的原因。
对了,
还有一个很重要的事情,要在位置格式的地方,选择用户自定义方式,输入如下参数:
中央经线:视当地经度确定;
投影比例:1;
东西偏差:500000;
南北偏差:0
在这里面中央经线的确定很重要,根据实际所在地不同而有所和差异。至于这些参数的实际意义么,还是另外开题写吧,那也是不少的内容呢。
范文四:七参数坐标转换
实验一 七参数坐标的转换
Matlab 程序编码如下: clear
format long
B=[1,0,0,0,-36.307,519081.665,3516958.654; 0,1,0,36.307,0,-3516958.654,519081.665; 0,0,1,-519081.665,3516958.654,0,36.307; 1,0,0,0,-32.894,518582.319,3516256.471; 0,1,0,32.894,0,-3516256.471,518582.319; 0,0,1,-518582.319,3516256.471,0,32.894; 1,0,0,0,-36.821,518397.422,3516573.822; 0,1,0,36.821,0,-3516573.822,518397.422; 0,0,1,-518397.422,3516573.822,0,36.821] l1=[3516958.654; 519081.655; 36.307;
3516256.471; 518582.319; 32.894;
3516573.822; 518397.422; 36.821]; l2=[3516911.078; 519030.457; 34.479;
3516208.905; 518531.123; 31.070;
3516526.254; 518346.230; 34.997];
l=l2-l1; NBB=B'*B; W=B'*l;
x=inv(NBB)*W V=B*x-l; r=2;
a=sqrt(V'*V/r)
g=[3516952.028,3516556.925,3516569.509,3516587.828,3516304.881,3516580.081,3516360.719; 519078.188,519081.738,518700.605,518636.550,518841.764,518699.370,518402.782; 36.393,31.720,31.097,30.246,29.427,31.056,34.964]; for i=1:1:7
c=[x(1,1);x(2,1);x(3,1)]; d=[1,x(6,1),-x(5,1); -x(6,1),1,x(4,1); x(5,1),-x(4,1),1]; e=[g(1:3,i)];
f=c+(1+x(7,1))*d*e end
范文五:GPS坐标系统与参数转换
GPS坐标系统与参数转换
GPS坐标值的输出格式有大地坐标系的B、L、H值和通过“用户自定义”方式设置有关参数后输出的高斯投影的平面直角坐标系的公里网格坐标(x、y)值。我们要转换到BJ-54或C-80坐标系统,可以利用其“用户自定义”方式设置DX、DY、DZ、Da、Df五个参数即可实现。
DX、DY、DZ的计算:在手持式GPS接收机应用的区域内(该区域不宜过大),从当地测绘部门收集3,5个已知点的BJ-54或C-80坐标系统的坐标值,然后在对应的点位上读取WGS-84坐标系的坐标值;将收集到的坐标值根据不同的坐标系按相应转换公式转换为空间坐标系的坐标值,
X54、Y54、Z54)或C-80 (X84、Y84、用WGS-84的X84、Y84、Z84分别减去相对应的BJ-54(
Z84)的坐标值,即可得出DX、DY、DZ的值;实地对参数检验,实测值与已知值比较,在机器标称精度范围内即可,否则查因重算。
投影分带主要有两种:6度带(每隔经差6度分为一带)和3度带(每隔经差3度分为一带)。每个带在平面上的投影如右下图所示,每个带设立一个平面直角坐标系,以轴子午线L为纵坐标0轴,以子午线与赤道的交点作为坐标系原点,赤道作为横坐标轴。这样,在轴子午线以东的点,其横坐标y是正值,轴子午线以西的点,其横坐标y是负值。为了避免在使用成果时可能搞错符号,把纵坐标轴向西平移500公里,这样得到的横坐标y+500的值就永远为正值。
山西的经度介于111?,115?之间,多数处于6?带的第19带内,只有东北部的天镇、灵丘、广灵处于20带内。其带数可从地形图公里网的横坐标上直接查找得到。
(中央)子午线经度的计算:
六度带横坐标的表示方法:带号+横坐标值
如:地形图上的横坐标为19365,其中19即为带号;365为横坐标值。
六度带中央子午线经度的计算:中央子午线经度= 带号×6-3
如:地形图上的横坐标为19365,其所处的六度带的中央子午线经度为:19×6?-3?=111?,同样计算20带的中央子午线为117?。
三度带中央子午线经度的计算:中央子午线经度= 带号×3
投影带带号的计算
六度投影带带号的计算:带号=INT(INT子午线经度?6)+1
三度投影带带号的计算:带号=INT,(INT子午线经度-1.5) ?3 , +1 如:某点处于110 ?的子午线经度上,其:
六度带带号=INT(110?6)+1=19
三度带带号=INT, (110-1.5) ?3 , +1 =37
我省采用统一的参数值:
DX:-5.4
DY:-113.5
DZ:-40.6
DA:-108.0
DF:0.0000005
位置格式:User Grid
距离/速度:公制
高度/垂直速度:米
中央经线:据分带位置确定,如111.00.00E、114.00.00E SCALE 比例系数: 1.0000000
FALSEE 东西偏移:500000.00
FALSEN 南北偏移: 0.0
地图基准:User
手持GPS使用注意事项
?参数必须正确,并做必要的误差实测,过大不宜使用。
?输入坐标时一定要正确,不要颠倒横、纵坐标的输入位置:一般来讲,第一行为横坐标,
第二行为纵坐标;如地形图坐标为19365,3985,则输入为0365000,3985000。 ?定点时,必须保证信号稳定,并应尽量避免在有障碍物的地方读数,特别是测量的距离较
短时,更应注意。必须在有电磁干扰的地方定点时,应采取差位定点的办法,即把定点位置平移
几米,并在草图中注明,在正式绘图时复位。每个定位点都应做好记录。 ?工作期间不要清除航迹,以便在数据万一丢失时,通过它来重新取得数据。 ?测面积时,为保证精度,应保持匀速前进,不要过快。
?注意节电,事先必须充好或更换,以免半途而废。