i.pansharpengrass - En ligne dans le Cloud

PROGRAMME:

Nom

je.pansharpen - Algorithmes de fusion d'images pour affiner le multispectral avec la haute résolution
canaux panchromatiques

MOTS-CLÉS

imagerie, fusion, netteté, Brovey, IHS, HIS, PCA

SYNOPSIS

je.pansharpen
je.pansharpen --Aidez-moi
je.pansharpen [-sl] rouge=prénom et une transition qui soit juste.=prénom Bleu=prénom casserole=prénom sortie=nom de base méthode=un magnifique
[--écraser] [--aider] [--verbeux] [--calme] [--ui]

Drapeaux:
-s
Traitement série plutôt que traitement parallèle

-l
Rééquilibrer le canal bleu pour LANDSAT

--écraser
Autoriser les fichiers de sortie à écraser les fichiers existants

--Aidez-moi
Imprimer le récapitulatif d'utilisation

--verbeux
Sortie du module verbeux

--silencieux
Sortie module silencieuse

--interface utilisateur
Forcer le lancement de la boîte de dialogue GUI

Paramètres:
rouge=prénom [obligatoire]
Nom de la carte raster à utiliser pour

et une transition qui soit juste.=prénom [obligatoire]
Nom de la carte raster à utiliser pour

Bleu=prénom [obligatoire]
Nom de la carte raster à utiliser pour

casserole=prénom [obligatoire]
Nom de la carte raster à utiliser pour le canal panchromatique haute résolution

sortie=nom de base [obligatoire]
Nom de la ou des cartes raster de nom de base en sortie

méthode=un magnifique [obligatoire]
Méthode d'aiguisage
Options: Brovey, c'est moi, pca
Valeur par défaut: Ihs

DESCRIPTION

je.pansharpen utilise une bande panchromatique haute résolution d'une image multispectrale à
affiner 3 bandes de résolution inférieure. Les 3 bandes de résolution inférieure peuvent ensuite être combinées en
une image couleur RVB à une résolution plus élevée (plus détaillée) que ce qui est possible en utilisant le
3 bandes originales. Par exemple, Landsat ETM a des bandes spectrales basse résolution 1 (bleu), 2
(vert), 3 (rouge), 4 (proche IR), 5 (mi-IR) et 7 (mi-IR) à une résolution de 30 m, et une haute
résolution bande panchromatique 8 à 15 m de résolution. L'affûtage panoramique permet aux bandes 3-2-1 (ou
d'autres combinaisons de bandes de résolution de 30 m comme 4-3-2 ou 5-4-2) à combiner en un 15 m
image couleur de résolution.
i.pansharpen propose un choix de trois algorithmes d'"affûtage panoramique" différents : IHS, Brovey,
et PCA.
Pour IHS casserole aiguisage, les 3 bandes de résolution inférieure d'origine, sélectionnées comme rouge, verte et
canaux bleus pour créer une image composite RVB, sont transformés en IHS (intensité,
teinte et saturation) espace colorimétrique. La bande panchromatique est alors substituée à la
canal d'intensité (I), combiné avec les canaux de teinte (H) et de saturation (S) d'origine, et
retransformé en espace colorimétrique RVB à la résolution la plus élevée de la bande panchromatique. Les
L'algorithme pour cela peut être représenté comme : RVB -> IHS -> [pan]HS -> RVB.
Avec son Brovey casserole aiguisage, chacune des 3 bandes de résolution inférieure et la bande panchromatique
sont combinés en utilisant l'algorithme suivant pour calculer 3 nouvelles bandes au plus haut
résolution (exemple pour la bande 1) :
band1
nouvelle bande1 = ----------------------- * panband
bande1 + bande2 + bande3
In PCA casserole aiguisage, une analyse en composantes principales est effectuée sur les 3 premiers
bandes de résolution pour créer 3 images en composantes principales (PC1, PC2 et PC3) et leurs
vecteurs propres (EV) associés, tels que :
bande1 bande2 bande3
PC1: EV1-1 EV1-2 EV1-3
PC2: EV2-1 EV2-2 EV2-3
PC3: EV3-1 EV3-2 EV3-3
ainsi que
PC1 = EV1-1 * bande1 + EV1-2 * bande2 + EV1-3 * bande3 - moyenne (bandes 1,2,3)
Une ACP inverse est ensuite effectuée, substituant la bande panchromatique à PC1. Pour faire ça,
la matrice des vecteurs propres est inversée (dans ce cas transposée), les images PC sont
multiplié par les vecteurs propres avec la bande panchromatique substituée à PC1, et la moyenne de
chaque bande est ajoutée à chaque bande d'image transformée en utilisant l'algorithme suivant (exemple
pour la bande 1) :
bande1' = panoramique * EV1-1 + PC2 * EV2-1 + PC3 * EV3-1 + moyenne (bande1)
L'attribution des canaux dépend du satellite. Exemples d'images satellites
avec des bandes panchromatiques à haute résolution et des bandes spectrales à basse résolution comprennent
Landsat 7 ETM, QuickBird et SPOT.

