Վինսենտի բանաձև

Կաղապար:Չթարգմանված Վինսենտի բանաձևը դա երկու կապված իտերատիվ մեթոդներ են որոնք օգտագործվում են երկրագիտության մեջ երկու կետերի միջև եղած հեռավորությունը հաշվելու համար սֆերոիդի մակերեսին, նրանք ստեղծվել են Թադեոս Վինսենտի կողմից (1975 թ)։ Նրանք հիմնված էին այն ենթադրության վրա, որ Երկրի ձևը նման է ձգված սֆերոիդի, հետևաբար դա ավելի ձշգրիտ է քան մեծ տարածություն ունեցող շրջանակը որը անվանվում է գնդաձև Երկիր։

Առաջին (ուղիղ)մեթոդը համակարգում է գործունեության վայրի կետը, որը ցույց է տալիս հեռավորությունը և ազիմուտը ուրիշ կետից։ Երկրորդ (հակադարձ) մեթոդը համակարգում է աշխարհագրական հեռավորությունը և ազիմուտը երկու դիտարկված կետերի միջև։ Նրանք արդեն լայնորեն օգտագործվում են գեոդեզիայի մեջ որովհետև նրանք հավասար են 0.5 mm (0.020″) Երկրի էլիպսոյիդում.

Ընդհանուր տեղեկություններ

Վինսենտի նպատակն էր արտահայտել առկա ալգորիթմի համար ուղիղ և հակառակ գեոդեսիկ խնդիրը նվազեցնելով ծրագրային երկարությունը։ Նրա չհրատակված հոդավածը (1975 թ)նշում է Wang 720 desk հաշվիչի օգտագործումը, որը ուներ ընդամենը մի քանի կիլոբայթ հիշողություն։ Ձեռք բերելով լավ երկարություն ունեցող գծեր, լուծումը օգտագործում է Legendre (1806), Bessel (1825)և Helmert (1880) դասական լուծումներ հիմնվելով լրացուցիչ ոլորտի վրա։ (Վինսենտեն ապավինում էր որ այդ մեթոդի ձևակերպումը ստեղծվել է Ռեյնսֆորդի կողմից, 1955 թ) Legendre ցույց է տալիս որ էլիպսակերպ գեոդեսիկը կարող է արտացոլվել մեծ շրջանում լրացուցիչ ոլորտը արտացոլումով աշխարհագրական լայնությունը նվազեցնելով ազիմուտի լայնությունով և կարգավորումով։ Լայնությունը էլիպսաձևի վրա և հեռավորությունը գեոդեսիկի երկայնքով հետո վերցված է ոլորտի երկարությունից սֆերայի վրա և arc երկարությունից մեց շրջանի երկայնքով պարզ ինտեգրալներով։ Bessel և Helmert-ը տրված է արագ համամիտող շարքի կողմից այս ինտեգրալների համար որոնք թույլ են տալիս գեոդեսիկին լինել հաշվարկված կամայական ճշգրտությամբ։

Որպեսզի նվազեցնել ծրագրի չափը, Վինսենտեն վերցրել է այս տեսակները, վերանայել է նրանց օգտագործման առաջին ժամկետը, և բեռնել է նրանց հրամանին ƒ³.: Դա հանգեցրեց կոմպակտ արտահայտություններին, որոնք վերաբերվում էին լայն և հեռավոր ինտեգրալներին։ Արտահայտւթյունները վեռցված էին Horner ձևից, քանի որ դա թույլ էր տալիս պոլինոմիալներին ճիշտ վերլուծվել, օգտագործելով միայն ժամանակավոր գրանցում։ Վերջապես, պարզ իտերատիվ տեխնիկան օգտագործվել է թաքնված խնդիրներ լուծելու համար; չնայած նրանք դանդաղ են (իսկ հակառակ դեպքում, այն մեթոդը երբեմն չի զուգադիմել), նրանք հանգեցնում են նվազագույն չափի ավելացման ծածկագրի։

Նշագրման սիստեմ

Սահմանվում է հետևյալ նշագրման սիստեմ։

a	կիսա -գլխավոր էլիպսակերպի երկարություն (շառավիղը հասարակածում);	(6378137.0 մետր in WGS-84)
ƒ	շտկող հաստոցներ էլիփսակերպում;	(1/298.257223563 in WGS-84)
b = (1 - ƒ) a	կիսա -գլխավոր էլիպսակերպի երկարություն (շառավիղը բևռներում);
φ1, φ2	բևեռների լայնությունը;
U1 = arctan[(1 − ƒ) tan φ1], U2 = arctan[(1 − ƒ) tan φ2]	reduced latitude (լայնության օժանդակ ոլորտ վրա)
L = L2 - L1	լայնության տարբերությունը երկու կետերում;
λ1, λ2	երկայնությունը կետերի օժանդակ ոլորտում;
α1, α2	փոխանցված ազիմուտ կետերում;
α	ազիմուտ հասարակածում;
s	էլիփսակերպ տարածություն երկու կետերի միջև;
σ	arc երկարությունը կետերի միջև օժանդակ ոլորտում;

