PROJ gyorstalpaló¶
A PROJ egy szabványos Unix szűrő, mely földrajzi hosszúság és szélesség koordinátákat derékszögű koordinátákká alakítja (és vissza), ez egy C API a szoftver fejlesztőknek arra, hogy a koordináta transzformációt a saját programjukba illesszék.
PROJ¶
Megjegyzés
Mit fogunk tanulni ? Ebben a fejezetben a PROJ eszköz használatát tanuljuk meg.
Mi a PROJ?¶
A PROJ és az invproj oda illetve vissza transzformációt hajt végre földrajzi és derékszögű koordináták között, a kiválasztható vetületek széles tartományán.
PROJ összegzés¶
proj [ -bcCeEfiIlmorsStTvVwW [ args ] ] [ +args ] file[s]
invproj [ -bcCeEfiIlmorsStTwW [ args ] ] [ +args ] file[s]
Ezeket a parancsokat részletesen leírják PROJ dokumentáció legújabb verziójában.
PROJ példa¶
A következő szkript
proj +proj=utm +lon_0=112w +ellps=clrk66
-r <<EOF
45d15'33.1" 111.5W
45d15.551666667N -111d30
+45.25919444444 111d30'000w
EOF
UTM vetületbe számítást hajt végre a 112Ny hosszúság közlében levő standard UTM közép meridiánnal. A példa földrajzi értékei azonosak és példaként szolgálnak a DMS adatbevitel különböző formáira. Az x-y output adatok három sorban jelennek meg:
460769.27 5011648.45
cs2cs¶
Megjegyzés
Mit fogunk tanulni ? Ebben a fejezetben a cs2cs eszköz használatát tanuljuk meg.
Mi az a cs2cs?¶
A cs2cs transzformációt hajt végre a forrás és a cél koordináta-rendszer között az input pontok halmazán. A koordináta-rendszer transzformáció tartalmazhat egy eltolást a vetületi és földrajzi koordináták között mint ahogy a dátum eltolást.
cs2cs összefoglaló¶
cs2cs [ -eEfIlrstvwW [ args ] ] [ +opts[=arg] ] [ +to [+opts[=arg]] ] file[s]
Ezeket a parancsokat részletesen leírják cs2cs dokumentációjának legújabb verziójában.
cs2cs példa¶
A következő szkript
cs2cs +proj=latlong +datum=NAD83
+to +proj=utm +zone=10 +datum=NAD27
-r <<EOF
45d15'33.1" 111.5W
45d15.551666667N -111d30
+45.25919444444 111d30'000w
EOF
az input NAD83 földrajzi koordinátákat UTM vetületbe transzformálja a kiválasztott 10-es zónába. A példa földrajzi értékei példaként szolgálnak a DMS adatbevitel különböző formáira. Az x-y output adatok három sorban jelennek meg:
1402285.99 5076292.42 0.000
geod¶
Megjegyzés
Mit fogunk tanulni ? Ebben a fejezetben a geod eszköz használatát tanuljuk meg.
Mi a geod?¶
geod (közvetlen) és invgeod (inverz) geodéziai („nagy kör”) számításokat hajt végre a szélesség, a hosszúság és a végpont vissza azimutjának meghatározásához egy adott kezdőpont szélesség, hosszúság, azimut, és távolság (közvetlen) vagy az oda és vissza azimut és a távolság, a kezdő és végpont között, szélesség és hosszúság (inverz). Az eredményekben a kerekítési hiba |f| < 1/50, ahol a f a lapultság.
invgeod nem érhető el minden platformon; ebben az esetben használja a geod-t a -l kapcsolóval.
geod összegzés¶
geod +ellps=<ellipse> [ -afFIlptwW [ args ] ] [ +args ] file[s]
invgeod +ellps=<ellipse> [ -afFIlptwW [ args ] ] [ +args ] file[s]
Ezeket a parancsokat részletesen leírják PROJ dokumentáció legújabb verziójában.
geod példa¶
A következő szkript kiszámítja a geodéziai azimutot és a távolságot U.S történelmi mérföldben Bostonból Portlandbe:
geod +ellps=clrk66 <<EOF -I +units=us-mi
42d15'N 71d07'W 45d31'N 123d41'W
EOF
Mely a következő eredményt adja:
-66d31'50.141" 75d39'13.083" 2587.504
ahol az első két érték az azimut Bostonból Portlandbe, a vissza azimut Portlandból Bostonba, melyet a távolság követ.
Egy példa a geodéziai alkalmazásra, használjuk Boston helyzetét Potland helyzetének meghatározására az azimut és a távolság ismeretében:
geod +ellps=clrk66 <<EOF +units=us-mi
42d15'N 71d07'W -66d31'50.141" 2587.504
EOF
melynek eredménye:
45d31'0.003"N 123d40'59.985"W 75d39'13.094"
Megjegyzés
a távolság érték pontatlansága megfelel Portland helyzeti pontosságának.