これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、MAC OS オンライン エミュレーターなど、複数の無料オンライン ワークステーションのいずれかを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド i.pansharpengrass です。
プログラム：
NAME
i.パンシャープン - 高解像度でマルチスペクトルを鮮明にする画像融合アルゴリズム
パンクロマティック チャンネル
KEYWORDS
画像、融合、シャープ化、Brovey、IHS、HIS、PCA
SYNOPSIS
i.パンシャープン
i.パンシャープン - 助けて
i.パンシャープン [-sl] 赤=名 （緑字）=名 青=名 パン=名 出力=ベース名 方法=string
[-上書きする] [-助けます] [-詳細] [-静かな] [-ui]
フラグ：
-s
並列処理ではなく逐次処理
-l
LANDSAT の青チャンネルのバランスを調整
-上書き
出力ファイルが既存のファイルを上書きできるようにする
- 助けて
使用状況の概要を印刷する
-詳細
冗長モジュール出力
- 静かな
静かなモジュール出力
--ui
GUIダイアログを強制的に起動する
パラメーター：
赤=名 [必要]
使用するラスターマップの名前
（緑字）=名 [必要]
使用するラスターマップの名前
青=名 [必要]
使用するラスターマップの名前
パン=名 [必要]
高解像度パンクロマティック チャネルに使用されるラスター マップの名前
出力=ベース名 [必要]
出力ベース名ラスター マップの名前
方法=string [必要]
パンシャープニングの方法
オプション： ブロービー、 ええ、 PCA
デフォルト： 私は
DESCRIPTION
i.パンシャープン マルチスペクトル画像の高解像度パンクロマティック バンドを使用して、
3 つの低解像度バンドをシャープにします。 次に、解像度の低い 3 つのバンドを結合して、
を使用して可能な解像度よりも高い (より詳細な) 解像度の RGB カラー画像
オリジナル3バンド。 たとえば、Landsat ETM には低解像度のスペクトル バンド 1 (青)、2 があります。
(緑)、3 (赤)、4 (近赤外線)、5 (中赤外線)、および 7 (中赤外線) を 30m の解像度で、および高い
解像度 8m の解像度でパンクロマティック バンド 15。 パンシャープニングにより、バンド 3-2-1 (または
30-4-3 や 2-5-4 などの 2m 解像度バンドの他の組み合わせ) を 15m に組み合わせる
解像度のカラー画像。
i.pansharpen では、IHS、Brovey、
そしてPCA。
IHS パン シャープ、赤、緑、およびとして選択された元の 3 つの低解像度バンド
RGB 合成画像を作成するための青チャンネルは、IHS (強度、
色相、彩度) 色空間。 次に、パンクロマティック バンドが代わりに使用されます。
元の色相 (H) および彩度 (S) チャンネルと組み合わせた強度チャンネル (I)、および
パンクロマティック バンドのより高い解像度で RGB 色空間に変換されます。 の
このアルゴリズムは、RGB -> IHS -> [pan]HS -> RGB として表すことができます。
とともに ブロベイ パン シャープ、3 つの低解像度バンドとパンクロマティック バンドのそれぞれ
次のアルゴリズムを使用して結合され、より高い位置にある 3 つの新しいバンドが計算されます。
解像度 (バンド 1 の例):
band1
new band1 = ----------------------- * パンバンド
バンド1 + バンド2 + バンド3
In PCA パン シャープ、主成分分析は元の 3 下位で実行されます
3 つの主成分画像 (PC1、PC2、および PC3) を作成するための解像度バンドとその
関連する固有ベクトル (EV)、次のように:
バンド1 バンド2 バンド3
PC1: EV1-1 EV1-2 EV1-3
PC2: EV2-1 EV2-2 EV2-3
PC3: EV3-1 EV3-2 EV3-3
と
PC1 = EV1-1 * バンド 1 + EV1-2 * バンド 2 + EV1-3 * バンド 3 - 平均 (バンド 1,2,3)
次に、逆 PCA が実行され、PC1 がパンクロマティック バンドに置き換えられます。 これをする、
固有ベクトル行列が反転され (この場合は転置)、PC 画像は次のようになります。
PC1 の代わりにパンクロマティック バンドを使用して固有ベクトルを乗算し、
各バンドは、次のアルゴリズムを使用して、変換された各イメージ バンドに追加されます (例
バンド 1 の場合):
band1' = pan * EV1-1 + PC2 * EV2-1 + PC3 * EV3-1 + mean(band1)
チャンネルの割り当ては、衛星によって異なります。 衛星画像の例
高解像度のパンクロマティック バンドと、低解像度のスペクトル バンドを含む
Landsat 7 ETM、QuickBird、および SPOT。
注意事項
このモジュールは現在、8 ビット イメージでのみ機能します。
このコマンドは、計算領域を一時的に高解像度に変更します。
シャープニング計算中のパンクロマティック バンド、その後前の領域を復元
設定。 現在のリージョン座標 (および null 値) が考慮されます。 高
解像度のパンクロマティック画像は、置換前のバンドに一致するヒストグラムです。
置換 (つまり、IHS シャープニングの強度チャンネル、選択された低解像度帯域
Brovey シャープニングを使用した各カラー チャネル、および PCA シャープニング用の PC1 画像)。
デフォルトでは、コマンドは最大 3 つのコアを使用して並列処理を採用しようとします。
同時に。 -s フラグは並列処理を無効にしますが、最適化された
ディスク I/O を削減するための r.mapcalc 式。
XNUMX つのパンシャープン出力チャンネルは、 d.rgb or r.コンポジット. 色
オプションで最適化することができます i.colors.enhance. 結果のカラー画像は
すべてのケースでより高い解像度になるため、3 つのパン シャープニング アルゴリズムは条件が異なります。
スペクトル応答の。
例
米国コロラド州ボールダーの Landsat 画像のパン シャープニング:
# 30m での R、G、B コンポジット
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
# IHS アルゴリズムによる i.pansharpen
i.pansharpen red=p034r032_7dt20010924_z13_30 green=p034r032_7dt20010924_z13_20
blue=p034r032_7dt20010924_z13_10 pan=p034r032_7dp20010924_z13_80
出力=ihs321 メソッド=ihs
# 15mで表示
g.region raster=ihs321_blue -p
d.rgb b=ihs321_blue g=ihs321_green r=ihs321_red
結果について
R, G, B コンポジット of ランドサット at 30m R, G, B コンポジット of ブロベイ シャープ 画像 at 15m
R, G, B コンポジット of IHS シャープ 画像 at 15m R, G, B コンポジット of PCA シャープ 画像 at 15m」
例: LANDSAT ETM+ (Landsat 7)、ノースカロライナ州のサンプル データセット:
# 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
# IHS アルゴリズムによる i.pansharpen
i.pansharpen red=lsat7_2002_30@PERMANENT \
green=lsat7_2002_20 blue=lsat7_2002_10 \
パン=lsat7_2002_80 メソッド=ihs \
出力=lsat7_2002_ihs
# 14.25mで表示
g.region raster=lsat7_2002_ihs_red -p
d。消去
d.rgb r=lsat7_2002_ihs_red g=lsat7_2002_ihs_green b=lsat7_2002_ihs_blue
# 前後の比較 (「Advanced」での RGB サポート):
g.gui.マップスワイプ
# オプションでカラーバランス:
i.colors.enhance r=lsat7_2002_ihs_red g=lsat7_2002_ihs_green b=lsat7_2002_ihs_blue
onworks.net サービスを使用してオンラインで i.pansharpengrass を使用する
