坐标反算

坐标反算
坐标反算

在现场工作中,以往我们都是已知某点的里程及边距,来计算出该点的坐标,但有时我们如果能在测得某点坐标后,计算出该点的里程和距线路中心的距离(在这里我姑且称之为坐标反算)的话,将会帮我们大大减轻野外工作量,提高我们的工作效率。例如:路基填了几层后要精确检查一下路基是否够宽,那么按照我们以往的做法,就是要先将线路中心线放出来,然后用尺拉一下路基宽度,与其在此高程的设计宽度作比较,这样做对高填方而言极不方便。或者是先按所测高程,计算的宽度放出路基边桩,再与所填边线作比较。以上两种方法现场工作量都比较大。较为简便的方法是,我们可以测一下已填路基边线上任一点的三维坐标,然后将其反算求出该点的里程,及其距中线的距离(即所填宽度),由计算出的里程,可算出该里程的路面设计高程,再有所测高程,可计算出该点的设计宽度,两宽度作比较即可。同样在桥面铺装施工时,我们也无须再像以往那样,先放出某点再测其高程,然后与设计高程比较计,算出该点铺装厚度,而可以沿桥面外边线随意布点,测其三维坐标,计算出其里程及到中心线的距离,便可由其里程及距中心距离,计算出该点的设计高程,与其测得高程作比较得出应铺厚度。这样便大大减轻了外业工作强度(由放出点后再测其高程,变为测任意点高程),而内业计算量与常规相当。另外在临时增加桥涵时,也常用到此方法来计算变更桥涵的中心里程(斜交或正交均可).如目前我标段就存在很多临时变更涵洞,按以往我们的方法是先估计该处大概里程,然后放出所估计里程的中心桩,再用皮尺量出所要增加涵洞处与该中心桩的距离,以此来推算出涵洞的中心里程,这一过程即繁琐又不准确。而目前我们采用的方法是用全站仪测得跨路基现有水沟两端的沟底坐标,计算出其与路基的夹角,按所测坐标及此夹角就可以准确、快速地反算出水沟中心所对应的线路中心里程了。我们在日常测量工作中的很多方面,也会用到这一方法来减轻野外工作量。在目前我标段的S334分离式立交桥的架设过程中,也同样用到了此方法.支座安装好后,对支座中心位置检及高程查无误后开始梁板架设,但是尽管测量控制放样符合规范要求,可是因为其它方面的各种原因可能会使梁板出现偏位高程也可能会出现偏差,那么对现在这种问题该如何检查呢?其实方法是一样的,首先我们可以用全站仪测得架设好后梁的边板外边缘任一点的三维坐标,由此坐标反算出该点所对应的中心里程和距中心的距离,就可以和设计图纸上的距离作比较来检查其是否存在偏位,该点的设计高程也可以由反算所得的中心里程和距中心的边距算出,与所测得的实际高程作一下比较也就可以了.那么通过以上讨论问题归结到了一点,那就是如何在测得任一点坐标后,计算出其所对应的线路中心里程,及其到线路中心的距离(或是斜交的长度)呢?解决此类问题,对目前一些测量软件来说早已不成问题,但是在现场工作中我们用的更多、更方便的还是计算器,那么能否用我们常用的4800或5800计算器编程,来计算此类问

题呢?对此我做了一下尝试,取得了不错的效果,现作一简要介绍:(此法需分线元计算,其计算原理如下:)

一:直线

如图:设OP为线路中心线,C为中心线外任一点,已知起点O的坐标、OP方位角αf、和角A,测得C点坐标后,αJ 和LOC便可计算出来.

由三角关系可得:∠K =αJ-αf∠C=180-∠K-(180-∠A)

LHC=sin∠K*LOC/sin(180-∠A)

LOH= sin∠C*LOC/sin(180-∠A)

便可求得A点所对应的线路H点里程(O点里程加上LOH),及HC的长度(正交时为其边距)。

二:圆曲线

如图:圆曲线可借助于圆心O、与圆曲线起点A、测点C以及切点B、构成的三角形,利用三角关系,来求得测点里程及其到中心线的长度或距离(正交时)。在正交时,如OD可直接由测点D的坐标,及圆心O的坐标求出OD的长度,再减去半径,即可得其到线路中心距离ED。测点里程可由圆曲线起点里程(即A点里程),加上角a所对的圆弧的长度(即弧AE的长度)求得。圆曲线计算原理与直线相似,但计算过程较为复杂(程序并不繁琐),在此不再赘述。

三:缓和曲线

1:此法原理是:首先确定一与待求点里程相近的一里程,输入“缓和曲线”正算程序,求出该里程的中心坐标及其切线方位角。

2:将所得坐标及方位角,作为起点坐标及方位角,输入到上述第一条中的“直线反算”程序中,以求出待求点,在该渐近直线中的里程及支距。

3:再将所得新里程输入到“缓和曲线”程序中去,以求出新的坐标及方位角。

4:再将所得坐标及方位角,作为“直线反算”中的起点坐标输入到“直线反算”程序中去。

如此反复渐近,直到所得切线,与待求点与其中心线点连线夹角,约等于已知的夹角为止(在程序中可设一判断值,实践得<0.12即可),所得里程及距离,即为要测点的里程和到中线的距离。(或是斜交的长度)。

附:4800计算程序

一:直线Lbl1:{XY A}:X“(CDX)X”:Y“(CDY)Y”:Q“(QDLC)Q”:N“(QDX)X1”:E“(QDY)Y1”:

F“(QDFWJ)F”:A“(JIAJAO)A”:Pol(X-N,Y-E):Fixm:J<0=>J=J+360:GoTo2: ≠>J=J:

Lbl2:G=J-F:B=√((Y-E)2+(X-N)2):

