مسارات تقسيم الأرض

مقطع مستوٍ لشكل بيضاوي

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

المسألة العكسية

المسألة العكسية لقطاعات الأرض هي: بمعلومية نقطتين،P1{\displaystyle P_{1}}وP2{\displaystyle P_{2}}أوجد الطول على سطح القطع الناقص المرجعي.s12{\displaystyle s_{12}}، من القوس القصير لمقطع كروي منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}وكذلك إيجاد سمت المغادرة والوصول (الزاوية من الشمال الحقيقي) لهذا المنحنى،α1{\displaystyle \alpha _{1}}وα2{\displaystyle \alpha _{2}}يوضح الشكل الموجود على اليمين الرموز المستخدمة هنا. ليكنPك{\displaystyle P_{k}}لها خط عرض جيوديسي ϕك{\displaystyle \phi _{k}}وخط الطولλك{\displaystyle \lambda _{k}}( k = 1, 2). يُفضل حل هذه المسألة باستخدام الهندسة التحليلية في إحداثيات ديكارتية مركزية الأرض وثابتة حولها (ECEF).R1=هـجهـF(P1){\displaystyle R_{1}=\mathrm {ECEF} (P_{1})}وR2=هـجهـF(P2){\displaystyle R_{2}=\mathrm {ECEF} (P_{2})}لتكن إحداثيات ECEF للنقطتين، محسوبة باستخدام التحويل الجيوديسي إلى ECEF الذي تمت مناقشته هنا .

يوضح هذا الترميز المستخدم في المسائل الجيوديسية التي نوقشت هنا.

مستوى المقطع

لتحديد مستوى القطع، حدد أي نقطة ثالثةR0{\displaystyle R_{0}}ليس على الخط منR1{\displaystyle R_{1}}لR2{\displaystyle R_{2}}اختيارR0{\displaystyle R_{0}}أن يكون على السطح عموديًا عندP1{\displaystyle P_{1}}سيحدد القسم العادي عندP1{\displaystyle P_{1}}. لوR0{\displaystyle R_{0}}إذا كانت نقطة الأصل هي نقطة الأرض، فإن مقطع الأرض هو القطع الناقص العظيم. (ستكون نقطة الأصل على خط مستقيم واحد مع نقطتين متقابلتين، لذا يجب استخدام نقطة مختلفة في هذه الحالة). بما أن هناك عددًا لا نهائيًا من الخيارات لـR0{\displaystyle R_{0}}المسألة المذكورة أعلاه هي في الواقع فئة من المسائل (مسألة لكل مستوى). لنفترضR0{\displaystyle R_{0}}لنفترض أن لدينا معطى. لوضع معادلة المستوى في الصورة القياسية،لx+مy+نz=د{\displaystyle lx+my+nz=d}، أينل2+م2+ن2=1{\displaystyle l^{2}+m^{2}+n^{2}=1}، يتطلب مكونات متجه الوحدة ،شمال^=(ل،م،ن){\displaystyle \mathbf {\hat {N}} =(l,m,n)}، عمودي على مستوى المقطع. يمكن حساب هذه المكونات كما يلي: المتجه منR0{\displaystyle R_{0}}لR1{\displaystyle R_{1}}يكونV0=R1-R0{\displaystyle \mathbf {V_{0}} =\mathbf {R_{1}} -\mathbf {R_{0}} }والمتجه منR1{\displaystyle R_{1}}لR2{\displaystyle R_{2}}يكونV1=R2-R1{\displaystyle \mathbf {V_{1}} =\mathbf {R_{2}} -\mathbf {R_{1}} }. لذلك،شمال^=uنأنات(V0×V1){\displaystyle \mathbf {\hat {N}} =\mathrm {unit} (\mathbf {V_{0}} \times \mathbf {V_{1}} )})، أينuنأنات(V){\displaystyle \mathrm {unit} (\mathbf {V} )}هو متجه الوحدة في اتجاهV{\displaystyle \mathbf {V} }. اصطلاح التوجيه المستخدم هنا هو أنشمال^{\displaystyle \mathbf {\hat {N}} }تشير إلى يسار المسار. إذا لم يكن الأمر كذلك، فأعد تعريفها.V0=-V0{\displaystyle \mathbf {V_{0}} =-\mathbf {V_{0}} }وأخيرًا، يمكن حساب المعامل d للمستوى باستخدام الضرب النقطي لـشمال^{\displaystyle \mathbf {\hat {N}} }بمتجه من نقطة الأصل إلى أي نقطة على المستوى، مثلR1{\displaystyle R_{1}}، أيد=شمال^R1{\displaystyle d=\mathbf {\hat {N}} \cdot \mathbf {R_{1}} }وبالتالي، فإن معادلة المستوى (في الصورة المتجهة) هيشمال^R=د{\displaystyle \mathbf {\hat {N}} \cdot \mathbf {R} =d}، أينR{\displaystyle \mathbf {R} }هو متجه الموضع لـ(x،y،z){\displaystyle (x,y,z)}.

السمت

يكشف فحص تحويل ENU إلى ECEF أن إحداثيات ECEF لمتجه وحدة يشير إلى الشرق عند أي نقطة على القطع الناقص هي:هـ^=(-الخطيئةλ،كوسλ،0){\displaystyle \mathbf {\hat {e}} =(-\sin \lambda ,\cos \lambda ,0)}متجه الوحدة الذي يشير إلى الشمال هون^=(-الخطيئةϕكوسλ،-الخطيئةϕالخطيئةλ،كوسϕ){\displaystyle \mathbf {\hat {n}} =(-\sin \phi \cos \lambda ,-\sin \phi \sin \lambda ,\cos \phi )}والمتجه الوحدوي الذي يشير إلى الأعلى هوu^=(كوسϕكوسλ،كوسϕالخطيئةλ،الخطيئةϕ){\displaystyle \mathbf {\hat {u}} =(\cos \phi \cos \lambda ,\cos \phi \sin \lambda ,\sin \phi )}المتجه المماس للمسار هو: ت=شمال^×u^{\displaystyle \mathbf {t} =\mathbf {\hat {N}} \times \mathbf {\hat {u}} }إذن، الجزء الشرقي منت{\displaystyle \mathbf {t} }يكونتهـ^{\displaystyle \mathbf {t} \cdot \mathbf {\hat {e}} }والمكون الشمالي هوتن^{\displaystyle \mathbf {t} \cdot \mathbf {\hat {n}} }لذلك، يمكن الحصول على السمت من دالة الظل العكسي ذات الوسيطين .α=atan2(تهـ^،تن^){\displaystyle \alpha =\operatorname {atan2} (\mathbf {t} \cdot \mathbf {\hat {e}} ,\mathbf {t} \cdot \mathbf {\hat {n}} )}استخدم هذه الطريقة في كليهماP1{\displaystyle P_{1}}وP2{\displaystyle P_{2}}للحصول علىα1{\displaystyle \alpha _{1}}وα2{\displaystyle \alpha _{2}}.

