Làm thế nào để pansharpen LandSAT 8 trong GRASS?


9

Đã sử dụng QGIS trong một thời gian khá lâu để lập bản đồ, tôi đã 'buộc' phải đào sâu vào GRASS vì việc pansharpening không được thực thi trực tiếp trong QGIS. Sau khi hiểu ít nhiều về cấu trúc, tôi đã thực hiện một số thử nghiệm với i.fusion.brovey. Tôi đã có thể tạo các trình quét với độ phân giải 15m, nhưng khi tải vào QGIS, tôi không thể đạt được / phục hồi độ tương phản màu đẹp như hiện tại trong các trình quét 30m ban đầu.

Trong một bài đăng trước đó, một người nào đó đã xác nhận điều đó và đề xuất sử dụng một chức năng khác cho Landsat: 'i.his.rgb' (sau khi sử dụng i.landsat.toar, mà tôi không sử dụng). Anh ta bắt đầu chuyển đổi tối thiểu / tối đa thành 1-255 trước khi sử dụng i.his.rgb ... Các trình quét gốc của tôi có phạm vi từ 0-65535. Nó không thành công với tôi ...

Tôi đã bỏ việc sau vài ngày vật lộn để chỉ thực hiện một pansharp, tôi cảm thấy bế tắc ... Tại sao chỉ nên áp dụng i.pansharp không hoạt động cho Landsat như đề xuất của người khác? Ai đó có thể đề xuất một quy trình làm việc đơn giản hơn để thực hiện pansharpening hình ảnh Landsat 8 không? Mong nhận được một lời khuyên, cảm ơn trước ...

Câu trả lời:


7

Hãy xem i.landsat.rgb - Thực hiện tự động cân bằng màu sắc cho hình ảnh ĐẤT, có thể trước khi chạy pansharpening. Bạn cũng có thể xem xét để chuyển đổi số kỹ thuật số của các kênh riêng lẻ sang độ rọi hoặc phản xạ trên cùng của khí quyển với i.landsat.toar . Xem thêm http://grasswiki.osgeo.org/wiki/LANDSAT

BTW: có phạm vi 0-65535 cho dữ liệu LandSAT ban đầu là số lẻ vì chúng được phân phối dưới dạng kênh 8 bit (0-255). Vì vậy, kiểm tra dữ liệu đầu vào của bạn đầu tiên.

Cập nhật: 10/2013

Để biết quy trình đơn giản hóa với GRASS GIS 7 hiện tại, hãy xem trang nhập mô tả hình ảnh ở đâyXử lý dữ liệu Landsat 8 trong GRASS GIS 7: Vật liệu tổng hợp RGB và mài sắc