L“(BZJL)L”=B*sinG/sin(180-AbsA)▲AbsA>AbsG=>U=1:GoTo3:≠>=U=-1:

Lbl3:C“(CDLC)C”=Q+U*√(L2+B2-Asb2BLcos(AbsA-AbsG))▲GoTo1

二:圆曲线

L bl1:{XYJ}:X”(CDX)X”:Y”(CDY)Y””(YQDLC)Q”:F”(YQDFWJ)F”:N”(YQDX)X1”:E”(Y QDY)Y1”:G=N+R”(BANJING)R”*cos(F+90)=E+R*sin(F+90):B=√((X-G)2+(Y-D)2):

A=sin-1(AbsR*sin(90+J“(JIAJIAO)J”)/B):AbsJ=90=>L“(BZJL)L”=AbsR-B▲

Goto2: ≠>L“(BZJL)L”=AbsR*sin(A+AbsJ-90)/sinA▲

Lbl2:H=√((X-N)2+(Y-E)2):O=cos-1((R2+B2-H2)/2AbsRB):AbsJ≤90=>C“(CDLC)C”=Q+∏Abs R*(O-Abs(AbsJ+A-90))/180▲

Goto1:≠>C“(CDLC)C”=Q+∏AbsR*(O+Abs(AbsJ+A-90))/180▲Goto1

3:缓和曲线

Lbl1:{WENO}:W”(HCDLC)W””(HQDLC)Q”:F”(HQDQXFWJ)F””(JHHQXC)D”:R”(Q

DYBJ)R””(ZDYBJ)P”:R>9998=>C=PD=0:Goto3:≠>>9998=>C=RD

=D:Goto3:≠>C=Abs(DRP/(R-D)):L=C/R:Lbl3:R>=>I=L+Abs(W-Q):W>Q=>U=1:Goto4:

≠>U=-1:Goto4?≠>R<=>L=L-Abs(W-Q):W>Q=>U=-1:Goto4: ≠>U=1

bl4:H“(HQDX)X”:K“(HQDY)Y”:T=F-UL2/2C*180/∏:V=UI2/2C*180/∏:A=(I-L)-(I5-L5)/40C2 +(I9-L9)/3456C4:

B=(I3-L3)/6C-(I7-L7)/336C3:X=H+UAcosT-BsinT:Y=K+UAsinT+BcosT:N”(CDX)X”:E”(CDY) Y”:O”(JIAJAO)J”ol(N-X,E-Y):Fixm:J<0=>J=J+360:Goto5:≠>J=J:Lbl5:G=J-T-V:

B=√((E-Y)2+(N-X)2):

M=BsinG/sin(180-AbsO):AbsO>AbsG=>U=1:Goto6:≠>U=-1:

Lbl6:W=W+U*√(M2+B2-Abs2Bmcos(AbsO-AbsG)):

Abs(AbsO-AbsG)<0.12=>M“(BZJL)L” ▲W“(CDLC)C” ▲

Goto1:≠>Goto3

坐标正算反算公式讲解(借鉴材料)

一 方位角: 在高斯直角坐标系中,由坐标纵轴方向的北端起,顺时针量到直线间的夹角,称为该直线的坐标方位角,常简称方位角,用a 表示。 1、第一象限的方位角 Y X 第一象限第二象限 第三象限 第四象限 o A a 图1 2、第二象限的方位角 Y X 第一象限 第二象限第三象限 第四象限 o A a 图2

3、第三象限的方位角 Y X 第一象限 第二象限 第三象限 第四象限 o A a 图3 4、第四象限的方位角 Y X 第一象限 第二象限 第三象限 第四象限 o A a 图4 方位角计算公式:

x =a -1 tan A Y O Y -A X O X - 方位角的计算器计算程序:Pol(X A -X O ,Y A -Y O ) 直线OA 方位角度值赋予给计算器的字母J ,0≤J <360。 直线段OA 的距离值赋予给计算器的字母I,I >0 直线OA 与直线AO 的方位角关系: 1、 当直线OA 的方位角≤180°时,其反方位角等于a+180°。 2、 当直线 OA 的方位角>180°时,其反方位角等于a-180°。 二 方位角的推算 (一)几个基本公式 1、坐标方位角的推算

或: 注意:若计算出的方位角>360°,则减去360°;若为负值,则加上360°。 例题:方位角的推算 已知:α12=30°,各观测角β如图,求各边坐标方位角α23、α34、α 45 、α51。 13 图5

解:α23= α12-β2+180°=30°-130°+180°=80° α34= α23-β3+180°=80°-65°+180°=195° α45=α34-β4+180°=195°-128°+180°=247° α51=α45-β5+180°=247°-122°+180°=305° α12=α51-β1+180°=305°-95°+180°=30°(检查) 三坐标正算 一、直线段的坐标计算 o B D A C E a a p 图6 设起点O的坐标(X O,Y O),直线OP的方位角为F op,求A、C、E点的坐标 1、设直线段OA长度为L,则A点坐标为 X A=X O+L×Cos(F op)

高斯投影坐标正反算VB程序

高斯投影坐标正反算 学院: 班级: 学号: 姓名: 课程名称: 指导老师:

实验目的: 1.了解高斯投影坐标正反算的基本思想; 2.学会编写高斯正反算程序,加深了解。 实验原理: 高斯投影正算公式中应满足的三个条件: 1. 中央子午线投影后为直线; 2. 中央子午线投影后长度不变; 3. 投影具有正形性质,即正形投影条件。 高斯投影反算公式中应满足的三个条件: 1. x坐标轴投影成中央子午线,是投影的对称轴; 2. x轴上的长度投影保持不变; 3. 正形投影条件,即高斯面上的角度投影到椭球面上后角度没 有变形,仍然相等。 操作工具: 计算机中的VB6.0 代码: Dim a As Double, b As Double, x As Double, y As Double, y_# Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#,