مقطع بيضاوي

إن تقاطع المستوى مع القطع الناقص (وهو تقاطع غير بديهي) ينتج عنه قطع ناقص. لذلك، فإن طول القوس،s12{\displaystyle s_{12}}، على مسار القسم منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}هو تكامل إهليلجي يمكن حسابه بدقة مطلوبة باستخدام متسلسلة مقطوعة أو تكامل عددي . قبل ذلك، يجب تعريف القطع الناقص وحساب حدود التكامل. ليكن القطع الناقص معطى بالعلاقة التالية:x2أ2+y2أ2+z2ب2=1{\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{a^{2}}}+{\frac {z^{2}}{b^{2}}}=1}ودعص=ل2+م2{\displaystyle p={\sqrt {l^{2}+m^{2}}}}. لوص=0{\displaystyle p=0}إذن، المقطع عبارة عن دائرة أفقية نصف قطرهاأ1-د2ب2{\textstyle a{\sqrt {1-{\frac {d^{2}}{b^{2}}}}}}والتي ليس لها حل إذا|د|>ب{\displaystyle |d|>b}.

لوص>0{\displaystyle p>0}ثم أظهر جيلبرتسون [ 1 ] أن إحداثيات ECEF لمركز القطع الناقص هيRج=دج(لأ2،مأ2،نب2){\textstyle {R_{c}}={\frac {d}{C}}(la^{2},ma^{2},nb^{2})}، أينج=أ2ص2+ب2ن2{\displaystyle C=a^{2}p^{2}+b^{2}n^{2}}،

المحور شبه الرئيسي هوأ*=أ1-د2ج{\textstyle a^{*}=a{\sqrt {1-{\frac {d^{2}}{C}}}}}، في الاتجاهأنا^*=(مص،-لص،0){\textstyle \mathbf {{\hat {i}}^{*}} =\left({\frac {m}{p}},{\frac {-l}{p}},0\right)}والمحور شبه الأصغر هوب*=بجأ*{\textstyle b^{*}={\frac {b}{\sqrt {C}}}a^{*}}، في الاتجاهج^*=(لنص،منص،-ص){\textstyle \mathbf {{\hat {j}}^{*}} =\left({\frac {ln}{p}},{\frac {mn}{p}},-p\right)}والتي ليس لها حل إذا|د|>ج{\displaystyle |d|>{\sqrt {C}}}.

طول القوس

تُقدّم الورقة البحثية المشار إليها أعلاه اشتقاقًا لصيغة طول القوس التي تتضمن الزاوية المركزية وقوىهـ2{\displaystyle e^{2}}لحساب طول القوس بدقة تصل إلى المليمتر، حيثهـ2=1-(ب*أ*)2{\textstyle e^{2}=1-\left({\frac {b^{*}}{a^{*}}}\right)^{2}}يمكن إعادة ترتيب صيغة طول القوس هذه ووضعها في الشكل التالي: s12=s(θ2)-s(θ1){\displaystyle s_{12}=s(\theta _{2})-s(\theta _{1})}، أين s(θ)=ب*(ج0θ+ج2الخطيئة(2θ)+ج4الخطيئة(4θ)+ج6الخطيئة(6θ)){\displaystyle s(\theta )=b^{*}({C_{0}}\theta +{C_{2}}\sin(2\theta )+{C_{4}}\sin(4\theta )+{C_{6}}\sin(6\theta ))}والمعاملات هي

ج0=1.0+هـ2(1/4+13هـ2/64+45هـ4/256+2577هـ6/16384){\displaystyle C_{0}=1.0+e^{2}(1/4+13e^{2}/64+45e^{4}/256+2577e^{6}/16384)}
ج2=هـ2(1/8+3هـ2/32+95هـ4/1024+385هـ6/4096){\displaystyle C_{2}=e^{2}(1/8+3e^{2}/32+95e^{4}/1024+385e^{6}/4096)}
ج4=-هـ4(1/256+5هـ2/1024+19هـ4/16384){\displaystyle C_{4}=-e^{4}(1/256+5e^{2}/1024+19e^{4}/16384)}
ج6=-هـ6(15/3072+35هـ2/4096){\displaystyle C_{6}=-e^{6}(15/3072+35e^{2}/4096)}

لحساب الزاوية المركزية، دعP{\displaystyle P}أي نقطة على القطع الناقص وR=هـجهـF(P){\displaystyle R=\mathrm {ECEF} (P)}. ثمV=R-Rج{\displaystyle \mathbf {V} =\mathbf {R} -\mathbf {R_{c}} }هو متجه من مركز القطع الناقص إلى النقطة. الزاوية المركزيةθ{\displaystyle \theta }هي الزاوية من المحور شبه الرئيسي إلىV{\displaystyle \mathbf {V} }تأجيرV^=uنأنات(V){\displaystyle \mathbf {\hat {V}} =\mathrm {unit} (\mathbf {V} )}لديناθ=atan2(V^ج^*،V^أنا^*){\displaystyle \theta =\operatorname {atan2} (\mathbf {\hat {V}} \cdot \mathbf {{\hat {j}}^{*}} ,\mathbf {\hat {V}} \cdot \mathbf {{\hat {i}}^{*}} )}وبهذه الطريقة نحصل علىθ1{\displaystyle \theta _{1}}وθ2{\displaystyle \theta _{2}}.

