工程测量程序集

1、测量程序集 survey program 重庆臻品公路工程建设有限公司测量部 e-mail: 目 录1：zxzbjs:直线坐标计算 32：fy：点位转点（1） 33：fy：点位转点（2） 34：zbjs:坐标计算 35：bjjs:边长、方位角计算 46：zbfs坐标反算 47：hyzbjs:缓和曲线坐标计算（包括圆曲线） 48：mjjs：面积计算 59：yoxzbjs：圆曲线坐标计算 510：zxlcjs：直线里程计算 511： yqxlcjs:园曲线点的里程计算 612：sqxgcjs:竖曲线高程计算 613： hyzxjdjs:缓

2、和与直线交点坐标计算 614：hylcjs：缓圆曲线里程计算 715：fgpjgcjs：方格网平均高程计算 816:scqxjg:四次曲线加宽 817：scqxbz直线段四次曲线，边向中的法线 818:scqxzb:直线段四次曲线，中线里程的 919：zxjdzblc直线交点坐标、里程 920：zxjdzblc-1直线交点坐标、里程 1021:sjcl视距测量 （1） 1022:sjcl视距测量（2） 101：zxzbjs:直线坐标计算kxyf:lb1 0:s:m=x+(s-k)cosfn=y+(s-k)sinflb1 1:g:g0=goto 0: rec(g,f+90):a=m+gcos(f

3、+90) b=n+gsin(f+90) c=m-gcos(f+90) d=n-gsin(f+90) goto1说明：依次输入直线里程k，坐标x、y，方位角f。输入里程s可得中线点坐标m、n。输入距中线宽g可得右边点坐标a、b，左边点坐标c、d。g输入-1返回输入新的里程s。2：fy：点位转点（1）(4500) x=g+r*cosf l1y=h+r*sinf l2fixm l3pol(x-a,y-b) w=ww0=w=w+360 说明: g、h为被转点x、y坐标.r为转点距离。f为转点的方位角。a、b为测站点x、y坐标。3：fy：点位转点（2） lb1 1:x=tx1+r*cosfy=uy1r*

4、sinfab：rf：s=（x-a）2（y-b）2）w=sin-1(y-b)/s:x-aw=180-w: ww=w+360: wgoto 1 说明：同。4：zbjs:坐标计算 ab：lb1 1:sf：x=a+scosfy=b+ssinfgoto 1 说明: 输入坐标a、b,边长s,方位角f,得坐标x、y。5：bjjs:边长、方位角计算 ab:lb1 1:xy:s=(x-a)2+(y-b)2) f=sin-1(y-b)/s):x-af=180-f: ff=f+360: fgoto1 说明:输入起点坐标a、b, 终点坐标x、y得边长s,方位角f。6：zbfs坐标反算 pol(a-b,c-d) w=w

5、w0=w=w+360 说明：a、b已知点x坐标，c、d已知点y坐标。显示距离、方位角。 7：hyzbjs:缓和曲线坐标计算（包括圆曲线） kxyfjirh：p=hh/24r-h4/2688r3：q=h/2-h3/240rr：b=90h/r：c=rh：t=（r+p）tan（i/2）+q:l=r（i-2b）/180+2h lb1 0:s:sgoto3：s=s-k:ss=-s: sh=goto1: m=s-s5/40cc+s9/3456c4:n=s3/6c-s7/336c3+s11/42240c5:goto2lb1 1:sl-h=goto5: o=180(s-h)/ r+b:m=rsino+q:n=

6、r(1-coso)+p:lb1 2:n=nj:u=x+mcosf-nsinfv=y+msinf+ncosf lb1 4:g:ggoto 0: sh=a=f+j(b+180(s-h)/ r)+90: sa=f+90jss/c+90: m=u+gcosa n=v+gsinaw=u-gcosa z=v-gsina goto 4lb1 5:d=x+2tcos(i/2)cos(f+ij/2):e=y+2tsin(f+ij/2) cos(i/2):o=f+ij+180:o360=o=o-360: s=l-s:m=s-s5/40cc+s9/3456c4:n=s3/6c-s7/336c3+s11/42240c

7、5:n=-nj:u=d+mcoso-nsino v=e+msino+ncoso lb1 6: g:ggoto 0: a=o-90jss/c+90:m=u-gcosa n=v-gsina w=u+gcosa z=v+gsina goto 6:lb1 3说明：依次输入直缓点里程x、 y坐标，方位角f，曲线偏转方向j左转-1右转1、偏角i、半径r、缓和曲线长h。输入里程s可得中线坐标u、v，输入距中线宽g可得右边点坐标m、n左边w、z。g为-1返回s输入新的里程。 8：mjjs：面积计算lb1 2：s=0：ab：e=a：f=b：lb1 0：cd:c=-99=goto1: s=s+ad-bc:a=c:

8、b=d:goto 0:lb1 1:s=(s+af-be)/2goto2 说明：ab起点坐标，cd依次各点坐标，最后输入为-99结束，s为面积。 9：yoxzbjs：圆曲线坐标计算 kxyfhr：lb1 0：s：o=90（s-k）/r:e=2rsino：rec（e，f+ho）：m=x+in=y+jlb1 1：g：ggoto 0: rec(g,f+2ho+90):a=m+ib=n+jc=m-id=n-jgoto 1说明：依次输入圆曲线起点里程k坐标xy方位角f曲线方向h左偏为-1右偏为1半径r。以下步骤同直线计算。 10：zxlcjs：直线里程计算kxyf：lb1 0：mn：pol（m-x，n-y

9、）：jj=j+360: l=k+icos(j-f) s=isin (j-f) goto 0 说明:输入直线起点里程k坐标xy方位角f。输入所求点坐标mn，l为所求点里程s为该点到中线距离，s为正该点在中线右,s为负在中线左边。 11： yqxlcjs:园曲线点的里程计算kxyfhr：c=f+90h：rec（r，c）：a=x+i：b=y+j：lb1 0mn：pol（m-a，n-b）：jj=j+360: g=c-180:gg=g+360: hgoto2: jj=j+360: l=r(j-g) /180:d=k+ls=i-r:s=-sgoto 0 lb1 2:gg=g+360: l =r(g-j)

10、/180:d=k+ls=i-rgoto 0说明：输入园曲线起点里程k，坐标xy，方位角f，方向h左转为负，右转为正，半径r。输入所求点坐标mn，显示里程d，距中距离s，s为正点在中线右边，s为负点在中线左边。 12：sqxgcjs:竖曲线高程计算： badtrzlb1 0: c:cz1=h=a+（b-a）/t*（c-d）-（c-d）2/（2r）h=a+（b-a）/t*（c-d）+（c-d）2/（2r）goto 0说明：“a”起点（或终点）高程；“c”需计算点桩号；“t”竖曲线的切线长；“d”起点或终点的桩号；“r”竖曲线的半径；“z”曲线凸凹判别符z1时为凸曲线，zh=goto 2:dgoto

11、 3:a=s-s5/90r2h2:b=f+30gs2/rh:p=x+acosb:q=y+asinb:goto1:lb1 2:a=h-h3/90r2:b=f+30hg/r:c=x+acosb:e=y+asinb:p=c:q=e:s=0:l=h:lb1 4:d=abs(q-n)cosi-(p-m)sini:s=s+d:dgoto3:a=2rsin(90s/r):b=f+90g(h+s)/ r:p=c+acosb: q=e+sinb:goto4:lb1 3:k+sz+zlpqgoto 0说明：输入直缓点里程k、坐标xy、方位角f、曲线方向g左转-1右转1、半经r、缓和曲线长h、正反算z正算1反算-1

12、。输入直线点起点坐标mn、方位角i，显示交点里程ksz+zl、坐标pq。继续输入新的直线起点坐标mn、方卫角i。 14：hylcjs：缓圆曲线里程计算 kxyfgrhz:c=hr:lb1 0:mn:pol(m-x,n-y):jj=j+360: d=j-f:p=icosd:q=isind:g0=q=-q: l=i: lb1 :a=l+qll/2c+7l5/120cc+ql6/24c3+ql10+/240c5+0.0005806l13/c6+17ql14/40320c7+139l9/24192c4-p: b=1+ql/c+7l4/24cc+ql5/4c3+139l8/2688c4+ql9/24c5