m4#, m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val(Text1.Text) min1 = Val(Text2.Text) sec1 = Val(Text3.Text) deg2 = Val(Text4.Text) min2 = Val(Text5.Text) sec2 = Val(Text6.Text) l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val(Text9.Text) k1 = ((l_ * 180 / 3.14159 + 3) / 6) k2 = (l_ * 180 / 3.14159 / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else If dh = 3 Then

高斯投影正反算公式 新

高斯投影坐标正反算 一、相关概念 大地坐标系由大地基准面和地图投影确定,由地图投影到特定椭圆柱面后在南北两极剪开展开而成,是对地球表面的逼近,各国或地区有各自的大地基准面,我国目前主要采用的基准面为:基准面,为GPS基准面,17届国际大地测量协会上推荐,椭圆柱长半轴a=6378137m,短半轴b=; 2.西安80坐标系,1975年国际大地测量协会上推荐,椭圆柱长半轴a=6378140m,短半轴b=; 3.北京54坐标系,参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=; 通常所说的高斯投影有三种,即投影后: a)角度不变(正角投影),投影后经线和纬线仍然垂直; b)长度不变; c)面积不变; 大地坐标一般采用高斯正角投影,即在地球球心放一点光源,地图投影到过与中央经线相切的椭圆柱面上而成;可分带投影,按中央经线经度值分带,有每6度一带或每3度一带两种(起始带中央经线经度为均为3度,即:6度带1带位置0-6度,3度带1带位置度),即所谓的高斯-克吕格投影。

图表11高斯投影和分带 地球某点经度(L)为过该点和地球自转轴的半圆与子午线所在半圆夹角,东半球为东经,西半球为西经;地球某点纬度(B)为所在水平面法线与赤道圆面的线面角。 正算是已知大地坐标(L,B),求解高斯平面坐标(X,Y),为确保Y值为正,Y增加500公里;反算则是由高斯平面坐标(X,Y)求解大地坐标(L,B)。 二、计算模型: 地球椭球面由椭圆绕地球自转轴旋转180度而成。 图表 1 椭圆 椭圆长半轴a,椭圆短半轴b, 椭圆方程:

(1) 图表2椭球面 椭球面方程: y2 a2+ x2 b2 + z2 a2 =1 /*************************************** 与网上充斥的将函数关系先展开为泰勒级数,再依据投影规则确定各参数不同,本文直接依据空间立体三角函数关系得出结果。 *****/ (一)正算 由图表1,

坐标正反算计算公式

坐标正反算公式

一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值,待定参数主要为GPS 网中点的坐标;同时,利用基线解算时随基线向量一同输出的基线向量的方差阵,形成平差的随机模型,最终形成平差完整的数学模型。随后对所形成的数学模型进行求解,根据平差结果来确定观测值中是否存在粗差,数学模型是否有需要改进的部分,若存在问题,则采用相应的方法进行处理并重新进行求解;若未发现问题,则输出最终结果,并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值,在平差过程中引入会使GPS 网的尺度、方向和位置发生变化的外部起算数据,从而实现GPS 网成果由基线解算时GPS 卫星星历所采用的参照系(WGS84 )到特定参照系的转换,得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。 二、南方GPS数据处理软件的平差方式

三维约束平差是指在基线解算后,WGS84坐标系下的三维平差,在三维平差中是不需要当地平面直角坐标系下的已知点坐标,当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标(如果只有经纬度时,可采用COORD4.1软件进行转换,本站免费提供)进行三维约束平差,即可得到与已知点相匹配的WGS84坐标。 一般情况下,在“已知点坐标录入”窗口中,我们都没有输入WGS8坐标,而只输入当地坐标系下的已知坐标,此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图:

如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。

坐标正反算定义及公式

坐标正反算定义及公式 Corporation standardization office #QS8QHH-HHGX8Q8-GNHHJ8

第六章→第三节→导线测量内业计算 导线计算的目的是要计算出导线点的坐标,计算导线测量的精度是否满足要求。首先要查实起算点的坐标、起始边的方位角,校核外业观测资料,确保外业资料的计算正确、合格无误。 一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离,计算待定点的坐标,称为坐标正算。如图6-6 所示,点的坐标可由下式计算: 式中、为两导线点坐标之差,称为坐标增量,即: 【例题6-1】已知点A坐标,=1000、=1000、方位角=35°17'36.5",两点水平距离=200.416,计算点的坐标?

35o17'36.5"=1163.580 35o17'36.5"=1115.793 2、坐标反算 已知两点的坐标,计算两点的水平距离与坐标方位角,称为坐标反算。可知,由下式计算水平距离与坐标方位角。 (6-3) (6-4) 式中反正切函数的值域是-90°~+90°,而坐标方位角为0°~360°,因此坐标方位角的值,可根据、的正负号所在象限,将反正切角值换算为坐标方位角。 【例题6-2】=3712232.528、=523620.436、=3712227.860、=523611.598,计算坐标方位角计算坐标方位角 、水平距离。

=62°09'29.4"+180°=242°09'29.4" 注意:一直线有两个方向,存在两个方位角,式中:、的计算是过A点坐标纵轴至直线的坐标方位角,若所求坐标方位角为,则应是A点坐标减点坐标。 坐标正算与反算,可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算,普通科学电子计算器的类型比较多,操作方法不相同,下面介绍一种方法。 【例题6-3】坐标反算,已知=2365.16、=1181.77、 =1771.03、=1719.24,试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键[=]等于纵坐标增量,按储存键[], 键入1719.24-1181.77按等号键[=]等于横坐标增量,按[]键输入,按[]显示横坐标增量,按[]键输入,按第二功能键[2ndF],再按[]键,屏显为距离,再按[]键,屏显为方位角。 【例题6-4】坐标正算,已知坐标方位角=294°42'51", =200.40,试计算纵坐标增量横坐标增量。