من ناحية أخرى، من الممكن استخدام صيغ قوس خط الزوال في الحالة الأكثر عمومية بشرط استخدام معلمات القطع الناقص بدلاً من معلمات الشكل الكروي. إحدى هذه المتسلسلات سريعة التقارب معطاة في المتسلسلات بدلالة خط العرض البارامتري . إذا استخدمناε{\displaystyle \varepsilon }للدلالة على انحراف الشكل الكروي، أيε2=1-(بأ)2{\textstyle \varepsilon ^{2}=1-\left({\frac {b}{a}}\right)^{2}}، ثمهـ8{\displaystyle e^{8}}ε8{\displaystyle \varepsilon ^{8}}1.8 × 10⁻⁹ . وبالمثل ، فإن التسطح الثالث للقطع الناقص محدود بالقيمة المقابلة للشكل الكروي، وبالنسبة للشكل الكروي لدينان3{\displaystyle n^{3}}4.4 × 10⁻⁹ ، ون4{\displaystyle n^{4}}7.3 × 10⁻¹² . لذلك ، قد يكفي تجاهل الحدود التي تتجاوزب6{\displaystyle B_{6}}في سلسلة خطوط العرض البارامترية. للتطبيقs(β)=أ*+ب*2(ب0β+ب2الخطيئة2β+ب4الخطيئة4β+ب6الخطيئة6β){\textstyle s(\beta )={\frac {a^{*}+b^{*}}{2}}(B_{0}\beta +B_{2}\sin 2\beta +B_{4}\sin 4\beta +B_{6}\sin 6\beta )}يتطلب السياق الحالي تحويل الزاوية المركزية إلى الزاوية البارامترية باستخدامβ=لون برونزي-1(لون برونزيθ/((1-و)){\displaystyle \beta =\tan ^{-1}\left(\tan \theta /((1-f)\right)}وباستخدام تقنية التسطيح الثالث للقطع الناقص. أيًا كانت الطريقة المستخدمة، يجب توخي الحذر عند استخدامها.θ1{\displaystyle \theta _{1}}وθ2{\displaystyle \theta _{2}}أوβ1{\displaystyle \beta _{1}}وβ2{\displaystyle \beta _{2}}لضمان استخدام القوس الأقصر الذي يربط النقطتين.

مشكلة مباشرة

تم طرح المسألة المباشرةP1{\displaystyle {P_{1}}}المسافةs12{\displaystyle s_{12}}، واتجاه المغادرةα1{\displaystyle \alpha _{1}}، يجدP2{\displaystyle {P_{2}}}وسمت الوصولα2{\displaystyle \alpha _{2}}.

مستوى المقطع

يعتمد حل هذه المشكلة على اختيارV0{\displaystyle \mathbf {V_{0}} }أي على نوع القسم. لاحظ ذلكV0{\displaystyle \mathbf {V_{0}} }يجب ألا يكون داخل span{ن^1،هـ^1{\displaystyle \mathbf {\hat {n}} _{1},\mathbf {{\hat {e}}_{1}} }} (وإلا ستكون الطائرة مماسية للأرض عندP1{\displaystyle {P_{1}}}(لذا لن ينتج أي مسار). بعد اتخاذ هذا الخيار، وبالنظر إلى التوجيه، اتبع الخطوات التالية. أنشئ متجه المماس عندP1{\displaystyle {P_{1}}}،ت^1=ن^1كوسα1+هـ^1الخطيئةα1{\displaystyle \mathbf {\hat {t}} _{1}=\mathbf {\hat {n}} _{1}\cos {\alpha _{1}}+\mathbf {{\hat {e}}_{1}} \sin {\alpha _{1}}}، أينن^1{\displaystyle \mathbf {\hat {n}} _{1}}وهـ^1{\displaystyle \mathbf {{\hat {e}}_{1}} }هي متجهات وحدة تشير إلى الشمال والشرق (على التوالي) عندP1{\displaystyle {P_{1}}}المتجه العموديشمال^=uنأنات(V0×ت^1{\displaystyle \mathbf {\hat {N}} =\mathrm {unit} (\mathbf {V_{0}} \times \mathbf {\hat {t}} _{1}})، بالإضافة إلىP1{\displaystyle \mathbf {P_{1}} }يُحدد المستوى. بعبارة أخرى، يحل المماس محل الوتر لأن الوجهة غير معروفة.

تحديد نقطة الوصول

هذه مسألة ثنائية الأبعاد في الفضاء الممتد {أنا^*،ج^*{\displaystyle \mathbf {\hat {i}} ^{*},\mathbf {\hat {j}} ^{*}}، والتي سيتم حلها باستخدام صيغة طول القوس المذكورة أعلاه. إذا كان طول القوس،s12{\displaystyle s_{12}}إذا كانت المعطيات معطى، فإن المسألة تكمن في إيجاد التغير المقابل في الزاوية المركزية θ12{\displaystyle \theta _{12}}، لهذا السببθ2=θ1+θ12{\displaystyle \theta _{2}=\theta _{1}+\theta _{12}}ويمكن حساب الموضع. بافتراض أن لدينا متسلسلة تعطيs=s(θ){\displaystyle s=s(\theta )}إذن ما نسعى إليه الآن هوθ2=s-1(s1+s12){\displaystyle \theta _{2}=s^{-1}(s_{1}+s_{12})}يمكن إيجاد معكوس متسلسلة طول قوس الزاوية المركزية المذكورة أعلاه في الصفحة 8أ من كتاب راب، المجلد 1، [ 2 ] الذي ينسب الفضل فيه إلى غانشين. [ 3 ] وثمة بديل لاستخدام المتسلسلة المعكوسة وهو استخدام طريقة نيوتن للتقريبات المتتالية لـθ12{\displaystyle \theta _{12}}تُقدّم مسألة خط الزوال العكسي للقطع الناقص معكوس متسلسلة أطوال أقواس بيسل بدلالة الزاوية البارامترية. قبل استخدام المتسلسلة العكسية، يجب استخدام متسلسلة الزاوية البارامترية لحساب طول القوس من نصف المحور الأكبر إلىP1{\displaystyle P_{1}}،s1=s(β1)=أ*+ب*2(ب0β1+ب2الخطيئة2β1+ب4الخطيئة4β1+ب6الخطيئة6β1){\textstyle s_{1}=s(\beta _{1})={\frac {a^{*}+b^{*}}{2}}(B_{0}\beta _{1}+B_{2}\sin 2\beta _{1}+B_{4}\sin 4\beta _{1}+B_{6}\sin 6\beta _{1})}. مرة واحدةs1{\displaystyle s_{1}}من المعروف أنه يمكن تطبيق الصيغة العكسية للحصول علىβ2=β(s1+s12)=μ2+ب2الخطيئة2μ2+ب4الخطيئة4μ2+ب6الخطيئة6μ2{\displaystyle \beta _{2}=\beta (s_{1}+s_{12})=\mu _{2}+B'_{2}\sin 2\mu _{2}+B'_{4}\sin 4\mu _{2}+B'_{6}\sin 6\mu _{2}}، أينμ2=2(s1+s12)/(ب0(أ*+ب*)){\displaystyle \mu _{2}=2(s_{1}+s_{12})/(B_{0}(a^{*}+b^{*}))}الإحداثيات المستطيلة في مستوى المقطع هيx2=أ*كوسβ2،y2=ب*الخطيئةβ2{\displaystyle x_{2}=a^{*}\cos \beta _{2},y_{2}=b^{*}\sin \beta _{2}}لذا، يمكن حساب متجه ECEF باستخدامV2=Rج+(x2أنا^*+y2ج^*){\displaystyle \mathbf {V_{2}} =\mathbf {R_{c}} +(x_{2}\mathbf {{\hat {i}}^{*}} +y_{2}\mathbf {{\hat {j}}^{*}} )}وأخيرًا، احسب الإحداثيات الجغرافية عبرP2=جيهـo(V2){\displaystyle P_{2}=\mathrm {Geo} (V_{2})}باستخدام خوارزمية Bowring لعام 1985، [ 4 ] أو الخوارزمية الموجودة هنا .

السمت

يمكن الحصول على السمت بنفس طريقة حل المسألة غير المباشرة:ت2=شمال^×u^2{\displaystyle \mathbf {t_{2}} =\mathbf {\hat {N}} \times \mathbf {{\hat {u}}_{2}} }وα2=atan2(ت2هـ^2،ت2ن^2){\displaystyle {\alpha _{2}}=\operatorname {atan2} (\mathbf {t_{2}} \cdot \mathbf {{\hat {e}}_{2}} ,\mathbf {t_{2}} \cdot \mathbf {{\hat {n}}_{2}} )}.

أمثلة

يوضح الانحراف الجيوديسي لمختلف المقاطع التي تربط نيويورك بباريس

القطع الناقص العظيم

القطع الناقص العظيم هو المنحنى الناتج عن تقاطع القطع الناقص مع مستوى يمر بمركزه. لذلك، لاستخدام الطريقة المذكورة أعلاه، ما عليك سوى أنR0{\displaystyle R_{0}}يكون الأصل، بحيثV0=R1{\displaystyle \mathbf {V_{0}} =\mathbf {R_{1}} }(متجه الموضع لـR1{\displaystyle R_{1}}تتجنب هذه الطريقة الصيغ المعقدة، والتي قد تكون غامضة أحيانًا، لحساب المثلثات الكروية ، وتوفر بديلًا لصيغ بورينغ. [ 5 ] يُعرف أقصر مسار بين نقطتين على سطح كروي بالمسار الجيوديسي. تُشتق هذه المسارات باستخدام الهندسة التفاضلية . خط الاستواء وخطوط الطول هي أشكال بيضاوية كبيرة تُعتبر أيضًا مسارات جيوديسية [ أ ] . يبلغ أقصى فرق في الطول بين شكل بيضاوي كبير والمسار الجيوديسي المقابل له، والذي يبلغ طوله 5000 ميل بحري، حوالي 10.5 متر. قد يصل الانحراف الجانبي بينهما إلى 3.7 ميل بحري. يكون المقطع العمودي الذي يربط النقطتين أقرب إلى المسار الجيوديسي منه إلى الشكل البيضاوي الكبير، إلا إذا لامس المسار خط الاستواء.

على سطح الإهليلج WGS84 ، النتائج الخاصة بالقوس الإهليلجي الكبير من نيويورك،ϕ1{\displaystyle \phi _{1}}= 40.64130°،λ1{\displaystyle \lambda _{1}}= -73.77810° إلى باريس،ϕ2{\displaystyle \phi _{2}}= 49.00970°،λ2{\displaystyle \lambda _{2}}= 2.54800° هي:

α1{\displaystyle \alpha _{1}}= 53.596810°،α2{\displaystyle \alpha _{2}}= 111.537138° وs12{\displaystyle s_{12}}= 5849159.753 (م) = 3158.293603 (نانومتر). الأرقام المقابلة للخط الجيوديسي هي:

α1{\displaystyle \alpha _{1}}= 53.511007°،α2{\displaystyle \alpha _{2}}= 111.626714° وs12{\displaystyle s_{12}}= 5849157.543 (م) = 3158.292410 (نانومتر).

لتوضيح اعتماد المسألة المباشرة على نوع المقطع، لنفترض أن سمت المغادرة ومسافة الرحلة هما نفس قيم المسار الجيوديسي المذكور أعلاه، ولنستخدم القطع الناقص الكبير لتعريف المسألة المباشرة. في هذه الحالة، تكون نقطة الوصول هي ϕ2{\displaystyle \phi _{2}}= 49.073057°،λ2{\displaystyle \lambda _{2}}= 2.586154°، أي ما يعادل 4.1 ميل بحري تقريبًا من نقطة الوصول في باريس المحددة أعلاه. وبالطبع، باستخدام سمت المغادرة والمسافة من القطع الناقص الكبير في المسألة غير المباشرة، سيتم تحديد موقع الوجهة بدقة.ϕ2{\displaystyle \phi _{2}}= 49.00970°،λ2{\displaystyle \lambda _{2}}= 2.54800°، وسمت الوصولα2{\displaystyle \alpha _{2}}= 111.537138°.

يوضح الانحراف الجيوديسي لمختلف المقاطع التي تربط سيدني ببانكوك

الأقسام الطبيعية

قسم عادي فيP1{\displaystyle P_{1}}يتم تحديده عن طريق السماحV0=u^1{\displaystyle \mathbf {V_{0}} =\mathbf {\hat {u}} _{1}}(العمودي على السطح عندP1{\displaystyle P_{1}}ينتج مقطع عمودي آخر، يُعرف باسم المقطع العمودي المتبادل، من استخدام العمودي على السطح عندP2{\displaystyle P_{2}}ما لم تقع النقطتان على نفس خط العرض أو خط الطول، فإن المقطع العمودي المتبادل سيمثل مسارًا مختلفًا عن المقطع العمودي الأصلي. يوفر هذا النهج بديلًا لنهج آخرين، مثل نهج بورينغ. [ 7 ] وقد تناول بحث ديكين وشيبارد وروس أهمية المقاطع العمودية في المسح، بالإضافة إلى مناقشة معنى مصطلح "الخط" في هذا السياق. [ 8 ]

على سطح القطع الناقص WGS84، نتائج المقطع الطبيعي من نيويورك،ϕ1{\displaystyle \phi _{1}}= 40.64130°،λ1{\displaystyle \lambda _{1}}= -73.77810° إلى باريس،ϕ2{\displaystyle \phi _{2}}= 49.00970°،λ2{\displaystyle \lambda _{2}}= 2.54800° هي:

α1{\displaystyle \alpha _{1}}= 53.521396°،α2{\displaystyle \alpha _{2}}= 111.612516° وs12{\displaystyle s_{12}}= 5849157.595 (م) = 3158.292438 (نانومتر). نتائج المقطع العمودي المقلوب من باريس إلى نيويورك هي:

α1{\displaystyle \alpha _{1}}= 53.509422°،α2{\displaystyle \alpha _{2}}= 111.624483° وs12{\displaystyle s_{12}}= 5849157.545 (م) = 3158.292411 (ميل بحري). [※ هذا خطأ. بالنسبة للمقطع العادي المقلوب، تقع نقطة الانطلاق، باريس، شرق نقطة الوصول، نيويورك. لذلك، يجب أن يتجاوز السمت 180 درجة، وبالتالي فإن سمت الانطلاق البالغ 53.509422 درجة وسمت الوصول البالغ 111.624483 درجة خطأ. يجب أن يكون سمت نقطة الانطلاق 291.624483 درجة، على الرغم من أن المسافة صحيحة]. يبلغ أقصى فرق في الطول بين المقطع العادي والمسار الجيوديسي المقابل الذي يبلغ طوله 5000 ميل بحري حوالي 6.0 أمتار. قد يصل الانحراف الجانبي بينهما إلى 2.8 ميل بحري.

لتوضيح اعتماد المسألة المباشرة على نوع المقطع، لنفترض أن سمت المغادرة ومسافة الرحلة هما نفس قيم الجيوديسية المذكورة أعلاه، ولنستخدم متجه السطح العمودي عند NY لتعريف المسألة المباشرة. في هذه الحالة، تكون نقطة الوصول هي ϕ2{\displaystyle \phi _{2}}= 49.017378°،λ2{\displaystyle \lambda _{2}}= 2.552626°، أي ما يقارب نصف ميل بحري من نقطة الوصول المحددة أعلاه. بالطبع، باستخدام سمت المغادرة والمسافة من المقطع العمودي، ستحدد المسألة غير المباشرة موقع الوجهة في باريس بدقة. من المفترض استخدام المسألة المباشرة عندما تكون نقطة الوصول غير معروفة، ومع ذلك، من الممكن استخدام أي متجه.V0{\displaystyle \mathbf {V_{0}} }يُرضي المرء. على سبيل المثال، باستخدام متجه السطح العمودي في باريس،u^2{\displaystyle \mathbf {\hat {u}} _{2}}، مما ينتج عنه نقطة وصولϕ2{\displaystyle \phi _{2}}= 49.007778°،λ2{\displaystyle \lambda _{2}}= 2.546842°، أي ما يعادل 1/8 ميل بحري تقريبًا من نقطة الوصول المحددة أعلاه. باستخدام العمودي على السطح في ريكيافيك (مع الاستمرار في استخدام سمت المغادرة ومسافة الرحلة للخط الجيوديسي إلى باريس)، ستصل إلى مسافة 347 ميلًا بحريًا تقريبًا من باريس، بينما يجعلك العمودي في زيورخ على مسافة 5.5 ميلًا بحريًا.

أدى البحث عن مقطع أقرب إلى الخط الجيوديسي إلى المثالين التاليين.

يوضح كيف يتغير الانحراف الجيوديسي مع السمت للمقاطع التي تبدأ عند خط عرض 20 درجة.

المقطع الطبيعي المتوسط

المقطع الطبيعي المتوسط ​​منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}يتم تحديده عن طريق السماحV0=0.5(u^1+u^2){\displaystyle \mathbf {V_{0}} =0.5(\mathbf {\hat {u}} _{1}+\mathbf {\hat {u}} _{2})}هذا تقريب جيد للخط الجيوديسي منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}في مجال الطيران أو الإبحار، يبلغ أقصى فرق في الطول بين المقطع العمودي المتوسط ​​والمسار الجيوديسي المقابل له، والذي يبلغ طوله 5000 ميل بحري، حوالي 0.5 متر. أما الانحراف الجانبي بينهما فلا يتجاوز 0.8 ميل بحري. بالنسبة للمسارات التي يبلغ طولها 1000 ميل بحري، يكون خطأ الطول أقل من مليمتر واحد، ويبلغ أسوأ انحراف جانبي حوالي 4.4 متر. وبمتابعة المثال من نيويورك إلى باريس على نظام الإحداثيات WGS84، نحصل على النتائج التالية للمقطع العمودي المتوسط:

α1{\displaystyle \alpha _{1}}= 53.515409°،α2{\displaystyle \alpha _{2}}= 111.618500° وs12{\displaystyle s_{12}}= 5849157.560 (م) = 3158.292419 (نانومتر).

يوضح الانحراف الجيوديسي لمقاطع عادية مختلفة بطول 5000 نانومتر من خط الاستواء.

المقطع الطبيعي عند نقطة المنتصف

المقطع العمودي لنقطة المنتصف منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}يتم تحديده عن طريق السماحV0{\displaystyle \mathbf {V_{0}} }= العمودي على السطح عند نقطة منتصف الخط الجيوديسي منP1{\displaystyle P_{1}}لP2{\displaystyle P_{2}}هذا المسار أقرب قليلاً إلى الخط الجيوديسي من المقطع العمودي المتوسط. يبلغ أقصى فرق في الطول بين المقطع العمودي عند نقطة المنتصف والخط الجيوديسي المقابل له، والذي يبلغ طوله 5000 ميل بحري، حوالي 0.3 متر. أما أسوأ انحراف جانبي بينهما فيبلغ حوالي 0.3 ميل بحري.

إنهاء المثال من نيويورك إلى باريس على نظام WGS84 يعطي النتائج التالية للقسم العمودي لنقطة المنتصف الجيوديسية: α1{\displaystyle \alpha _{1}}= 53.506207°،α2{\displaystyle \alpha _{2}}= 111.627697° وs12{\displaystyle s_{12}}= 5849157.545 (م) = 3158.292411 (نانومتر).

مناقشة

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

يُعدّ الرسم البياني الأول نموذجياً لحالات خطوط العرض المتوسطة حيث يُمثّل القطع الناقص الكبير حالةً شاذة. ويُعتبر المقطع الطبيعي المرتبط بأبعد نقطة عن خط الاستواء خياراً مناسباً لهذه الحالات.

المثال الثاني أطول وهو نموذجي لحالات عبور خط الاستواء، حيث يتجاوز القطع الناقص الكبير المقاطع العادية. ومع ذلك، ينحرف المقطعان العاديان على جانبين متقابلين من الخط الجيوديسي، مما يجعل المقطع العادي المتوسط ​​خيارًا مناسبًا هنا.

يوضح الرسم البياني الثالث كيف تتغير الانحرافات الجيوديسية بتغير السمت الجيوديسي الأولي انطلاقًا من خط عرض 20 درجة شمالًا. يبلغ أسوأ انحراف في المقاطع العادية التي يبلغ طولها 5000 ميل بحري حوالي 2.8 ميل بحري، ويحدث عند سمت جيوديسي أولي قدره 132 درجة من خط عرض 18 درجة شمالًا (48 درجة جنوبًا).

يُظهر المخطط الرابع شكل المخطط الثالث عند الابتعاد عن خط الاستواء. عند خط الاستواء، تزداد التناظرات نظرًا لأن المقاطع عند سمتي 90° و270° هي أيضًا خطوط جيوديسية. ونتيجةً لذلك، يُظهر المخطط الرابع 7 خطوط متميزة فقط من أصل 24 خطًا، تفصل بينها مسافة 15 درجة. تحديدًا، تتطابق الخطوط عند سمتي 15 و75 و195 و255، وكذلك الخطوط عند سمتي 105 و165 و285 و345 على الجانب الآخر باعتبارها الأقرب (باستثناء الخطوط الجيوديسية). أما الخطوط المتطابقة التالية الأبعد عن الخطوط الجيوديسية الأربعة، فتقع عند سمتي 30 و60 و210 و240 على جانب، وعند سمتي 120 و150 و300 و330 على الجانب الآخر. تقع الخطوط الخارجية عند سمت 45 و225 على جانب، و135 و315 على الجانب الآخر. ومع تحرك نقطة الانطلاق شمالًا، لم تعد الخطوط عند سمت 90 و270 خطوطًا جيوديسية، وتتباعد الخطوط المتطابقة الأخرى وتنتشر حتى خط عرض 18 درجة، حيث يبلغ الانحراف أقصى حد له. بعد هذه النقطة، تتقلص الانحرافات على شكل مروحة يابانية مع تقدم نقطة البداية شمالًا. وبذلك، عند خط عرض 84 درجة، يبلغ أقصى انحراف للمقاطع العادية حوالي 0.25 ميل بحري.

يُعد المقطع الطبيعي عند نقطة المنتصف خيارًا جيدًا دائمًا (تقريبًا).

التقاطعات

لنفترض وجود مستويين مقطعيين:شمال^1R=د1{\displaystyle \mathbf {\hat {N}} _{1}\cdot \mathbf {R} =d_{1}}، وشمال^2R=د2{\displaystyle \mathbf {\hat {N}} _{2}\cdot \mathbf {R} =d_{2}}بافتراض أن المستويين غير متوازيين، فإن خط التقاطع يقع على كلا المستويين. وبالتالي فهو متعامد مع كلا العموديين، أي في اتجاهشمال3=شمال^1×شمال^2{\displaystyle \mathbf {N_{3}} =\mathbf {\hat {N}} _{1}\times \mathbf {\hat {N}} _{2}}(لا يوجد سبب للتطبيع)شمال3{\displaystyle \mathbf {N_{3}} }).

منذشمال^1{\displaystyle \mathbf {\hat {N}} _{1}}وشمال^2{\displaystyle \mathbf {\hat {N}} _{2}}ليست على خط مستقيم واحدشمال^1{\displaystyle \mathbf {\hat {N}} _{1}}،شمال^2{\displaystyle \mathbf {\hat {N}} _{2}}،شمال3{\displaystyle \mathbf {N_{3}} }يشكل أساسًا لـR3{\displaystyle \mathbb {R} ^{3}}لذلك، توجد ثوابتج1{\displaystyle C_{1}}وج2{\displaystyle C_{2}}بحيث يكون خط تقاطع المستويين معطى بالعلاقة التاليةR=ج1شمال^1+ج2شمال^2+تشمال3{\displaystyle R=C_{1}\mathbf {\hat {N}} _{1}+C_{2}\mathbf {\hat {N}} _{2}+t\mathbf {N_{3}} }، حيث t هو مُعامل مستقل.

بما أن هذا الخط يقع على كلا مستويي القطع، فإنه يحقق الشرطين التاليين: ج1+ج2(شمال^1شمال^2)=د1{\displaystyle C_{1}+C_{2}(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})=d_{1}}، و ج1(شمال^1شمال^2)+ج2=د2{\displaystyle C_{1}(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})+C_{2}=d_{2}}.

حل هذه المعادلات لـج1{\displaystyle {C_{1}}}وج2{\displaystyle {C_{2}}}أعطِ ج1[1-(شمال^1شمال^2)2]=د1-د2(شمال^1شمال^2){\displaystyle C_{1}[1-(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})^{2}]=d_{1}-d_{2}(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})}، و ج2[1-(شمال^1شمال^2)2]=د2-د1(شمال^1شمال^2){\displaystyle C_{2}[1-(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})^{2}]=d_{2}-d_{1}(\mathbf {\hat {N}} _{1}\cdot \mathbf {\hat {N}} _{2})}.

