Tính toán diện tích hình tròn trong QGIS là không chính xác

8

Tôi đã tạo một vòng tròn có bán kính 3 km, tính diện tích bằng máy tính trường và nhận được kết quả khác với diện tích thực.

Kết quả QGIS là 28273078.8 trong khi diện tích được tính toán thủ công là 28274333.9 cũng là con số tương tự tôi nhận được khi tính diện tích trong ArcMap.

Phép chiếu "đang hoạt động" bị vô hiệu hóa và CRS dự án được đặt thành utm32

Bất kỳ ý tưởng những gì đang gây ra điều này, làm thế nào để có được tính toán diện tích chính xác hơn?

qgis  arcgis-desktop  area 
— A. Jousef
nguồn
3
Làm thế nào bạn tạo ra vòng tròn, cụ thể là bạn đã sử dụng bao nhiêu đỉnh để xấp xỉ vòng tròn?
— Jake
Tôi đã thử hai phương pháp, phương pháp đầu tiên là đệm một vòng tròn từ một điểm có 99 đoạn gần đúng và cách thứ hai là tạo một vòng tròn trong ArcMap và sử dụng nó trong QGIS. Cả hai phương pháp đều tạo ra cùng một vấn đề.
— A. Jousef
Nó phụ thuộc vào bán kính, nhưng tôi thích sử dụng ít nhất 181 đỉnh nếu khu vực này quan trọng đối với tôi.
— Vince
Xem câu trả lời của tôi trong gis.stackexchange.com/questions/168168/ cho độ chính xác của khu vực theo mật độ đỉnh.
— Vince
Vince đã cho bạn câu trả lời, nhưng tôi vẫn bị thu hút bởi kết quả mà QGIS đã cho bạn. Với 99 đỉnh và bán kính 3 km, kết quả sẽ là 28255356.207 m2. Bạn có thể kiểm tra xem khoảng cách từ tâm vòng tròn đến một trong các đỉnh là 3 km không? Có lẽ các đỉnh thực sự xa hơn một chút để có được xấp xỉ tốt hơn.
— Jake

Câu trả lời:

15

Tất cả các phiên bản gần đây của gian lận ArcGIS - loại hình dạng cơ bản đều biết khi nào nó được tạo thành một vòng tròn và báo cáo một khu vực vòng tròn hoàn hảo, ngay cả khi số lượng đỉnh sẽ biến nó thành một hình bát giác.

Tôi vừa đánh dấu một số mã kiểm tra bằng API ArcSDE 'C' và nhận được các kết quả sau (tất cả được tạo dưới dạng một vòng tròn từ điểm {0,100000}trong 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_PERFECTđược tạo bằng cách sử dụng SE_shape_get_area()và Area_ACTUALđược tạo bằng cách chuyển đổi đa giác thành hình dạng đường đơn giản và trở lại đa giác trước khi gọi hàm diện tích. (Tôi đã phải quay lại API 9.3.1sp2 để có được một khu vực THỰC TẾ từ chức năng khu vực của hình dạng vòng tròn.)

Tôi không chắc chắn những gì QGIS đang làm, vì khu vực được báo cáo trong câu hỏi ( 28273078.8 m2) và khu vực được báo cáo trong các nhận xét ( 28273147.55 m2) không chính xác tương ứng với 99 và 397 đỉnh (hoặc 98, 100, 384, 385, 386, 387, 396 hoặc 398 cho vấn đề đó). Có thể có một phép tính trắc địa liên quan (mặc dù tôi không có cách dễ dàng để đánh giá rằng nếu không có dữ liệu tọa độ chính xác), hoặc giá trị có thể đã được tính cho độ nổi 32 bit.

Tôi đoán bài học ở đây là: Nếu bạn thực sự quan tâm đến khu vực hình học, bạn nên đảm bảo các vòng tròn chứa rất nhiều đỉnh.

— Hoàng tử
nguồn
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookie và Chính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.