Setiap kali saya mencoba untuk mempertajam komposit beberapa gambar Landsat dalam GRASS menggunakan i.pansharpen
, i.fusion.brovey
atau metode penajaman IHS, output akan memiliki beberapa atau semua karakteristik berikut:
- warna komposit berada dalam rona yang berbeda dibandingkan dengan komposit yang tidak diasah
- tingkat kecerahan kacau
- seluruh komposit menjadi serba putih / hitam-semua (saat menggunakan gambar pra-diproses untuk reflektansi atmosfer-atas atau koreksi reflektansi permukaan di
i.landsat.toar
)
Saya juga sudah mencoba semua yang berikut; tetapi warna / kecerahannya tetap sama atau berubah menjadi lebih buruk:
- Diterapkan
i.landsat.rgb
, sebelum dan sesudah proses penajaman panci - Dimainkan dengan tanda
-f
atau-p
masuki.landsat.rgb
- Mencoba
r.colors
mengedit tabel warna ke grey / grey255 / grey.eq - Sudah mencoba
i.pansharpen
menggunakan semua metode Brovey / IHS / PCA - Dimainkan dengan
-l
benderai.pansharpen
untuk menyeimbangkan saluran biru
The pengguna GRASS GIS tidak menjelaskan tentang cara melakukan pan-sharpening dan warna-balancing, tapi aku tidak bisa mencari cara untuk menggabungkan kedua proses dalam alur kerja bersamaan. Saya menduga bahwa ini adalah karena pemahaman saya yang buruk tentang tabel warna, histogram warna, dll dalam GRASS ..
Jadi, bisakah seseorang menjelaskan kepada saya - bagaimana Anda mengatasi masalah penyeimbangan warna ketika berhadapan dengan gambar Landsat setelah pemrosesan gambar dalam GRASS? Bisakah Anda berbagi dengan saya alur kerja / metode favorit Anda?
Terima kasih banyak atas umpan baliknya!
i.pansharpen
menangani semua jenis format. Silakan lihat "tiket" terkait: Tiket # 2048: i.pansharpen terbatas pada citra 8-bit . Di sisi lain, saya mungkin tidak mengerti hal-hal dan, dengan demikian, tidak dapat menggunakan i.pansharpen dengan benar ...?Saya telah mencari tinggi dan rendah, dan saya pikir saya telah menemukan akar masalah saya. Saya yakin saya punya solusi untuk mereka sekarang - tapi agak berantakan. Saya yakin ada cara yang lebih baik untuk menyelesaikannya. Berbagilah jika Anda tahu cara yang lebih mudah!
AKAR MASALAH:
i.landsat.toar
dalam floating point . Saya menyadari bahwa ketika saya menggunakan raster floating point dalam metode pan-sharpening apa pun, warnanya akan berantakan. Algoritme itu entah bagaimana memilih raster dalam bentuk integer asli.i.pansharpen
dani.fusion.brovey
modul akan mengacaukan warnanya . Saya belum cukup memahami algoritma yang mereka gunakan dalam modul tersebut - tetapi entah bagaimana tabel warna akan terpengaruh, dan merusak gambar pan-sharpening yang dihasilkan.LARUTAN:
i.landsat.toar
dari float ke int, menggunakanr.recode
.i.rgb.his
dani.his.rgb
). Saya terus menggunakani.pansharpen
dani.fusion.brovey
.WORKFLOW:
Gunakan
r.info
dengan-r
bendera untuk mendapatkan nilai DN min dan maks dari setiap pita raster yang telah diprosesi.landsat.toar
. Sebagai contoh:Seperti yang dapat kita lihat, nilainya antara 0-1, yang sangat berbeda dari yang asli (yaitu antara 0-255). Itu menjelaskan mengapa output dari pan-sharpening ternyata kosong, karena rentang nilai yang digunakan sangat rendah (di bawah 1).
Konversikan band raster itu menggunakan
r.recode
. Gunakan nilai minimum dan maksimum yang diperoleh dari langkah 1 untuk mengkonversi ke rentang baru 0-255. Cuplikan kode contoh:Kami dapat memeriksa nilai yang dikonversi baru dengan
r.info
:Nilai ada di 0-255: sekarang dapat digunakan untuk proses mengasah pan.
Terapkan tabel warna skala abu-abu dengan pita yang dikonversi
r.colors
.Sejauh ini, saya mendapatkan hasil terbaik menggunakan tabel warna abu-abu - komposit yang ditajamkan pan cocok dengan komposit asli. Alternatif lain adalah untuk menyamakan tabel warna abu-abu dengan
color=grey.eq
atau menggunakan-e
bendera dengancolor=grey
. Atau kita bisa menggunakani.landsat.rgb
modul, bukanr.colors
..Ulangi langkah 1-3 dengan pita raster lain yang ingin kami gunakan sebagai komposit, termasuk pan raster (band 8). Penggunaan skrip akan sangat dihargai di sini.
Kemudian gunakan raster yang diproses sebagai input dalam metode pan-sharpening IHS . Misalnya, ketika membuat komposit dari band 7,4,2:
Ini akan menghasilkan 3 lapisan: lapisan rona
HUE
, lapisan intensitasINT
, dan juga lapisan saturasiSAT
. Kami kemudian akan mengganti layer intensitasINT
dengan band pan rasterNEWBAND8
dii.his.rgb
:Menghasilkan saluran merah pada COMP742_red , COMP742_green , COMP742_blue kemudian dapat digabungkan menggunakan
d.rgb
ataur.composite
..SAMPEL SEBELUM & SETELAH:
Sebelum mengasah pan:
Setelah pan-sharpening:
Mungkin sulit untuk membedakannya , ketika melihat dari gambar sekecil itu. Namun yang penting adalah warna pan-sharpened image cocok dengan komposit dari aslinya. Misi selesai!
CATATAN LAINNYA:
r.recode
menggunakan pita termal (pita 6) .i.landsat.toar
mengeluarkan pita termal ini dalam nilai suhu Kelvin (tidak ada hubungannya dengan nilai DN). Jagar.recode
rutinitas pada band multi-spektral dan pan normal (band 1-5,7,8).i.landsat.toar
, tetapi komposit yang dihasilkan terlihat sangat salah , biasanya karena ketidakcocokan tabel warna sebelum dan sesudah proses pan-sharpening. Saya menerapkanr.colors RASTER color=grey
pada band raster asli sebelum pan-sharpening, dan ke saluran yang dihasilkan setelah pan-sharpening untuk memastikan pencocokan warna.color=grey255
, output darii.landsat.rgb
dalamcolor=grey.eq
. Tidak heran mereka berdua terlihat berbeda!Semoga ini bermanfaat bagi seseorang: butuh waktu berhari-hari untuk menemukan apa yang salah ..
sumber
r.mapcalc
bersama dengan fungsinya yang terintegrasiint()
. Saya pikirr.recode
ini tidak dimaksudkan untuk digunakan dalam konteks aplikasi Anda.i.pansharpen
adalah bahwa ia mengharapkan peta raster 8-bit sebagai input, yaitu nilai input harus berkisar dari 0 hingga 255.i.pansharpen
hanya di grass7_trunk. Itu adalah versi pengembangan ... :-)