bearing static method
Initial bearing from the first point to the second, degrees clockwise from north.
Implementation
static double bearing(double lat1, double lng1, double lat2, double lng2) {
final double y = sin(_rad(lng2 - lng1)) * cos(_rad(lat2));
final double x = cos(_rad(lat1)) * sin(_rad(lat2)) - sin(_rad(lat1)) * cos(_rad(lat2)) * cos(_rad(lng2 - lng1));
return (atan2(y, x) * 180 / pi + 360) % 360;
}