Dove posso trovare i dati GIS per testare le operazioni di coordinate GIS?


13

Sto cercando di trovare alcuni dati che possono essere utilizzati per testare le operazioni di coordinate GIS da un CRS a un altro CRS. Finora tutto ciò che ho trovato è il set di dati Gold e anche vari punti di esempio nell'EPSG G7-2 . I dati sull'oro sono OK per testare alcune proiezioni di WGE e i pochi punti nel documento G7-2 sono buoni per testare singole operazioni. Ciò di cui ho bisogno ora sono alcuni punti che possono testare operazioni coordinate coordinate esplicite e concatenate come EPSG: da 4277 a EPSG: 4230 per esempio, di cui avevo trovato un campione in http://www.epsg.org/Exchange /1065.pdf . Ci deve essere una buona fonte là fuori con più di un punto e CRS che usano più riferimenti. Vola, scimmie di Google, vola!

MODIFICARE:

Ho anche trovato qualcosa chiamato GIGS


2
C'è un motivo esplicito per cui non puoi simulare i dati da solo (ovvero perché non puoi creare il tuo gold standard e vedere come il cambiamento dei sistemi di riferimento influisce su questo)?
Andy W,

Uno dei problemi con i test è che se trovo i dati del test è come se fossi l'autorità su ciò che è corretto, il che è tutt'altro che vero. Questo sarà il piano di backup, ma preferirei i dati di un esperto.
Dandy,

Andy W: può essere una buona idea rielaborare il tuo commento in una risposta. Se non succede nulla, questa è sicuramente la risposta corretta.
Dandy,

Non riesco a trovare nulla di veramente utile su Geospatial Integrity of Geoscience Applications (GIGS). Davvero non mi piace quello che ho fatto in osr_ct_test.py . Mi dà solo un'idea di dove cercare possibili problemi, ma ho appena iniziato con 0, 0 quindi i punti spesso non hanno molto senso.
Kurt Schwehr,

Risposte:


5

Per aggiungere al commento / risposta di Andy W:

Potresti (dovresti) generare un set di dati di test con controllo incrociato tramite più API di riproiezione di terze parti diverse, dandoti più sicurezza che non stai solo riproducendo, ad esempio, un bug proj.4.

Inoltre, è possibile creare varie identità, errore in virgola mobile modulo, che è possibile testare in modo casuale - ad esempio, il round trip di CRS A -> CRS B -> CRS A dovrebbe fornire lo stesso punto con cui si è iniziato.


Proj4 e GCTP sono gli unici due che conosco, ce ne sono altri là fuori (non basati su Proj4)?
Dandy,

Potrebbe essere difficile scoprire se usano Proj4 / GCTP rinominati sotto il cofano, ma ESRI spedisce un motore di riproiezione e così fa FME. (Alcuni suggerimenti su Google che ESRI offre più di un motore di riproiezione: forums.esri.com/Thread.asp?c=9&f=1193&t=147790 ) FME ti consente di scegliere il motore di riproiezione: docs.safe.com/fme/2009 /html/FFT/func_reproject.htm
Dan S.

Il sito FME menziona GTRANS che potrebbe essere earth-info.nga.mil/GandG/geotrans Non sono sicuro su cosa si basi, se si basa su qualcosa ma le persone sembrano confrontarlo con Proj4, quindi suppongo che sia un candidato anche.
Dandy,

Penso che sia questo GTRANS, in realtà: translate.google.com/translate?hl=it&sl=sv&u=http://… Non che sia necessariamente un link utile. ;)
Dan S.

Data la mancanza di un set di dati così completo, ho creato il mio utilizzando Proj4 e il database EPSG. Ho testato ogni CRS nel database contro ogni altro CRS. Ho prodotto un rapporto e ho confrontato le differenze, se le differenze sono grandi cerco quella singola operazione.
Dandy,

4

Noi (software sicuro) offriamo un set di dati di esempio da utilizzare in formazione ed esercitazioni. Esistono diversi tipi di set di dati, formati e sistemi di coordinate (principalmente Texas TX83)

Non sono sicuro che sia esattamente quello che stai cercando, ma sei libero di usarlo se ti aiuta.


Ho bisogno di qualcosa in cui il punto n. 1 nel set di dati A in CRS Y rappresenti la stessa posizione fisica del punto n. 1 nel set di dati B in CRS Z. Ciò mi consente di verificare se le operazioni di coordinate utilizzate per trasformare le coordinate tra questi due CRS sono corrette. Analizzerò un po 'i tuoi dati e vedrò se riesco a trovare qualcosa del genere.
Dandy,


1

In generale, PROJ.4 e la sua progenie sono abbastanza buoni, puoi testare le coordinate sul web usando spatialreference.org - la mappa web riporta le coordinate di input e output, sebbene le coordinate di input siano limitate alla geografia.

