Come determinare i parametri di proiezione durante la personalizzazione di una proiezione


9

Sto cercando di personalizzare una proiezione Albers e una Hotine Oblique Mercator (HOM) per minimizzare la distorsione nella regione che sto analizzando. La regione si estende da circa 51 a 62 gradi di latitudine, coprendo un'area delle dimensioni dell'Ucraina. La regione è orientata NW - SE.

Voglio assicurarmi di utilizzare i metodi corretti per determinare i due parametri di proiezione: lat / long del centro di proiezione e azimut della linea centrale . Sto usando ArcMap v10. Ecco la procedura che ho seguito finora:

  1. Crea un singolo poligono che definisce la regione di analisi (creando, generalmente, uno scafo convesso attorno all'estensione dei bacini idrografici che coprono la regione). Questo poligono è l'area per la quale sto personalizzando la proiezione.
  2. Proiettato il poligono su Geographic / NAD 83.
  3. Utilizzato gli strumenti di Jeff Jenness per la grafica e le forme
    ( http://www.jennessent.com/arcgis/shapes_graphics.htm ) per determinare il centro di massa del poligono sullo sferoide GRS80. Le coordinate risultanti sono quelle che ho usato per il parametro "centro di proiezione".
  4. Per determinare l'azimut della linea centrale, per prima cosa ho proiettato il poligono su una proiezione equidistante azimutale, specificando il centro di proiezione alle coordinate determinate nel passaggio 3.
  5. Quindi ho disegnato una polilinea (nella proiezione equidistante azimutale), scattata al punto centrale della proiezione, che rappresenta la tendenza direzionale del poligono della regione. Per ottenere l'azimut al centro della proiezione, ho usato gli strumenti per la grafica e le forme di Jeff Jenness per determinare l'azimut iniziale della curva geodetica nel punto centrale.
  6. Per la proiezione di Albers sto usando la longitudine per il centro di proiezione, come determinato al punto 3. Sto anche usando il fantastico foglio di calcolo creato da Bill Huber ( http://forums.esri.com/Attachments/34278.xls ) per determinare dove posizionare i paralleli standard per ridurre al minimo la distorsione della scala all'interno della regione poligonale.

Se necessario, sto usando la versione di HOM di ArcMap che utilizza una linea centrale definita da un punto al centro della proiezione e dal suo angolo di azimut. ESRI chiama questo "Hotine_Oblique_Mercator_Azimuth_Center". In EPSG, credo che questo sia Oblique Mercator, Hotine Variant B, codice metodo EPSG 9815.

Spero che ci siano alcuni esperti di proiezione che possano dirmi se la procedura sopra descritta, in particolare i passaggi 3 e 4 , è un modo corretto per determinare i parametri di proiezione necessari. Sono sulla buona strada? È corretto determinare il centro sullo sferoide e l'angolo del geodetico dal punto centrale (anziché un centro geometrico "2d" e azimut)?

Spero che la descrizione del problema fosse chiara. Non vedo l'ora di ricevere risposte, consigli, discussioni, ecc.!


In qualche modo correlato, concentrandosi su questioni più strategiche: gis.stackexchange.com/questions/2769/…
whuber

Risposte:


11

L'approccio descritto nella domanda mostra un'attenzione eccezionale nella selezione delle proiezioni per una determinata area di studio. Questa risposta mira solo a stabilire una connessione più diretta tra l'obiettivo (di ridurre al minimo la distorsione) e i passi che sono stati e possono essere adottati, in modo che possiamo essere certi che tale approccio avrà successo (sia qui che nelle applicazioni future).

Tipo di distorsione

Aiuta a inquadrare il problema in modo un po 'più chiaro e quantitativo. Quando diciamo "distorsione", possiamo riferirci a diverse cose correlate ma diverse:

  • In ogni punto in cui la proiezione è liscia (cioè non fa parte di una "piega" o unione di due diverse proiezioni e non è sul suo confine o una "lacrima"), c'è una distorsione di scala che generalmente varia con il cuscinetto lontano dal punto. Ci saranno due direzioni opposte in cui la distorsione è maggiore. La distorsione sarà almeno nelle direzioni perpendicolari. Queste sono chiamate le direzioni principali . Possiamo riassumere la distorsione della scala in termini di distorsioni nelle direzioni principali.

  • La distorsione nell'area è il prodotto delle principali distorsioni della scala.

  • Anche le direzioni e gli angoli possono essere distorti. Una proiezione è conforme quando due percorsi sulla terra che si incontrano in un angolo sono mappati su linee garantite per incontrarsi con lo stesso angolo: i progetti conformi preservano gli angoli. Altrimenti, ci sarà una distorsione degli angoli. Questo può essere misurato.

Anche se vorremmo minimizzare tutte queste distorsioni, in pratica ciò non è mai possibile: tutte le proiezioni sono compromessi. Quindi una delle prime cose da fare è dare la priorità: che tipo di distorsione deve essere controllata?

Misurare la distorsione generale

Queste distorsioni variano da punto a punto e, in ciascun punto, spesso variano in base alla direzione. In alcuni casi prevediamo di eseguire calcoli che coprono l'intera regione di interesse: per loro, una buona misura della distorsione complessiva è il valore medio su tutti i punti, in tutte le direzioni. In altri casi è più importante mantenere le distorsioni entro i limiti stabiliti, qualunque cosa accada. Per loro, una misura più appropriata della distorsione generale è la gamma di distorsioni riscontrate in tutta la regione, tenendo conto di tutte le direzioni possibili. Queste due misure possono essere sostanzialmente diverse, quindi è necessario un pensiero per decidere quale sia la migliore.

La scelta di una proiezione è un problema di ottimizzazione

Una volta che abbiamo selezionato un modo per misurare la distorsione e per esprimere il suo valore per l'intera regione di interesse, il problema diventa relativamente semplice: selezionare una proiezione tra quelle supportate dal proprio software e trovare parametri ammissibili per quella proiezione (come la sua centrale meridiano, fattore di scala e così via) che minimizzano la misura complessiva della distorsione.

Nell'applicazione, questo non è facile da eseguire, perché ci sono molte proiezioni possibili, ognuna in genere ha molti parametri che possono essere impostati e se le distorsioni medie sulla regione devono essere ridotte al minimo, dobbiamo anche calcolare quelle medie (quali importi eseguire un'integrazione bidimensionale o tridimensionale ogni volta che si varia un parametro di proiezione). In pratica, quindi, le persone di solito impiegano l'euristica per ottenere una soluzione ottimale approssimativa:

  • Identificare una classe di proiezioni adatta per l'attività. Ad esempio , se la corretta valutazione degli angoli sarà importante, limitarsi a proiezioni conformi (come il HOM). Quando il calcolo di aree o densità è importante, limitarsi alle proiezioni di aree uguali (come gli Albers). Quando è importante mappare i meridiani su linee parallele su e giù, scegliere una proiezione cilindrica. Ecc. Ecc.

  • All'interno di quella classe, concentrati su un piccolo numero noto - attraverso l'esperienza - per essere appropriato per la propria regione di interesse. Questa scelta viene in genere fatta in base all'aspetto della proiezione che potrebbe essere necessario (per la HOM, si tratta di un aspetto "obliquo" o ruotato) e alle dimensioni della regione (in tutto il mondo, un emisfero, un continente o uno più piccolo ). Più grande è la regione, maggiore è la distorsione che devi sopportare. Con regioni di dimensioni più piccole o più piccole, un'attenta selezione di una proiezione diventa sempre meno importante, perché le distorsioni non sono così grandi.

  • Questo ci porta alla domanda attuale: dopo aver selezionato alcune proiezioni, come scegliere i loro parametri? È qui che emergono i primi sforzi per inquadrarlo come problema di ottimizzazione. Selezionare i parametri per ridurre al minimo la misura di distorsione generale scelta. Questo è spesso fatto per tentativi ed errori, usando valori di partenza ragionevolmente intuitivi.

Applicazione pratica

Esaminiamo i passaggi della domanda da questa prospettiva.

1) ( Definizione della regione di interesse. ) È una semplificazione usare lo scafo convesso. Non c'è nulla di importante in questo, ma perché non usare esattamente la regione di interesse? Il GIS può gestirlo.

2 e 3) ( Trovare un centro di proiezione. ) Questo è un ottimo modo per ottenere una stima iniziale del centro, ma - anticipando le fasi successive in cui varieremo i parametri di proiezione - non è necessario essere pignoli su questo. Qualsiasi tipo di centro "bulbo oculare" andrà bene per cominciare.

