Xác định lại giá trị nốt bằng 0 trong QGIS?

Tôi có một hình ảnh mà tôi đã cắt bằng đa giác AOI. Vùng đường viền được đặt thành pixel value = 0 và tôi có giá trị nốt là -32768. Tôi muốn xác định lại giá trị nốt bằng không. Tôi đã từng làm điều này với ArcMap một cách dễ dàng trong máy tính raster, nhưng trong QGIS tôi không thấy làm thế nào tôi có thể làm điều đó. Có ý kiến ​​gì không?

Trong QGIS, bạn có thể sử dụng Máy tính Raster với phép tính sau:

("your_raster" != -32768) * "your_raster"

Với phép tính này, nếu giá trị ô là -32768, bạn sẽ nhận được 0 trong ô đó và nếu nó khác với -32768, ô sẽ giữ giá trị mà nó có.

Giải pháp đơn giản nhất một bước và, IMO, đáng tin cậy nhất để phân loại lại NoData thành 0 là sử dụng công cụ Reclassify Grid Values ​​(SAGA) trong hộp công cụ xử lý.

Sau khi chọn raster sẽ được phân loại lại, chỉ cần cuộn xuống dưới cùng của hộp thoại, đảm bảo replace no data valueshộp được chọn. Giá trị mặc định là 0 (nhưng dĩ nhiên bạn có thể thay thế thứ khác). Sau đó bỏ chọn replace other valueshộp và để tất cả các trường khác một mình (vì vậy tất cả những gì bạn đang làm là đảm bảo một hộp được chọn và bỏ chọn một hộp khác - không yêu cầu cài đặt nào khác).

Tôi tìm thấy một giải pháp cho việc này. Đầu tiên tôi chuyển đổi hình ảnh thành một lớp vrt (gdalbuildvrt) với -srcnodata 0, sau đó tôi dịch nó trở lại thành một tệp tif với gdal_translate -a_nodata. Làm việc tốt :)

Tôi đoán bạn có thể làm điều đó với máy tính raster .

Nếu không có sự trợ giúp của plugin GRASS và r.null bạn sẽ làm điều đó một cách dễ dàng. Tùy chọn đó có nghĩa là trước tiên bạn sẽ tạo một bộ bản đồ GRASS để nhập raster của bạn vào.

Có cùng một vấn đề, cuối cùng tôi đã sử dụng Python trực tiếp - bạn có thể phải điều chỉnh numpy.where cho mục đích cụ thể của mình. Trong trường hợp bên dưới, các giá trị pixel được giữ nguyên nếu chúng là> = 0, tất cả các pixel khác - trong trường hợp này chỉ là các giá trị có giá trị không có dữ liệu - được đặt thành "0"

import gdal, gdalconst, numpy
maskfile = gdal.Open('C:\Users\max\Desktop\Baltic2.tif', gdalconst.GA_Update)
maskraster = maskfile.ReadAsArray()
maskraster = numpy.where((maskraster >= 0), maskraster, 0 ) 
maskband = maskfile.GetRasterBand(1)
maskband.WriteArray( maskraster )
maskband.FlushCache()

Tôi đã thử sử dụng cả plugin rastercalc và Raster Calculator và không hoạt động - Lỗi cú pháp. Công cụ Warp đã làm gì và chỉ định các giá trị nốt và nguồn đích (xem câu trả lời của phần dưới cho câu hỏi tương tự Làm thế nào để loại trừ các giá trị bị thiếu khỏi lớp raster? )

Điều làm việc cho tôi là xuất tệp raster dưới dạng tệp .asc, sau đó mở trong trình soạn thảo văn bản và sử dụng find và thay thế để thay thế giá trị NoData cũ bằng giá trị mà tôi muốn.

Trong trường hợp của tôi, đó là hình ảnh vệ tinh tôi cần để chuyển đổi thành 8 bit và nó không xuất ra giá trị dữ liệu nào khi sử dụng -scale trong gdal_translate.

Tôi đã có thể giải quyết vấn đề bằng cách xuất raster (lớp nhấp chuột phải -> xuất -> lưu lớp dưới dạng ...) vào một tệp mới dưới dạng dữ liệu được hiển thị (hộp kiểm ở đầu cửa sổ hộp thoại).

Đối với kịch bản python, bạn có thể sử dụng dòng sau trước WriteArray.

maskband.SetNoDataValue(0)