Thống kê Zonal hoạt động chính xác như thế nào?


10

Tôi hơi bối rối về cách các công cụ thống kê khu vực hoạt động trong QGIS (công cụ 'Thống kê khu vực' trong công cụ Raster từ hộp công cụ thuật toán địa lý của QGIS hoặc công cụ 'Thống kê raster từ đa giác' trong Công cụ SAGA Vector <-> Raster ).

Về cơ bản tôi có một lớp bản đồ raster và một số đa giác nhỏ trong một lớp riêng biệt. Tôi muốn gán cho mỗi đa giác giá trị của raster bên dưới. Trong trường hợp một đa giác chồng lên nhiều ô raster, tôi chỉ muốn có mức tối đa.

Tuy nhiên, đầu ra không có ý nghĩa gì và tôi không thể hiểu được lý do của những kết quả đó (rõ ràng không phải là mức tối đa của tế bào raster bên dưới):

nhập mô tả hình ảnh ở đây

Bất cứ ý tưởng nào tôi đang làm sai hoặc những gì tôi chưa xem xét? Các kết quả với API python trong qgis.analysis.QgsZonalStatistics()( theo câu hỏi này ở đây ) tốt hơn nhiều, nhưng nó chỉ cung cấp count, meansum... nhưng vì một số lý do không có maxhoặc min.


1
Chỉ cần thử điều này trong ArcMap bằng cách chuyển đổi Raster của tôi thành Điểm và có vẻ như công cụ Thống kê Zonal của QGIS đang thực hiện tương tự bên trong và chỉ sử dụng các trường hợp một điểm giao nhau với đa giác. Tuy nhiên, dường như không có tùy chọn đối sánh phù hợp trong QGIS (một cái gì đó có thể so sánh với các tùy chọn được cung cấp trong công cụ Spatial Join của ArcMap, ví dụ INTERSECT, WithIN_A_DISTANCE, CLOSEST, ...).
GeoEki

Như một manh mối để tìm ra điều này, tôi đã tính "số đếm" cho một lớp raster thô (khoảng một pixel nằm trong đa giác) với Thống kê Zonal của QGIS và có giá trị 0,45 - đáng ngạc nhiên không phải là số nguyên. Điều này cho thấy rằng QGIS bằng cách nào đó có thể làm trọng số của các giá trị pixel. Đây có thể là lý do tại sao tối đa / phút không tương ứng với những gì bạn mong đợi - có thể là việc thực hiện kém trọng số này khi tính toán tối đa / phút.
Jon

Xem gis.stackexchange.com/a/281753/18189 để biết mô tả về cách hoạt động của thuật toán thống kê khu vực QGIS.
dbaston

Câu trả lời:


4

Cập nhật vào tháng 4 năm 2018 và QGIS 2.18

Phiên bản hiện tại của Plugin Zonal Statistics có thể cung cấp:

Count: to count the number of pixels
Sum: to sum the pixel values
Mean: to get the mean of pixel values
Median: to get the median of pixel values
StDev: to get the standard deviation of pixel values
Min: to get the minimum of pixel values
Max: to get the maximum of pixel values
Range: to get the range (max - min) of pixel values
Minority: to get the less represented pixel value
Majority: to get the most represented pixel value
Variety: to count the number of distinct pixel values

Câu trả lời gốc

Nói tóm lại, cách tiếp cận này với QGIS, theo như tôi biết là không thể xem tối đa pixel raster. Chỉ tính, trung bình và tổng. Tôi đã cung cấp một số tài nguyên bao gồm điều này.

Đây là trang tài liệu cho plugin Zonal Statistics Plugin

Với plugin biểu tượng thống kê Zonal icon_zonal_statistic, bạn có thể phân tích kết quả của phân loại theo chủ đề. Nó cho phép bạn tính toán một số giá trị của các pixel của lớp raster với sự trợ giúp của lớp vectơ đa giác (xem hình_zonal_statistic). Bạn có thể tính tổng, giá trị trung bình và tổng số pixel nằm trong đa giác. Plugin tạo các cột đầu ra trong lớp vectơ với tiền tố do người dùng xác định.

Ngoài ra, Hướng dẫn lấy mẫu dữ liệu raster lấy mẫu này bằng cách sử dụng Điểm hoặc Đa giác cung cấp dữ liệu mẫu để đạt được kết quả thống kê khu vực. Tôi đã bao gồm ảnh chụp màn hình họ cung cấp như một ví dụ về quy trình và kết quả công việc của họ.

Kích hoạt các plugin thống kê Zonal. Đây là một plugin cốt lõi vì vậy nó đã được cài đặt. Xem Sử dụng Plugin để biết cách bật plugin lõi

nhập mô tả hình ảnh ở đây

nhập mô tả hình ảnh ở đây

nhập mô tả hình ảnh ở đây

Khi quá trình xử lý kết thúc, chọn lớp hạt. Sử dụng công cụ Xác định và nhấp vào bất kỳ đa giác quận. Bạn sẽ thấy ba thuộc tính mới được thêm vào lớp: ZS_count, ZS_mean và ZS_sum. Các thuộc tính này chứa số lượng pixel raster, giá trị trung bình của các giá trị pixel raster và tổng các giá trị pixel raster tương ứng. Vì chúng tôi quan tâm đến nhiệt độ trung bình, trường ZS_mean sẽ là trường được sử dụng.

nhập mô tả hình ảnh ở đây


2

Kể từ QGIS 3.0 (có thể sớm hơn, tôi không chắc chắn), các số liệu thống kê sau có sẵn:

  • Đếm
  • Tổng
  • Nghĩa là
  • Trung bình
  • Độ lệch chuẩn
  • Tối thiểu
  • Tối đa
  • Phạm vi
  • Dân tộc thiểu số
  • Đa số (chế độ)
  • Đa dạng
  • Phương sai
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 cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.