4 e 5) ( Scelta dell'aspetto. ) Per la proiezione HOM, il problema riguarda come orientarlo. Ricordiamo che la proiezione Mercator standard, nel suo aspetto equatoriale, mappa accuratamente l'Equatore e la sua vicinanza, ma poi aumenta la sua distorsione esponenzialmente con la distanza dall'Equatore. HOM utilizza essenzialmente la stessa proiezione, ma sposta l'Equatore sulla regione di interesse e la ruota. Lo scopo è quello di posizionare la regione equatoriale a bassa distorsione sulla maggior parte della regione di interesse. A causa della crescita esponenziale della distorsione dall'Equatore, minimizzare la distorsione generale ci impone di prestare attenzione alle parti della nostra regione di interesse che si trovano più lontano dalla linea centrale. Pertanto, il nome di questo gioco è trovare una linea (una geodetica sferica) che transita nella regione in modo tale che (a) la maggior parte dell'area sia il più vicino possibile a quella linea (ciò riduce al minimo la distorsione media) o ( b) le parti della regione più lontane da quella linea siano il più vicino possibile (ciò minimizza la massima distorsione).

Un ottimo modo per eseguire questa procedura per tentativi ed errori è quello di indovinare una soluzione e quindi esplorarla rapidamente con un'applicazione interattiva Tissot Indicatrix. (Fare riferimento a questo esempio sul nostro sito. Per i calcoli necessari, consultare https://gis.stackexchange.com/a/5075 .) L'esplorazione si concentra in genere sui punti in cui la proiezione avrà la maggiore distorsione. Il TI non misurerà solo i vari tipi di distorsione - scala, area, angolo, rilevamento - ma rappresenterà anche graficamente tale distorsione. L'immagine vale più di mille parole (e una mezza dozzina di numeri).

6) ( Scelta dei parametri ) Questo passaggio è molto ben fatto: la domanda descrive un modo quantitativo per valutare la distorsione nella proiezione Albers (Conic Equal Area). Con il foglio di calcolo in mano, è semplice regolare i due parallelismi in modo tale da ridurre al minimo la massima distorsione. È un po 'più difficile regolarli per ridurre al minimo la distorsione media in tutta la regione, quindi ciò viene fatto raramente.

