def get_bounding_box(self, maxdist):
"""
Bounding box containing all points, enlarged by the maximum distance
and the maximum rupture projection radius (upper limit).
"""
lon = self.location.x
lat = self.location.y
maxradius = self._get_max_rupture_projection_radius()
a1 = (maxdist + maxradius) * KM_TO_DEGREES
a2 = angular_distance(maxdist + maxradius, lat)
return lon - a2, lat - a1, lon + a2, lat + a1
def test(self):
aae(angular_distance(km=1000, lat=80), 51.7897747)
aae(angular_distance(km=1000, lat=88), 257.68853)
def test(self):
aae(angular_distance(km=1000, lat=80), 51.7897747)
aae(angular_distance(km=1000, lat=88), 257.68853)
