Python获得两个GPS点轴承

Question

尝试获得两个gps坐标度，但得到错误的结果。 我认为问题是蟒蛇数学工作在光芒四射，但当转换到度数也有错误的结果。

好的结果是在30-33左右。 使用此代码，我得到0,75，当将rad转换为deg时，得到43。 有我的代码：

import math

lat1 = 47.357818 
lon1 = 18.985417
lat2 = 47.369282 
lon2 = 18.996625
pi = 3.1415

X = math.cos(lat2) * math.sin(lon1-lon2)
Y = math.cos(lat1) * math.sin(lat2) - math.sin(lat1) * math.cos(lat2) * math.cos(lon2-lon1)
B = math.atan2(X,Y)

print B

请帮帮我，转换的好方法是什么？ ps：抱歉我的英文。