تحويل هيلمرت

يمكن وصف التحويل من إطار مرجعي 1 إلى إطار مرجعي 2 بثلاث عمليات إزاحة Δx، Δy، Δz، وثلاث عمليات دوران Rx، Ry، Rz، ومعامل قياس μ.

تحويل هيلمرت ( نسبةً إلى فريدريك روبرت هيلمرت ، 1843-1917 ) هو طريقة لتحديد معلمات التشابهات الهندسية ، ويُستخدم في علم المساحة لحساب التحويلات بين المراجع . يُعرف تحويل هيلمرت أيضًا باسم التحويل ذي المعلمات السبع .

تعريف

يمكن التعبير عن تحويل هيلمرت على النحو التالي:

Xتي=ج+μRX{\displaystyle X_{T}=C+\mu RX\,}

أين

  • X هي النقطة الأولية التي يمثلها متجه الإحداثيات
  • X T هو المتجه المحوّل نقطيًا

المعايير هي :

الاختلافات

تُعدّ تحويلة هيلمرت ثنائية الأبعاد حالةً خاصة. في هذه الحالة، لا نحتاج إلا إلى أربعة مُعاملات (إزاحتان، وتغيير في الحجم، ودوران). يُمكن تحديد هذه المُعاملات من نقطتين معروفتين؛ وإذا توفرت نقاط إضافية، يُمكن إجراء عمليات التحقق.

في بعض الأحيان يكون من الكافي استخدام التحويل ذي المعلمات الخمس ، والذي يتكون من ثلاث عمليات نقل، ودوران واحد فقط حول المحور Z، وتغيير واحد في المقياس.

قيود

يستخدم تحويل هيلمرت عامل قياس واحد فقط، لذا فهو غير مناسب لما يلي:

في هذه الحالات، يكون التحويل الأفيني الأكثر عمومية هو الأفضل.

طلب

تُستخدم تحويلة هيلمرت، من بين أمور أخرى، في علم المساحة لتحويل إحداثيات النقطة من نظام إحداثيات إلى آخر. وباستخدامها، يصبح من الممكن تحويل نقاط المسح الإقليمية إلى مواقع WGS84 التي يستخدمها نظام تحديد المواقع العالمي (GPS) .

على سبيل المثال، بدءًا من إحداثيات جاوس-كروجر ، x و y ، بالإضافة إلى الارتفاع h ، يتم تحويلها إلى قيم ثلاثية الأبعاد على مراحل:

  1. إلغاء إسقاط الخريطة : حساب خط العرض وخط الطول والارتفاع الإهليلجي ( W , L , H )
  2. التحويل من الإحداثيات الجيوديسية إلى الإحداثيات الجيومركزية : حساب قيم x و y و z بالنسبة إلى القطع الناقص المرجعي للمسح
  3. تحويل ذو 7 معلمات (حيث تتغير x و y و z دائمًا تقريبًا ببضع مئات من الأمتار على الأكثر، والمسافات ببضعة ملليمترات لكل كيلومتر).
  4. وبسبب هذا، يمكن مقارنة المواقع المقاسة أرضياً ببيانات نظام تحديد المواقع العالمي (GPS)؛ ويمكن بعد ذلك إدخال هذه البيانات في عملية المسح كنقاط جديدة - يتم تحويلها بالترتيب المعاكس.

تتألف الخطوة الثالثة من تطبيق مصفوفة الدوران ، والضرب بمعامل المقياس.μ=1+s{\displaystyle \mu =1+s}(بقيمة قريبة من 1) وإضافة عمليات الإزاحة الثلاث، c x ، c y ، c z .

يتم اشتقاق إحداثيات نظام مرجعي B من النظام المرجعي A بواسطة الصيغة التالية (اتفاقية تحويل متجه الموضع وتبسيط زوايا الدوران الصغيرة جدًا): [ 1 ]

[XYZ]ب=[جxجyجz]+(1+s×10-6)[1-رzرyرz1-رx-رyرx1][XYZ]أ{\displaystyle {\begin{bmatrix}X\\Y\\Z\end{bmatrix}}^{B}={\begin{bmatrix}c_{x}\\c_{y}\\c_{z}\end{bmatrix}}+(1+s\times 10^{-6})\cdot {\begin{bmatrix}1&-r_{z}&r_{y}\\r_{z}&1&-r_{x}\\-r_{y}&r_{x}&1\end{bmatrix}}\cdot {\begin{bmatrix}X\\Y\\Z\end{bmatrix}}^{A}}

أو لكل معلمة منفردة من معلمات الإحداثيات:

Xب=جx+(1+s×10-6)(Xأ-رzYأ+رyZأ)Yب=جy+(1+s×10-6)(رzXأ+Yأ-رxZأ)Zب=جz+(1+s×10-6)(-رyXأ+رxYأ+Zأ).{\displaystyle {\begin{aligned}X_{B}&=c_{x}+(1+s\times 10^{-6})\cdot (X_{A}-r_{z}\cdot Y_{A}+r_{y}\cdot Z_{A})\\Y_{B}&=c_{y}+(1+s\times 10^{-6})\cdot (r_{z}\cdot X_{A}+Y_{A}-r_{x}\cdot Z_{A})\\Z_{B}&=c_{z}+(1+s\times 10^{-6})\cdot (-r_{y}\cdot X_{A}+r_{x}\cdot Y_{A}+Z_{A}).\end{aligned}}}

أما بالنسبة للتحويل العكسي، فيتم ضرب كل عنصر في -1 .

تُحدد المعايير السبعة لكل منطقة بثلاث نقاط متطابقة أو أكثر من كلا النظامين. ولتحقيق التوافق بينهما، تُعالج الاختلافات الطفيفة (عادةً بضعة سنتيمترات فقط) باستخدام طريقة المربعات الصغرى ، أي تُزال بطريقة معقولة إحصائيًا.

المعايير القياسية

ملاحظة: زوايا الدوران الواردة في الجدول هي بالثواني القوسية ويجب تحويلها إلى راديان قبل استخدامها في الحساب.
رمز EPSGمنطقةبيانات المصدرالبيانات المستهدفةالدقة ( متر )ج × (متر)ج س (متر)ج ز (متر)s ( جزء في المليون )r x ( ثانية قوسية )r y ( ثانية قوسية )r z ( ثانية قوسية )
8048أسترالياGDA94 ( EPSG:4283 )GDA2020 ( EPSG:7844 )0.010.06155-0.01087-0.04019-0.009994-0.0394924-0.0327221-0.0328979
9690أسترالياWGS84 ( EPSG:4326 )GDA2020 ( EPSG:7844 )30.06155-0.01087-0.04019-0.009994-0.0394924-0.0327221-0.0328979
1618النمساMGI ( EPSG:4312 )WGS84 ( EPSG:4326 )1.5577.32690.129463.9192.42325.1371.4745.297
1776ألمانيا (الغربية)DHDN ( EPSG:4314 )ETRS89 ( EPSG:4258 )3598.173.7418.26.70.2020.045-2.455
1777ألمانيا (الغربية)DHDN ( EPSG:4314 )WGS84 ( EPSG:4326 )3598.173.7418.26.70.2020.045-2.455
15869ألمانيا (الشرقية)DHDN ( EPSG:4314 )WGS84 ( EPSG:4326 )2612.477440.22.55-0.0540.057-2.797
1641أيرلنداTM65 ( EPSG:4299 )WGS84 ( EPSG:4326 )1482.5-130.6564.68.15-1.042-0.214-0.631
1953أيرلنداTM75 ( EPSG:4300 )ETRS89 ( EPSG:4258 )1482.5-130.6564.68.15-1.042-0.214-0.631
1954أيرلنداTM75 ( EPSG:4300 )WGS84 ( EPSG:4326 )1482.5-130.6564.68.15-1.042-0.214-0.631
8689سلوفينياMGI 1901 ( EPSG:3906 )سلوفينيا 1996 ( EPSG:4765 )1476.08125.947417.819.896638-4.610862-2.38813711.942335
1314المملكة المتحدةOSGB36 ( EPSG:4247 )WGS84 ( EPSG:4326 )2446.448-125.157542.06-20.4890.150.2470.842
1315المملكة المتحدةOSGB36 ( EPSG:4247 )ED50 ( EPSG:4230 )2535.948-31.357665.16-21.6890.150.2470.998
1901الولايات المتحدةNAD83(HARN) ( EPSG:4152 )WGS84 ( EPSG:4326 )1-0.9911.90720.512901.250 33 × 10 −74.6785 × 10 −85.6529 × 10 −8

