دوال للتعامل مع الإحداثيات الجغرافية

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 — قيمة Tuple مكوّنة من رقمين، أو Point.
• [(a, b), (c, d) ...] — رؤوس المضلع. نوع البيانات — Array أو Ring. تُمثَّل كل قمة بزوج من الإحداثيات (a, b). يجب تحديد الرؤوس بترتيب مع عقارب الساعة أو عكس اتجاه عقارب الساعة. الحد الأدنى لعدد الرؤوس هو 3.
• تدعم الدالة أيضًا المضلع ذي الفتحات (الأجزاء المقتطعة). نوع البيانات — Polygon. مرِّر إما Polygon بالكامل باعتباره المُدخل الثاني، أو مرِّر الحلقة الخارجية أولًا ثم كل فتحة كمدخل إضافي منفصل.
• تدعم الدالة أيضًا 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 لاستبعاد الحبيبات غير ذات الصلة.