坐标反算正算计算公式

坐标反算正算计算公式 一、坐标正算 根据A点的坐标X A、Y A和直线AB的水平距离D AB与坐标方位角O AB,推算B点的坐标X B、Y B,为坐标正算,其计算公式为: X B = X A + AX AB Y B = X A + AY AB(1-18 ) 二式中,AX AB与AY AB分别称为A?B的纵、横坐标增量,其计算公式为: AX AB = X B—X A = D AB COS O AB AY AB = Y B—Y A = D AB sin O AB(1-19 ) 注意,AX AB和AY AB均有正、负,其符号取决于直线AB的坐标方位角所在的象限。 二、坐标反算 根据A、B两点的坐标X A、Y A和X B、Y B,推算直线AB的水平距离D AB与坐标方位角 OCAB , 为坐标反算。其计算公式为: (1-20 ) 注意,由(1-20 )式计算OCAB时往往得到的是象限角的数值,必须先根据AX AB、AY AB的正、负号,确定直线AB所在的象限,再将象限角换算为坐标方位角。 三角函数内容规律 三角函数看似很多,很复杂,但只要掌握了三角函数的本质及内部规律就会发现 三角函数各个公式之间有强大的联系。而掌握三角函数的内部规律及本质也是学好三 角函数的关键所在. 1、三角函数本质: 三角函数的本质来源于定义,如右图: 根据右图,有 sin 0 =y/ R; cos 0 =x/R; tan 0 =y/x; cot 0 =x/y。 深刻理解了这一点,下面所有的三角公式都可以从这里出发推导出来,比如以推导 si n( A+B) = si nAcosB+cosAs inB 为例: 推导: 首先画单位圆交X轴于C,D,在单位圆上有任意A,B点。角AOD为a,BO D为B,旋转AOB使0B与0D重合,形成新A'OD。 A(cos a ,sin a ),B(cos 3 ,sin 3 ),A'(cos( - BM,sin( 诩)) OA'=OA=OB=OD=1,D(1,0) [cos( a- 3 >1]A2+[sin( a- 3 )]A2=(cos a cos 3 )A2+(sin a-sin 3 )A2 和差化积及积化和差用还原法结合上面公式可推出(换(a+b)/2与(a-b)/2 ) [1] (1-21 )

5800-9860计算器坐标正反算通用程序

5800-9860计算器坐标正反算通用程序 1. 坐标正算主程序(命名为ZBZS) 第1行:Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B 第2行:Prog “A” 第3行:”X=”:N+Zcos(F+B)◢ 第4行:”Y=”:E+Zsin(F+B)◢ 第5行:”F=”:F◢ 第6行:Goto 0 K——计算点的里程 BIAN——计算点到中桩的距离(左负右正) α——取前右夹角为正 2. 坐标反算桩号和偏距主程序(命名为ZBFS) 第1行:”X1=”? C:”Y1=”?D: ”K1=”?K 第2行:Lbl 0:Prog “A” 第3行:Pol(C-N,D-E) 第4行:List Ans[1]→I 第5行:List Ans[2]→J 第6行:Icos(F-J)→S:K+S→K 第7行:Abs(S)>0.0001=>Goto 0 第8行:”K1=”:K◢ 第9行:”BIAN=”:Isin(J-F)→Z◢ X1——取样点的X坐标 Y1——取样点的Y坐标 K1——输入时为计算起始点(在线路内即可),输出时为反算点的桩号 Z——偏距(左负右正) 3. 计算坐标子程序(命名为XYF) 为了简洁,本程序由数据库直接调用,上述中的正反算主程序不直接调用此程序第1行:K-A→S:(Q-P)÷L→I 第2行:N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N 第3行:E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E 第4行:F+S(2P+SI)×90÷π→F 4. 数据库(命名为A)

第1行:K≤175.191=>Stop 第2行: 175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/240→Q:70.417→L: K≤A+L =>GoTo 1(第一缓和曲线,圆半径为240) 第3行:245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q: 72.915→L: K≤A+L =>Goto 1(第圆曲线,半径为240) 第4行:318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104 →L: K≤A+L =>Goto 1(第二缓和曲线,圆半径为240) 第5行:373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:- 1/180→Q:67.222→L: K≤A+L=>Goto 1:Stop(第一缓和曲线,圆半径为180) 第6行:Lbl 1:Prog “XYF” A——曲线段起点的里程 N——曲线段起点的x坐标 E——曲线段起点的y坐标 F——曲线段起点的坐标方位角 P——曲线段起点的曲率(左负右正) Q——曲线段终点的曲率(左负右正) L——曲线段长度(尽量使用长度,为计算断链方便) 说明: (1)在9860中,程序中所有公式和部分函数结果均存储在List Ans列表数组中,要想多次调用最好随公式取出结果,并赋给变量。 (2)正算主程序可以计算一般边桩的坐标,如要计算类似涵洞端墙的坐标需增加第二偏距和转角两个变量。 (3)程序规定,左偏曲线曲率(半径倒数)输入负值,右偏曲线曲率输入正值,直线上点曲率输入0,例如直线段,线元起点和终点均输入0,第一缓和曲线分别输入0和圆半径的倒数,圆曲线均输入半径倒数,第二缓和曲线分别输入圆半径倒数和0,卵形曲线分别输 入对应圆半径的倒数 (4)若是从大里程向小里程的反方向计算,则曲率取正方向时的负值,方位角减去(或加上)180度。 (5)有多个匝道的项目,可随时更改正反算主程序中的红色字体部分来调用其它线路的数据 (6)反算桩号偏差为1mm (7)可以计算任意线型的任意点坐标

计算坐标与坐标方位角基本公式

二 计算坐标与坐标方位角的基本公式 控制测量的主要目的是通过测量和计算求出控制点的坐标,控制点的坐标是根据边长及方位角计算出来的。下面介绍计算坐标与坐标方位角的基本公式,这些公式是矿山测量工中最基本最常用的公式。 一、坐标正算和坐标反算公式 1.坐标正算 根据已知点的坐标和已知点到待定点的坐标方位角、边长计算待定点的坐标,这种计算在测量中称为坐标正算。 如图5—5所示,已知A 点的坐标为A x 、A y ,A 到B 的边长和坐标方位角分别为AB S 和AB α,则待定点B 的坐标为 AB A B AB A B y y y x x x ?+=?+= } (5—1) 式中 AB x ? 、AB y ?——坐标增量。 由图5—5可知 AB AB AB AB AB AB S y S x ααsin cos =?=? } (5—2) 式中 AB S ——水平边长; AB α——坐标方位角。 将式(5-2)代入式(5-1),则有 AB AB A B AB AB A B S y y S x x ααsin cos +=+= }