هذه مجموعات معلمات قياسية للتحويل ذي السبع معلمات (أو تحويل البيانات) بين مرجعين. وللتحويل في الاتجاه المعاكس، يجب حساب معلمات التحويل العكسي أو تطبيق التحويل العكسي (كما هو موضح في الورقة البحثية "حول التحويلات الجيوديسية" [ 2 ] ). تُوصف الإزاحات c<sub> x </sub> و c<sub> y</sub> و c<sub> z </sub> أحيانًا بـ t<sub> x</sub> و t <sub>y </sub> و t <sub> z</sub> ، أو d<sub> x </sub> و dy و d <sub> z</sub> . كما تُوصف الدورانات r<sub> x</sub> و r<sub> y</sub> و r<sub> z</sub> أحيانًا بـω{\displaystyle \omega }،ϕ{\displaystyle \phi }وκ{\displaystyle \kappa }في المملكة المتحدة ، ينصبّ الاهتمام الرئيسي على التحويل بين نظام الإحداثيات OSGB36، الذي تستخدمه هيئة المساحة البريطانية لتحديد إحداثيات الشبكة على خرائطها من طراز Landranger وExplorer، ونظام WGS84 المستخدم في تقنية GPS. أما نظام إحداثيات Gauss-Krüger المستخدم في ألمانيا، فيعتمد عادةً على إهليلج بيسل . ومن أنظمة الإحداثيات الأخرى ذات الأهمية نظام ED50 (نظام الإحداثيات الأوروبي 1950)، الذي يستند إلى إهليلج هايفورد . وقد شكّل ED50 جزءًا من أساسيات إحداثيات حلف شمال الأطلسي (الناتو) حتى ثمانينيات القرن الماضي، وتُعرَّف العديد من أنظمة الإحداثيات الوطنية القائمة على Gauss-Krüger باستخدام ED50.

لا يمتلك سطح الأرض شكلاً بيضاوياً مثالياً، بل يُوصف بأنه جيود . وبدلاً من ذلك، يُوصف جيود الأرض بعدة أشكال بيضاوية. وبحسب الموقع الفعلي، يُستخدم "الشكل البيضاوي الأنسب محلياً" لأغراض المسح ورسم الخرائط. وتُعطي مجموعة المعايير القياسية دقة تبلغ حوالي7  أمتار لتحويل OSGB36/WGS84. هذه الدقة غير كافية لأغراض المسح، ولذا تُكمل هيئة المساحة البريطانية هذه النتائج باستخدام جدول بحث يتضمن تحويلات إضافية للوصول إلى دقة 1 سم .

تقدير المعلمات

إذا كانت معاملات التحويل غير معروفة، فيمكن حسابها باستخدام نقاط مرجعية (أي نقاط معروفة إحداثياتها قبل التحويل وبعده). وبما أنه يجب تحديد سبعة معاملات إجمالاً (ثلاثة تحويلات إزاحة، ومعامل تغيير حجم، وثلاثة تحويلات دوران)، فيجب معرفة نقطتين على الأقل وإحداثي نقطة ثالثة (على سبيل المثال، الإحداثي Z). ينتج عن ذلك نظام من سبع معادلات وسبعة مجاهيل، يمكن حله.

بالنسبة للتحويلات بين إسقاطات الخرائط المطابقة بالقرب من نقطة عشوائية، يمكن حساب معلمات تحويل هيلمرت بدقة من مصفوفة جاكوبيان لدالة التحويل.

من الناحية العملية، يُفضّل استخدام عدد أكبر من النقاط. فمن خلال هذا التوافق، يتم الحصول على دقة أعلى، ويصبح التقييم الإحصائي للنتائج ممكنًا. في هذه الحالة، يتم تعديل الحساب باستخدام طريقة المربعات الصغرى الغاوسية .

يتم الحصول على قيمة عددية لدقة معلمات التحويل عن طريق حساب القيم عند النقاط المرجعية، وترجيح النتائج بالنسبة إلى مركز النقاط.

على الرغم من أن هذه الطريقة دقيقة رياضياً، إلا أنها تعتمد كلياً على دقة المعايير المستخدمة. عملياً، تُحسب هذه المعايير من خلال تضمين ثلاث نقاط معروفة على الأقل في الشبكات. مع ذلك، ستؤثر دقة هذه النقاط على معايير التحويل اللاحقة، إذ ستتضمن هذه النقاط أخطاء رصد. لذا، فإن التحويل "الواقعي" لن يكون سوى أفضل تقدير، ويجب أن يتضمن مقياساً إحصائياً لجودته.

انظر أيضاً

مراجع

  1. "المعادلات المستخدمة لتحويلات البيانات" . هيئة معلومات الأراضي في نيوزيلندا (Toitū Te Whenua) . تم الاطلاع عليه بتاريخ 30 يونيو 2022 .
  2. حول التحويلات الجيوديسية، بو-غونار ريت، 2009 https://www.lantmateriet.se/contentassets/4a728c7e9f0145569edd5eb81fececa7/rapport_reit_eng.pdf