Sommario

Definendo la scelta della proiezione come un problema di ottimizzazione, stabiliamo criteri pratici per fare quella scelta con saggezza e difesa. La procedura può essere efficacemente eseguita per tentativi ed errori, il che implica che non è necessaria particolare attenzione per la selezione iniziale dei parametri: esperienza e intuizione sono di solito sufficienti per iniziare bene, e quindi strumenti interattivi come un'app Tissot Indicatrix e un software associato per le distorsioni di calcolo possono aiutare a finire il lavoro.


Grazie. Ho trovato v. Piccoli dettagli pratici su come scegliere il "centro di proiezione" (centro grafico? Centro sferico? Deve essere esattamente al centro? Come determinare le implicazioni di dove si trova?) E come scegliere il linea centrale azimut (quale sarebbe un flusso di lavoro in ArcMap per determinare correttamente l'azimut? come valutare la distorsione associata al posizionamento? ci sono regole rigide da seguire per queste decisioni? come fanno le altre persone?). Penso che tu abbia risposto alla mia domanda con "questo è spesso fatto per tentativi ed errori ..."
fbiles,

1
Definizione della regione di interesse - Questo è un buon suggerimento, usando "esattamente la regione di interesse". Trasforma l'area del progetto dall'aspetto di una bara all'aspetto di un'aragosta ... il che è bello. Ho usato lo scafo convesso perché l'intera metà occidentale dell'area del progetto è un arcipelago. Volevo essere sicuro che il confine del progetto catturasse tutte le piccole isole periferiche e l'area acquatica in mezzo.
fbiles,

Grazie per la spiegazione dell'uso di uno scafo convesso. In effetti, la mia esperienza è stata che molte analisi spaziali alla fine si estendono in qualche modo oltre i confini originali dell'area di studio (se non altro per aiutare a evitare effetti collaterali nelle analisi statistiche), quindi trovo che di solito sia una buona idea definire la regione di interesse come un buffer attorno all'area di studio originale.
whuber

1

Siamo spiacenti, sto pubblicando questo sotto "Risposta". Non sono sicuro che sia appropriato (è troppo lungo per i commenti). Sono nuovo di questo sito ... forse avrei dovuto iniziare una domanda correlata sulla valutazione della distorsione? Ma stavo lavorando a un'idea come risultato di questo post la scorsa settimana, per valutare la distorsione della scala associata alla scelta di diversi valori di centro di proiezione, azimut e fattore di scala per HOM. Ho deciso di pubblicare qui l'idea perché 1) forse sarà uno strumento utile che può essere utilizzato per aiutare a rispondere a parti della domanda originale e 2) Speravo in un feedback sul fatto che questo sembra un approccio ragionevole.

Utilizzando lo stesso concetto del foglio di calcolo creato per valutare la distorsione della scala di Albers, crea un foglio di calcolo riempito con le equazioni di Snyder per HOM (formula ellissoide, "alternativa B", pagina 74 di "Proiezioni della mappa - Un manuale di lavoro"). L'utente inserisce i parametri ellissoidi scelti (a ed e) e i parametri di proiezione "personalizzati" (lat / long del centro di proiezione, azimut della linea centrale, fattore di scala e falso est / nord). Le restanti costanti di proiezione vengono quindi calcolate automaticamente. Il foglio di calcolo contiene anche celle per ogni coppia lat / long (con incrementi di mezzo grado o qualunque incremento desiderato) nell'area di proiezione. Il fattore di scala e le coordinate rettificate su ciascuna coppia lat / long vengono calcolate automaticamente quando si modifica uno qualsiasi dei parametri di proiezione. Adesso, il fattore di scala può essere valutato numericamente 1) calcolando una media complessiva e un intervallo di distorsione della scala nella regione di proiezione, e 2) le coordinate dei punti e i relativi fattori di scala associati possono essere facilmente importati in ArcMap per creare un'immagine visiva di come la scala la distorsione è distribuita. Ovviamente i risultati sono solo un esempio e varieranno a seconda di quante posizioni lat / long vengono valutate, masembra un approccio ragionevole?


+1 È davvero un approccio ragionevole: automatizza ciò che si fa quando si crea una serie di indicatori Tissot in una regione per valutare una potenziale scelta di proiezione.
whuber
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.