坐标法测图的优化程序（转帖）

作者：张明睦 出处：《北京测绘》1999年第1期
一、问题提出
在实施坐标法施测大比例尺地形图时，常用是视距极坐标法观测，支距法现场展绘，观测数据按常规的顺序输入，展点以图幅的左下角为零点，数据在图上成1:1000的比例。由于CASIOfx--3600P计算器的程序步的极限， 则就要优化坐标法则图的程序。
二、优化的CASIOfx--3600P计算器程序:

MODE O INV PCL P1
Kout 4 Kin 1 Kout 5 Kin 2 Kout 6 Kin 3
ENT 1 × ENT 1 Kin - 3 C ENT 1 INV Min
Sin × MR COS × Kin + 3 MR tan ÷ 0.5
== INV P→R ENT 1 == Kin + 1 INV X→Y
Kin + 2 Kout 1 INV HLT Kout 2 INV HLT
Kout 3 INV HLT INV RTN MODE .
三、单元分配:
常量:A0=(XA-X0)÷K → Kin 4
B0=(YA-Y0)÷K → kin 5
HO=(HA+i) → kin 6
变量:AP=(XP-XO)÷K → kin 1
Bp=(Yp-Y0)÷k → Kin 2
Hp=H0+dh-L →Kin 3
Z→Min
四、操作应用:
仪器设置在测站点上，水平度盘配置后视点的方位值，程序编入计算器。把常量A0，B0，H0分别存入计算器的Kin4，Kin5，Kin6。启动该程序P1，即显示视线高H0，依次输入视距S，截尺L，天顶距Z，即显示(D÷K)图上成1:1000的平距值， 输入方位值α，显示Ap,Bp,Hp(即以图幅左下角为零点，在图上成1:1000 的坐标值及高程)。程序自动返回。
五、算例
测图比例设为1:500，则K=500÷1000=0.5
测站点:XA=46638.737 HA=323.17
YA=39215.655 i=1.50
图幅左下角零点坐标:X0=46400
Y0=39100
先计算存入:(XA-X0)÷k=477.474 → Kin 4
(YA-Y0)÷k=231.31 → Kin 5
HA+i=324.67 → Kin 6
再启动程序P1，显示HA+i(324.67)，输入视距(158.0)，截尺(1.0),天顶距(96°58′)显示成1:1000的平距值D÷K(311.36),输入方位值(185°02′)显示成1: 1000的AP(167.324),BP(203.994)及高程HP(342.69);程序自动返回。
六、结论:
该程序中，仅改变K值(即0.5)，就可以应用于不同比例尺的坐标法测图。改变符号(即Kin+3)中的“+”为“-”。就可以应用设置不同值的天顶距的仪器。 应用方便、直观，提高外业现场展点的速度，使坐标法则更加完善，逐步向数字化测图系统进入。