Пеленг между двумя точками

Я использовал пакет geopy, который отлично справляется, однако некоторые из полученных мной результатов несовместимы или имеют относительно большое смещение, я подозреваю, что проблема заключается в с моим расчетом пеленга:

def gb(x,y,center_x,center_y):
dx=x-center_x
dy=y-center_y
if ((dy>=0)and((dx>0)or(dx<0))):
    return math.degrees(math.atan2(dy,dx))
elif (dy<=0)and((dx>0)or (dx<0)):
    return (math.degrees(math.atan2(dy,dx))+360)
else:
    return (math.degrees(math.atan2(dy,dx))+360)%360

Мне нужно вычислить пеленг, st center_x и center_y являются точкой поворота. затем я использую геопию для обратного проектирования координат GPS:

latlon = VincentyDistance(miles=dist).destination(Point(lat1, lon1), bearing)

Кто-нибудь может указать мне, что я делаю не так?