Die QGIS-Flächenberechnung unterscheidet sich, wenn die CRS-Transformation im laufenden Betrieb aktiviert ist

Wenn ich QGIS öffne, die Ebene hinzufüge und die Flächen des Shapefiles über den Feldrechner berechne, erhalte ich eine andere Fläche als beim Öffnen von QGIS und aktiviere "Sofortige CRS-Transformation aktivieren" und berechne die Fläche. Dies gilt trotz der Sicherstellung, dass das Projekt und die Ebene dasselbe Koordinatensystem (dieselbe EPSG-Nummer) haben. Was mache ich falsch?

Ich habe ein Shapefile mit Flächenberechnungen, die mit ArcGIS erstellt wurden (nicht ich, die Daten wurden mir übergeben und ich habe keine Ahnung, für welches CRS die Fläche mit ArcGIS berechnet wurde). Die Shapefile-Schicht CRS ist EPSG: 21781 (Schweiz). Wenn ich in QGIS die OTF-Einstellungen nicht ändere und das Projekt CRS als EPSG: 4326 (WGS84) belasse, erhalte ich denselben Wert wie der ArcGIS-Bereichswert. Wenn ich jedoch den OTF ändere, bevor ich die Ebene zu EPSG: 21781 hinzufüge, erhalte ich unterschiedliche Flächenwerte. Soweit ich weiß, deutet dies darauf hin, dass ArcGIS Area mit dem CRS EPSG: 4326 berechnet wurde.

Erster Workflow:

1. Öffnen Sie QGIS
2. Projekt CRS: EPSG 4326
3. Ebene hinzufügen
4. Projekt CRS passt sich automatisch an und ist jetzt EPSG 21781
5. $ Fläche mit Feldrechner berechnen

Zweiter Workflow:

1. Öffnen Sie QGIS
2. Projekt CRS: EPSG 4326
3. Schalten Sie OTF ein und setzen Sie das Projekt CRS auf EPSG 21781
4. Ebene hinzufügen
5. $ Fläche mit Feldrechner berechnen

Schritt 5 des ersten und zweiten Workflows erzeugt NICHT den gleichen Bereich.

EDIT - Haftungsausschluss: Ich möchte die Leser auf die Diskussion mit ChrisW unten verweisen. Es kann sein, dass das Abrufen eines Bereichs, der auf einem OTF-CRS basiert, doch kein Fehler ist. Das heißt, zumindest in Arcgis wird es auch verwendet, um die Geoverarbeitung von zwei Schichten aus verschiedenen CRS zu ermöglichen.

Um auf das obige Problem einzugehen. Wie AndreJ vorgeschlagen und gezeigt hat - dies ist wahrscheinlich ein Fehler in der aktuellen Version von qgis. Es sollte jedoch beachtet werden, dass das Problem nicht das falsche Gebiet ist, sondern dass sich die On-the-Fly-Transformation ohnehin auf Gebietsberechnungen auswirkt.

Der Zweck der spontanen Transformation / Projektion besteht darin, Daten aus verschiedenen Quellen und mit verschiedenen CRS auszurichten. Das ist hauptsächlich für Anzeigezwecke. EG arcmap führt automatisch eine On-the-Fly-Projektion durch, in jedem Fall stimmt ein Layer-CRS nicht mit dem Datenrahmen-CRS überein.

Arcmap bietet auch die Möglichkeit, Daten zu bearbeiten, während sie im laufenden Betrieb projiziert werden. Beachten Sie jedoch auch Folgendes: ( Quelle )

Es ist jedoch wichtig zu beachten, dass bestimmte Bearbeitungsvorgänge abhängig von den verwendeten Koordinatensystemen zu unerwarteten Ausrichtungs- oder Genauigkeitsproblemen führen können.

Bestimmte Bearbeitungsvorgänge, die Probleme verursachen können, umfassen das Ändern der Formen von Features, das Einrasten an der Kante oder Grenze von Features oder das Erweitern und Zuschneiden von Features. Diese Probleme treten eher auf, wenn sich die von Ihnen bearbeiteten Features nahe am Rand oder außerhalb des Verwendungsbereichs des Koordinatensystems befinden

Das heißt: Die spontane Transformation ist weniger genau als nur die Projektierung der Daten auf ein anderes CRS (was auch eigene Probleme mit sich bringt).

Es ist jedoch nicht verwunderlich, dass auf der Grundlage einer On-the-Fly-Transformation eine falsche Fläche berechnet wird. Es ist jedoch überraschend, dass die Tatsache, dass On-the-Fly aktiviert wurde, die Berechnung der Geometrie in irgendeiner Weise beeinflusst, die dies sollte auf den Daten basieren. Daher spielt es keine Rolle, ob die On-the-Fly-Transformation auf demselben oder einem anderen CRS basiert. Die Flächenberechnung sollte jedes Mal identisch sein.

Um praktischer zu sein, wenn Sie das Gebiet berechnen möchten, verwenden Sie es nicht im laufenden Betrieb. Wenn Sie das falsche CRS haben, projizieren Sie Ihre Daten.

Ich kann bestätigen, dass es ein Fehler zu sein scheint.

Erstellen Sie eine CSV-Datei mit folgendem Inhalt:

E N
600000 200000
700000 200000
700000 300000
600000 300000

Importieren Sie es mit EPSG: 21781 als begrenzten Text, aktivieren Sie das Einrasten und zeichnen Sie ein Polygon-Shapefile auf die vier Punkte.

Ohne OTF ergibt sich ein Ergebnis von $area/1000000.010000 m² (was offensichtlich richtig ist).

Drehen OTF auf und die gleiche EPSG Auswahl: 21781, erhalten Sie 9988,2338 m².

Die Wahl eines anderen CRS wie EPSG: 4326 liefert 9990,5339 m², da die Berechnung auf einem anderen Ellipsoid (WGS84 anstelle von Bessel) erfolgt.

Vector --> Geometry Tools --> Export/Add Geometry Columns scheint korrekte Werte zu liefern.

Der Fehler hat bereits einige Tickets: https://issues.qgis.org/issues/10966 und https://issues.qgis.org/issues/12473