Учитывая широту и долготу, а также расстояние, я хочу найти ограничивающий прямоугольник
-
22-07-2019 - |
Вопрос
Учитывая широту, долготу и расстояние, я хочу найти ограничивающую рамку, где расстояния меньше заданного расстояния.
Этот вопрос был задан здесь: Как рассчитать ограничивающий прямоугольник для заданного местоположения широты и ширины?
Я не хочу, чтобы это было частично точным, поэтому я изменил и упростил его до
def boundingBox(latitudeInDegrees, longitudeInDegrees, halfSideInKm):
lat = math.radians(latitudeInDegrees)
lon = math.radians(longitudeInDegrees)
halfSide = 1000*halfSideInKm
RADIUS_OF_EARTH = 6371
# Radius of the parallel at given latitude
pradius = radius*math.cos(lat)
latMin = lat - halfSide/radius
latMax = lat + halfSide/radius
lonMin = lon - halfSide/pradius
lonMax = lon + halfSide/pradius
rad2deg = math.degrees
return (rad2deg(latMin), rad2deg(lonMin), rad2deg(latMax), rad2deg(lonMax))
Но я не могу понять, как это работает, в частности, эта строка не имеет смысла для меня halfSide = 1000 * halfSideInKm
Решение
Эта линия переводит ограничивающие рамки из километров в метры.
Другие советы
Этот код не совсем работает, он перемещается между KM и M.
Исправлен код, добавлены имена в стиле PEP8 и добавлен простой объект box:
class BoundingBox(object):
def __init__(self, *args, **kwargs):
self.lat_min = None
self.lon_min = None
self.lat_max = None
self.lon_max = None
def get_bounding_box(latitude_in_degrees, longitude_in_degrees, half_side_in_miles):
assert half_side_in_miles > 0
assert latitude_in_degrees >= -90.0 and latitude_in_degrees <= 90.0
assert longitude_in_degrees >= -180.0 and longitude_in_degrees <= 180.0
half_side_in_km = half_side_in_miles * 1.609344
lat = math.radians(latitude_in_degrees)
lon = math.radians(longitude_in_degrees)
radius = 6371
# Radius of the parallel at given latitude
parallel_radius = radius*math.cos(lat)
lat_min = lat - half_side_in_km/radius
lat_max = lat + half_side_in_km/radius
lon_min = lon - half_side_in_km/parallel_radius
lon_max = lon + half_side_in_km/parallel_radius
rad2deg = math.degrees
box = BoundingBox()
box.lat_min = rad2deg(lat_min)
box.lon_min = rad2deg(lon_min)
box.lat_max = rad2deg(lat_max)
box.lon_max = rad2deg(lon_max)
return (box)
Не связан с StackOverflow