13、+0.00754785l12/c6+17ql13/2880c7:t=l-a/b:u=t-l:l=t:absu0.0005=goto1: lh=goto2: w=l-l5/40cc+l9/3456c4:e=l3/6c-l7/336c3+l11/42240c5:o=90ll/c: rec(ww+ee),f+og/3):x+iy+jpol(m-x-i,n-y-j):jj=j+360:sin(j-f-og)i=-i: k+zls=igoto 0:lb1 2:a=h-h3/40rr+h5/3456r4:b=hh/6r-h4/336r3:b=bg:pol(a,b):rec(i,f+j):a=x=i:b=y

14、+j:d=90h/r:rec(r,f+gd+90g):p=a+i:q=b+j:o=f+gd+90g+180:o360=o=o-360:pol(m-p,n-q):jj=j+360:t=j:v=i:rec(r,t):p+iq+jggoto4: tt=t+360: l=(t-o) r/180:k+zl+zhs=r-vgoto 0:lb1 4:oo=o+360: l=(o-t) r+180:k+zl+zhs=v-rgoto 0 说明：输入直缓点里程k、坐标xy、方位角f、曲线方向g左转-1右转1、半经r、缓和曲线长h、正反算z正算1反算-1。输入所求点坐标mn，首先显示线路中点坐标x+i、y+j或p+

15、i、q+j，显示里程k+zl或k+zl+zh，最后显示距中点的距离s，s为正点在前进方向右边、负在左边。继续输入新的所求点坐标mn 15：fgpjgcjs：方格网平均高程计算 lb1 0:n:s=0:lb1 :a:a=-99=goto 2: s=s+a:goto 1:lb1 2: b:b=-99=goto 3: s=s+2b:goto 2lb1 3: c:c=-99=goto 4: s=s+3c:goto 3:lb1 4:d:d=-99=goto 5: s=s+4d:goto: lb1 5:p=s/4ngoto 0 说明：输入方格个数n，依次输入角点a、边点b、拐点c、中点d每类点输入完后均输

16、入-99，显示平均高度p。 16:scqxjg:四次曲线加宽 klb：lb1 0:s :x=s-k:y=4(x/l)3-3(x/l)4:y=bygoto 0 说明：输入起点里程k、曲线长l、宽b，得输入里程加宽值y 17：scqxbz直线段四次曲线，边向中的法线 kxyflbdz：zd=d+b: lb1 0:s:a=s-k:n=b(4 (a/l)3-3(a/l)4):zn=-n: n=n+d:rec(a,f):p=x+iq=y+ju=x+acosf-nsinfv=y+asinf+ncosfw=b(12a2/l3-12a3/l4):zw=-w: g=tan-1w:c=wn+slb1 2:rec(

17、c-s,f)p=p+iq=q+jpol(p-u,q-v):jgoto 0 说明：输入起点里程k、坐标x、y、方位角f、曲线长l、加宽值b、宽度d，由窄到宽或宽到窄z（1）。输入所求点里程，显示中点坐标p、q，法线方位角j，其中输入b、d带符号，前进方向右正、右负。 18:scqxzb:直线段四次曲线，中线里程的 kxyfbdz:zd=d+b: lb1 0:s:a=s-k:q=a:lb1 1:c=b2(48q5/l6-84q6/l736q7/l8)+q-a:m=db(12q2/l3-12q3l4) e=b2(240q4/l6-504q5/l7+252q6/l8)+1:n=db(24q/l3-36

18、q2/l4):zc=-c: ze=-e: m=c+m:n=e+n:p=q-m/n:abs(q-p)goto 2: q=p:goto 1:lb1 2:rec(a,f):u=x+iv=y+jm=p+kn=b(4(p/l)3-3(p/l)4):zn=-n: n=n+dpol(m-s,n):j=j+f:jj=j+360: j360=j=j-360: c=j:rec(i,j):u=u+iv=v+jt=c+180:t360=t=t-360: tgoto 0 说明： 输入同15，显示里程s中点坐标u、v，法线交点里程m、交点至中线距离n、交点坐标u、v法线方位交t。 19：zxjdzblc直线交点坐标、里程 kabf:lbl 0:cdg:c=c:d=d:x=(atanf-ctang-b+d)/(tanf-tang) y=(x-c)tang+dl=k+(x-a)2+(y-b)2) goto 0 说明：输入直线起点里程k、坐标a、b、 方位角f。再输入所求直线点坐标c、d