Ten artykuł od 2024-04 wymaga zweryfikowania podanych informacji.

Należy podać wiarygodne źródła w formie przypisów bibliograficznych.
Część lub nawet wszystkie informacje w artykule mogą być nieprawdziwe. Jako pozbawione źródeł mogą zostać zakwestionowane i usunięte.
Sprawdź w źródłach: Encyklopedia PWN • Google Books • Google Scholar • Federacja Bibliotek Cyfrowych • BazHum • BazTech • RCIN • Internet Archive (texts / inlibrary)
Dokładniejsze informacje o tym, co należy poprawić, być może znajdują się w dyskusji tego artykułu.
Po wyeliminowaniu niedoskonałości należy usunąć szablon {{Dopracować}} z tego artykułu.

Algorytm Hirvonena – algorytm służący do transformacji współrzędnych ortokartezjańskich (prostokątnych) x, y, z na współrzędne geodezyjne B, L, h. Jest to proces iteracyjny. W wyniku 3-4-krotnego powtarzania procedury można przeliczyć współrzędne na poziomie dokładności 1 cm.

wychodzimy ze wzorów na x,y,z:

${\displaystyle X=(N+h)\cos B\cos L,}$

${\displaystyle Y=(N+h)\cos B\sin L,}$

${\displaystyle Z=[N(1-e^{2})+h]\sin B.}$

${\displaystyle {\frac {Y}{X}}={\frac {(N+h)\cos B\sin L}{(N+h)\cos B\cos L}}={\text{tg }}L,}$

${\displaystyle L={\text{arctg }}{\frac {Y}{X}}.}$

${\displaystyle {\frac {Z}{\sqrt {X^{2}+Y^{2}}}}={\frac {[N(1-e^{2})+h]\sin B}{(N+H)\cos B}}={\frac {N+Ne^{2}+h}{N+h}}{\text{tg }}B=\left({\frac {N+h}{N+h}}-{\frac {Ne^{2}}{N+h}}\right){\text{tg }}B,}$

${\displaystyle {\frac {Z}{\sqrt {X^{2}+Y^{2}}}}=\left(1-e^{2}{\frac {N}{N+h}}\right){\text{tg }}B,}$

${\displaystyle {\text{tg }}B={\frac {Z}{\sqrt {X^{2}+Y^{2}}}}\left(1-e^{2}{\frac {N}{N+h}}\right)^{-1}.}$

Ostatecznie:

${\displaystyle {\text{tg }}B_{(k+1)}={\frac {Z}{\sqrt {X^{2}+Y^{2}}}}\left(1-e^{2}{\frac {N_{k}}{N_{k}+h_{k}}}\right)^{-1}}$

przy czym:

${\displaystyle {\text{tg }}B^{(k=0)}={\frac {Z}{\sqrt {X^{2}+Y^{2}}}}(1-e^{2})^{-1}.}$

${\displaystyle {\sqrt {X^{2}+Y^{2}}}=(N+h)\cos B,}$

${\displaystyle {\sqrt {X^{2}+Y^{2}}}\cos B=N+h,}$

${\displaystyle h={\frac {\sqrt {X^{2}+Y^{2}}}{\cos B}}-N.}$

• obliczenie długości geodezyjnej L,
• obliczenie szerokości geodezyjnej, B dla k = 0,
• Na podstawie tego B₀ liczymy N₀ i h₀. Następnie mając dane N₀ i h₀ liczymy na ich podstawie B₁. Na podstawie B₁ liczymy N₁ i h₁ itd. Zazwyczaj wystarczają 3 iteracje. Proces jest powtarzany do uzyskania zadowalającej dokładności.