عرّف " الزاوية ثنائية السطوحν{\displaystyle \nu }، بواسطةكوسν=شمال^1شمال^2{\displaystyle \cos \nu ={\mathbf {\hat {N}} _{1}}\cdot {\mathbf {\hat {N}} _{2}}}. ثمج1=(د1-د2كوسν)الخطيئة2ν{\displaystyle C_{1}={\frac {(d_{1}-d_{2}\cos \nu )}{\sin ^{2}\nu }}}، وج2=(د2-د1كوسν)الخطيئة2ν{\displaystyle C_{2}={\frac {(d_{2}-d_{1}\cos \nu )}{\sin ^{2}\nu }}}.

على خط التقاطع لديناR=R0+تشمال3{\displaystyle \mathbf {R} =\mathbf {R_{0}} +t\mathbf {N_{3}} }، أينR0=ج1شمال^1+ج2شمال^2{\displaystyle \mathbf {R_{0}} =C_{1}\mathbf {\hat {N}} _{1}+C_{2}\mathbf {\hat {N}} _{2}}. لذلك:x=x0+تل3{\displaystyle x=x_{0}+tl_{3}}،y=y0+تم3{\displaystyle y=y_{0}+tm_{3}}، وz=z0+تن3{\displaystyle z=z_{0}+tn_{3}}، أين x0=ج1ل1+ج2ل2{\displaystyle x_{0}=C_{1}l_{1}+C_{2}l_{2}}،y0=ج1م1+ج2م2{\displaystyle y_{0}=C_{1}m_{1}+C_{2}m_{2}}، وz0=ج1ن1+ج2ن2{\displaystyle z_{0}=C_{1}n_{1}+C_{2}n_{2}}، شمال^أنا=(لأنا،مأنا،نأنا){\displaystyle \mathbf {\hat {N}} _{i}=(l_{i},m_{i},n_{i})}، لـ i=1،2، وشمال3=(ل3،م3،ن3){\displaystyle \mathbf {N_{3}} =(l_{3},m_{3},n_{3})}.