Հակադարձ խնդիրներ

Վերցնելով երկու կետերի կոորդինատները (φ₁, L₁) և (φ₂, L₂), Հակադարձ խնդիրը կգտնի ազիմուտները α₁, α₂ և էլիփսակերպ տարածությունը s.

Հաշվենք U₁, U₂ և L, սահմանելով սկզբնական արժեքը λ = L. Այնուհետև իտերատիվ գնահատումը հետևյալ հավասարման համար, մինչև λ :

$\sin \sigma =\sqrt{(\cos U_2\sin \lambda {)}^2+(\cos U_1\sin U_2-\sin U_1\cos U_2\cos \lambda {)}^2}$

$\cos \sigma =\sin U_1\sin U_2+\cos U_1\cos U_2\cos \lambda$

$\sigma =\arctan \frac{\sin \sigma }{\cos \sigma }$

$\sin \alpha =\frac{\cos U_1\cos U_2\sin \lambda }{\sin \sigma }$

${\cos }^2\alpha =1-{\sin }^2\alpha$

$\cos (2{\sigma }_m)=\cos \sigma -\frac{2\sin U_1\sin U_2}{{\cos }^2\alpha }$

$C=\frac{f}{16}{\cos }^2\alpha [4+f(4-3{\cos }^2\alpha )]$

$\lambda =L+(1-C)f\sin \alpha \{\sigma +C\sin \sigma [\cos (2{\sigma }_m)+C\cos \sigma (-1+2{\cos }^2(2{\sigma }_m))]\}$

$u^2={\cos }^2\alpha \frac{a^2-b^2}{b^2}$

$A=1+\frac{u^2}{16384}\{4096+u^2[-768+u^2(320-175u^2)]\}$

$B=\frac{u^2}{1024}\{256+u^2[-128+u^2(74-47u^2)]\}$

$\mathrm{\Delta }\sigma =B\sin \sigma \{\cos (2{\sigma }_m)+\frac{1}{4}B[\cos \sigma (-1+2{\cos }^2(2{\sigma }_m))-\frac{1}{6}B\cos (2{\sigma }_m)(-3+4{\sin }^2\sigma )(-3+4{\cos }^2(2{\sigma }_m))\left]\right\}$

$s=bA(\sigma -\mathrm{\Delta }\sigma )$

${\alpha }_1=\arctan \left(\frac{\cos U_2\sin \lambda }{\cos U_1\sin U_2-\sin U_1\cos U_2\cos \lambda }\right)$

${\alpha }_2=\arctan \left(\frac{\cos U_1\sin \lambda }{-\sin U_1\cos U_2+\cos U_1\sin U_2\cos \lambda }\right)$

Երկու մոտ տրամագծորեն հակառակ կետերի միջև, իտերատիվ բանաձևը չի կարող զուգամիտել։

Ուղիղ խնդիր

Սկսում ենք հաշվարկման հետևյալ տարբերակը։

$\tan U_1=(1-f)\tan {\varphi }_1$

${\sigma }_1=\arctan \left(\frac{\tan U_1}{\cos {\alpha }_1}\right)$

$\sin \alpha =\cos U_1\sin {\alpha }_1;{\cos }^2\alpha =(1-\sin \alpha )(1+\sin \alpha )$

$u^2={\cos }^2\alpha \frac{a^2-b^2}{b^2}$

$A=1+\frac{u^2}{16384}\{4096+u^2[-768+u^2(320-175u^2)]\}$

$B=\frac{u^2}{1024}\{256+u^2[-128+u^2(74-47u^2)]\}$

Այնուհետև, օգտագործում ենք սկզբնական արժեքը $\sigma =\frac{s}{bA}$:

$2{\sigma }_m=2{\sigma }_1+\sigma$

$\mathrm{\Delta }\sigma =B\sin \sigma \{\cos (2{\sigma }_m)+\frac{1}{4}B[\cos \sigma (-1+2{\cos }^2(2{\sigma }_m))-\frac{1}{6}B\cos (2{\sigma }_m)(-3+4{\sin }^2\sigma )(-3+4{\cos }^2(2{\sigma }_m))\left]\right\}$

$\sigma =\frac{s}{bA}+\mathrm{\Delta }\sigma$

Մի անգամ σ-ն օգտագործվեց որպես ճշտության գնահատող։

${\varphi }_2=\arctan \left(\frac{\sin U_1\cos \sigma +\cos U_1\sin \sigma \cos {\alpha }_1}{(1-f)\sqrt{{\sin }^2\alpha +(\sin U_1\sin \sigma -\cos U_1\cos \sigma \cos {\alpha }_1{)}^2}}\right)$

$\lambda =\arctan \left(\frac{\sin \sigma \sin {\alpha }_1}{\cos U_1\cos \sigma -\sin U_1\sin \sigma \cos {\alpha }_1}\right)$

