Функции для работы с географическими координатами

greatCircleDistance(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — долгота первой точки в градусах. Диапазон: [-180°, 180°].
• lat1Deg — широта первой точки в градусах. Диапазон: [-90°, 90°].
• lon2Deg — долгота второй точки в градусах. Диапазон: [-180°, 180°].
• lat2Deg — широта второй точки в градусах. Диапазон: [-90°, 90°].

SELECT greatCircleDistance(55.755831, 37.617673, -55.755831, -37.617673) AS greatCircleDistance

┌─greatCircleDistance─┐
│            14128352 │
└─────────────────────┘

geoDistance(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — долгота первой точки в градусах. Диапазон: [-180°, 180°].
• lat1Deg — широта первой точки в градусах. Диапазон: [-90°, 90°].
• lon2Deg — долгота второй точки в градусах. Диапазон: [-180°, 180°].
• lat2Deg — широта второй точки в градусах. Диапазон: [-90°, 90°].

SELECT geoDistance(38.8976, -77.0366, 39.9496, -75.1503) AS geoDistance

┌─geoDistance─┐
│   212458.73 │
└─────────────┘

greatCircleAngle(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — долгота первой точки в градусах.
• lat1Deg — широта первой точки в градусах.
• lon2Deg — долгота второй точки в градусах.
• lat2Deg — широта второй точки в градусах.

SELECT greatCircleAngle(0, 0, 45, 0) AS arc

┌─arc─┐
│  45 │
└─────┘

geoToUTM(longitude, latitude[, zone])

• longitude — Долгота в градусах. Диапазон: [-180°, 180°]. Float32/Float64.
• latitude — Широта в градусах. Диапазон: [-80°, 84°]. Float32/Float64.
• zone — Необязательно. Принудительно использовать эту зону UTM для проекции вместо автоматического выбора. Диапазон: [1, 60]. (U)Int*.

SELECT geoToUTM(2.294497, 48.858222) AS utm; -- Eiffel Tower

(448251.5978370684,5411935.125629659,31,'U')

UTMToGeo(easting, northing, zone, is_north)

• easting — восточное смещение в метрах (включая ложное восточное смещение 500000 м). (U)Int*/Float*.
• northing — северное смещение в метрах (включая ложное северное смещение 10000000 м для южного полушария). (U)Int*/Float*.
• zone — номер зоны UTM. Диапазон: [1, 60]. (U)Int*.
• is_north — полушарие: 1 для северного полушария, 0 для южного. (U)Int*.

SELECT UTMToGeo(448251.6, 5411935.13, 31, 1) AS coord;

(2.2944970289079203,48.85822204127082)

geoToMGRS(longitude, latitude[, precision])

• longitude — Долгота в градусах. Диапазон: [-180°, 180°]. Float32/Float64.
• latitude — Широта в градусах. Диапазон: [-80°, 84°]. Float32/Float64.
• precision — Необязательно. Количество цифр для значений восточного и северного смещения. По умолчанию: 5. Диапазон: [0, 5]. (U)Int*.

SELECT geoToMGRS(2.294497, 48.858222) AS mgrs, geoToMGRS(2.294497, 48.858222, 3) AS mgrs_100m;

┌─mgrs────────────┬─mgrs_100m───┐
│ 31UDQ4825111935 │ 31UDQ482119 │
└─────────────────┴─────────────┘

MGRSToGeo(mgrs)

• mgrs — строка координат MGRS для декодирования. String/FixedString.

SELECT MGRSToGeo('31UDQ4825111935') AS coord;

(2.294495618908297,48.85822536113692)

pointInEllipses(x, y, x₀, y₀, a₀, b₀,...,xₙ, yₙ, aₙ, bₙ)

• x, y — координаты точки на плоскости.
• xᵢ, yᵢ — координаты центра i-го эллипса.
• aᵢ, bᵢ — оси i-го эллипса в единицах координат x, y.

SELECT pointInEllipses(10., 10., 10., 9.1, 1., 0.9999)

┌─pointInEllipses(10., 10., 10., 9.1, 1., 0.9999)─┐
│                                               1 │
└─────────────────────────────────────────────────┘

pointInPolygon((x, y), [(a, b), (c, d) ...], ...)

• (x, y) — Координаты точки на плоскости. Тип данных — Tuple — кортеж из двух чисел или Point.
• [(a, b), (c, d) ...] — Вершины полигона. Тип данных — Array или Ring. Каждая вершина задаётся парой координат (a, b). Вершины следует указывать по часовой стрелке или против часовой стрелки. Минимальное число вершин — 3.
• Функция поддерживает полигоны с отверстиями (вырезанными областями). Тип данных — Polygon. Либо передайте весь Polygon как второй аргумент, либо сначала передайте внешний Ring, а затем каждое отверстие как отдельный дополнительный аргумент.
• Функция также поддерживает MultiPolygon. Тип данных — MultiPolygon. Либо передайте весь MultiPolygon как второй аргумент, либо перечислите каждый входящий в него полигон как отдельный аргумент.
• Аргумент полигона также может быть столбцом Geometry, содержащим значения полигональных типов (Ring, Polygon или MultiPolygon).

SELECT pointInPolygon((3., 3.), [(6, 0), (8, 4), (5, 8), (0, 2)]) AS res

┌─res─┐
│   1 │
└─────┘

CREATE TABLE poly (id UInt32, shape Polygon) ENGINE = Memory;
INSERT INTO poly VALUES (1, [[(0, 0), (10, 0), (10, 10), (0, 10)], [(4, 4), (6, 4), (6, 6), (4, 6)]]);
SELECT id, pointInPolygon((2., 2.), shape) AS res FROM poly;

┌─id─┬─res─┐
│  1 │   1 │
└────┴─────┘

Примечание • Вы можете установить validate_polygons = 0, чтобы отключить проверку геометрии. • pointInPolygon исходит из того, что каждый полигон сформирован корректно. Если входной полигон самопересекается, содержит кольца в неправильном порядке или перекрывающиеся рёбра, результаты становятся ненадёжными — особенно для точек, которые лежат точно на ребре, в вершине или внутри самопересечения, где невозможно однозначно определить, что считать “внутри”, а что — “снаружи”. • Когда аргумент полигона является константой, а точка задаётся с использованием индексируемых столбцов ключа (например, pointInPolygon((x, y), constant_polygon) в таблице, где x, y входят в PRIMARY KEY или покрываются индексом minmax), ClickHouse может использовать и первичный ключ, и индексы пропуска данных minmax, чтобы отсечь нерелевантные гранулы.