(5—3) 当A 点的坐标A x 、A y 和边长AB S 及其坐标方位角AB α为已知时,就可以用上述公式计算出待定点B 的坐标。式(5—2)是计算坐标增量的基本公式,式(5—3)是计算坐标的基本公式,称为坐标正算公式。 从图5—5可以看出AB x ?是边长AB S 在x 轴上的投影长度, AB y ?是边长AB S 在 y 轴上的投影长度,边长是有向线段,是在 实地由A 量到B 得到的正值。而公式中的坐标方位角可以从0°到360°变化,根据三角函数定义,坐标方位角的正弦值和余弦值就有正负两种 情况,其正负符号取决于坐标方位角所在的象限,如图5—6所示。从式(5—2)知,由于三角函数值的正负决定了坐标增量的正负,其符号归纳成表5—3。

测绘技术之坐标反算与正算

5.3坐标反算 坐标反算,就是根据直线两个端点的已知坐标,计算直线的边长和坐标方位角的工作。如图5.3所示,若A、B为两已知点,其坐标分别为(XA,YA)和(XB,YB),根据三角函数,可以得出直线的边长和坐标方位角计算公式: tgα=△YAB/△XAB=(YB-YA)/(XB-XA) αAB =tg-1 (△YAB/△XAB)= tg-1 ((YB-YA)/(XB-XA)) /td> DAB=△YAB/sin αAB=XAB/cos αAB 或 (5.6) DAB=√(△X2+△Y2) 应当注意,按公式(5.5)用计算器计算时显示的反正切函数值在-90°~+90°之间,而坐标方位角范围是0°~360°,所以按(5.5)式反算方位角时,要根据ΔX、ΔY的正负符号确定直线AB 所在的象限,从而得出正确的坐标方位角。如使用fx140等类型的计算器,可使用功能转换键 INV 和极坐标与直角坐标换算键P→R以及x←→y键直接计算求得方位角。按键顺序为: ΔX INV R→P ΔY =显示D X←→y 显示α。 例5.2 已知B点坐标为(1536.86 ,837.54),A点坐标为(1429.55,772.73),求距离DBA和坐标方位角αBA。 解:先计算出坐标增量: ΔXBA=1429.55-1536.86=-107.31 ΔYBA=772.73-837.54=-64.81 直接用计算器计算: 按-107.31 INV P→R -64.81 =显示125.36(距离DBA); 按 x←→y 显示211°07′53″(坐标方位角αBA)。 5.2 坐标正算 坐标正算,就是根据直线的边长、坐标方位角和一个端点的坐标,计算直线另一个端点的坐标的工作。如图5.3所示,设直线AB的边长DAB和一个端点A的坐标XA、YA为已知,则直线另一个端点B 的坐标为: XB=XA+ΔXAB (5.1) YB=YA+ΔYAB (5.2) 式中,ΔXAB、ΔYAB称为坐标增量,也就是直线两端点A、B的坐标值之差。由图5.3中,根据三角函数,可写出坐标增量的计算公式为: 图5.3 ΔXAB=DAB·cosαAB (5.3)

坐标正反算

一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离,计算待定点的坐标,称为坐标正算。如图6-6 所示,点的坐标可由下式计算: 式中、为两导线点坐标之差,称为坐标增量,即: 【例题6-1】已知点A坐标,=1000、=1000、方位角 =35°17'36.5",两点水平距离=200.416,计算点的坐标? 35o17'36.5"=1163.580 35o17'36.5"=1115.793 2、坐标反算 已知两点的坐标,计算两点的水平距离与坐标方位角,称为坐标反算。如图6-6可知,由下式计算水平距离与坐标方位角。 (6-3) (6-4)式中反正切函数的值域是-90°~+90°,而坐标方位角为0°~360°,因此坐标方位角的值,可根据、的正负号所在象限,将反正切角值换算为坐标方位角。

【例题6-2】=3712232.528、=523620.436、 =3712227.860、=523611.598,计算坐标方位角计算坐标方位角 、水平距离。 =62°09'29.4"+180°=242°09'29.4" 注意:一直线有两个方向,存在两个方位角,式中:、 的计算是过A点坐标纵轴至直线的坐标方位角,若所求坐标方位角为,则应是A点坐标减点坐标。 坐标正算与反算,可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算,普通科学电子计算器的类型比较多,操作方法不相同,下面介绍一种方法。 【例题6-3】坐标反算,已知=2365.16、=1181.77、 =1771.03、=1719.24,试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键[=]等于纵坐标增量,按储存键[], 键入1719.24-1181.77按等号键[=]等于横坐标增量,按[]键输入,按[]显示横坐标增量,按[]键输入,按第二功能键[2ndF],

高斯投影正反算公式83

