5800坐标正反算
5800坐标反算：
曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序程序功能
四、算例
某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组
成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度
LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：
S0 X0 Y0 F0 LS R0 RN Q
500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0 769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1 806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1 919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1 999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0
1、正算
（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元
要素）
S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590
S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837
S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084
S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270
S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358
S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446
S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091
S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379
S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668
S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802
S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642
S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168
2、反算
X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164
X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136
X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137 X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663 X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578
X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049
X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814
X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656
X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937
X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710
X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694
进行坐标反算:
输入2 X0= 612961.9828 Y0=513124.0300 QDZH=4535.940 FWJ=160.699092 LS=2266.37 QDR=11005 ZDR=11005 Q=1
三、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的
值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。
X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
QDZH ？线元起点里程
FWJ ？线元起点切线方位角
LS ？线元长度
QDR ？线元起点曲率半径
ZDR ？线元止点曲率半径
Q ？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）
S ？正算时所求点的里程
Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=×××正算时，计算得出的所求点的X坐标
YS=×××正算时，计算得出的所求点的Y坐标
S=×××反算时，计算得出的所求点的里程
Z=×××反算时，计算得出的所求点的边距
坐标正算反算公式讲解(借鉴材料)
一 方位角： 在高斯直角坐标系中，由坐标纵轴方向的北端起，顺时针量到直线间的夹角，称为该直线的坐标方位角，常简称方位角，用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)
FX5800P计算器坐标正反算程序Word文档
(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。) 卡西欧fx5800p计算器坐标正反算程序
一、程序功能 本程序由 6 个主程序、 5 个次子程序及 5 个参数子程序组成。主要用于公路测量中坐标正反算，设计任意点高程及横坡计算 , 桥涵放样，路基开挖口及填方坡脚线放样。程序坐标计算适应于任何线型 . 二、源程序 1. 主程序 1 ：一般放样反算程序(① 正算坐标、放样点至置仪点方位角及距离；② 反算桩号及距中距离 ) 程序名 :1ZD-XY Lb1 0:Norm 2 F=1 ： ( 正反算判别， F=1 正算， F=2 反算 , 也可以改 F 前加？，改 F 为变量 ) Z[1]=90 （与路线右边夹角） Prog ＂ THB ＂： F=1=>Goto 1:F=2=>Goto 2 Lb1 1: Ｆ ix 3: ＂Ｘ = ＂： Locate 6,4, Ｘ◢ ＂Ｙ＝＂： Locate 6,4, Ｙ◢ Ｐ rog ＂3JS”:Goto 0: Lb1 2:Fix 3: ＂ＫＭ＝＂： Locate 6,4, Ｚ◢ ＂Ｄ＝＂： Locate 6,4, Ｄ◢ Ｇ oto 0 ２．主程序２：高程序横坡程序 ( 设计任意点高程及横坡 ) 程序名： 2GC LbI 0:Norm 2 “KM”?Z:?D: Prog”H”:Fix 3:” H=”:Locate 6,4,H◢ “ I=”: Locate 6,4,I◢ Goto 0 3. 主程序 3 ：极坐放样计算程序 ( 计算放样点至置仪点方位角及距离 ) 程序名： 3JS X ： Y ： 1268 ．123→K( 置仪点 X 坐标 ) 2243 ．545→L （置仪点 Y 坐标，都是手工输入 , 也可以建导线点数据库子程序 , 个人认为太麻烦） Y-L→E ： X-K→F ： Pol(F,E):IF J<0:Then J+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:( 不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢ 来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢( 不习惯小数点后四位为角度显示的，也可以用命令 J◢DMS◢ 来直接显示 ) Fix 3:” S=”:Locate 6,4,I◢ 4 ．主程序 4 ：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序 3 中输入置仪点坐标后计算放样点至置仪点方位角及距离 ) 程序名： 4JH-XY LbI 0:Norm 2 90→Z[1]( 涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量 )
高斯投影坐标正反算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，
卡西欧5800计算器坐标正反算程序
M = (1.0/Re-1.0/Rs)/Ls; x=∫｛cos(Ta + L/Rs + 0.5*M *L*L)，0，L｝; y=∫｛sin(Ta +L/Rs + 0.5*M *L*L)，0，L｝; a(i)= Ta +L/Rs + 0.5*M *L*L Rs:缓和曲线起点半径 Re:缓和曲线止点半径 Rs，Re （NE坐标系下，右偏为正，左偏为负） Ta:缓和曲线起点的真北方位角 Ls:不完整缓和曲线长度。 此公式为缓和曲线在坐标系下任意位置的通用积分公式，能完全适应缓和曲线左偏、右偏、Rs >Re 、Rs
5800曲线道路全线坐标正反算(已验证)
CASIO fx-5800P线元法坐标正反算程序V2.0 说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。 主程序名：ABCYT 第1步Deg：Fix 3：10→DimZ 第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂K0+＂?S：Prog＂ABCYTSJ＂：If P=0：Then 10^(45)→P：IfEnd：If R=0：Then 10^(45)→R：IfEnd 第3步1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If W=1：Then Goto 1：Else Goto2：IfEnd 第4步Lbl 1：＂E＂?Z 90=N：Abs(S-O)→W：Prog ＂ABCYTZ＂ 第5步Cls：＂F=＂：Locate 3,1,F°：＂X=＂：Locate 3,2,X：＂Y=＂：Locate 3,3,Y◢ 第6步Prog＂ABCYTSQX＂：Cls：＂H=＂：Locate 3,1,H◢ 第7步1→W：90→N：Goto 3 第8步Lbl 2：?X：?Y：X→I：Y→J：Prog＂ABCYTF＂：O+W→S 第9步Cls：＂K=＂：Locate 3,1,S：＂E=＂：Locate 3,2,Z◢ 第10步2→W：Goto 3 正算子程序名：ABCYTZ 第1步0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L 第2步1-L→F：1-K→M 第3步U+W×(A×cos(G+Q×E×K×W×(C+K×W×D))+B×cos(G+Q×E×L×W×(C+L×W× D))+B×cos(G+Q×E×F×W×(C+F×W×D))+A×cos(G+Q×E×M×W×(C+M×W×D)))→X 第4步V+W×(A×sin(G+Q×E×K×W×(C+K×W×D))+B×sin(G+Q×E×L×W×(C+L×W× D))+B×sin(G+Q×E×F×W×(C+F×W×D))+A×sin(G+Q×E×M×W×(C+M×W×D)))→Y 第5步G+Q×E×W×(C+W×D)→F：F+ N→Z[1] 第6步X+Z×cos（Z[1]）→X：Y+Z×sin（Z[1]）→Y 反算子程序名：ABCYTF 第1步Lbl 2：（S-O）→W：0→Z：Prog ＂ABCYTZ＂：F-90→Z[9]：(J-Y)×cos(Z[9])-(I-X)×sin(Z[9])
坐标正反算计算公式
坐标正反算公式
一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现GPS 网成果由基线解算时GPS 卫星星历所采用的参照系（WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。 二、南方GPS数据处理软件的平差方式
三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图：
如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。
卡西欧5800公路坐标正反算程序教程文件
卡西欧5800公路坐标正反算程序
目录 一、坐标正算基本公式 (02) 二、坐标反算原理 (04) 三、高程数据库录入变换 (05) 四、计算器程序 (07) 01、ZBZS(坐标正算) (07) 02、ZBFS（坐标反算） (08) 03、GCJF(高程积分) (09) 04、PJFY（坡脚放样） (10) 05、JFCX（积分程序） (11) 06、ZBFY（坐标放样） (11) 07、DT（递推） (12) 08、HP（横坡） (13) 09、LK（路宽） (14) 10、SJK1（平面数据库） (14) 11、SJK2（纵面数据库） (14) 12、SJK3（左路宽度数据库） (15) 13、SJK4（右路宽度数据库） (15) 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库） (16) 16、SJK7（左上边坡数据库） (17) 17、SJK8（右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序
一、正算所涉及的计算公式 X R i d X αβB d Y d l d βI 图表 1 在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标 点。设： A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为 ρA ,A 点的里程为L A ,B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的 里程为L I 。I 点的切线角为β。 由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。故此任 意点的曲率为; C L R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为; A B A i A B A i L L L L --?-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。设：
坐标正反算定义及公式
坐标正反算定义及公式 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，试计算纵坐标增量横坐标增量。
5800坐标程序
颜色文字Casio5800交点法与线元法（积木法）匝道坐标正反算放样程序（XUFENG 2011.2.14） 本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《Casio5800交点法程序》与《线元法（积木法）匝道坐标正反算放样程序》，根据歪哥意见“需要的自行修改结合XY框架自己修改为数据库反算程序等”，本人不才，采用最笨的办法将两个程序综合了一下，使之能既能进行交点法正反算，又能进行线元法正反算。在此特别感谢大歪哥！ 将程序发上来，愿与大家一同交流学习欢迎大家吐口水，只要能进步就行！ 程序由一个主程序ZBZFS和8个子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标正反算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并能对单一线元进行坐标正反算。 1主程序名：ZBZFS（功能：进入计算主程序） 65→Dimz↙ Deg:Fix 3↙ ＂1.JD ZFS 2. ZHADAO ZFS＂? I: I→Z[61]: ＂1.ZHONG SHU JS 2. JS＂? I ↙ If I=1: Then Goto1: Else Goto2:IfEnd↙ LbI 1 :If Z[61]=1: Then Prog＂JDYS＂:Else Cls:＂K0＂?A:＂KN＂?L :＂X0＂?U :＂Y0＂?V :＂F0＂?W :＂R0＂?P :＂RN＂?Q:＂ZX:-1,+1,0＂?G:IfEnd↙ LbI 2 :Prog＂JS＂ 2子程序名：JS（功能：选择正算或反算模式） Cls:＂XC＂?H:＂YC＂?Z↙
坐标反算正算计算公式
坐标反算正算计算公式 一、坐标正算 根据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公路坐标正反算程序
坐标正反算程序 程序变量说明: D:选线。E、F:置仪点坐X、Y。K:待算点桩号。B:待算点距离。Q:待算点与线路交角（以顺时针方向）L:线元起点至待算点长度。0:待算点方位角。S:线元长度。U:线元起点半径。V:线元终点半径。 X、Y:正算结果。Z、计算序列。C、I、J T:内部运算。A、G、H、M、N、P、R、W未用变量。 ZB- Deg:Norm 2 Z[1]T Z:”两个空格ZS=1 Or FS=2Z If Z M 1 And Z M2:Then Cls:Locate 5,2,”NO 10r 2”：Locate 6,3 !Stop!”： Stop:Else Z Z[1]:lfEnd Z[10]T D:四个空格XL T 1 T 2T3九个空格XL N=N”D If D M 1 And D M 2 And D M 3:Then Cls:Locate 5,2,'XL r 1 T2T3”:Locate 5,3 !Stop!”： Stop: Else D T Z[9]: IfEnd If D=Z[10]:The n Cls:Do:Locate 3,1,"[EXE] GO!GO!':Locate 8,2,9-I ntg(Z/7) Locate 4,3,"10S T DAT":Locate 12,3,D:Isz Z:Z=70=>Goto 1:Lpwhile Getkey M 47:Goto 2:IfE nd Lbl 1 D=1=>Prog” ”:D=2=>Prog” ”:D=3=>Prog” ”:D=X=>Prog‘ D T Z[7]:For 1 T Z To Z[7]-1 Cls: Locate 9,2,Z[7]-Z Z[10Z+11]-Z[10Z+1]T L Prog" TQSJ:Next Lbl 2 Z[9] T Z[10] Cls:Z[1]=2=>Prog”FS' Z[2]T E: (X)=”E:E T Z[2]:Z[3] T F：”(Y)仝F:F T Z[3]:Prog ”ZS'
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。
卡西欧5800公路坐标正反算程序
目录 一、坐标正算基本公式 (02) 二、坐标反算原理 (04) 三、高程数据库录入变换 (05) 四、计算器程序 (07) 01、ZBZS(坐标正算) (07) 02、ZBFS（坐标反算） (08) 03、GCJF(高程积分) (09) 04、PJFY（坡脚放样） (10) 05、JFCX（积分程序） (11) 06、ZBFY（坐标放样） (11) 07、DT（递推） (12) 08、HP（横坡） (13) 09、LK（路宽） (14) 10、SJK1（平面数据库） (14) 11、SJK2（纵面数据库） (14) 12、SJK3（左路宽度数据库） (15) 13、SJK4（右路宽度数据库） (15) 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库） (16) 16、SJK7（左上边坡数据库） (17) 17、SJK8（右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序
一、正算所涉及的计算公式 X R i d X αβB d Y d l d βI 图表 1 在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标点。设： A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为ρ A ,A 点的里程为L A , B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的里 程为L I 。I 点的切线角为β。 由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。故此任 意点的曲率为; C L R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为; A B A i A B A i L L L L --?-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。设：
5800P编程计算器坐标正反算程序
5800P编程计算器坐标正反算程序 一，主程序文件名（TYQXJS） LbI 4:“1.SZ=> XY’’:’’2.XY=>SZ’’︰?N:?S:Prog“SUB0” 1÷P→C：(P－R) ÷(2HPR) →D：180÷π→E：N=1=> Goto1:Goto2 LbI 1: ？Z:Abs(s－O)→W︰Prog“SUB1”︰“XS=”:X?“YS=”:Y? F－90→F:”FS=”:F?DMS? Goto 4 LbI 2:?X:?Y:X→I:Y→J:Prog”SUB2”:O＋W→S﹕“S=”S?“Z=”：Z? Goto 4 二:子程序（文件名SUB1） 0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1－L→F:1－K→M:U＋W(Acos(G＋QEKW(C＋KWD))＋Bcos(G＋QELW(C＋QELW(C＋LWD))＋Bcos(G＋QEFW(C＋FWD))＋Acos(G＋QEMW(C＋MWD)))→X:V＋W(Asin(G＋QEKW(C＋KWD))＋Bsin(G ＋QELW(C＋LWD))＋Bsin(G＋QEFW(C＋FWD))＋Asin(G ＋QEMW(C＋QEMW(C＋QEMW(C＋MWD)))→Y:G＋QEW(C＋WD)＋90→F:X＋Zcos(F)→X:Y＋Zsin(F)→Y 三:反算子程序(文件名SUB2)
G－90→T:Abs((Y－V)Cos(T)－(X－U)Sin(T)→W:0→Z:LbI 0:Prog“SUB1”:T＋QEW(C＋WD)→L:(J－Y)cos(L)－(I－X)sin(L)→Z:IfAbs(Z)＜0.000001:ThenGoto1:Else W＋Z→W:Goto0 LbI 1:0→Z:Prog“SUB1”:(J－Y)÷sin(F)→Z: 四:数据库(文件名) LbI 1:If S＜止点里程:Then起点方位角→G:起点里程→O:起点X坐标→X:起点Y坐标→Y:起点半径→P:止点半径→R:曲线长→H:左转(-1)右转(1)直线(0)→Q:Return:IfEnd 数据库说明: 每增加一段曲线就增加以上一句上面的语句，当半径为无穷大时以10的45次方代替，左转输入负1，右转输入正1 五:程序输入及显示说明 1:输入部分 1，SZ=>XY 2，XY=>SZ N? 选择计算方式输入1表示由里程边距计算坐标(正算)输入2表示由坐标计算里程及边距(反算) S? 正算时输入所求点的里程，反算时为所求点的近似里程(注意反算时的近似里程不能跨曲线，如果跨曲线需要进行两次反算，并以第一的里程结果作为近似里程，否则，反
测绘技术之坐标反算与正算
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)
卡西欧fx-5800P直线圆曲线坐标正反算程序
[精] fx 5800 直线圆曲线坐标正反算程序 (2010-07-01 21:50:11) 标签： 杂谈 直线计算程序 0→I:0→J:”X 0”?D:”Y0”?E:”X1”?B:”Y1”?C:Pol(B-D,C-E):J→A:If A<0: Then A+360→A:Else A→A:IfEnd:Lbl 0:?O:?S:If O≠0:Then Goto 1:IfEnd: D+Scos(A) →X:”X=”:X◢ E+Ssin(A) →Y:”Y=”:Y◢ Goto 0:Lbl 1:D+Scos(A)+Ocos(A+90) →X:”X=”:X◢ E+Ssin(A)+Osin(A+90) →Y:”Y=”:Y◢ Goto 0 输入程序时注意区别字母O与数字0 程序运行时符号说明 X0? Y0?分别输入直线起点的XY坐标值 X1? Y1?分别输入直线终点的XY坐标值 O? 输入边桩与中桩的距离(左边桩为负值,右边桩为正值),如计算中桩坐标输入0 S? 输入所求点到直线起点的距离 圆曲线计算程序 0→I:0→J: ”X0”?C:”Y0”?D:”X1”?E:”Y1”?F:?R:”L:-1 R:1”?N:”ZY”?W:Lbl 0:”LN”?T: T-W→O:Pol(E-C,F-D):I→S:J→A: If A<0:Then A+360→A:Else A→A:IfEnd: sin-1(S÷(2R)) →K:2∏RK÷180→L:180O÷(2∏R) →G:(2R)sin(G) →H:C+Hcos(A-KN+GN) →X:”X=”:X◢ D+Hsin(A-KN+GN) →Y:”Y=”:Y◢ R-0.5√(4R2-H2) →Q:”Q=”:Q◢ “S=”?V:If V=0:Then Goto 0:IfEnd:”L:-90 R:90”?U:X+Vcos(A-KN+2GN+U) →X:”X=”:X◢ Y+Vsin(A-KN+2GN+U) →Y:”Y=”:Y◢ Goto 0 输入程序时注意区别字母O与数字0 程序运行时符号说明 X0? Y0?分别输入直线起点的XY坐标值 X1? Y1?分别输入直线终点的XY坐标值 R? 输入圆曲线半径 L:-1 R:1?圆曲线向左转弯时输入-1,向右转弯时输入1 ZY? 输入起点桩号
坐标正反算
一、坐标正算与坐标反算 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］，