MANUAL GEO Lib
MANUAL-GEO-lib
User Manual: Pdf
Open the PDF directly: View PDF 
.
Page Count: 4
!"#$%&'()(*$+#($,)-%)'$
./ &01234561730$
-879$:7;2<2=$6301<709$<$9>1$3?$?50617309$<7@>4$13$4><:$A718$47??>2>01$6332470<1>$9=91>@9$
<04$706:54>9$9>B>2<:$517:717>9$13$63@C51>$285@;9$<04$47??>2>01$D704$3?$:70>9$6300>6170E$
C37019/$F970E$71$A718$18>$G,%$:7;2<2=$79$<$59>?5:$133:$70$18>$6301>H1$3?$18>$95;I>61/$'>:3A$
79$;27>?$4>9627C1730$3?$><68$?5061730$13E>18>2$A718$18>$C<2<@>1>29$59>4/$
J/ K30B>29730$;>1A>>0$6332470<1>$9=91>@9$$
+FLK&#L$L),"$
%%)J"K"+$
M"NK(&O-&#L$
K30B>219$:<1P:30P8>7E81$Q%%)R$6332470<1>9$13$"<218SK>01>2>4$$"<218S+7H>4$
Q"K"+R$6332470<1>9T$5970E$U!NVW$>::7C9374$
,)-%)'$
N&L-)X&N$
ecef = lla2ecef(llh);
&LOF-$
O)(),"-"(N$
llh = YH.$<22<=$6301<7070E$:<1T$:30ET$8>7E81$:36<1730$Q2<4T$2<4T$59>2$
50719R$
#F-OF-$
O)(),"-"(N$
ecef = YH.$<22<=$A718$HT=TZ$Q59>2$50719R$
+FLK&#L$L),"$
"K"+J%%)$
M"NK(&O-&#L$
K30B>219$ 6<21>97<0$ QHT=TZR$ 6332470<1>9$ 3?$ <$ 2>?>2>06>$ C3701$ Q"K"+R$ 13$ :<1T$
:30T$8>7E81$$Q%%)R$6332470<1>9$70$18>$U!N$VW$9=91>@$
,)-%)'$
N&L-)X&N$
llh = ecef2lla(ecef);
&LOF-$
O)(),"-"(N$
ecef = YH.$<22<2$HT=TZ$Q59>2$50719R$
#F-OF-$
O)(),"-"(N$
llh = YH.$<22<=$:<1T$:30ET$8>7E81$:36<1730$Q2<4T$2<4T$59>2$50719R$
Y/ !2><1$6726:>$<04$285@;$:70>9$63@C51<1730$
Y/./ M72>61$C23;:>@$93:51730$
+FLK&#L$L),"$
#(-#[("KG#L$
M"NK(&O-&#L$
-2<B>:9$ <$ E7B>0$ 4791<06>$ <:30E$ <$ E7B>0$ <Z7@518$ 91<2170E$ <1$ <$ E7B>0$ 70717<:$
C3701T$5970E$18>$U!NSVW$"<218$>::7C9374/$$
(>15209$18>$>04C3701$<04$?70<:$<Z7@51/$
,)-%)'$
N&L-)X&N$
[lat2,lon2,a21] = orto_reckon(lat1,lon1,s,a12);
&LOF-$
O)(),"-"(N$
lat1 = 70717<:$:<17154>$Q4>E2>>9R$
lon1 = 70717<:$:30E7154>$Q4>E2>>9R
s = distance (meters)
a12 = initial azimuth (degrees)
#F-OF-$
O)(),"-"(N$
lat2, lon2 = >04C3701$Q4>E2>>9R$
a21 = 2>B>29>$<Z7@518$Q4>E2>>9RT$<1$?70<:$C3701$?<670E$;<6D$13A<24$
70717<:$C3701
$
+FLK&#L$L),"$
%#X#[("KG#L$
M"NK(&O-&#L$
-2<B>:9$ <$ E7B>0$ 4791<06>$ ?3::3A70E$ <$ 285@;$ :70>T$ 7/>T$ A718$ <$ 63091<01$
<Z7@518/$
(>15209$18>$>04C3701$<04$?70<:$<Z7@51/$
,)-%)'$
N&L-)X&N$
[lat2,lon2,a21] = loxo_reckon(lat1,lon1,s,a12)
&LOF-$
O)(),"-"(N$
lat1 = 7071<:$:<17154>$Q4>E2>>9R$
lon1 = 70717<:$:30E7154>$Q4>E2>>9R$
s = 4791<06>$Q@>1>29R$
a12 = 70717<:$<Z7@518$Q4>E2>>9R
#F-OF-$
O)(),"-"(N$
lat2, lon2 = >04C3701$Q4>E2>>9R$
a21 = 2>B>29>$<Z7@518$Q4>E2>>9RT$<1$?70<:$C3701$?<670E$;<6D$13A<24$
70717<:$C3701
+FLK&#L$L),"$
\&LK"L-*[("KG#L$
M"NK(&O-&#L$
-2<B>:9$ <$ E7B>0$ 4791<06>$ <:30E$ <$ E7B>0$ <Z7@518$ 91<2170E$ <1$ <$ E7B>0$ 70717<:$
C3701T$5970E$18>$U!NSVW$"<218$>::7C9374/$
$(>15209$ 18>$ >04C3701$ <04$ ?70<:$ <Z7@51 within <$ ?>A$ @7::7@>1>29$ 3?$
<6652<6=$$5970E$\706>01=]9$<:E32718@$
,)-%)'$
N&L-)X&N$
[lat2,lon2,a21] = vincenty_reckon(lat1,lon1,s,a12)
&LOF-$
O)(),"-"(N$
lat1 = 7071<:$:<17154>$Q4>E2>>9R$
lon1 = 70717<:$:30E7154>$Q4>E2>>9R
s = 4791<06>$Q@>1>29R
a12 = 70717<:$<Z7@518$Q4>E2>>9R
#F-OF-$
O)(),"-"(N$
lat2, lon2 = >04C3701$Q4>E2>>9R
a21 = 2>B>29>$<Z7@518$Q4>E2>>9RT$<1$?70<:$C3701$?<670E$;<6D$13A<24$
70717<:$C3701
$
Y/J/ &0B>29>$C23;:>@$93:51730$
+FLK&#L$L),"$
#(-#[M&N-)^&$
M"NK(&O-&#L$
K3@C51>9$ 18>$ 4791<06>$ ;>1A>>0$ 1A3$ C37019$ <995@70E$ <$ C>2?>61:=$
9C8>276<:$><218$<04$_<B>2970>$<CC23H7@<1730/$
K3@C51>9$<:93$?32A<24$<Z7@518T$<04$;<6DA<24$<Z7@518/$
,)-%)'$
N&L-)X&N$
s = orto_distazi(lat1,lon1,lat2,lon2);
[s,a12] = orto_distazi(lat1,lon1,lat2,lon2);
[s,a12,a21] = orto_distazi(lat1,lon1,lat2,lon2);
&LOF-$
O)(),"-"(N$
lat1 = !"#M"-&K$:<17154>$3?$?7291$C3701$Q4>E2>>9R
lon1 = :30E7154>$3?$?7291$C3701$Q4>E2>>9R
lat2, lon2 = 9>6304$C3701$Q4>E2>>9R
(70C519$@<=$;>$96<:<29T$B>61329T$32$@<1276>9R
#F-OF-$
O)(),"-"(N$
s = 4791<06>$70$@>1>29
a12 = <Z7@518$70$4>E2>>9$?23@$?7291$C3701$13$9>6304$C3701$Q?32A<24R
a21 = <Z7@518$70$4>E2>>9$?23@$9>6304$C3701$13$?7291$C3701$Q;<6DA<24R
Q)Z7@5189$<2>$70$4>E2>>9$6:36DA79>$?23@$03218/R
+FLK&#L$L),"$
%#X#[M&N-)^&$
M"NK(&O-&#L$
K3@C51>9$ 18>$ 4791<06>$ ;>1A>>0$ 1A3$ C37019$ 18235E8$ <$ :3H3423@76$
Q285@;R$ :70>$ 70$ <$ 9C8>2/$K3@C51>9$ ?32A<24$ <Z7@518T$ <04$ ;<6DA<24$
<Z7@518/$-8>$<Z7@518$2>@<709$63091<01$93$?32A<24$<04$;<6DA<24$@<168/$
,)-%)'$
N&L-)X&N$
s = loxo_distazi(lat1,lon1,lat2,lon2);
[s,a12] = loxo_distazi(lat1,lon1,lat2,lon2);
[s,a12,a21] = loxo_distazi(lat1,lon1,lat2,lon2);
&LOF-$
O)(),"-"(N$
lat1 = !"#M"-&K$:<17154>$3?$?7291$C3701$Q4>E2>>9R
lon1 = :30E7154>$3?$?7291$C3701$Q4>E2>>9R$
lat2, lon2 = s>6304$C3701$Q4>E2>>9R
Q70C519$@<=$;>$96<:<29T$B>61329T$32$@<1276>9R
#F-OF-$
O)(),"-"(N$
s = 4791<06>$70$@>1>29
a12 = aZ7@518$70$4>E2>>9$?23@$?7291$C3701$13 second point
Q?32A<24R
a21 = aZ7@518$70$4>E2>>9$?23@$9>6304$C3701$13$?7291$C3701$Q;<6DA<24R
Q)Z7@5189$<2>$70$4>E2>>9$6:36DA79>$?23@$03218/R
$
$
+FLK&#L$L),"$
\&LK"L-*[M&N-)^&$
M"NK(&O-&#L$
K3@C51>$ 18>$ 4791<06>$ ;>1A>>0$ 1A3$ C37019$ A71870$ <$ ?>A$ @7::7@>1>29$ 3?$
<6652<6=$ 5970E$ 18>$ U!NSVW$ "<218$ >::7C9374/$ $ K3@C51>$ ?32A<24$ <Z7@518T$
<04$ 63@C51>$ ;<6DA<24$ <Z7@518T$ <::$ 5970E$ <$ B>61327Z>4$ B>29730$ 3?$
\706>01=]9$<:E32718@$
,)-%)'$
N&L-)X&N$
s = vincenty_distazi(lat1,lon1,lat2,lon2);
[s,a12] = vincenty_distazi(lat1,lon1,lat2,lon2);
[s,a12,a21] =
vincenty_distazi(lat1,lon1,lat2,lon2);
&LOF-$
O)(),"-"(N$
lat1 = !"#M"-&K$:<17154>$3?$?7291$C3701$Q4>E2>>9R
lon1 = :30E7154>$3?$?7291$C3701$Q4>E2>>9R
lat2, lon2 = 9>6304$C3701$Q4>E2>>9R$
Q70C519$@<=$;>$96<:<29T$B>61329T$32$@<1276>9R
#F-OF-$
O)(),"-"(N$
s = 4791<06>$70$@>1>29$
a12 = <Z7@518$70$4>E2>>9$?23@$?7291$C3701$13$second point
Q?32A<24R
a21 = <Z7@518$70$4>E2>>9$?23@$9>6304$C3701$13$?7291$C3701$Q;<6DA<24R
Q)Z7@5189$<2>$70$4>E2>>9$6:36DA79>$?23@$03218/R
$
$ $
W/ F17:717>9$
$
+FLK&#L$L),"$
N"XJM"K$
M"NK(&O-&#L$
K30B>219$?23@$9>H$4>E2>>9$13$4>67@<:$4>E2>>9$$
,)-%)'$
N&L-)X&N$
x = sex2dec(rfs);
&LOF-$
O)(),"-"(N$
-U#$+#(,)-N$)M,&--"M$
$+#(,)-$.$
+32@<1$$]L`YabJa]]JJ/``c]$$
+32@<1$$]U`JVbY`]]dW/``c]$
+#(,)-$J$
+32@<1$$W`JadaL$P$44@@99$?32$:<17154>$
+32@<1$```d`J`U$P$44@@99$?32$:30E7154>$
#F-OF-$
O)(),"-"(N$
He$$2>95:1$70$4>67@<:$4>E2>>9$$
+FLK&#L$L),"$
M&N-[O#&L-[-#[%&L"$
M"NK(&O-&#L$
(>15209$18>$4791<06>$70$?>>19$;>1A>>0$<$C3701$<04$<$:70>$4>?70>4$;=$1A3$
C37019$$
,)-%)'$
N&L-)X&N$
dist = dist_point_to_line(Lonp1, Latp1, Lonr1,
Latr1, Lonr2, Latr2);
&LOF-$
O)(),"-"(N$
lonp1, Latp1e$%30E7154>$<04$:<17154>$3?$18>$C3701$$
Lonr1, Latr1, Lonr2, Latr2e$%30E7154>$<04$:<17154>$3?$18>$
C37019$2.$<04$2J 18<1$18>=$4>?70>$18>$:70>$
#F-OF-$
O)(),"-"(N$
dist = 70$?>>19$
+FLK&#L$L),"$
M&N-[+%)-$
M"NK(&O-&#L$
(>15209$18>$4791<06>$70$D7:3@>1>29$;>1A>>0$1A3$C37019/ &1$<995@>9$]?:<1$
><218]$97@C:7?76<1730/$ )C23He$
.$$:<17154>$4>E2>>$f$..`$D@$$
.$$:30E7154$4>E2>>9f$...$D@$
,)-%)'$
N&L-)X&N$
dist = dist_flat(lon1, lat1, lon2, lat2);
&LOF-$
O)(),"-"(N$
lon1, lat1 = ?7291$C3701$:30E7154>$<04$:<17154>
lon2, lat2 =$$9>6304$C3701$:30E7154>$<04$:<17154>$
#F-OF-$
O)(),"-"(N$
dist = 70$D@$
$
$
