坐標(biāo)方位角計算公式

坐標(biāo)方位角計算公式(通用)用極坐標(biāo)法放樣必須計算出測站點(儀器點)到放樣點得距離和方位角，才能進(jìn)行放樣。原計算公式為：S12=sqr((x2-x1)2+(y2-y1)2)=sqr(△x221+△y221)A12=arcsin((y2-y1)/S12)S12為測站點１至放樣點２的距離；A12為測站點１至放樣點２的坐標(biāo)方位角。x1，y1為測站點坐標(biāo)；x2，y2為放樣點坐標(biāo)。按公式A12=arcsin((y2-y1)/S12)計算出的方位角都要進(jìn)行象限判斷后加常數(shù)才是真正的方位角。新計算公式為：A12=arccos(△x21/S12)*sgn(△y21)+360°式中sgn()為取符號函數(shù)，改公式只需加上條件（A12>360°,A12=A12-360°）就可以計算出坐標(biāo)方位角，不需要進(jìn)行象限判斷。我的這個公式要更好一些，計算結(jié)果就是正確結(jié)果：

SGN是正負(fù)號的函數(shù)。括號內(nèi)的數(shù)字大于零SGN（）就是+號，反之就是-號。===================================函數(shù)開始==================================='jiaodu10(x,splitStr)函數(shù)將60進(jìn)制度轉(zhuǎn)換為10進(jìn)制度格式.x為度數(shù)，splitStr為分隔符號，'如x為43%67%367，則splitStr為"%",參數(shù)要用雙引號括起來，jiaodu10("x","%")Functionjiaodu10(x,splitStr)IfInStr(1,x,splitStr)ThenDimss=Split(x,splitStr)jiaodu10=s(0)+s(1)/60+s(2)/3600Elsejiaodu10="錯誤"EndIfEndFunction'--------------------------------------------------------------------------------'jiaodu60(x,splitStr)函數(shù)將10進(jìn)制度轉(zhuǎn)換為60進(jìn)制度格式，splitStr分隔表示'x為數(shù)字，可以不用雙引號括起來，參數(shù)splitStr要用雙引號括起來iaodu10(12.31313,"-")Functionjiaodu60(x,splitStr)Dimfen,miaoFen=Round((fen-Int(fen))*60,0)Ifmiao>=60Thenmiao=miao-60fen=fen+1EndIfjiaodu60=Int(x)&splitStr&Int(fen)&splitStr&miaoEndFunction'--------------------------------------------------------------------------------'juli(待算點縱坐標(biāo)x，待算點橫坐標(biāo)y，測站點縱坐標(biāo)m，測站點縱坐標(biāo)n)用于計算距離。Functionjuli(x,y,m,n)juli=Math.Spr((x-m)^2+(y-n)^2)EndFunction'--------------------------------------------------------------------------------'jiaodu(x,y,m,n)計算角度Functionjiaodu(x,y,m,n)Dimdx,dy,a,jdu10dx=x-mdy=y-ma=Math.Abs(Math.Atn(dy/dx)*180/3.14159265)jdu10=0If(dx>0)ThenIf(dy>0)Thenjdu10=aElsejdu10=360-aEndIfElseIf(dy>0)Thenjdu10=180-aElsejdu10=180+aEndIfEndIfjiaodu=jiaodu60(jdu