i.pansharpenggrass - Online di Cloud

PROGRAM:

NAMA

i.pansharpen - Algoritma fusi gambar untuk mempertajam multispektral dengan resolusi tinggi
saluran pankromatik

KEYWORDS

citra, fusi, pertajam, Brovey, IHS, HIS, PCA

RINGKASAN

i.pansharpen
i.pansharpen --membantu
i.pansharpen [-sl] merah=nama hijau=nama biru=nama panci=nama keluaran=nama dasar metode=string
[--menimpa] [--membantu] [--bertele-tele] [--tenang] [--ui]

Bendera:
-s
Pemrosesan serial daripada pemrosesan paralel

-l
Seimbangkan kembali saluran biru untuk LANDSAT

--timpa
Izinkan file keluaran untuk menimpa file yang ada

--membantu
Cetak ringkasan penggunaan

--bertele-tele
Keluaran modul verbose

--diam
Keluaran modul yang tenang

--ui
Paksa meluncurkan dialog GUI

Parameter:
merah=nama [yg dibutuhkan]
Nama peta raster yang akan digunakan

hijau=nama [yg dibutuhkan]
Nama peta raster yang akan digunakan

biru=nama [yg dibutuhkan]
Nama peta raster yang akan digunakan

panci=nama [yg dibutuhkan]
Nama peta raster yang akan digunakan untuk saluran pankromatik resolusi tinggi

keluaran=nama dasar [yg dibutuhkan]
Nama untuk peta raster nama dasar keluaran

metode=string [yg dibutuhkan]
Metode untuk mengasah panci
Pilihan: broy, ih, pc
Default: ihs

DESKRIPSI

i.pansharpen menggunakan pita pankromatik resolusi tinggi dari gambar multispektral untuk
mempertajam 3 pita resolusi lebih rendah. 3 pita resolusi lebih rendah kemudian dapat digabungkan menjadi
gambar warna RGB pada resolusi yang lebih tinggi (lebih detail) daripada yang dimungkinkan menggunakan
asli 3 band. Misalnya, Landsat ETM memiliki pita spektrum resolusi rendah 1 (biru), 2
(hijau), 3 (merah), 4 (IR dekat), 5 (IR tengah), dan 7 (IR tengah) pada resolusi 30m, dan tinggi
resolusi pankromatik band 8 pada resolusi 15m. Penajaman panci memungkinkan pita 3-2-1 (atau
kombinasi lain dari pita resolusi 30m seperti 4-3-2 atau 5-4-2) untuk digabungkan menjadi 15m
resolusi gambar warna.
i.pansharpen menawarkan tiga pilihan algoritma 'pan sharpening' yang berbeda: IHS, Brovey,
dan PCA.
Untuk IHS panci mengasah, 3 pita resolusi rendah asli, dipilih sebagai merah, hijau dan
saluran biru untuk membuat gambar komposit RGB, diubah menjadi IHS (intensitas,
hue, dan saturasi) ruang warna. Pita pankromatik kemudian diganti dengan
saluran intensitas (I), dikombinasikan dengan saluran warna asli (H) dan saturasi (S), dan
ditransformasikan kembali ke ruang warna RGB pada resolusi pita pankromatik yang lebih tinggi. NS
algoritma untuk ini dapat direpresentasikan sebagai: RGB -> IHS -> [pan]HS -> RGB.
Dengan Brovey panci mengasah, masing-masing dari 3 pita resolusi lebih rendah dan pita pankromatik
digabungkan menggunakan algoritma berikut untuk menghitung 3 pita baru di yang lebih tinggi
resolusi (contoh untuk band 1):
band1
band1 baru = ----------------------- * panband
band1 + band2 + band3
In PCA panci mengasah, analisis komponen utama dilakukan pada 3 yang lebih rendah asli
pita resolusi untuk membuat 3 gambar komponen utama (PC1, PC2, dan PC3) dan
terkait eigenvectors (EV), sehingga:
band1 band2 band3
PC1: EV1-1 EV1-2 EV1-3
PC2: EV2-1 EV2-2 EV2-3
PC3: EV3-1 EV3-2 EV3-3
dan
PC1 = EV1-1 * band1 + EV1-2 * band2 + EV1-3 * band3 - mean(band 1,2,3)
PCA terbalik kemudian dilakukan, menggantikan pita pankromatik untuk PC1. Untuk melakukan ini,
matriks vektor eigen dibalik (dalam hal ini ditransposisikan), gambar PC adalah
dikalikan dengan vektor eigen dengan pita pankromatik yang diganti dengan PC1, dan rata-rata
setiap band ditambahkan ke setiap band gambar yang diubah menggunakan algoritma berikut (contoh
untuk kelompok 1):
band1' = pan * EV1-1 + PC2 * EV2-1 + PC3 * EV3-1 + mean(band1)
Penugasan saluran tergantung pada satelit. Contoh citra satelit
dengan pita pankromatik resolusi tinggi, dan pita spektral resolusi lebih rendah termasuk
Landsat 7 ETM, QuickBird, dan SPOT.