$C=\frac{f}{16}{\cos }^2\alpha [4+f(4-3{\cos }^2\alpha )]$

$L=\lambda -(1-C)f\sin \alpha \{\sigma +C\sin \sigma [\cos (2{\sigma }_m)+C\cos \sigma (-1+2{\cos }^2(2{\sigma }_m))]\}$

${\alpha }_2=\arctan \left(\frac{\sin \alpha }{-\sin U_1\sin \sigma +\cos U_1\cos \sigma \cos {\alpha }_1}\right)$

Եթե նախնական կոտը գտնվում Հյուսիսային կամ Հարավային բևեռում, ապա առաջին հավասարումը անորոշ է։ Եթե նախնական ազիմուտը ուղղված է Արևմուտք կամ Արևելք, ապա երկրորդ հավասարումը անորոշ է։ Եթե երկարժեք atan2 ֆունկցիայի տեսակն է օգտագործվում, ապա այս արժեքները սովորաբար գործածվում են ճիշտ։

Վինսենթի ձևափոխությունը

Այս նամակում, ուղղված Survey Review-ին 1976 թ, Վինսեթը առաջարկեց փոխարինել իր արտահայտությունները A-ով և B-ով, մի պարզ բանաձևով, օգտագործելով Հելմերթի ընդարձակ պարամետրը՝ k₁։

$A=\frac{1+\frac{1}{4}(k_1{)}^2}{1-k_1}$

$B=k_1(1-\frac{3}{8}(k_1{)}^2)$

որտեղ Կաղապար:Pad$k_1=\frac{\sqrt{(1+u^2)}-1}{\sqrt{(1+u^2)}+1}$

Մոտավոր տրամագծորեն հակառակ կետերը

Ինչպես վերը նշվեց, հակառակ պրոբլեմի ինտերատիվ լուծումը դառնում է զուգադիմված կամ դանդաղորեն զուգադիմում է մոտավոր տրամագծորեն հակառակ կետերին։ Դանդաղ զուգադիպման օրինակ է հանդիսանում (φ₁, L₁) = (0°, 0°) և (φ₂, L₂) = (0.5°, 179.5°) WGS84 էլիփսոիդի համար։ Սա պահանջում է մոտավորապես 130 իտերատիաներ, 1 mm արդյունքի հասնելուն համար։ Կախված նրանից, թե ինչպես է կարգավորված մեթոդը կիրառվում ալգորիթմը կարող է վերադարձնել ճիշտ արդյունքի, (19936288.579 m), ոչ ճիշտ արդյունքի կամ սխալ ինդիկատորի։ Ոչ ճիշտ արդյունքի օրինակը ներկայացվում է NGS online utility-ի կողմից, որը բերում է տարածության մոտ 5 km, որը շատ երկար է։ Այս դեպքերում Վինսենթը առաջարկում է զուգադիպման արագացնող մեթոդը (Rapp, 1973)։

Զուգատիպման կարգավորվող մեթոդի ձախողումը (φ₁, L₁) = (0°, 0°) է և (φ₂, L₂) = (0.5°, 179.7°), WGS84 էլիթսոիդի համար։ Մի ոչ տպագրված հաշվետվության մեջ Վինսենթը (1975b) ներկայացրեց իտերատիվ սխեմայի այլընտրանքային տարբերակ՝ այսպիսի դեպքերը կարգավորելու համար։ Սա բերում է ճիշտ արդյունքի 19944127.421 m, մոտավորապես 60 իտերատիայի; այնուամենայնիվ ուրիշ դեպքերում հազարավոր իտերացիաներ կարող են գոյություն ունենալ։

Նյութոնի մեթոդը հաջողությամբ օգտագործվեց արագ զուգադիպումների ժամանակ բոլոր զույգերի ներմուծած միավորների համար (Կարնի, 2011).

Տեղեկագրություն

• Կաղապար:Cite journal
• Կաղապար:Cite book
• Կաղապար:Cite tech report
• Կաղապար:Cite journal
• Կաղապար:Cite journal
• Կաղապար:Cite tech report
• Կաղապար:Cite journal
• Կաղապար:Cite journal
• Կաղապար:Cite tech report
• Կաղապար:Cite book

Արտաքին հղումներ

• Online calculators from Geoscience Australia:
  • Vincenty Direct (destination point)
  • Vincenty Inverse (distance between points)
• Calculators from the U.S. National Geodetic Survey:
  • Online and downloadable PC-executable calculation utilities, including forward (direct) and inverse problems, in both two and three dimensions (accessed 2011-08-01).
• Online calculators with JavaScript source code by Chris Veness (Creative Commons Attribution license)։
  • Vincenty Direct (destination point)
  • Vincenty Inverse (distance between points)
• GeographicLib provides a utility Geod (with MIT/X11 licensed source code) for solving direct and inverse geodesic problems. Compared to Vincenty, this is about 1000 times more accurate (error = 15 nm) and the inverse solution is complete. Here is an online version of Geod.