Bagaimana cara pansharpen LANDSAT 8 di GRASS?

9

Setelah menggunakan QGIS untuk beberapa waktu untuk pemetaan, saya 'dipaksa' untuk menggali GRASS karena pansharpening tidak secara langsung dapat dieksekusi di QGIS. Setelah kurang lebih memahami strukturnya, saya melakukan beberapa percobaan dengan i.fusion.brovey. Saya dapat menghasilkan raster dengan resolusi 15m, tetapi ketika memuat ke QGIS saya tidak dapat mencapai / memulihkan kontras warna yang bagus seperti yang ada dalam raster 30m asli.

Dalam posting sebelumnya seseorang mengkonfirmasi itu dan mengusulkan untuk menggunakan fungsi lain untuk landsat: 'i.his.rgb' (setelah menggunakan i.landsat.toar, yang tidak saya gunakan). Dia mulai mengonversi min / maks ke 1-255 sebelum menggunakan i.his.rgb ... Raster asli saya memiliki kisaran 0-65535. Itu tidak berhasil bagi saya ...

Saya berhenti setelah beberapa hari berjuang untuk hanya melakukan pansharp, saya merasa macet ... Mengapa hanya menerapkan i.pansharp tidak bekerja untuk Landsat seperti yang disarankan oleh orang lain? Bisakah seseorang menyarankan alur kerja yang lebih mudah untuk mengeksekusi pansharpening dari gambar Landsat 8? Menantikan beberapa saran, terima kasih sebelumnya ...

Arnoud Cuppens
sumber

Jawaban:

7

Lihatlah i.landsat.rgb - Melakukan penyeimbangan warna secara otomatis untuk gambar LANDSAT, mungkin sebelum menjalankan pansharpening. Anda juga dapat mempertimbangkan untuk mengubah angka digital dari masing-masing saluran menjadi cahaya atmosfer atau reflektansi dengan i.landsat.toar . Lihat juga http://grasswiki.osgeo.org/wiki/LANDSAT

BTW: memiliki rentang 0-65535 untuk data LANDSAT asli adalah aneh karena dikirimkan sebagai saluran 8 bit (0-255). Jadi periksa data input Anda terlebih dahulu.

Pembaruan: 10/2013

Untuk prosedur yang disederhanakan dengan GRASS GIS 7 saat ini, lihat halaman masukkan deskripsi gambar di siniMemproses data Landsat 8 dalam GRASS GIS 7: komposit RGB dan penajaman panci

markN
sumber
Saya awalnya menemukan artikel yang Anda tautkan dari neteler.org, tetapi petunjuknya tidak bekerja sama sekali (karena format floating point raster yang dihasilkan i.landsat.toar, jadi saya akhirnya mengikuti saran dari sini: gis.stackexchange.com / a / 39239/9857
Juan
3
Perhatikan bahwa saran itu diberikan bertahun-tahun yang lalu dan tidak mencerminkan kode GIS GRASS terbaru. BTW: modul pansharpening yang benar-benar hebat adalah tambahan baru: i.fusion.hpf - Sekering Panchromatic resolusi tinggi dan data Multi-Spektral resolusi tinggi berdasarkan pada teknik Penambahan Filter Pass Tinggi, lihat grass.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html
markusN
4

Terima kasih, Michal dan Markus. Akhirnya saya bisa melakukan pansharpening dengan fungsi yang ditunjukkan brov. Setelah dilakukan di GRASS, saya mengekspor raster RGB ke gdal GTiff dengan tipe data Uint16. Markus, saya telah mengunduh dan menggunakan beberapa gambar Landsat 8 (masing-masing dari 11 band), memperoleh hasil visual yang bagus, baik dalam warna alami / salah. Sepengetahuan saya, tiff landsat asli datang sebagai Uint16 dengan nilai mulai dari 0 hingga 65535 (atau apakah saya kehilangan sesuatu?) ... Terima kasih atas tautan GRASS LANDSAT ...

Arnoud Cuppens
sumber
3
Ya, gambar Landsat 8 disimpan sebagai tipe Uint16.
webrian
3

FIHS atau Brovery Transform pan-sharpening cukup sederhana jika Anda mengikuti formula dari artikel seperti ini: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Dengan asumsi Anda menggunakan GRASS GIS dan pansharpening Landsat 8.

Atur wilayah Anda menggunakan pan 8 tif - penting untuk memiliki resolusi 15m yang akan ditetapkan jika Anda melakukan ini.

Untuk Brovery tanpa bobot, masukkan ini ke dalam Kalkulator Peta masing-masing untuk warna biru, hijau dan merah:

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

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

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

Jalankan r.colors untuk masing-masing output di atas dan atur histogram equalisasi dan tabel warna menjadi abu-abu.

Jalankan r.composite, pilih band yang sesuai (merah = 4, hijau = 3, biru = 2)

Lihat gambar.

Jelas proses ini dapat dilakukan dengan lebih banyak pertimbangan untuk koreksi atmosfer, peregangan warna dan bobot pita seperti yang dibahas dalam artikel di atas.

Ben
sumber
2

Anda juga dapat mencoba Monteverdi untuk melakukan pansharpening. Alat ini dapat menggunakan GeoTIFF biasa sebagai input (namun, pita multispektral yang terpisah harus digabung menjadi satu raster multiband), dan berkinerja cukup baik.

Michal Mackiewicz
sumber
2

Setelah bergumul dengan masalah ini, saya menemukan solusi sederhana yang memungkinkan Anda melakukan pan-sharpen di QGIS (2.18.21 adalah versi saya) menggunakan SCP. Pembaruan untuk SCP yang menggabungkan metode ini dijelaskan di sini https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Ini mudah dilakukan dan berjalan seperti pesona! Saya memeriksa hasil akhir dan ukuran piksel pada resolusi 15 m - cocok dengan band 8 panchromatic. Jika ada yang penasaran - ini menggunakan metode Brovey Transform seperti yang dijelaskan dalam manual (lihat https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Langkah 1 untuk alat pansharpen SCP Langkah 3-4 untuk pansharpen SCP

Mark Thompson
sumber
1

Selamat datang di situs Janos. Idealnya jawaban seperti milik Anda harus memiliki deskripsi mengapa itu adalah metode yang baik. ..

Untuk pansharpening gambar Landsat 8 metode yang lebih mudah, http://www.geosage.com/highview/download.html . Saran saya, coba saja. Hanya satu klik ....Transformator spektral

János Breznai
sumber
-2

Untuk pansharpenig saya menggunakan transformator spektral GUI untuk Landsat 8 gambar masukkan deskripsi gambar di sini http://www.geosage.com/highview/download.html Sangat bagus. BR Janos

János Breznai
sumber
2
Selamat datang di situs Janos. Idealnya jawaban seperti milik Anda harus memiliki deskripsi mengapa itu adalah metode yang baik. Jawaban hanya tautan cenderung menjadi tidak berfungsi jika tautan hilang.
Aaron