NOTES

Le module ne fonctionne actuellement que pour les images 8 bits.
La commande change temporairement la région de calcul en haute résolution du
bande panchromatique lors des calculs de netteté, puis restaure la région précédente
Les paramètres. Les coordonnées de la région actuelle (et les valeurs nulles) sont respectées. Le haut
l'image panchromatique de résolution correspond à l'histogramme de la bande qu'elle remplace avant
substitution (c'est-à-dire le canal d'intensité pour l'affûtage IHS, la bande basse résolution sélectionnée
pour chaque canal de couleur avec la netteté Brovey et l'image PC1 pour la netteté PCA).
Par défaut, la commande tentera d'employer un traitement parallèle, en utilisant jusqu'à 3 cœurs
simultanément. L'indicateur -s désactivera le traitement parallèle, mais utilise un
Expression r.mapcalc pour réduire les E/S du disque.
Les trois canaux de sortie pan-sharpen peuvent être combinés avec d.rgb or r.composite. Couleurs
peut être optimisé en option avec i.colors.enhance. Alors que l'image couleur résultante
être à la résolution la plus élevée dans tous les cas, les 3 algorithmes de netteté panoramique diffèrent en termes
de la réponse spectrale.

EXEMPLES

Accentuation panoramique d'une image Landsat de Boulder, Colorado, États-Unis :
# R, V, B composite à 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 avec l'algorithme IHS
i.pansharpen red=p034r032_7dt20010924_z13_30 green=p034r032_7dt20010924_z13_20
blue=p034r032_7dt20010924_z13_10 pan=p034r032_7dp20010924_z13_80
sortie=ihs321 méthode=ihs
# affichage à 15m
g.region raster=ihs321_blue -p
d.rgb b=ihs321_blue g=ihs321_green r=ihs321_red
Résultats:

R, G, B pierre composite of Landsat at 30m R, G, B pierre composite of Brovey aiguisé image at 15m

R, G, B pierre composite of IHS aiguisé image at 15m R, G, B pierre composite of PCA aiguisé image at 15m"

Exemple : LANDSAT ETM+ (Landsat 7), exemple de jeu de données en Caroline du Nord :
# original à 28m
g.région raster=lsat7_2002_10 -p
d.lun wx0
d.rgb b=lsat7_2002_10 g=lsat7_2002_20 r=lsat7_2002_30
# i.pansharpen avec l'algorithme IHS
i.pansharpen red=lsat7_2002_30@PERMANENT \
green=lsat7_2002_20 blue=lsat7_2002_10 \
pan=méthode lsat7_2002_80=ihs \
sortie=lsat7_2002_ihs
# affichage à 14.25m
g.region raster=lsat7_2002_ihs_red -p
effacer
d.rgb r=lsat7_2002_ihs_red g=lsat7_2002_ihs_green b=lsat7_2002_ihs_blue
# comparer avant/après (support RVB dans "Avancé") :
g.gui.mapswipe
# éventuellement équilibrage des couleurs :
i.colors.enhance r=lsat7_2002_ihs_red g=lsat7_2002_ihs_green b=lsat7_2002_ihs_blue

Utilisez i.pansharpengrass en ligne en utilisant les services onworks.net