§8.3高斯投影坐标正反算公式 任何一种投影①坐标对应关系是最主要的;②如果是正形投影,除了满足正形投影的条件外(C-R 偏微分方程),还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式: B,l ? x,y 高斯投影必须满足以下三个条件: ①中央子午线投影后为直线;②中央子午线投影后长度不变;③投影具有正形性质,即正形投影条件。 由第一条件知中央子午线东西两侧的投影必然对称于中央子午线,即(8-10)式中,x 为l 的偶函数,y 为l 的奇函数;0330'≤l ,即20/1/≈''''ρl ,如展开为l 的级数,收敛。 +++=++++=553316644220l m l m l m y l m l m l m m x (8-33) 式中 ,,10m m 是待定系数,它们都是纬度B 的函数。 由第三个条件知: q y l x l y q x ??-=????=??, (8-33)式分别对l 和q 求偏导数并代入上式 ----=++++++=+++553315 63424 42204 52 3164253l dq dm l dq dm l dq dm l m l m l m l dq dm l dq dm dq dm l m l m m (8-34) 上两式两边相等,其必要充分条件是同次幂l 前的系数应相等,即

dq dm m dq dm m dq dm m 231 20 13121? =? -== (8-35) (8-35)是一种递推公式,只要确定了 0m 就可依次确定其余各系数。 由第二条件知:位于中央子午线上的点,投影后的纵坐标x 应等于投影前从赤道量至该点的子午线弧长X ,即(8-33)式第一式中,当0=l 时有: 0m X x == (8-36) 顾及(对于中央子午线) B V M r M B N dq dB M dB dX cos cos 2 ==== 得: B V c B N r dq dB dB dX dq dX dq dm m cos cos 01===?===(8-37,38) B B N dq dB dB dm dq dm m cos sin 2 2121112=?-=?-= (8-39) 依次求得6543,,,m m m m 并代入(8-33)式,得到高斯投影正算公式

卡西欧计算器坐标的正反算

可以算任意斜交涵洞轴线的坐标,增加T为斜交角度,规定T为涵轴右侧方向与“线路前进方向切线”之间的夹角,当涵轴与线路正交时,T=90,其他操作与原程序一样; 1. 正算子程序(SUB1) [color=Red]A=0.26:B=0.74:K=0.02:L=0.82:F=1-L: M=1-K:X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LW D))+Bcos(G+57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD))):Y=V+W(Asin(G+57.2958QKW(1/ P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW (1/P+MWD))):F=G+57.2958QW(1/P+ WD)+90:X=X+Zcos(F-90+T):Y=Y+Zsin(F-90+T) 2. 反算子程序(SUB2) W=Abs((Y-V)cos(G-90)-(X-U)sin(G-90)):Z=0:Lbl 0:Prog "SUB1":L=(G-90)+5 7.2958QW(1/P+ WD):Z=(J-Y)cosL-(I-X)sinL:AbsZ<1E-6=>Goto1:≠>W=W+Z:Goto 0Δ←┘ Lbl 1:Z=0:Prog "SUB1":Z=(J-Y)÷sinF 二.增设数据库程序(SJK主程序) Lb1 4:"1.SZ => XY":"2.XY => SZ":{NS}:S∠下一线元起点里程=>O =本线元起点里程:U=本线元起点X:V=本线元起点Y:G=本线元起算方位角:H =本线元长度:P=起点曲率半径:R=终点曲率半径:Q=0或1、-1:Prog“TYQXJS”:Goto0Δ←┘(第一线元数据要素) S∠下一线元起点里程=>O=本线元起点里程:U=本线元起点X:V=本线元起点Y:G=本线元起算方位角:H=本线元长度:P=起点曲率半径:R=终点曲率半径:Q=0或1、-1:Goto0Δ←┘(第二线元数据要素)

坐标正算与反算

一、坐标正算 根据A点的坐标X A、Y A和直线AB的水平距离D AB与坐标方位角αAB,推算B点的坐标X B、Y B,为坐标正算,其计算公式为: X B=X A + ΔX AB Y B=X A+ ΔY AB(1-18) 二式中,ΔX AB与ΔY AB分别称为A~B的纵、横坐标增量,其计算公式为: ΔX AB=X B-X A=D AB · cosαAB ΔY AB=Y B-Y A=D AB · sinαAB(1-19) 注意,ΔX AB和ΔY AB均有正、负,其符号取决于直线AB的坐标方位角所在的象限。 二、坐标反算 根据A、B两点的坐标X A、Y A和X B、Y B,推算直线AB的水平距离D AB与坐标方位角αAB,为坐标反算。其计算公式为: (1-20) (1-21) 注意,由(1-20)式计算αAB时往往得到的是象限角的数值,必须先根据ΔX AB、ΔY AB的正、负号,确定直线AB所在的象限,再将象限角换算为坐标方位角。

三角函数内容规律 三角函数看似很多,很复杂,但只要掌握了三角函数的本质及内部规律就会发现三角函数各个公式之间有强大的联系。而掌握三角函数的内部规律及本质也是学好三角函数的关键所在. 1、三角函数本质: 三角函数的本质来源于定义,如右图: 根据右图,有 sinθ=y/ R; cosθ=x/R; tanθ=y/x; cotθ=x/y。 深刻理解了这一点,下面所有的三角公式都可以从这里出发推导出来,比如以推导 sin(A+B) = sinAcosB+cosAsinB 为例: 推导: 首先画单位圆交X轴于C,D,在单位圆上有任意A,B点。角AOD为α,BO D为β,旋转AOB使OB与OD重合,形成新A'OD。 A(cosα,sinα),B(cosβ,sinβ),A'(cos(α-β),sin(α-β)) OA'=OA=OB=OD=1,D(1,0) ∴[cos(α-β)-1]^2+[sin(α-β)]^2=(cosα-cosβ)^2+(sinα-sinβ)^2 和差化积及积化和差用还原法结合上面公式可推出(换(a+b)/2与(a-b)/2) [1] 两角和公式 sin(A+B) = sinAcosB+cosAsinB sin(A-B) = sinAcosB-cosAsinB

坐标正反算定义及公式-1

第六章→第三节→导线测量内业计算 导线计算的目的是要计算出导线点的坐标,计算导线测量的精度是否满足要求。首先要查实起算点的坐标、起始边的方位角,校核外业观测资料,确保外业资料的计算正确、合格无误。 一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离,计算待定点的坐标,称为坐标正算。如图6-6 所示,点的坐标可由下式计算: 式中、为两导线点坐标之差,称为坐标增量,即:

【例题6-1】已知点A坐标,=1000、=1000、方位角 =35°17'36.5",两点水平距离=200.416,计算点的坐标? 35o17'36.5"=1163.580 35o17'36.5"=1115.793 2、坐标反算 已知两点的坐标,计算两点的水平距离与坐标方位角,称为坐标反算。如图6-6 可知,由下式计算水平距离与坐标方位角。 (6-3) (6-4) 式中反正切函数的值域是-90°~+90°,而坐标方位角为0°~360°,因此坐标方位角的值,可根据、的正负号所在象限,将反正切角值换算为坐标方位角。 【例题6-2】=3712232.528、=523620.436、 =3712227.860、=523611.598,计算坐标方位角计算坐标方位角 、水平距离。

=62°09'29.4"+180°=242°09'29.4" 注意:一直线有两个方向,存在两个方位角,式中:、 的计算是过A点坐标纵轴至直线的坐标方位角,若所求坐标方位角为,则应是A点坐标减点坐标。 坐标正算与反算,可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算,普通科学电子计算器的类型比较多,操作方法不相同,下面介绍一种方法。 【例题6-3】坐标反算,已知=2365.16、=1181.77、 =1771.03、=1719.24,试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键[=]等于纵坐标增量,按储存键[], 键入1719.24-1181.77按等号键[=]等于横坐标增量,按[]键输入,按[]显示横坐标增量,按[]键输入,按第二功能键[2ndF],再按[]键,屏显为距离,再按[]键,屏显为方位角。 【例题6-4】坐标正算,已知坐标方位角=294°42'51", =200.40,试计算纵坐标增量横坐标增量。 键入294.4251,转换为以度为单位按[DEG],按[]键输入,

线路坐标正算及反算程序

一.程序清单: 1.主程序(TYQXJS) "K0="?O:"X0="?U:"Y0="?V:"FWJ="?G:"LS="?H:"+1,0,-1="?Q:If 0= Q:Then 1e45→P:1e45→R:Goto 0:Else "R1="?P:"R2="?R:Goto 0:←┘ Lbl 0:(P-R)÷(2HPR) →D:"Z=1,F=2"? N:If 1=N:Then Goto 1:Else Goto 2:←┘ Lbl 1:"SK="?S:"JZ="?Z:Abs(S-O) →W:Prog "SUB1":F-90→F:"X=":X◢"Y=":Y◢"FWJ=":F◢Goto 1◢ Lbl 2:"CX="?X:"CY="?X→I:Y→J:Prog "SUB2":"O+W→S:Cls:Locate 1,1,"SK"←┘Locate 6,1,S←┘Locate 1,3,"JZ" ←┘Locate 6,3,Z←┘Goto 2◢ 2. 正算子程序(SUB1) 0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M: U+W(Acos(G+57.29577951QKW(1/P+KWD))+Bcos(G+57.29577951QLW(1/P+LWD))+Bco s(G+57.29577951QFW(1/P+FWD))+Acos(G+57.29577951QMW(1/P+MWD))) →X: V+W(Asin(G+57.29577951QKW(1/P+KWD))+Bsin(G+57.29577951QLW(1/P+LWD))+Bsin (G+57.29577951QFW(1/P+FWD))+Asin(G+57.29577951QMW(1/P+MWD))) →Y: G+57.29577951QW(1/P+WD)+90→F:X+ZcosF→X:Y+ZsinF→Y:Return 3. 反算子程序(SUB2) G-90→T:Abs((Y-V)cosT-(X-U)sinT)→W:0→Z:Lbl 0:Prog "SUB1":T+Q57.295779 51W(1/P+WD) →L:(J-Y)cosL-(I-X)sinL→Z:If AbsZ<1E-6:Then Goto1:Else W+Z→W:Goto 0←┘ Lbl 1:0→Z:Prog "SUB1":(J-Y)÷sinF→Z:Return 二.增设数据库程序(SJK) Lbo0:{S}:If S∠下一线元起点里程Then O=起点里程:U=起点X:V=起点Y:G=起算方位角:H=元长度:P=起点半径:R=终点半径:Q=0或1、-1:Prog“TYQX

坐标正算反算公式讲解..

在高斯直角坐标系中,由坐标纵轴方向的北端起,顺时针量到 直线间的夹角,称为该直线的坐标方位角,常简称方位角,用a表示。 1、第一象限的方位角 X A 第四象限第一象限 a 0 Y 第三象限第二象限 图1 2、第二象限的方位角 X 第四象限第一象限 a 0 Y A 第三象限第二象限 图2

3、第三象限的方位角 4、第四象限的方位角 X 第四象限 第一象限 1 - 'o a 第三象限 第二象限 X 第四象限 第一象限 A 第三象限 第二象限

方位角计算公式: A 图4

M AO 方位角的计算器计算程序:P O I(X A-X O,Y A-Y。) 直线OA方位角度值赋予给计算器的字母J, 0< JV 360。 直线段OA的距离值赋予给计算器的字母1,1 >0 直线OA与直线AO的方位角关系: 当直线OA的方位角W 180°时,其反方位角等于a+180°。2、当直线OA的方位角〉180°时,其反方位角等于a-180 °。 二方位角的推算 (一)几个基本公式 1、坐标方位角的推算 a = a + P ± d on 前后左I 80 ziYAo *A a = tan - X.-X) 1、

或: a 前 =a 后 一 3右± 180 注意:若计算出的方位角>360°,则减去360°;若为负值,则加上 360° 例题:方位角的推算 45 、 a 51 o a 23= a 12- P 2+180° =30° -130 +180° =80已知:a 12=30 ,各观测角P 如图,求各边坐标方位角 ot 23 、 a 34、 a 解:

坐标正算 直线段的坐标计算 E 、 X a V B 厂 D ” 丿a A ” \ ' 、C 设起点0的坐标(X o ,Y ),直线OP 的方位角为F oP ,求A 、C 、E 点 的坐标 X A =X O +L X Cos(F op ) Y A 二Y o +L X Sin(F op ) 设直线段OB 长度为L OB ,直线段BC 长度为L BC ,则C 点坐标为 X B =X O +L OB X Cos(F op ) p 3+180° =80° -65 ° + 180° = 195° B4+180° = 195° -128 + 180° =247° p5+180° =247° -122 + 180° =305° p 1 + 180° =305 ° -95 + 180° =30° Ot 45= a 34- Ot 51= a 45- ot 12= a 51- (检 查) Ot 34 = a 23- 1、 设直线段OA 长度为L,则A 点坐标为 2、

高斯投影坐标正反算公式

高斯投影坐标正反算公式 未知2010-04-03 10:47:15 本站 §高斯投影坐标正反算公式 任何一种投影①坐标对应关系是最主要的;②如果是正形投影,除了满足正形投影的条件外( C-R 偏微分方程),还有它本身的特殊条件。 1.1 高斯投影坐标正算公式: B, x,y 高斯投影必须满足以下三个条件: ①中央子午线投影后为直线;②中央子午线投影后长度不变;③投影具有正形性质,即正形投影条件。 由第一条件知中央子午线东西两侧的投影必然对称于中央子午线,即 (8-10) 式中, x 为 的偶函数, y 为的奇函数;,即,如展开为的级数,收敛。 ( 8-33 ) 式中是待定系数,它们都是纬度 B 的函数。 由第三个条件知: (8-33) 式分别对和 q 求偏导数并代入上式

(8-34) 上两式两边相等,其必要充分条件是同次幂前的系数应相等,即 (8-35) (8-35) 是一种递推公式,只要确定了就可依次确定其余各系数。 由第二条件知 : 位于中央子午线上的点,投影后的纵坐标 x 应等于投影前从赤道量至该点的子午线弧长 X ,即 (8-33) 式第一式中,当时有: (8-36) 顾及 ( 对于中央子午线 )

得: (8-37,38) (8-39) 依次求得并代入 (8-33) 式,得到高斯投影正算公式 (8-42) 1.2 高斯投影坐标反算公式 x,y B, 投影方程: (8-43)

满足以下三个条件:

①x 坐标轴投影后为中央子午线是投影的对称轴;② x 坐标轴投影后长度不变;③投影具有正形性质,即正形投影条件。 高斯投影坐标反算公式推导要复杂些。 ①由 x 求底点纬度 ( 垂足纬度 ), 对应的有底点处的等量纬度,求 x,y 与 的关系式,仿照 (8-10) 式有, 由于 y 和椭球半径相比较小 (1/16.37) ,可将展开为 y 的幂级数;又由于是对称投影, q 必是 y 的偶函数,必是 y 的奇函数。 (8-45) 是待定系数,它们都是 x 的函数 . 由第三条件知: , , (8-21)

坐标正反算vb测量程序

坐标正反算vb测量程序 河北工程大学测绘1001 刘长君 以下是登陆界面和计算界面的截图: 以下就是对应的代码: 《登陆界面的代码》 Private Sub Command1_Click() If = "刘长君" And = "0" Then Unload Me Else Dim a As Integer a = MsgBox("密码输入有误,请重新输入", 1, "提示框")

If a = 1 Then = "" End If End If End Sub Private Sub Command2_Click() End End Sub 《计算界面的代码》 Private Sub Command1_Click(Index As Integer) Dim a1%, a2%, a3%, a4!, a5!, pi! Const conpi = If IsNumeric(Text1(1)) = False Or IsNumeric(Text2(2)) = False Or IsNumeric(Text3(3)) = False Or IsNumeric(Text4(4)) = False Or IsNumeric(Text5(5)) = False Or IsNumeric(Text6(6)) = False Then MsgBox "输入有误" Text1(1).SetFocus Else If Val(Trim(Text5(5).Text)) > 60 Or Val(Trim(Text5(5).Text)) < 0 Or Val(Trim(Text6(6).Text)) > 60 Or Val(Trim(Text6(6).Text)) < 0 Then MsgBox "方位角输入有误,请重新输入" Text4(4).SetFocus Text4(4).Text = "" Text5(5).Text = "" Text6(6).Text = "" End If a1 = Val(Text4(4).Text) a2 = Val(Text5(5).Text) a3 = Val(Text6(6).Text) a4 = a1 + a2 / 60 + a3 / 3600 a5 = a4 * conpi / 180 = Val(Text3(3).Text) * Cos(a5) = Val(Text3(3).Text) * Sin(a5) = Val(Text1(1).Text) + Val = Val(Text2(2).Text) + Val = Format((Val * 10000 + / 10000, "####.0000") = Format((Val * 10000 + / 10000, "####.0000") = Format((Val * 10000 + / 10000, "####.0000") = Format((Val * 10000 + / 10000, "####.0000") End If End Sub

相关文档
最新文档