أكثر

تقريب إحداثيات GPS لنقطة من n مواقع متوقعة


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

بعد رسم هذه الإحداثيات على الخريطة ، لاحظت ما يلي: مرتين أو ثلاث مرات من أصل عشر مرات ، أعطى جهاز GPS الخاص بي مجموعة خاطئة من الإحداثيات والتي تبين أنها هادئة بعيدًا عن النقطة X. أيضًا بسبب حركة المرور ، فأنا أحيانًا غير قادر على الوقوف بالقرب من النقطة x ، وبالتالي في هذه الحالة أيضًا ، تكون الإحداثيات التي تم الحصول عليها بعيدة عن النقطة X.

المشكلة: من بين مجموعات الإحداثيات الـ 16 التي تم الحصول عليها ، ما هي العملية التي أستخدمها لتضييق نطاقها إلى مجموعة واحدة من الإحداثيات تكون قريبة جدًا من نقطة انتباهي (النقطة X)؟


طريقة واحدة للتعامل مع هذه المشكلة المثيرة للاهتمام هي النظر إليها على أنها مقدر قوي لمركز توزيع نقطة ثنائية المتغير. الحل (المعروف) هو تقشر الهياكل المحدبة حتى لا يتبقى شيء. النقطه الوسطى لآخر بدن غير فارغ يحدد موقع المركز.

(هذا مرتبط بـ bagplot. لمزيد من المعلومات ، ابحث في الويب عن "محدب تقشير البدن الخارجي متعدد المتغيرات.")

تظهر نتيجة 16 نقطة مصورة كمثلث مركزي في هذه الخريطة. تظهر المضلعات الثلاثة المحيطة الهياكل المحدبة المتتالية. تمت إزالة النقاط الخمس النائية (30٪ من الإجمالي!) في أول خطوتين.


تم حساب المثال فيص. يتم تنفيذ الخوارزمية نفسها في الكتلة الوسطى ، "التقشير المحدب". يستخدم المدمج فيتشولروتين ، والذي يقوم بإرجاع فهارس النقاط على بدن السفينة. تتم إزالة هذه النقاط عن طريق تعبير الفهرسة السلبيxy [-hull،]. يتكرر هذا حتى تتم إزالة النقاط الأخيرة. في الخطوة الأخيرة ، يتم حساب النقطه الوسطى عن طريق حساب متوسط ​​الإحداثيات.

لاحظ أنه في كثير من الحالات ، لا يكون إسقاط البيانات ضروريًا: لن تتغير الهياكل المحدبة ما لم تمتد الميزات الأصلية على خط الطول المضاد (+/- 180 درجة خط الطول) ، إما قطب ، أو واسعة جدًا لدرجة أن انحناء الأجزاء بينها سوف احداث فرق. (حتى ذلك الحين ، لن يكون الانحناء مصدر قلق كبير ، لأن التقشير سيظل يتقارب إلى نقطة مركزية).

# # مشروع البيانات. # dy <- c (8،7،5،10،7،17،19،19،21،22،22،22،24،24،26،26) dx <- c (66،67،66،89) ، 89،79،78،76،75،81،78،77،75،80،77،83) خط العرض <- (28.702 + dy / 1e5) / 180 * pi lon <- (77.103 + DX / 1e5) / 180 * pi y <- dy x <- cos (mean (lat)) * dx # # تقشير محدب. # xy <- cbind (x، y) while (TRUE) {hull <- chull (xy) if (length (hull) 


شاهد الفيديو: كيفية تصحيح نقطة مرجعيه تم رصدها بــ GPS أونلاين post processing static base point gps data (شهر اكتوبر 2021).