CATATAN

Modul saat ini hanya berfungsi untuk gambar 8-bit.
Perintah untuk sementara mengubah wilayah komputasi ke resolusi tinggi dari
pita pankromatik selama perhitungan penajaman, lalu mengembalikan wilayah sebelumnya
pengaturan. Koordinat wilayah saat ini (dan nilai nol) dihormati. yang tinggi
resolusi citra pankromatik adalah histogram yang dicocokkan dengan pita yang digantikannya sebelum
substitusi (yaitu, saluran intensitas untuk penajaman IHS, pita beresolusi rendah dipilih
untuk setiap saluran warna dengan penajaman Brovey, dan gambar PC1 untuk penajaman PCA).
Secara default, perintah akan mencoba menggunakan pemrosesan paralel, menggunakan hingga 3 core
serentak. Bendera -s akan menonaktifkan pemrosesan paralel, tetapi menggunakan yang dioptimalkan
ekspresi r.mapcalc untuk mengurangi I/O disk.
Tiga saluran keluaran pan-sharpened dapat digabungkan dengan d.rgb or r.komposit. warna
dapat dioptimalkan secara opsional dengan i.warna.meningkatkan. Sedangkan gambar berwarna yang dihasilkan akan
berada pada resolusi yang lebih tinggi dalam semua kasus, algoritma penajaman 3 panci berbeda dalam hal
dari respon spektral.

CONTOH

Penajaman pan gambar Landsat dari Boulder, Colorado, AS:
# R, G, B komposit pada 30m
g.region raster=p034r032_7dt20010924_z13_10 -p
d.rgb b=p034r032_7dt20010924_z13_10 g=lp034r032_7dt20010924_z13_20
r=p034r032_7dt20010924_z13_30
#i.pansharpen dengan algoritma IHS
i.pansharpen red=p034r032_7dt20010924_z13_30 green=p034r032_7dt20010924_z13_20
blue=p034r032_7dt20010924_z13_10 pan=p034r032_7dp20010924_z13_80
keluaran=ihs321 metode=ihs
# tampilan pada 15m
g.wilayah raster=ihs321_blue -p
d.rgb b=ihs321_biru g=ihs321_hijau r=ihs321_merah
hasil:

R, G, B gabungan of landat at 30m R, G, B gabungan of Brovey dipertajam gambar at 15m

R, G, B gabungan of IHS dipertajam gambar at 15m R, G, B gabungan of PCA dipertajam gambar at 15m"

Contoh: LANDSAT ETM+ (Landsat 7), kumpulan data sampel Carolina Utara:
# asli pada 28m
g.region raster=lsat7_2002_10 -p
d.mon wx0
d.rgb b=lsat7_2002_10 g=lsat7_2002_20 r=lsat7_2002_30
#i.pansharpen dengan algoritma IHS
i.pansharpen red=lsat7_2002_30@PERMANENT \
green=lsat7_2002_20 blue=lsat7_2002_10 \
pan=lsat7_2002_80 metode=ihs \
keluaran=lsat7_2002_ihs
# tampilan pada 14.25m
g.region raster=lsat7_2002_ihs_red -p
d.hapus
d.rgb r=lsat7_2002_ihs_red g=lsat7_2002_ihs_green b=lsat7_2002_ihs_blue
# bandingkan sebelum/sesudah (dukungan RGB di "Lanjutan"):
g.gui.mapswipe
# penyeimbangan warna opsional:
i.warna.tingkatkan r=lsat7_2002_ihs_red g=lsat7_2002_ihs_hijau b=lsat7_2002_ihs_biru

Gunakan i.pansharpenggrass online menggunakan layanan onworks.net