VB6.0物探線平面坐標(biāo)推算程序設(shè)計(jì)

王寶江

1 問題的提出

地質(zhì)部門的測量任務(wù)中有許多是關(guān)于物探線的測量定線問題,按測量學(xué)分類屬于直線的測設(shè)問題。該類問題放樣點(diǎn)、線數(shù)據(jù)量很大,如何提取每條線的平面放樣數(shù)據(jù)就成為擺在地質(zhì)測量人員面前急需解決的問題。

現(xiàn)在地質(zhì)及測量部門經(jīng)常使用的成圖軟件有MAPGIS和加載在CAD基礎(chǔ)上的南方CASS成圖系統(tǒng),都不具備按線的方向和間距提取大量坐標(biāo)數(shù)據(jù)的功能,為此本人用VB6.0軟件編制了一個(gè)推算物探線平面坐標(biāo)數(shù)據(jù)的應(yīng)用程序,推算后生成的數(shù)據(jù)文件為南方CASS數(shù)據(jù)的.DAT格式,可以在南方CASS成圖軟件中直接展繪。

2 處理思路

物探線一般要求5 m,10 m的點(diǎn)距,長度可達(dá)幾千米,假定一條1 km的物探線點(diǎn)距5 m,則要測設(shè)201個(gè)點(diǎn)位。線距一般有50 m,100 m不等。而直線由起點(diǎn)和終點(diǎn)唯一確定直線的方向和距離。有了線的起點(diǎn)坐標(biāo)和方位角就可以按下式計(jì)算間隔點(diǎn)的坐標(biāo):間隔點(diǎn)的 X(或 Y)坐標(biāo)=起始點(diǎn) X(或 Y)坐標(biāo)+X(或 Y)方向的坐標(biāo)增量×(間隔點(diǎn)號-起始點(diǎn)號);

現(xiàn)在的關(guān)鍵是按點(diǎn)距D和線的方位角α要求算出X方向和Y方向的坐標(biāo)增量ΔX,ΔY。

ΔX=D ×cos(α);ΔY=D ×sin(α)。

注意:數(shù)學(xué)坐標(biāo)格式為:X,Y與南方CASS格式中坐標(biāo)格式相同,但在實(shí)際應(yīng)用中的測量坐標(biāo)格式為北N,東E格式。

程序結(jié)構(gòu)示意圖見圖1。

3 具體實(shí)現(xiàn)

程序主界面見圖2。

點(diǎn)擊以上復(fù)選框,彈出如圖3所示表單。

以下為線上點(diǎn)坐標(biāo)增量計(jì)算程序的主要代碼:

Private Sub Command1-Click()

Dim c

Dim cc

Dem“以下為將測量的方位角:度、分、秒(北東格式)裝換為VB6.0要求的弧度格式。”

cc=(Int(fwj.Text)*3600+Int((fwj.Text-Int(fwj.Text))*100)*60+(fwj.Text*100-Int(fwj.Text*100))*100)/206264.81

Dem“以下為測量方位角轉(zhuǎn)換為數(shù)學(xué)方位角。”

c=3.14159265358979/2-cc

Dem“進(jìn)行坐標(biāo)增量計(jì)算,如Xzl對應(yīng)于線上點(diǎn)的X增量。”

Form1.Xzl.Text=Cos(c)*Form2.djg

Form1.Yzl.Text=Sin(c)*Form2.djg

Form1.Xjg.Text=Cos((c+3.14159265358979/2))*Form2.xjg

Form1.Yjg.Text=Sin((c+3.14159265358979/2))*Form2.xjg

Form1.jszl.Value=0

End Sub

主表單的程序相對簡單,故不贅述。

輸入全部數(shù)據(jù)的示例見圖4。

以上示例中線號選擇了按偶數(shù)增加。

輸入必要的數(shù)據(jù)后先點(diǎn)擊“推算坐標(biāo)”再點(diǎn)擊“顯示結(jié)果”。本軟件用Windows操作系統(tǒng)下普遍使用的文本編輯器NOTPAD.exe打開一個(gè)名稱為Temp.dat的數(shù)據(jù)文件。

Temp.dat會在下一次推算數(shù)據(jù)時(shí)被刷新,故推算數(shù)據(jù)在確信正確無誤后,應(yīng)及時(shí)另存且不可與Temp.dat重名。

4 推算數(shù)據(jù)的檢核

本軟件默認(rèn)安裝于D:chy下。專用于物探線平面坐標(biāo)的推算,推算后生成的數(shù)據(jù)文件為南方CASS數(shù)據(jù)的.DAT格式,可以在南方CASS成圖軟件中直接展繪。

現(xiàn)取推算后數(shù)據(jù)的展繪結(jié)果,用CASS的距離量測工具檢核:

1)兩點(diǎn)間實(shí)地距離=495.000 m,圖上距離=495.000 mm,方位角=30°10′29.86″。2)兩點(diǎn)間實(shí)地距離 =50.000 m,圖上距離=50.000 mm,方位角 =300°10′28.90″。3)兩 點(diǎn)間實(shí)地距 離 =200.000 m,圖上距離=200.000 mm,方位角=300°10′30.31″。

經(jīng)過實(shí)例檢測:推算精度距離誤差在2 000 m內(nèi)遠(yuǎn)低于毫米級,角度誤差低于秒級。實(shí)際使用中由于舍入到毫米即可滿足需要,故推算數(shù)據(jù)保留到毫米,這使得在CAD中的查詢方向有微小的誤差,實(shí)際應(yīng)用中低于秒級可忽略不計(jì)。

[1] 姬根杰.Excel VBA編程在導(dǎo)線測量中的應(yīng)用[J].山西建筑,2008,34(1):367-368.