python 計算方位角實(shí)例(根據(jù)兩點(diǎn)的坐標(biāo)計算)

更新時間：2020年01月17日 11:48:24 作者：JohnieLi

今天小編就為大家分享一篇python 計算方位角實(shí)例(根據(jù)兩點(diǎn)的坐標(biāo)計算)，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧

知道兩點(diǎn)坐標(biāo)，怎么計算兩點(diǎn)方向的方位角？

答：首先計算坐標(biāo)增量dx，dy(兩個對應(yīng)坐標(biāo)分量相減，終點(diǎn)的減始點(diǎn)的)。

若dx,dy中有一個為零時，根據(jù)另一個的正負(fù)決定方位角（0，90，180，270這四個中的一個，可畫坐標(biāo)軸圖分析，但不要畫為數(shù)學(xué)坐標(biāo)哦）。

基本思路：

若dx,dy都不為零；則

計算a=arcatn(|dy/dx|)（這好像叫象限角）

當(dāng)dx>0dy>0時方位角=a;

當(dāng)dx<0dy>0時方位角=180-a;

當(dāng)dx<0dy<0時方位角=180+a; 負(fù)范圍為a-pi

當(dāng)dx>0dy<0時方位角=360-a; 負(fù)范圍為-a

還有一種方法，使用 atan2來計算方位角，范圍為-pi,pi

atan2(y,x)所表達(dá)的意思是坐標(biāo)原點(diǎn)為起點(diǎn)，指向(x,y)的射線在坐標(biāo)平面上與x軸正方向之間的角的角度。

結(jié)果為正表示從 X 軸逆時針旋轉(zhuǎn)的角度，結(jié)果為負(fù)表示從 X 軸順時針旋轉(zhuǎn)的角度。

atan 和 atan2 都是求反正切函數(shù)，如：有兩個點(diǎn) point(x1,y1), 和 point(x2,y2);

那么這兩個點(diǎn)形成的斜率的角度計算方法分別是：

float angle = atan( (y2-y1)/(x2-x1) );
float angle = atan2( y2-y1, x2-x1 );

atan 和 atan2 區(qū)別：

1：參數(shù)的填寫方式不同；

2：atan2 的優(yōu)點(diǎn)在于如果 x2-x1等于0 依然可以計算，但是atan函數(shù)就會導(dǎo)致程序出錯；

3：atan2(a,b)的取值范圍介于 -pi 到 pi 之間（不包括 -pi），而atan(a/b)的取值范圍介于-pi/2到pi/2之間（不包括±pi/2)。

另外要注意的是，函數(shù)atan2(y,x)中參數(shù)的順序是倒置的，atan2(y,x)計算的值相當(dāng)于點(diǎn)(x,y)的角度值。

atan2(y, x)是4象限反正切，它的取值不僅取決于正切值y/x，還取決于點(diǎn) (x, y) 落入哪個象限：

當(dāng)點(diǎn)(x, y) 落入第一象限時，atan2(y, x)的范圍是 0 ~ pi/2;
當(dāng)點(diǎn)(x, y) 落入第二象限時，atan2(y, x)的范圍是 pi/2 ~ pi;
當(dāng)點(diǎn)(x, y) 落入第三象限時，atan2(y, x)的范圍是 －pi～－pi/2;
當(dāng)點(diǎn)(x, y) 落入第四象限時，atan2(y, x)的范圍是 -pi/2～0.


而 atan(y/x) 僅僅根據(jù)正切值為y/x求出對應(yīng)的角度 （可以看作僅僅是2象限反正切）：

當(dāng) y/x > 0 時，atan(y/x)取值范圍是 0 ~ pi/2；
當(dāng) y/x < 0 時，atan(y/x)取值范圍是 -pi/2～0.

如果要實(shí)現(xiàn)方位角的計算，代碼如下：

# 計算方位角函數(shù)
def azimuthAngle( x1, y1, x2, y2):
  angle = 0.0;
  dx = x2 - x1
  dy = y2 - y1
  if x2 == x1:
    angle = math.pi / 2.0
    if y2 == y1 :
      angle = 0.0
    elif y2 < y1 :
      angle = 3.0 * math.pi / 2.0
  elif x2 > x1 and y2 > y1:
    angle = math.atan(dx / dy)
  elif x2 > x1 and y2 < y1 :
    angle = math.pi / 2 + math.atan(-dy / dx)
  elif x2 < x1 and y2 < y1 :
    angle = math.pi + math.atan(dx / dy)
  elif x2 < x1 and y2 > y1 :
    angle = 3.0 * math.pi / 2.0 + math.atan(dy / -dx)
  return (angle * 180 / math.pi)

math中關(guān)于三角函數(shù)常用的操作：

import math
math.acos(x)  # 返回 x 的反余弦 弧度值。  
math.asin(x)  # 返回 x 的反正弦 弧度值。  
math.degrees(x)  # 將 弧度 轉(zhuǎn)換為 角度, 如 degrees(math.pi/2) ， 返回90.0  
math.radians(x)  # 將 角度 轉(zhuǎn)換為 弧度
注意負(fù)數(shù)角度的轉(zhuǎn)換。

以上這篇python 計算方位角實(shí)例(根據(jù)兩點(diǎn)的坐標(biāo)計算)就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。

您可能感興趣的文章: