Я использую метод "buffer" из API GEOS для создания круга, основанного на точке и радиусе в GeoDjango, в соответствии с этим ответом: GeoDjango: как создать круг, основанный на точке и радиусе
Как заметил один из комментаторов, вам нужно будет преобразовать радиус в км в градусы... но как? Я предполагаю, что это быстрая формула, но она греческая для меня. Я был бы благодарен за любые указатели.
Ответ может зависеть от местоположения на земле. Около экватора 1 км будет равняться примерно 0,008 градусам (1 км /40 000 км * 360 градусов) широты и долготы, но около полюсов 1 км будет равняться примерно 0,008 градусам по широте, но может быть много много градусов долготы. (Если вы на расстоянии 1 км от полюса, расстояние в 1 км "запад" приблизит вас к 57 градусам долготы на запад.)
Но, если API просто хочет получить градусы по большому кругу в качестве измерения, возможно, этого будет достаточно (n км /40 000 км * 360 градусов). По крайней мере, 40 000 км "достаточно для меня":) более точные цифры доступны.
1 градус лата составляет около 110,567 км на экваторе до 111,699 км на полюсах.
Этот link может вас заинтересовать. Преобразование основано на окружности деления земли на 360 градусов.
1 град. длиной также около 111 км на экваторе и → 0, когда вы идете на полюса.
Буфер в API-интерфейсе GEOS django создаст буфер, используя те единицы, которые использует ваша текущая система координат.
Если вы сохраняете все в 4326 (в латах/длинных градусах), вам придется найти какой-то сложный способ преобразования KM в градусы. Но теперь ваш буфер будет сильно искажен, чем больше север вы идете.
Лучшим решением является перепроектирование вашей геометрии в проекцию, которая поддерживает область, и часто такая проекция может отслеживать единицы измерения в метрах.
Например, если вы создаете буферные области в Северной Америке, вы можете использовать эту проекцию, которая использует метры http://spatialreference.org/ref/sr-org/7314/
Вот пример того, как это сделать с использованием Django GEOS API:
from django.contrib.gis.geos import Point
# Defines a point in lat/long
p = Point(-70, 50)
# This projection defines lat/long coordinate system
p.srid = 4326
# Transform into the 7314 projection using the OGC WKT format to define that projection
p.transform('PROJCS["NA Lambert Azimuthal Equal Area",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["false_easting",0.0],PARAMETER["false_northing",0.0],PARAMETER["longitude_of_center",-100.0],PARAMETER["latitude_of_center",45.0],UNIT["meter",1.0]]')
# Creates a buffered polygon of 1000 meters in radius
poly = p.buffer(1000)
расстояние = степень * диаметр PI */360
Используя приблизительную оценку радиуса Земли: 6378 км
1 Степень @на экваторе или вдоль долготы: 111.317 км
1 Степень вдоль широты в градуированном градуите: (градусы * PI * диаметр /360) * cos (широта)
Итак, 60 широт: 55 659 м в направлении долготы.
Я использую этот конвертер
Формулы для расчета градусов в километрах: L = π * R * a/180