لإيجاد نقطة تقاطع هذا الخط مع الأرض، قم بتعويض معادلات الخط فيx2أ2+y2أ2+z2ب2=1{\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{a^{2}}}+{\frac {z^{2}}{b^{2}}}=1}للحصول على أت2+2بت+ج=0{\displaystyle At^{2}+2Bt+C=0}، أينأ=ل32+م32+أ2ب2ن32{\displaystyle A=l_{3}^{2}+m_{3}^{2}+{\frac {a^{2}}{b^{2}}}n_{3}^{2}}، ب=x0ل3+y0م3+أ2ب2z0ن3{\displaystyle B=x_{0}l_{3}+y_{0}m_{3}+{\frac {a^{2}}{b^{2}}}z_{0}n_{3}}، ج=x02+y02+أ2ب2z02-أ2{\displaystyle C=x_{0}^{2}+y_{0}^{2}+{\frac {a^{2}}{b^{2}}}z_{0}^{2}-a^{2}}.

لذلك، يتقاطع الخط مع الأرض عندت=-ب±ب2-أجأ{\displaystyle t={\frac {-B\pm {\sqrt {{B}^{2}-AC}}}{A}}}. لوب2<أج{\displaystyle B^{2}<AC}إذاً، لا يوجد تقاطع.ب2=أج{\displaystyle B^{2}=AC}ثم يكون الخط مماساً للأرض عندت=-ب/أ{\displaystyle t=-B/A}(أي أن الأجزاء تتقاطع عند تلك النقطة الواحدة).

