李香清

(1.蘭州市勘察測(cè)繪研究院，甘肅蘭州 730030； 2.蘭州市地理信息中心，甘肅蘭州 730030)

基于對(duì)話框編寫實(shí)現(xiàn)坐標(biāo)方位角正算與反算的程序

李香清1，2?

(1.蘭州市勘察測(cè)繪研究院，甘肅蘭州 730030； 2.蘭州市地理信息中心，甘肅蘭州 730030)

通過Visual C++6.0來設(shè)計(jì)坐標(biāo)正反算的程序，基于對(duì)話框來編寫，特別是在控制網(wǎng)平差計(jì)算中，用來坐標(biāo)方位角的推算和平面坐標(biāo)的正反算，在導(dǎo)線測(cè)量中，通過坐標(biāo)的正反算來確?？刂凭W(wǎng)的精度。該程序可以廣泛應(yīng)用于施工放樣等相關(guān)工作中，并且用實(shí)際的測(cè)量工作驗(yàn)證了其準(zhǔn)確性和實(shí)用性。

坐標(biāo)正反算；Visual C++；坐標(biāo)方位角；導(dǎo)線測(cè)量

1 引 言

在實(shí)際測(cè)量工作時(shí)，比如在市政工程施工測(cè)量中，經(jīng)常會(huì)遇到根據(jù)已知導(dǎo)線控制點(diǎn)，利用經(jīng)緯儀、鋼尺測(cè)設(shè)待定點(diǎn)的實(shí)際問題，解決此類問題往往需要計(jì)算坐標(biāo)方位角或距離，如何能夠更快、更準(zhǔn)確的利用程序?qū)崿F(xiàn)坐標(biāo)方位角的正算與反算是一個(gè)亟待研究的問題。

2 相關(guān)概念

2.1 方位角

由直線一端的基本方向起，順時(shí)針方向至該直線的水平角度稱為該直線的方位角。方位角的取值范圍是0°～360°。下面介紹相關(guān)的三種方位角:

真方位角:由真北方向起算的方位角，用A表示。坐標(biāo)方位角:由坐標(biāo)北方向起算的方位角，用a表示。磁方位角:由磁北方向起算的方位角，用Am表示。

2.2 方位角之間的相互換算

由于三個(gè)指北的標(biāo)準(zhǔn)方向并不重合，所以一條直線的三種方位角并不相等，它們之間存在著一定的換算關(guān)系。如圖1所示，一條直線的真方位角A、磁北方位角Am、坐標(biāo)方位角a之間有如下的關(guān)系式，式中的δ為磁偏角，γ為子午線收斂角。

圖1 三種方位角的關(guān)系

關(guān)系表達(dá)式如下:

2.3 正、反坐標(biāo)方位角

一條直線的坐標(biāo)方位角，由于起始點(diǎn)的不同而存在著兩個(gè)值。如圖2所示，P1、P2為直線P1P2的兩端點(diǎn)，a12表示P1P2方向的坐標(biāo)方位角，a21表示P2P1方向的坐標(biāo)方位角。a12和a21互為正、反坐標(biāo)方位角。若以a12為正坐標(biāo)方位角，則稱a21為反坐標(biāo)方位角。由于在同一高斯平面直角坐標(biāo)系內(nèi)各點(diǎn)處坐標(biāo)北方向均是平行的，所以一條直線的正、反坐標(biāo)方位角相差180°，即a12=a21±180°。

圖2 正、反坐標(biāo)方位角 圖3 坐標(biāo)正、反算示意圖

3 平面直角坐標(biāo)正、反算

如圖3所示，設(shè)A為已知點(diǎn)，B為未知點(diǎn)，當(dāng)A點(diǎn)坐標(biāo)(XA，YA)、A點(diǎn)至B點(diǎn)的水平距離SAB和坐標(biāo)方位角aAB均為已知時(shí)，則可求得B點(diǎn)的坐標(biāo)(XB，YB)，通常稱為坐標(biāo)正算問題。

由圖可知坐標(biāo)正算的公式為:

所以上面的式子可以寫成:

式中的△XAB和△XAB稱之為坐標(biāo)增量。

同理可以得出坐標(biāo)反算的公式，設(shè)A、B兩已知點(diǎn)的坐標(biāo)分別為(XA，YA)和(XB，YB)，則這條直線AB的坐標(biāo)方位角aAB和水平距離SAB分別為:

應(yīng)當(dāng)指出，在上面的式子中△XAB和△YAB應(yīng)該取絕對(duì)值，計(jì)算得到的為象限角R，象限角取值范圍為0°～90°。而測(cè)量工作中通常用坐標(biāo)方位角表示直線的方向，因此，計(jì)算出象限角R后，應(yīng)該將其轉(zhuǎn)換為坐標(biāo)方位角，其轉(zhuǎn)換方法表1。

象限角轉(zhuǎn)換成坐標(biāo)方位角 表1

4 Visual C++6.0編程設(shè)計(jì)

4.1 對(duì)話框界面的設(shè)計(jì)

圖4 對(duì)話框的界面設(shè)計(jì)

4.2 坐標(biāo)正算程序的核心代碼

void Czbzs::OnButtonEqual2()//響應(yīng)坐標(biāo)正算按鈕

{

//TODO:Add your control notification handler code here

UpdateData(true)；//獲得輸入數(shù)據(jù)

i=m_44； //坐標(biāo)方位角獲取

//為減少運(yùn)算，先計(jì)算幾種特殊情況

if(i==0｜｜i==360)

{

m_55=m_11+m_33；

m_66=m_22；

｝

else if(i==180)

{

m_55=m_11-m_33；

m_66=m_22；

｝

else if(i==90)

{

m_55=m_11；

m_66=m_22+m_33；

｝

else if(i==270)

{

m_55=m_11；

m_66=m_22-m_33；

｝

else if(i＞0＆＆i＜90)//第一象限角

{

j=(m_44?PI)/180；//角度轉(zhuǎn)換為弧度

m_55=m_11+m_33?cos(j)；//由表一和公式8、公式9綜合

m_66=m_22+m_33?sin(j)；

｝

else if(i＞90＆＆i＜180)//第二象限角

{

j=((m_44-90)?PI)/180；

m_55=m_11+m_33?(-cos(j))；

m_66=m_22+m_33?sin(j)；

｝

else if(i＞180＆＆i＜270)//第三象限角

{

j=((m_44-180)?PI)/180；

m_55=m_11+m_33?(-cos(j))；

m_66=m_22+m_33?(-sin(j))；

｝

else if(i＞270＆＆i＜360)//第四象限角

{

j=((m_44-270)?PI)/180；

m_55=m_11+m_33?cos(j)；

m_66=m_22+m_33?(-sin(j))；

｝

UpdateData(false)；//更新顯示

｝

4.3 坐標(biāo)反算程序的核心代碼

void Czbfs::OnButtonEqual()//響應(yīng)坐標(biāo)反算按鈕

{

//TODO:Add your control notification handler code here

UpdateData(true)；//獲得輸入數(shù)據(jù)

dx=m_3-m_1；

dy=m_4-m_2；

az1=dxdy_az(dx，dy)；//計(jì)算坐標(biāo)方位角

az=radio_to_angle(az1)；//弧度轉(zhuǎn)換為角度

s=dist(dx，dy)；//計(jì)算距離

m_5=m_9；

m_6=m_10；

m_7=az1；

m_8=s；

UpdateData(false)；//更新顯示

｝

double Czbfs::dxdy_az(double dx，double dy)//坐標(biāo)方位角計(jì)算

{

double R，az；

R=atan(dy/dx)；//公式10

if(dx==0)//特殊情況

{

if(dy＞0)

az=PI/2；

else

az=3?PI/2；

｝

else if(dx＞0)//由表一得出az與R關(guān)系

{

if(dy＞0)

az=R；

else

az=R+2?PI；

｝

else

az=PI+R；

return(az)；

｝

double Czbfs::radio_to_angle(double radio)//弧度轉(zhuǎn)換為角度

{

double a1，a2，angle；

int d，f，m；

a1=radio?180/PI；

d=int(a1)；//計(jì)算度

a2=(a1-d)?60；

f=int(a2)；//計(jì)算分

m=int((a2-f)?60)；//計(jì)算秒angle=d+f/100.0+m/10000.0；

return(angle)；

｝

double Czbfs::dist(double dx，double dy)//計(jì)算距離

{

return(sqrt(dx?dx+dy?dy))；

｝

5 運(yùn)行結(jié)果

將外業(yè)實(shí)測(cè)數(shù)據(jù)進(jìn)行坐標(biāo)方位角的推算和平面坐標(biāo)的正反算，通過坐標(biāo)的正反算來確?？刂凭W(wǎng)的精度。在坐標(biāo)反算主界面中輸入起點(diǎn):1點(diǎn)坐標(biāo)(4053345.321，486456.321)，終點(diǎn)2點(diǎn)坐標(biāo)(4053323.678，486123.443)，通過計(jì)算得出坐標(biāo)方位角(弧度)為4.64746249773597，1點(diǎn)到2點(diǎn)的距離為333.580848270678；在坐標(biāo)正算主界面中輸入起點(diǎn):1點(diǎn)坐標(biāo)(4053321.476，485123.443)，距離334.756，坐標(biāo)方位角(角度)213.342，通過計(jì)算得出終點(diǎn)坐標(biāo)為(4053041.81926841，484939.449271809)。計(jì)算結(jié)果如圖5所示:

圖5 坐標(biāo)正反算的驗(yàn)證結(jié)果

6 結(jié) 論

本文通過對(duì)坐標(biāo)正反算的驗(yàn)證得出該程序的可靠性與實(shí)用性，基于對(duì)話框編寫的該程序，界面直觀簡(jiǎn)潔，使計(jì)算方位角的工作變得簡(jiǎn)單易行，實(shí)現(xiàn)坐標(biāo)方位角的正反算更加方便快捷，可以廣泛應(yīng)用于工程測(cè)量的施工放樣及相關(guān)工作中，大幅降低了測(cè)繪工作量。

[1] 潘正風(fēng)，楊正堯，程效軍等.數(shù)字測(cè)圖原理與方法[M].武漢大學(xué)出版社，2004.

[2] 王育堅(jiān).Visual C++面向?qū)ο缶幊探坛蘙M].清華大學(xué)出版社，2007.

[3] 沈晶，劉海波，周長(zhǎng)建.Visual C++數(shù)字圖像處理典型案例詳解[M].機(jī)械工業(yè)出版社，2012.

[4] 武漢大學(xué)測(cè)繪學(xué)院測(cè)量平差學(xué)科組.誤差理論與測(cè)量平差基礎(chǔ)[M].武漢大學(xué)出版社，2003.

[5] 陳天華.面向?qū)ο蟪绦蛟O(shè)計(jì)與Visual C++6.0教程題解與實(shí)驗(yàn)指導(dǎo)[M].清華大學(xué)出版社，2006.

[6] 張正祿.工程測(cè)量學(xué)[M].武漢大學(xué)出版社，2005.

Based on the Dialog Box to W rite the Program to Realize Direct and Inverse Calculation of Geodetic Coordinates and Azimuth

Li Xiangqing1，2
(1.Lanzhou Istitute of Survey and Investigation，Lanzhou 730030，China；2.Lanzhou Geographic Information Center Gansu，Lanzhou 730030，China)

Through the Visual C++6 to design the inverse coordinate calculation program，based on the dialog box to write，especially the calculation in the control network adjustment，and to calculate the coordinate azimuth angle and plane coordinates calculation，in a traverse，through coordinate and ensure the precision of the control network.This program iswidely used in construction lofting and related work and verified its accuracy and reliability by using examples.

direct and inverse calculation of coordinates；visual C++；coordinate azimuth；traverse survey

1672-8262(2013)05-137-03

P209

B

2013—06—28

李香清(1978—)，女，工程師，現(xiàn)從事地理信息空間數(shù)據(jù)質(zhì)量檢查、GIS應(yīng)用研究及開發(fā)等。