Không thể tính tổng các raster vì các pixel nốt

Tôi có một bộ raster (.vrt) với dữ liệu độ ẩm đất hàng ngày. Tôi muốn tính tổng số pixel của tất cả các trình quét để có số đo theo tháng. Tuy nhiên, toàn bộ thế giới không được bảo hiểm mỗi ngày, điều này dẫn đến các giá trị nốt ở những nơi mà các biện pháp tồn tại cho những ngày khác.

Những gì tôi muốn làm, là tổng các giá trị của mỗi raster. Tuy nhiên, dường như mỗi khi một pixel nốt được tính tổng, pixel kết quả được phân loại trực tiếp thành nốt. Tôi muốn có điều ngược lại: bỏ qua tất cả các giá trị nốt và tóm tắt phần còn lại.

Tôi nghĩ ra 2 cách giải quyết vấn đề:

1. tổng hợp các raster bỏ qua các giá trị nốt
2. chuyển đổi các pixel nốt thành giá trị 0, sau đó tổng hợp tất cả các trình quét

Thật không may, tôi không thể tìm thấy bất kỳ công cụ để làm điều này.

Ai giúp tôi với?

Nếu bạn có quyền truy cập vào ArcGIS thì công cụ Thống kê di động có tùy chọn bỏ qua nốt mà bạn tính tổng \ mean \ min \ max các trình quét của bạn.

Đầu tiên, bạn có thể sử dụng gdal_calc.py để thay đổi tất cả -9999 thành 0 và đặt giá trị NoData thành 0.

Ví dụ:

gdal_calc.py -A input.tif --outfile=input_with_NoData.tif --calc="A+9999*(A==-9999)" --NoDataValue=0

Sau đó, bạn có thể bỏ qua giá trị NoData bằng cách sử dụng gdal_translate với tùy chọn -a_nodata theo sau là không có .

giá trị -a_nodata :

Gán một giá trị nốt được chỉ định cho các dải đầu ra. Bắt đầu với GDAL 1.8.0, có thể được đặt thành không để tránh đặt giá trị nốt cho tệp đầu ra nếu có tồn tại cho tệp nguồn

Thí dụ:

gdal_translate -a_nodata none input_with_NoData.tif output_without_NoData.tif

Trong R:

library(raster)
s <- stack('file1', 'file2', ...)
ss <- sum(s, na.rm=TRUE)

Công cụ có điều kiện trong Phân tích không gian là một cách dễ dàng để chuyển đổi giá trị null thành số không. Sau đó, các khoản tiền nên được một snap.

Tôi đã có cùng một vấn đề trước đây và tôi đã giải quyết nó.

Cũng giống như một lưu ý về việc hiển thị các giá trị nốt trong QGIS: các giá trị nốt luôn được hiển thị dưới dạng giá trị nốt trong cửa sổ bản đồ và khi sử dụng công cụ "thông tin đối tượng", bất kể chúng có thực sự có giá trị -9999, 0, NULL, v.v. Vì vậy, sau khi phân loại lại các giá trị gật đầu thành 0, nếu bạn đi vào cài đặt lớp, bạn sẽ tìm thấy trong Siêu dữ liệu mà gật đầu có giá trị 0.

Cuối cùng tôi đã tìm thấy công cụ hoàn hảo cho nhiệm vụ, GRASS r.series. Với cài đặt "hoạt động tổng hợp" "tối đa" và không được chọn "tuyên truyền NULL", công cụ sẽ tạo một Lớp phù hợp với yêu cầu của bạn.

Tôi nhớ có vấn đề này một thời gian trở lại. Như một giải pháp tôi nghĩ rằng tôi đã thực hiện +1 cho tất cả các ô trong raster, cộng các raster lại với nhau, sau đó thực hiện phép trừ cuối cùng (tuy nhiên có nhiều +1 bạn đã thêm). Nó khá thô nhưng nếu nó hoàn thành công việc ...!

Hoặc, sử dụng 'Phân loại lại các giá trị lưới' trong hộp công cụ Xử lý trong QGIS để chuyển đổi giá trị null của bạn thành 0.