لاحظ ذلكأ0{\displaystyle A\neq 0}منذشمال^1{\displaystyle \mathbf {\hat {N}} _{1}}وشمال^2{\displaystyle \mathbf {\hat {N}} _{2}}ليست على خط مستقيم واحد. بإدخال t في R=R0+تشمال3{\displaystyle \mathbf {R} =\mathbf {R_{0}} +t\mathbf {N_{3}} }، يوضح نقاط تقاطع أجزاء الأرض.

مثال

حدد نقطة تقاطع مقطع من نيويورك إلى باريس مع خط غرينتش. يمكن وصف مستوى خط غرينتش بالصيغة التالية:شمال^=(0،1،0){\displaystyle \mathbf {\hat {N}} =(0,1,0)}ود=0{\displaystyle d=0}والنتائج هي كالتالي:

التقاطعات
قسمخط العرض
الإهليلج العظيم49.634970°
طبيعي49.637377°
متوسط ​​طبيعي49.637568°
متبادل49.637759°
نقطة المنتصف49.637862°

خطوط العرض والطول القصوى

أقصى (أو أدنى) خط عرض هو النقطة التي يتقاطع فيها القطع الناقص مع خط عرضي في نقطة واحدة. ولإعداد المسألة، لنفترض شمال^1=(ل،م،ن){\displaystyle \mathbf {{\hat {N}}_{1}} =(l,m,n)}،د1=د{\displaystyle d_{1}=d} ليكن المقطع المعطى هو المستوى. الموازي هو شمال^2=(0،0،1){\displaystyle \mathbf {{\hat {N}}_{2}} =(0,0,1)}،د2=z0{\displaystyle d_{2}=z_{0}}، أينz0{\displaystyle z_{0}}يجب تحديدها بحيث تكون هناك نقطة تقاطع واحدة فقط. تطبيق طريقة التقاطع المذكورة أعلاه ينتج عنهشمال3=شمال^1×شمال^2=(م،-ل،0){\displaystyle \mathbf {N_{3}} =\mathbf {\hat {N}} _{1}\times \mathbf {\hat {N}} _{2}=(m,-l,0)}، شمال^1شمال^2=ن{\displaystyle {\mathbf {\hat {N}} _{1}}\cdot {\mathbf {\hat {N}} _{2}}=n}، ج1=1ص2(د-نz0){\textstyle C_{1}={\frac {1}{p^{2}}}(d-nz_{0})}، وج2=1ص2(z0-ند){\textstyle C_{2}={\frac {1}{p^{2}}}(z_{0}-nd)}، منذ1-ن2=ل2+م2=ص2{\displaystyle 1-n^{2}=l^{2}+m^{2}=p^{2}}تصبح المعادلات الخطية الناتجةx=x0+تم{\displaystyle x=x_{0}+tm}،y=y0-تل{\displaystyle y=y_{0}-tl}، وz=z0{\displaystyle z=z_{0}}، أين x0=ج1ل{\displaystyle x_{0}=C_{1}l}،y0=ج1م{\displaystyle y_{0}=C_{1}m}، وz0{\displaystyle z_{0}}سيتم تحديدها. معاملات المعادلة التربيعية الناتجة هي أ=م2+ل2=ص2{\displaystyle A=m^{2}+l^{2}=p^{2}}، ب=مx0-لy0=لمج1-لمج1=0{\displaystyle B=mx_{0}-ly_{0}=lmC_{1}-lmC_{1}=0}، ج=x02+y02+أ2ب2z02-أ2=ص2ج12+أ2ب2z02-أ2=1ص2(د-نz0)2+أ2ب2z02-أ2{\displaystyle C=x_{0}^{2}+y_{0}^{2}+{\frac {a^{2}}{b^{2}}}z_{0}^{2}-a^{2}=p^{2}C_{1}^{2}+{\frac {a^{2}}{b^{2}}}z_{0}^{2}-a^{2}={\frac {1}{p^{2}}}(d-nz_{0})^{2}+{\frac {a^{2}}{b^{2}}}z_{0}^{2}-a^{2}}لذلك، سينتج عن التقاطع حل واحد فقط إذاب2=أج{\displaystyle B^{2}=AC}لكن منذ ذلك الحينب=0{\displaystyle B=0}وأ>0{\displaystyle A>0}[ ب ] ، تصبح المعادلة الحرجةج=0{\displaystyle C=0}يمكن إعادة ترتيب هذه المعادلة ووضعها في الشكل التالي:هـz02-2Fz0+جي=0{\displaystyle Ez_{0}^{2}-2Fz_{0}+G=0}، أين هـ=أ2ب2ص2+ن2{\displaystyle E={\frac {a^{2}}{b^{2}}}p^{2}+n^{2}}، F=ند{\displaystyle F=nd}، و جي=د2-أ2ص2{\displaystyle G=d^{2}-a^{2}p^{2}}. لذلك،z0=F±F2-هـجيهـ{\textstyle z_{0}={\frac {F\pm {\sqrt {{F}^{2}-EG}}}{E}}}يُحدد هذا المقياس المسافة من نقطة الأصل إلى المستويات المتوازية المطلوبة. التعويضz0{\displaystyle z_{0}}داخلج1{\displaystyle C_{1}}يعطي القيم لـx0{\displaystyle x_{0}}وy0{\displaystyle y_{0}}تذكر أنت=-ب/أ=0{\displaystyle t=-B/A=0}لذاx=x0{\displaystyle x=x_{0}}،y=y0{\displaystyle y=y_{0}}تمثل هذه الإحداثيات المتبقية للتقاطعات. ويمكن بعد ذلك حساب الإحداثيات الجغرافية باستخدام تحويل ECEF_to_Geo.

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