Ban đầu tôi tình cờ thấy bài viết đó bạn đang liên kết từ neteler.org, nhưng các hướng dẫn hoàn toàn không hoạt động (vì định dạng dấu phẩy động được tạo bằng i.landsat.toar, vì vậy tôi đã kết thúc theo lời khuyên từ đây: gis.stackexchange.com / a / 39239/9857
Juan

3
Lưu ý rằng lời khuyên đó đã được đưa ra từ nhiều năm trước và không phản ánh mã GRASS GIS mới nhất. BTW: một mô-đun pansharpening thực sự tuyệt vời là addon mới: i.fusion.hpf - Hợp nhất dữ liệu Đa phổ độ phân giải cao và độ phân giải thấp dựa trên kỹ thuật Bổ sung bộ lọc High-Pass, xem grass.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html
markusN

4

Cảm ơn Michal và Markus. Cuối cùng tôi đã có thể thực hiện pansharpening với chức năng brov được chỉ định. Sau khi được thực hiện trong GRASS, tôi đã xuất các trình quét RGB sang gdal GTiff với kiểu dữ liệu Uint16. Markus, tôi đã tải xuống và sử dụng một số hình ảnh Landsat 8 (tiff của mỗi trong số 11 dải), thu được kết quả hình ảnh đẹp, cả hai màu tự nhiên / sai. Theo hiểu biết của tôi, các tiff Landsat ban đầu có tên là Uint16 với các giá trị nằm trong khoảng từ 0 đến 65535 (hoặc tôi có thiếu thứ gì không?) ... Cảm ơn vì liên kết GRASS LandSAT ...


3
Có, hình ảnh Landsat 8 được lưu trữ dưới dạng các loại Uint16.
webrian

3

Pan-sharpening FIHS hoặc Brovery Transform khá đơn giản nếu bạn làm theo các công thức từ một bài viết như thế này: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Giả sử bạn đang sử dụng GRASS GIS và pansharpening Landsat 8.

Đặt vùng của bạn bằng pan 8 tif - điều quan trọng là có độ phân giải 15m sẽ được đặt nếu bạn thực hiện việc này.

Đối với Brovery không có trọng lượng, hãy nhập phần này vào Máy tính bản đồ cho màu xanh lam, xanh lục và đỏ tương ứng:

band8 / (band2 + band3 + band4) / 3 * band2

band8 / (band2 + band3 + band4) / 3 * band3

band8 / (band2 + band3 + band4) / 3 * band4

Chạy r.colors cho mỗi đầu ra ở trên và đặt bảng cân bằng biểu đồ và bảng màu thành màu xám.

Chạy r.compozit, chọn các dải thích hợp (red = 4, green = 3, blue = 2)

Kiểm tra hình ảnh.

Rõ ràng quá trình này có thể được thực hiện với sự cân nhắc nhiều hơn cho hiệu chỉnh khí quyển, độ căng màu và trọng lượng dải như đã thảo luận trong bài viết trên.


2

Bạn cũng có thể thử Monteverdi để thực hiện pansharpening. Công cụ này có thể sử dụng các GeoTIFF thông thường làm đầu vào (tuy nhiên, các dải đa bán cầu riêng biệt phải được hợp nhất thành một raster multiband) và hoạt động khá tốt.


2

Sau khi vật lộn với vấn đề này, tôi đã tìm thấy một giải pháp đơn giản cho phép bạn làm sắc nét trong QGIS (2.18,21 là phiên bản của tôi) bằng SCP. Bản cập nhật cho SCP kết hợp phương pháp này được mô tả tại đây https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Điều này là dễ dàng để làm và nó chạy như một nét duyên dáng! Tôi đã kiểm tra đầu ra cuối cùng và các pixel đo ở độ phân giải 15 m - phù hợp với dải 8 panchromatic. Nếu bất kỳ ai cũng tò mò - nó sử dụng phương pháp Brovey Transform như được mô tả trong hướng dẫn sử dụng (xem https://media.readthedocs.org/pdf/semiautomatic classifyingmanual /latest / semiautomatic classificationmanual.pdf ).

Bước 1 cho công cụ pansharpen SCP Bước 3-4 cho SCP pansharpen


1

Chào mừng đến với trang web Janos. Các câu trả lời lý tưởng như của bạn nên có một mô tả về lý do tại sao chúng là phương pháp tốt. ..

Đối với pansharpening của Landsat 8 hình ảnh phương pháp dễ dàng hơn http://www.geosage.com/highview/doad.html . Lời khuyên của tôi, hãy thử nó. Chỉ một cú nhấp chuột ....Máy quang phổ


-2

Đối với pansharpenig tôi đã sử dụng biến áp phổ GUI cho hình ảnh Landsat 8 nhập mô tả hình ảnh ở đây http://www.geosage.com/highview/doad.html Rất tốt. BR Janos


2
Chào mừng đến với trang web Janos. Các câu trả lời lý tưởng như của bạn nên có một mô tả về lý do tại sao chúng là phương pháp tốt. Liên kết chỉ trả lời có xu hướng trở nên không còn nếu liên kết biến mất.
Aaron
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.