Precisione delle coordinate QGIS

In quale (decima) posizione decimale dei valori delle coordinate QGIS definisce un insieme di punti duplicato?

Pensavo che QGIS gestisse il 15 ° decimale; ma era la limitazione solo perché avevo lavorato principalmente su Shapefile.

Tuttavia, il limite non si applica quando scelgo un'altra fonte di dati, come un livello temporaneo di scratch.

Usando un dato fittizio come di seguito, ho fatto un rapido test con due strumenti, per trovare il valore più piccolo in grado di identificare le differenze di coordinate:

• Geoprocessing QGIS: Delete duplicate geometries
• Plugin MMQGIS: Delete Duplicate Geometries

L'output di QGIS Delete duplicate geometriesè lo stesso dei dati di input precedenti (tutti i 20 record sono stati conservati), quindi può significare che QGIS pensa che siano tutti distinti. Questo limite va oltre 1e-29 (o 1e-30) visto in questo piccolo test?

Proprio come un confronto, MMQGIS Delete Duplicate Geometriesprodotto di seguito. Sembra che MMQGIS stabilisca il 16 ° o il 17 ° decimale come limite.

[Modificare]

Temo di non essere stato chiaro quale fosse la parte centrale della mia domanda. Il mio obiettivo è comprendere la relazione tra i valori delle coordinate e i duplicati / sovrapposizioni, come indicato nel primo paragrafo della mia domanda. Speriamo che questo tipo di conoscenza possa aiutarci a controllare facilmente le funzioni sovrapposte modificando le espressioni dei calcolatori sul campo.

Tuttavia, il problema di fondo che avevo cercato di focalizzare era che le posizioni decimali in cui QGIS riconosce i punti (nodi) in quanto la sovrapposizione sembrano variare, a seconda dell'origine dati.

Se utilizziamo Shapefile per il nostro livello, le maniglie QGIS ~ 15 ° decimale e la differenza minore (al 16 ° o 17 °) non vengono riconosciute ... Anche MMQGIS è influenzato da questa limitazione (ai miei occhi). Se creiamo un livello temporalmente scratch o livelli DB questa limitazione va ben oltre il trentesimo? Questo cambiamento nel comportamento mi fa meravigliare.

Lo strumento QGIS Delete Duplicate Geometries (nel Geoalgorithms, Vector General Tools di QGIS 2.18) è uno script Python che opera sugli oggetti geometrici esportati / esposti dal livello Python di QGIS. E sotto questo livello è il livello GEOS.

L'operatore di uguaglianza geometrica del livello GEO accetta un parametro di tolleranza (in unità di distanza della mappa) quando si confrontano due geometrie. Sfortunatamente, questo parametro di tolleranza non è esposto dal layer Python di QGIS - quindi per Delete Duplicate Geometries , due geometrie sono considerate uguali se e solo se tutti i loro valori XY sono totalmente uguali. Ad esempio, POINT (1.000000 1.000000) non è uguale a POINT (1.00000001 1.00000001).

Una lettura interessante su questo argomento è questa domanda frequente: misurare l'accuratezza di latitudine e longitudine?

Non conosco gli strumenti progettati per identificare i duplicati, ma dal punto di vista computazionale è solo una questione di come i numeri sono gestiti dall'architettura di archiviazione delle entità e dal software che li gestisce. Verranno emessi due numeri diversi FLASEpurché la lingua rappresenti i numeri in modo tale da fare un po 'di differenza se il test è se non sono uguali. In molti casi, l'interesse del duplicato è trovare valori uguali, scoprire una sorta di errore di battitura durante la digitalizzazione, non quanto vicini siano.

In altri contesti, l'interesse dei duplicati è di identificare oggetti abbastanza vicini e in termini pratici corrispondenti.

Supponendo che QGIS utilizzi Python come linguaggio principale, la tua domanda può essere indirizzata ai documenti Python e avrai una lettura interessante su come un'aritmetica di base 2 differisce dai nostri libri di matematica elementari: Python: aritmetica a virgola mobile: problemi e limitazioni

Sono un geometra registrato in diversi stati. Quello che vedi con le diverse proiezioni sono linee rette / distanze che provengono dalle Griglie di Lambert o Mercator. È come posizionare un cono sulla terra, quindi la misurazione che stai prendendo viene posizionata sulle linee dei raggi sopra la curvatura della terra. Transverse Mercator è una proiezione cilindrica per gli stati dell'Est Ovest più lunghi come il Tennessee. Usando i sistemi GPS di grado di rilevamento spesso confrontiamo le informazioni con le proiezioni del mondo e persino le proiezioni di stato mentre si trovano sulle antenne. tutte le misurazioni si basano sui parametri integrati nei calcoli in ciascun punto misurato. L'uso dell'attrezzatura di grado GIS che ha un parametro del contatore secondario si adatta molto bene a tratti di uno o due miglia. Dopo circa 6 cifre decimali nelle misurazioni angolari è il massimo che ci si può aspettare dalla precisione. Indipendentemente dal numero di volte in cui si GPS un punto con una precisione di un metro, non si otterrà mai una "Media" o anche una "Modalità" che corrisponderà alle misurazioni di alta precisione. Per riassumere, quando lavoro con le informazioni Antenne o Satellite non mi lascio preoccupare un metro o giù di lì. Oh, ci sono calcoli sferici per le misure del pisolino della terra usando gli angoli theta che calcolano l'angolo a piombo nel campo rispetto all'angolo calcolato in quel punto al centro della terra. Non esiste una vera misurazione. Per dimostrarlo, quando ho iniziato a rilevare nel 1975 il nostro capo ha messo due chiodi in una strada collinare a un quarto di miglio di distanza e abbiamo dovuto incatenarlo e sparare alla distanza ogni giorno per un mese. Con le correzioni atmosferiche sulla catena e sul misuratore di distanza HP, abbiamo abbinato la distanza originale solo 3 volte al misuratore di distanza e 1 volta alla catena. Un esempio molto umiliante e resta con me per tutti questi anni. sono stato fischiato a una convention GIS una volta quando mi hanno annunciato come geometra (lui lui) e non avevo detto una parola, ma so quando indossare il mio cappello da ispettore e quando indossare il mio cappello GIS.