يتمثل النهج الأبسط في حساب نقاط نهاية المحورين الأصغر والأكبر للقطع الناقص باستخدامR=Rج±ب*ج^*{\displaystyle \mathbf {R} =\mathbf {R_{c}} \pm b^{*}\mathbf {\hat {j}} ^{*}}، و R=Rج±أ*أنا^*{\displaystyle \mathbf {R} =\mathbf {R_{c}} \pm a^{*}\mathbf {\hat {i}} ^{*}}ثم تحويلها إلى إحداثيات جغرافية. ومن الجدير بالذكر هنا أن خط تقاطع مستويين يتكون من مجموعة النقاط الثابتة، ومن ثم محور الدوران، لدوران إحداثي يرسم أحد المستويين على الآخر.

بالنسبة للمثال من نيويورك إلى باريس، فإن النتائج هي:

قسمنقطة المحور الثانوي 1النقطة الثانية على المحور الثانوينقطة المحور الرئيسي 1النقطة الثانية من المحور الرئيسي
الإهليلج العظيمϕ1{\displaystyle \phi _{1}}= 52.418061°،λ1{\displaystyle \lambda _{1}}= -25.123079°ϕ2{\displaystyle \phi _{2}}= -52.418061°،λ2{\displaystyle \lambda _{2}}= 154.876921°ϕ1{\displaystyle \phi _{1}}= 0.000000°،λ1{\displaystyle \lambda _{1}}= 64.876921°ϕ2{\displaystyle \phi _{2}}= 0.000000°،λ2{\displaystyle \lambda _{2}}= -115.123079°
طبيعيϕ1{\displaystyle \phi _{1}}= 52.433790°،λ1{\displaystyle \lambda _{1}}= -25.154863°ϕ2{\displaystyle \phi _{2}}= -52.739188°،λ2{\displaystyle \lambda _{2}}= 154.845137°ϕ1{\displaystyle \phi _{1}}= -0.093365°،λ1{\displaystyle \lambda _{1}}= 64.723898°ϕ2{\displaystyle \phi _{2}}= -0.093365°،λ2{\displaystyle \lambda _{2}}= -115.033623°
متوسط ​​طبيعيϕ1{\displaystyle \phi _{1}}= 52.435039°،λ1{\displaystyle \lambda _{1}}= -25.157380°ϕ2{\displaystyle \phi _{2}}= -52.764681°،λ2{\displaystyle \lambda _{2}}= 154.842620°ϕ1{\displaystyle \phi _{1}}= -0.100746°،λ1{\displaystyle \lambda _{1}}= 64.711732°ϕ2{\displaystyle \phi _{2}}= -0.100746°،λ2{\displaystyle \lambda _{2}}= -115.026491°
متبادلϕ1{\displaystyle \phi _{1}}= 52.436288°،λ1{\displaystyle \lambda _{1}}= -25.159896°ϕ2{\displaystyle \phi _{2}}= -52.790172°,λ2{\displaystyle \lambda _{2}}= 154.840104°ϕ1{\displaystyle \phi _{1}}= -0.108122°،λ1{\displaystyle \lambda _{1}}= 64.699565°ϕ2{\displaystyle \phi _{2}}= -0.108122°،λ2{\displaystyle \lambda _{2}}= -115.019357°
نقطة المنتصفϕ1{\displaystyle \phi _{1}}= 52.436959°،λ1{\displaystyle \lambda _{1}}= -25.161247°ϕ2{\displaystyle \phi _{2}}= -52.803863°،λ2{\displaystyle \lambda _{2}}= 154.838753°ϕ1{\displaystyle \phi _{1}}= -0.112082°،λ1{\displaystyle \lambda _{1}}= 64.693029°ϕ2{\displaystyle \phi _{2}}= -0.112082°،λ2{\displaystyle \lambda _{2}}= -115.015522°

