mapprojectgmt - Online nel cloud

PROGRAMMA:

NOME

mapproject - Esegui trasformazioni di mappe in avanti e inverse, conversioni di datum e geodesia

SINOSSI

mappaprogetto [ con tabelle ] parametri regione [ b|B|f|F|o|O[lon0/lat0]] [ [dx/dy]] [ c|i|p ] [
[datum]] [ [unità]] [ [x0/y0/][[+|-]unità][+|-]] [] [ riga.xy[/[+|-]unità][+]] [
[a|c|g|m]] [ [d|e ] [ ] [ [h]da[/a]] [ [livello]] [ [w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

Nota: Non è consentito alcuno spazio tra il flag dell'opzione e gli argomenti associati.

DESCRIZIONE

mappaprogetto legge (longitudine, latitudine) le posizioni da con tabelle [o input standard] e
calcola le coordinate (x,y) utilizzando la proiezione e le scale della mappa specificate. Facoltativamente,
può leggere le posizioni (x,y) e calcolare i valori (longitudine, latitudine) facendo l'inverso
trasformazione. Questo può essere usato per trasformare punti lineari (x,y) ottenuti dalla digitalizzazione
una mappa di proiezione nota alle coordinate geografiche. Può anche calcolare le distanze lungo
traccia, a un punto fisso o al massimo avvicinamento a una linea. Infine, può essere utilizzato per eseguire
varie conversioni di dati. Sono consentiti campi dati aggiuntivi dopo le prime 2 colonne
che deve avere (longitudine,latitudine) o (x,y). Vedi opzione -: su come leggere
(latitudine, longitudine).

OBBLIGATORIO ARGOMENTI

-Jparametri (Di Più ...)
Seleziona la proiezione della mappa.

-R[unità]xmin/Xmax/ymin/ymax[R] (Di Più ...)
Specificare la regione di interesse. Caso speciale per la proiezione UTM: If -C viene usato
ed -R non è dato allora la regione è impostata per coincidere con la data zona UTM così
in modo da preservare la soluzione ellissoidale completa (vedi RESTRIZIONI per ulteriori informazioni)
informazione).

OPTIONAL ARGOMENTI

tavolo Uno o più ASCII (o binari, vedi -bi[coli][Digitare]) file di tabelle dati contenenti a
numero di colonne di dati. Se non vengono fornite tabelle, leggiamo dallo standard input.

-Lab|Si|f|Fa|o|O[lon0/lat0]
-Af calcola l'azimut (in avanti) dal punto fisso long/lat a ciascun punto dati.
Usa il -Lontano per riportare l'azimut dai punti dati al punto fisso. Utilizzo -Ao ottenere
orientamenti (-90/90) anziché azimut (0/360). Maiuscolo F, B or O volere
convertire da latitudini geodetiche a geocentriche e stimare l'azimut delle geodetiche
(supponendo che l'ellissoide corrente non sia una sfera). Se non viene dato alcun punto fisso allora
calcoliamo l'azimut (o back-azimut) dal punto precedente.

-C[dx/dy]
Imposta il centro delle coordinate proiettate in modo che sia al centro di proiezione della mappa [Il valore predefinito è
angolo in basso a sinistra]. Facoltativamente, aggiungere offset nelle unità proiettate da aggiungere (o
sottratto quando -I è impostato) su (dalle) coordinate proiettate, come false
est e nord per particolari zone di proiezione [0/0]. L'unità utilizzata per il
offsets è l'unità di distanza del tracciato in vigore (vedi PROJ_LENGTH_UNIT) a meno che -F is
utilizzato, nel qual caso gli offset sono in metri.

-Dc|i|p
Sostituisci temporaneamente PROJ_LENGTH_UNIT e usa c (centimetro), i (pollice), o p (Punti)
Invece. Non può essere utilizzato con -F.

-E[datum]
Conversione da geodetica (lon, lat, altezza) a Earth Centered Earth Fixed (ECEF)
(x,y,z) coordinate (add -I per la conversione inversa). Aggiungi ID datum (vedi -Qd)
o dare ellissoide:dx,dy,dz where ellissoide potrebbe essere un ID ellissoide (vedi -Qe) o puoi
dato come a[,*inv_f*], dove a è il semiasse maggiore e inv_f è l'inverso?
appiattimento (0 se omesso). Se datum è - o non dato, assumiamo WGS-84.

-F[unità]
Ridimensionamento forzato 1:1, cioè output (o input, vedi -I) i dati sono nella proiezione effettiva
metri. Per specificare altre unità, aggiungere l'unità desiderata (vedere UNITÀ). Privo di -F,
l'uscita (o l'ingresso, vedi -I) sono nelle unità specificate da PROJ_LENGTH_UNIT (ma
vedere -D).

-G[x0/y0/][[+|-]unità][+|-]
Calcola le distanze lungo il binario or al punto opzionale impostato con -Gx0/a0. Aggiungere
l'unità di distanza (vedi UNITÀ), incluso c (Distanza cartesiana utilizzando l'input
coordinate) o C (Distanza cartesiana utilizzando le coordinate proiettate). Il C unità
richiede -R ed -J da impostare. Senza un punto fisso, calcoliamo il cumulo
distanze lungo il binario. Aggiungere - per ottenere la distanza incrementale tra successive
punti. Aggiungere + per specificare il 2° punto tramite due colonne aggiuntive nel file di input.

-I Eseguire la trasformazione inversa, ovvero ottenere (longitudine, latitudine) dai dati (x,y).

-Lriga.xy[/[+|-]unità][+]
Determinare la distanza più breve dai punti dati di input alle linee fornite in
il file multisegmento ASCII riga.xy. La distanza e le coordinate del
il punto più vicino verrà aggiunto all'output come tre nuove colonne. Aggiungi il
unità di distanza (vedi UNITÀ), incluso c (distanza cartesiana utilizzando le coordinate di input)
or C (Distanza cartesiana utilizzando le coordinate proiettate). Il C l'unità richiede -R ed
-J da impostare. Infine, append + per riportare l'id del segmento di linea e il frazionario
numero del punto invece di lon/lat del punto più vicino.

-N[a|c|g|m]
Converti da latitudini geodetiche (usando l'ellissoide corrente; vedi PROJ_ELLIPSOID)
a una delle quattro diverse latitudini ausiliarie (le longitudini non sono influenzate). Scegliere
da autalico, cinformale, geocentrico, e mlatitudini eridionali [geocentrico]. Utilizzo -I
convertire da latitudini ausiliarie a latitudini geodetiche.

-Q[d|e Elenca tutti i parametri di proiezione. Per elencare solo i dati, usa -Qd. Solo per elencare
ellissoidi, usa -Qe.

-S Sopprime i punti che cadono al di fuori della regione.

-Ns]da[/a]
Conversioni coordinate tra dati da ed a usando lo standard Molodensky
trasformazione. Utilizzo -Ns se la terza colonna di input ha un'altezza sopra l'ellissoide [Predefinito
assume altezza = 0, cioè sull'ellissoide]. Specificare i datum utilizzando l'ID datum (vedi
-Qd) o dare ellissoide:dx,dy,dz where ellissoide potrebbe essere un ID ellissoide (vedi -Qe) o puoi
dato come a[,*inv_f*], dove a è il semiasse maggiore e inv_f è l'inverso?
appiattimento (0 se omesso). Se datum è - o non dato, assumiamo WGS-84. -T può essere
usato insieme a -R -J per cambiare il datum prima della proiezione delle coordinate
(Inserisci -I per applicare la conversione datum dopo la proiezione inversa). Assicurati che
l'impostazione PROJ_ELLIPSOID è corretta per il tuo caso.

-V[livello] (Di Più ...)
Seleziona il livello di verbosità [c].

-W[w|h]
Stampa la larghezza e l'altezza della mappa sull'output standard. Nessun file di input viene letto. Per solo
stampa la larghezza o l'altezza, append w or h, rispettivamente. Le unità di
le dimensioni possono essere modificate tramite -D.

-bi[coli][T] (Di Più ...)
Seleziona input binario nativo. [Il valore predefinito è 2 colonne di input].

-bo[coli][Digitare] (Di Più ...)
Seleziona l'output binario nativo. [Il valore predefinito è lo stesso dell'input].

-d[io|o]nessun dato (Di Più ...)
Sostituisci le colonne di input uguali nessun dato con NaN e fai il contrario in uscita.

-f[io|o]colinfo (Di Più ...)
Specificare i tipi di dati delle colonne di input e/o output.

-g[a]x|y|d|X|Y|D|[col]z[+|-]divario[U] (Di Più ...)
Determina le lacune nei dati e le interruzioni di riga.

-h[io|o][n][+c][+d][+rosservazione][+rtitolo] (Di Più ...)
Salta o produce record di intestazione.

-icols[l] [sscala][ooffset][,...] (Di Più ...)
Seleziona le colonne di input (0 è la prima colonna).

-ocols[,...] (Di Più ...)
Seleziona le colonne di output (0 è la prima colonna).

-p[x|y|z]Azim/elevazione[/livello z][+wlon0/lat0[/z0]][+vx0/y0] (Di Più ...)
Seleziona la vista prospettica.

-S[cols][a|r] (Di Più ...)
Imposta la gestione dei record NaN.

-:[io|o] (Di Più ...)
Scambia la prima e la seconda colonna su input e/o output.

-^ or ad appena -
Stampa un breve messaggio sulla sintassi del comando, quindi esce (NOTA: su Windows
usa solo -).

-+ or ad appena +
Stampa un messaggio di utilizzo esteso (aiuto), inclusa la spiegazione di qualsiasi
opzione specifica del modulo (ma non le opzioni comuni GMT), quindi esce.

-? or no argomenti
Stampare un messaggio di utilizzo completo (aiuto), inclusa la spiegazione delle opzioni, quindi
esce.

--versione
Stampa la versione GMT ed esci.

--show-datadir
Stampa il percorso completo della directory di condivisione GMT ed esci.

UNITA '

Per l'unità di distanza sulla mappa, aggiungere unità d per grado d'arco, m per minuto d'arco, e s per arco
secondo, o e per contatore [Predefinito], f per piede, k per km, M per miglia legali, n per la nautica
miglio, e u per piede di indagine degli Stati Uniti. Per impostazione predefinita calcoliamo tali distanze utilizzando una sfera
approssimazione con cerchi massimi. anteponi - ad una distanza (o l'unità è nessuna distanza è
dato) per eseguire calcoli "Flat Earth" (più veloci ma meno precisi) o anteporre + a
eseguire calcoli geodetici esatti (più lenti ma più accurati).

ASCII FORMATO PRECISIONE

I formati di output ASCII dei dati numerici sono controllati da parametri nel tuo gmt.conf
file. La longitudine e la latitudine sono formattate in base a FORMAT_GEO_OUT, mentre altro
i valori sono formattati secondo FORMAT_FLOAT_OUT. Tieni presente che il formato in vigore può
portare a una perdita di precisione nell'output, che può portare a vari problemi a valle. Se
trovi che l'output non è scritto con sufficiente precisione, considera il passaggio a binario
produzione (-BO se disponibile) o specificare più decimali utilizzando l'impostazione FORMAT_FLOAT_OUT.

ESEMPI

Per convertire le coordinate UTM in metri in posizioni geografiche, dato un file utm.txt e
conoscendo la zona UTM (e la zona o l'emisfero), prova

gmt mapproject utm.txt -Ju+11/1:1 -C -I -F

Per trasformare un file con (longitudine,latitudine) in (x,y) posizioni in cm su un Mercator
griglia per una data scala di 0.5 cm per grado, corsa

gmt mapproject lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

Per trasformare diversi file binari a 2 colonne a doppia precisione con (latitudine, longitudine)
in (x,y) posizioni in pollici su una griglia di Mercatore trasversale (longitudine centrale 75W) per
scale = 1:500000 e sopprimere quei punti che cadrebbero al di fuori dell'area della mappa, eseguire

tracce gmt mapproject.* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

Per convertire le coordinate geodetiche (lon, lat, height) nel file old.dat dal NAD27
CONUS datum (Datum ID 131 che utilizza l'ellissoide Clarke-1866) a WGS 84, eseguire

gmt mapproject old.dat -Th131 > new.dat

Per calcolare la distanza più vicina (in km) tra ogni punto nel file di input quakes.dat
e i segmenti di linea dati nel file ASCII multisegmento coast.xy, run

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

RESTRIZIONI

La regione di input rettangolare impostata con -R sarà in generale mappato in un non rettangolare
griglia. Salvo che -C è impostato, il punto più a sinistra su questa griglia ha xvalue = 0.0, e il
il punto più basso avrà yvalue = 0.0. Quindi, prima di digitalizzare una mappa, esegui l'estremo
le coordinate della mappa attraverso mappaprogetto utilizzando la scala appropriata e vedere quali valori (x,y)
su cui sono mappati. Utilizzare questi valori durante l'impostazione per la digitalizzazione per avere il
la trasformazione inversa funziona correttamente o, in alternativa, utilizzare awk per ridimensionare e spostare il
(x,y) valori prima della trasformazione.

Per alcune proiezioni, può essere utilizzata una soluzione sferica nonostante l'utente abbia selezionato un
ellissoide. Ciò si verifica quando gli utenti -R l'impostazione implica una regione che supera il dominio
in cui valgono gli sviluppi in serie ellissoidali. Queste sono le condizioni: (1)
Conica conforme di Lambert (-JL) e Albers Equal-Area (JB) userà la soluzione sferica
quando la scala della mappa supera 1.0E7. (2) Mercatore trasversale (-JT) e UTM (-JU) sarà
usa la soluzione sferica quando il confine ovest o est è dato in -R è più che
10 gradi dal meridiano centrale, e (3) lo stesso per Cassini (-JC) ma con un limite di
solo 4 gradi.

ELLISSOIDI E SFEROIDI

GMT utilizzerà formule ellissoidali se sono implementate e l'utente ha selezionato un'opzione
ellissoide come forma di riferimento (vedi PROJ_ELLIPSOID). L'utente deve essere a conoscenza di alcuni
potenziali insidie: (1) Per alcune proiezioni, come Transverse Mercator, Albers e
Conica conforme di Lambert usiamo le espressioni ellissoidali quando le aree mappate sono
piccolo, e passare alle espressioni sferiche (e sostituendo l'ausiliario appropriato
latitudini) per mappe più grandi. Le formule ellissoidali sono usate come segue: (a) Trasversale
Mercatore: quando tutti i punti si trovano entro 10 gradi dal meridiano centrale, (b) proiezioni coniche
quando l'intervallo longitudinale è inferiore a 90 gradi, (c) proiezione di Cassini quando tutti i punti
sono entro 4 gradi dal meridiano centrale. (2) Quando stai cercando di abbinare alcuni dati storici
dati (es. coordinate ottenute con una certa proiezione e un certo riferimento
ellissoide) potresti scoprire che GMT fornisce risultati leggermente diversi. uno probabile
La fonte di questa discrepanza è che i calcoli più vecchi usavano spesso cifre meno significative. Per
Ad esempio, gli esempi di Snyder usano spesso l'ellissoide di Clarke 1866 (da lui definito come avente
un appiattimento f = 1/294.98). Da f otteniamo l'eccentricità al quadrato per essere 0.00676862818
(questo è ciò che usa GMT), mentre Snyder arrotonda e usa 0.00676866. Questa differenza può
danno discrepanze di diverse decine di cm. Se è necessario riprodurre le coordinate proiettate
con questa eccentricità leggermente diversa, dovresti specificare il tuo ellissoide con il
stessi parametri di Clarke 1866, ma con f = 1/294.97861076. Inoltre, tieni presente che più vecchio
i dati possono essere riferiti a dati diversi e, a meno che tu non sappia quale dato è stato utilizzato e
converti tutti i dati in un dato comune potresti riscontrare discrepanze da decine a centinaia di
metri. (3) Infine, tieni presente che PROJ_SCALE_FACTOR ha determinati valori predefiniti per alcuni
proiezioni, quindi potrebbe essere necessario sovrascrivere l'impostazione per abbinare i risultati prodotti
con altre impostazioni.

Usa mapprojectgmt online utilizzando i servizi onworks.net