Se vuoi ricontrollare l'accuratezza e non fidarti delle implementazioni software esistenti, puoi sempre provare a fare i calcoli direttamente - non è troppo oneroso per molte proiezioni. L'articolo di Wikipedia sulle trasformazioni di Helmert è utile, così come le Proiezioni della mappa di John Snyder: un manuale di lavoro ( PDF ).



0

La risposta di Dan S. è quello che stavo pensando originariamente, anche se pensarci un po 'di più su un simile approccio in due fasi non farebbe distinzione tra se si stesse verificando un errore nella prima trasformazione o nella seconda trasformazione.

Spererei ancora che esista un set di dati che desideri, anche se potrebbero esserci sempre delle limitazioni e qualunque sia il tuo obiettivo che potrebbe farti desiderare di generare dati in un modo specifico (oltre all'errore assoluto tra i punti potresti essere interessato a direzionale errore o errore nella distanza tra i punti o errore nella dimensione delle aree).

Che ne dici di questa soluzione, rubando parte dell'approccio di Dan S.;

Diciamo che hai i tuoi dati gold standard in CRS B. Quindi generi dati in CRS A che quando viene trasformato si allinea perfettamente con il tuo gold standard in CRS B (suppongo che tali trasformazioni non abbiano alcun errore stocastico). Quindi puoi trasformare i punti gold standard in CRS B in CRS A e saprai dove dovrebbero trovarsi.

Ciò elimina la possibilità che la trasformazione da CRS A -> CRS B sia la causa dell'errore e che qualsiasi errore sia attribuibile solo alla trasformazione CRS B -> CRS A.

MODIFICARE:

Purtroppo non ho trovato alcun set di dati che soddisfi le tue esigenze. La maggior parte dei documenti sulla precisione della geocodifica che ho nella mia biblioteca utilizzano le stazioni di monitoraggio aereo EPA. Questo documento utilizzava una più ampia varietà di fonti, di cui non credo soddisfacessero le vostre esigenze. Di quelli che potresti voler visitare la pagina web del National Geodetic Survey . Immagino che abbiano la migliore scommessa di avere tali informazioni (delle fonti che ho visto che è).

Buona fortuna, e se trovi qualcosa di nuovo con quella fonte.


Le operazioni di test e il loro inverso possono essere eseguite singolarmente su ciascuna operazione specifica. Ci sono alcuni dati là fuori per questo, che spesso accompagnano le equazioni. Il problema con il test si presenta davvero quando provo a convertire da CRS A a CRS D. Potrebbero esserci diversi modi per farlo, ma uno di questi è migliore. Potrebbe essere più veloce usare un'operazione A -> D ma il percorso corretto potrebbe essere A -> B -> C -> D o A -> C -> D o persino A -> Z -> D. Tutti dovrebbero farti arrivare a distanza di un cratere dalla luna, ma vorrei averlo il più vicino possibile a un cratere causato da una mucca che cade.
Dandy,

Non hai ancora bisogno di valutare da quale trasformazione specifica proviene l'errore? Il mio esempio può essere esteso a più CRS. Mi asterrò dal dare più suggerimenti, dato che sono un po 'fuori dalla mia zona di comfort e ovviamente sei più informato sull'argomento di me. Continuerò a controllare e vedere se alcuni set di dati soddisfano il tuo standard. Penso che alcuni documenti epidemiologici che valutano l'accuratezza della geocodifica che ho letto possano eventualmente soddisfare i tuoi requisiti. È un grande forse però.
Andy W,

La cosa principale è verificare se le singole operazioni vengono utilizzate correttamente anziché se sono esse stesse corrette.
Dandy,

0

Conosco 3 posti in cui è possibile trovare tali dati. I primi 2 sono correlati a Proj4:

1. I file di test di Proj4 : https://github.com/OSGeo/proj.4/tree/master/nad

I file che hanno "test" nei loro nomi sono script che eseguono comandi di trasformazione proj4. L'input è all'interno di questi file di test. I dati di output previsti si trovano nei file corrispondenti con "_out.dist" nei loro nomi. Alla fine dell'input e dell'output viene confrontato con il comando diff.

2. https://github.com/Beman/boost-trunk-git-svn/blob/master/libs/geometry/test_extensions/gis/projections/projections.cpp

Questi sono alcuni dati sintetici. Posso confermare che è corretto con Proj 4.8.0.

3. Il pacchetto di dati Gold di GeoTrans : http://earth-info.nga.mil/GandG/coordsys/Conversion_Software/index.html

Non l'ho mai provato, spero che ti possa essere utile!

Ps Un simpatico strumento web che uso e su cui faccio affidamento: http://cs2cs.mygeodata.eu/ (anche basato su proj4)

In bocca al lupo!

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.