def calc_focal_pts(self, lon0, lat0, a, b, theta):
# focal distance
f = np.sqrt(a**2 - b**2)
bearing = 90. - theta
# coordinates of focal points
lon1, lat1 = vincenty(lon0, lat0, bearing, f)
lon2, lat2 = vincenty(lon0, lat0, bearing + 180., f)
return lon1, lat1, lon2, lat2
def calc_focal_pts(self, lon0, lat0, a, b, theta):
# focal distance
f = np.sqrt(a**2 - b**2)
if isinstance( theta, u.Quantity):
bearing = 90. - (theta.to(u.degree)).value
else:
bearing = 90. - theta
# coordinates of focal points
lon1, lat1 = vincenty(lon0, lat0, bearing, f)
lon2, lat2 = vincenty(lon0, lat0, bearing + 180., f)
return lon1, lat1, lon2, lat2