انظر أيضاً

ملحوظات

  1. المسارات الاستوائية هي مسارات جيوديسية حتى نقطة معينة. على سبيل المثال، المسار الجيوديسي الذي يربط نقطتين تفصل بينهما زاوية 180° على خط الاستواء هو مسار زوال يمر فوق أحد القطبين، بينما يظل خط الاستواء نفسه قطعًا ناقصًا كبيرًا. في الواقع، يوجد عدد لا نهائي من القطوع الناقصة الكبيرة في هذه الحالة، اثنان منها فقط مساران جيوديسيان. بالنسبة للأقواس القصيرة، يتطابق المسار الجيوديسي مع القطع الناقص الكبير. إذن، عند أي نقطة يتغير؟ يحسب راب الإجابة لتكون 179° 23' 38.18182" [ 6 ] . عند هذه النقطة، يبدأ المسار الجيوديسي بالابتعاد عن خط الاستواء، وعند زاوية 180° يكون قد وصل إلى أحد القطبين.
  2. وإلا فإن المقطع عبارة عن خط عرضي، لذلك لا يوجد شيء لحله، لأن جميع خطوط العرض متساوية.

مراجع

  1. جيلبرتسون، تشارلز (ربيع 2012). "مسارات مقطع الأرض". الملاحة . 59 (1): 1-7 . doi : 10.1002/navi.2 .
  2. راب، آر إتش (1991)، الجيوديسيا الهندسية، الجزء الأول، جامعة ولاية أوهايو، hdl : 1811/24333
  3. ^ غانشين، في في (1969) [1967]. هندسة الأرض الإهليلجية. تمت الترجمة بواسطة ويليس، جي إم سانت لويس: مركز معلومات ومخططات الطيران. دوى:10.5281/زينودو.32854. OCLC 493553. الترجمة من الروسية لكتاب Геометрия земного эллипсоида (موسكو، 1967)
  4. بورينغ، بي آر (1985). "دقة معادلات خط العرض والارتفاع الجيوديسية". مراجعة المسح . 28 (218): 202-206 . رمز Bibcode : 1985SurRv..28..202B . doi : 10.1179/sre.1985.28.218.202 .
  5. بورينغ، بي آر (1984). "الحلول المباشرة والمعكوسة للخط الإهليلجي العظيم على الإهليلج المرجعي". النشرة الجيوديسية . 58 (1): 101-108 . Bibcode : 1984BGeod..58..101B . doi : 10.1007/BF02521760 . S2CID 123161737 . 
  6. راب، آر إتش (1993)، الجيوديسيا الهندسية، الجزء الثاني، جامعة ولاية أوهايو، hdl : 1811/24409
  7. بورينغ، بي آر (1971). "المقطع الطبيعي - الصيغ الأمامية والعكسية عند أي مسافة". مراجعة المسح . 21 (161): 131-136 . رمز Bibcode : 1971SurRv..21..131B . doi : 10.1179/sre.1971.21.161.131 .
  8. ديكين، ر. إي.؛ شيبارد، س. و.؛ روس، ر. (2011). "إعادة النظر في خط بلاك-آلان" (ملف PDF) . المؤتمر الإقليمي الرابع والعشرون لمسح ولاية فيكتوريا، شيبارتون، 1-3 أبريل 2011. مؤرشف من الأصل (ملف PDF) في 5 يناير 2012. تم الاطلاع عليه في 3 فبراير 2012 .

للمزيد من القراءة