地理座標を扱う関数

greatCircleDistance(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — 1 点目の経度 (度単位) 。範囲: [-180°, 180°]。
• lat1Deg — 1 点目の緯度 (度単位) 。範囲: [-90°, 90°]。
• lon2Deg — 2 点目の経度 (度単位) 。範囲: [-180°, 180°]。
• lat2Deg — 2 点目の緯度 (度単位) 。範囲: [-90°, 90°]。

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

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

geoDistance(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — 1 番目の地点の経度 (度) 。範囲: [-180°, 180°]。
• lat1Deg — 1 番目の地点の緯度 (度) 。範囲: [-90°, 90°]。
• lon2Deg — 2 番目の地点の経度 (度) 。範囲: [-180°, 180°]。
• lat2Deg — 2 番目の地点の緯度 (度) 。範囲: [-90°, 90°]。

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

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

greatCircleAngle(lon1Deg, lat1Deg, lon2Deg, lat2Deg)

• lon1Deg — 1 番目の点の経度 (度) 。
• lat1Deg — 1 番目の点の緯度 (度) 。
• lon2Deg — 2 番目の点の経度 (度) 。
• lat2Deg — 2 番目の点の緯度 (度) 。

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 m の仮東距を含む) 。(U)Int*/Float*。
• northing — 北距 (メートル単位。南半球では 10000000 m の仮北距を含む) 。(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 — 2 つの数値からなるタプル、または Point。
• [(a, b), (c, d) ...] — Polygon の頂点。データ型 — Array または Ring。各頂点は座標の組 (a, b) で表されます。頂点は時計回りまたは反時計回りの順序で指定する必要があります。頂点の最小数は 3 です。
• この関数は穴 (切り抜き部分) のある Polygon もサポートします。データ型 — Polygon。Polygon 全体を第 2 引数として渡すか、外側のリングを最初に渡し、その後に各穴を個別の追加引数として渡します。
• この関数は MultiPolygon もサポートします。データ型 — MultiPolygon。MultiPolygon 全体を第 2 引数として渡すか、各構成 Polygon をそれぞれ個別の引数として列挙します。
• Polygon 引数には、Polygon 形状の値 (Ring、Polygon、または MultiPolygon) を保持する Geometry カラムも使用できます。

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 は、すべてのポリゴンが正しい形式であることを前提としています。入力が自己交差している、リングの順序が誤っている、または辺同士が重なっている場合、結果は信頼できなくなります。特に、点が辺上や頂点上にちょうどある場合や、“内側” と “外側” の区別が定義できない自己交差部分の内部にある場合は、その傾向が顕著です。 • ポリゴン引数が定数で、点が索引付きのキーカラムで表されている場合 (たとえば、x, y が PRIMARY KEY の一部である、または minmax 索引でカバーされているテーブル上の pointInPolygon((x, y), constant_polygon) など) 、ClickHouse は主キーと minmax のデータスキップ索引の両方を使って、無関係なグラニュールを刈り込むことができます。