德州仪器TI nspire CAS工程测量程序

版权声明:本文为博主原创文章,遵循 CC 4.0 by-sa 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/hjpqindong/article/details/95944292

1、高斯坐标转工程坐标系:

Define a_gauss2gc()=

Prgm

setgeomode()

RequestStr "施工坐标系名=",tconame

0→stat

For i,1,dim(coname),1

 If coname[i]=tconame Then

 1→stat

 Request "待转换点X(N)=",tpx

 Request "待转换点Y(E)=",tpy

 xa[i]→txa

 ya[i]→tya

 xb[i]→txb

 yb[i]→tyb

 sa[i]→tsa

 R►Pθ(txb-txa,tyb-tya)→talpha

 EndIf

EndFor

If stat=1 Then

(tpx-txa).cos(talpha)+(tpy-tya).sin(talpha)+tsa→tstage

-(tpx-txa).sin(talpha)+(tpy-tya).cos(talpha)→toff

Disp "高斯坐标→施工坐标:",tconame

Disp "(桩号)x=",round(tstage,3)

Disp "(偏距)y=",round(roff,3)

Else

Disp "选择的坐标系未定义!"

EndIf

EndPrgm

 

 

2、工程坐标转高斯坐标系:

Define a_gc2gauss()=

Prgm

Setgeomode()

RequestStr  “施工坐标系名=”,tconame

0→stat

For i,1,dim(coname),1

 If  coname[i]=tconame Then

 1→stat

 Request “待转换点桩号x=”,tpx

 Request “待转换点偏距y=”,tpy

 xa[i]→txa

 ya[i]→tya

xb[i]→txb

yb[i]→tyb

sa[i]→tsa

R►Pθ (txb-txa,tyb-tya)→talpha

EndIf

EndFor

If stat=1 Then

txa+(tpx-tsa).cos(talpha)-tpy.sim(talpha)→tpn

tya+(tpx-tsa).sin(talpha)+tpy.cos(talpha)→tpe

Disp “施工坐标→高斯坐标 坐标系:”,tconame

Disp “X(N)=”,round(tpn,3)

Disp “Y(E)=”,round(tpe,3)

Else

Disp “选择的坐标系未定义!”

EndIf

EndPrgm

 

3、setgeomode子程序

Define setgeomode()=

Prgm

If getMode(2)≠2 Then

setMode(1,1)

setMode(2,2)

setMode(5,2)

setMode(7,1)

setMode(8,1)

EndIf

EndPrgm

 

4、工程坐标系转工程坐标系

Define a_gc2gc()=

Prgm

RequestStr “源施工坐标系名=”,tconame

0→stat

For i,1,dim(coname),1

 If coname[i]=tconame Then

 1→stat

 Request “待转换点桩号x=”,tpx

 Request “待转换点偏距y=”,tpy

 xa[i]→txa

 ya[i]→tya

 xb[i]→txb

 yb[i]→tyb

 sa[i]→tsa

 R►Pθ(txb-txa,tyb-tya)→talpha

 EndIf

EndFor

If stat=1 Then

txa+(tpx-tsa).cos(talpha)-tpy.sin(talpha)→tpn

tya+(tpx-tsa).sin(talpha)+tpy.cos(talpha)→tpe

RequestStr “目标坐标系名=”,tarconame

0→stat2

For i,1,dim(coname),1

 If coname[i]=tarconame Then

 1→stat2

 xa[i]→txa

 ya[i]→tya

 xb[i]→txb

 yb[i]→tyb

 sa[i]→tsa

 R►Pθ(txb-txa,tyb-tya)→talpha

 EndIf

EndFor

If stat2=1 Then

(tpn-txa).cos(talpha)+(tpe-tya).sin(talpha)+tsa→tstage

-(tpn-txa).sin(talpha)+(tpe-tya).cos(talpha)→toff

Disp “施工坐标系间转换:”,tconame & “→”& tarconame

Disp “(桩号)x=”,round(tstage,3)

Disp “(偏距) y=”,round(roff,3)

Else

Disp “选择的目标坐标系未定义!”

EndIf

Else

Disp “选择的源坐标系未定义!”

EndIf

EndPrgm

 

5、坐标反算

Define a_zbfs()=

Prgm

setgeomode()

Request “起点X=”,x1

Request “起点Y=”,y1

Request “终点X=”,x2

Request “终点Y=”,y2

Disp “距离=”,R►Pr(x2-x1,y2-y1)

R►Pθ (x2-x1,y2-y1)→α

If α<0 Then

α+360→α

EndIf

Disp “方位角=”,α►DMS

EndPrgm

 

6、坐标正算

Define a_zbzs()=

Prgm

setgeomode()

Request “起点X=”,x1

Request “起点Y=”,y1

Request “距离=”,s

RequestStr “方位角(D.MMSS)=”,str

expr(str)→α

fdms(α)→α

Disp “计算点X=”,round(x1+s.cos(α),3)

Disp “计算点Y=”,round(y1+s.sin(α),3)

EndPrgm

 

7、转换D.MMSS格式为十进制度

Define fdms(α)=

Func

Return intDiv(α,1)+(intDiv(remain(α,1).100,1))/60+(remain(α*100,1).100)/3600

EndFunc

 

 

展开阅读全文

没有更多推荐了,返回首页