Die Kreisflächenberechnung in QGIS ist ungenau

8

Ich habe einen Kreis mit einem Radius von 3 km erstellt, die Fläche mit dem Feldrechner berechnet und ein anderes Ergebnis als die reale Fläche erhalten.

Das QGIS-Ergebnis ist 28273078.8, während die manuell berechnete Fläche 28274333.9 ist. Dies ist auch die gleiche Zahl, die ich bei der Berechnung der Fläche in ArcMap erhalten habe.

Die Projektion "on the fly" ist deaktiviert und das Projekt-CRS ist auf utm32 gesetzt

Irgendwelche Ideen, was dies verursacht, wie man genauere Flächenberechnungen erhält?

qgis arcgis-desktop area A. Jousef
quelle
3
Wie haben Sie den Kreis erstellt, insbesondere wie viele Eckpunkte haben Sie verwendet, um den Kreis zu approximieren?
Jake
Ich habe zwei Methoden ausprobiert: Die erste bestand darin, einen Kreis von einem Punkt mit ungefähr 99 Segmenten zu puffern, und die zweite darin, einen Kreis in ArcMap zu erstellen und in QGIS zu verwenden. Beide Methoden führten zu demselben Problem.
A. Jousef
Es hängt vom Radius ab, aber ich bevorzuge es, mindestens 181 Eckpunkte zu verwenden, wenn der Bereich für mich wichtig ist.
Vince
Siehe meine Antwort in gis.stackexchange.com/questions/168168/… für die Flächengenauigkeit nach Scheitelpunktdichte.
Vince
Vince hat Ihnen bereits die Antwort gegeben, aber ich bin immer noch fasziniert von dem Ergebnis, das QGIS Ihnen gegeben hat. Mit 99 Eckpunkten und einem Radius von 3 km sollte das Ergebnis sein 28255356.207 m2. Können Sie überprüfen, ob der Abstand vom Kreismittelpunkt zu einem der Eckpunkte 3 km beträgt? Vielleicht sind die Eckpunkte tatsächlich etwas weiter entfernt, um eine bessere Annäherung zu erhalten.
Jake

Antworten:

15

Alle neueren Versionen des ArcGIS-Cheats - der zugrunde liegende Formtyp weiß, wann er als Kreis generiert wurde, und meldet eine perfekte Kreisfläche, selbst wenn die Scheitelpunktzahl ihn zu einem Achteck machen würde.

Ich habe gerade einen Testcode mit der ArcSDE-C-API erstellt und die folgenden Ergebnisse erhalten (alle als Kreis ab Punkt {0,100000}in generiert PCS_WGS_1984_UTM_32N):

ArcSDE 10.2 Dynamic ASCII Table Utility  Sat Jan 02 15:37:30 2016
------------------------------------------------------------------------
         1
              Radius: 3000.00
              nVerts: 9
               Shape: Area shape (1 part, 9 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 25455900.00000000
          pct_ACTUAL: 90.03182924
         2
              Radius: 3000.00
              nVerts: 13
               Shape: Area shape (1 part, 13 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 26999987.41125000
          pct_ACTUAL: 95.49292133
         3
              Radius: 3000.00
              nVerts: 19
               Shape: Area shape (1 part, 19 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 27703641.02895000
          pct_ACTUAL: 97.98158692
         4
              Radius: 3000.00
              nVerts: 37
               Shape: Area shape (1 part, 37 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28131014.34700000
          pct_ACTUAL: 99.49311083
         5
              Radius: 3000.00
              nVerts: 91
               Shape: Area shape (1 part, 91 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28251381.25285000
          pct_ACTUAL: 99.91882168
         6
              Radius: 3000.00
              nVerts: 99
               Shape: Area shape (1 part, 99 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28254960.01385000
          pct_ACTUAL: 99.93147896
         7
              Radius: 3000.00
              nVerts: 100
               Shape: Area shape (1 part, 100 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28255353.61820000
          pct_ACTUAL: 99.93287105
         8
              Radius: 3000.00
              nVerts: 181
               Shape: Area shape (1 part, 181 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28268601.72120000
          pct_ACTUAL: 99.97972663
         9
              Radius: 3000.00
              nVerts: 361
               Shape: Area shape (1 part, 361 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28272903.62000000
          pct_ACTUAL: 99.99494148
        10
              Radius: 3000.00
              nVerts: 385
               Shape: Area shape (1 part, 385 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273075.88680000
          pct_ACTUAL: 99.99555075
        11
              Radius: 3000.00
              nVerts: 386
               Shape: Area shape (1 part, 386 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273079.37910000
          pct_ACTUAL: 99.99556310
        12
              Radius: 3000.00
              nVerts: 387
               Shape: Area shape (1 part, 387 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273082.09195000
          pct_ACTUAL: 99.99557270
        13
              Radius: 3000.00
              nVerts: 396
               Shape: Area shape (1 part, 396 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273146.91000000
          pct_ACTUAL: 99.99580194
        14
              Radius: 3000.00
              nVerts: 397
               Shape: Area shape (1 part, 397 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273144.96360000
          pct_ACTUAL: 99.99579506
        15
              Radius: 3000.00
              nVerts: 398
               Shape: Area shape (1 part, 398 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273150.34970000
          pct_ACTUAL: 99.99581411
        16
              Radius: 3000.00
              nVerts: 721
               Shape: Area shape (1 part, 721 vertices)
          PiRSquared: 28274333.88230814
        Area_PERFECT: 28274333.88230814
         Area_ACTUAL: 28273981.39730000
          pct_ACTUAL: 99.99875334

16 rows read in 6.69 ms

Area_PERFECTwurde generiert mit SE_shape_get_area()und Area_ACTUALwurde generiert, indem das Polygon in eine einfache Linienform und zurück in ein Polygon konvertiert wurde, bevor die Flächenfunktion aufgerufen wurde. (Ich musste zur 9.3.1sp2-API zurückkehren, um einen IST-Bereich aus der Bereichsfunktion der Kreisform zu erhalten.)

Ich bin mir nicht ganz sicher, was QGIS tut, da der in der Frage ( 28273078.8 m2) und der in Kommentaren ( 28273147.55 m2) angegebene Bereich nicht genau 99 und 397 Eckpunkten (oder 98, 100, 384, 385, 386,) entspricht. 387, 396 oder 398). Möglicherweise handelt es sich um eine geodätische Berechnung (obwohl ich dies ohne die genauen Koordinatendaten nicht einfach auswerten kann), oder der Wert wurde möglicherweise auf einen 32-Bit-Float berechnet.

Ich denke, die Lektion hier lautet: Wenn Sie sich wirklich für den Geometriebereich interessieren, sollten Sie sicherstellen, dass Kreise viele Eckpunkte enthalten.

Vince
quelle