参考「地球観測におけるオープンソース地理空間ツールアプリケーション「マルチスペクトル土地被覆分類に関する17の事例-Chapter
今年の時間と、地球上の場所に応じて、衛星画像の
シーンは、いくつかの雲を含めることができます。このケーススタディ内で使用されるイメージが含まれています
いくつかの雲。雲が直接画質に影響を与え、深刻な妨げことができます
画像解析と分類。最小限の雲とのシーンを得るために、我々はコンバイン
シーンをダウンロードして、複合「無料のクラウド」を作成します。あなた自身を使用している場合
衛星画像は、別の場所のために取得し、それはあなたがスキップできるよりも雲無料です
次の処理ステップ。BQAデータセットは、ビットパックされたフォーマット2であり、含まれてい
水や雲に関する情報。これは、我々は情報をデコードする必要があることを意味します
例えば、クラウドマスクを得るためです。ビットパッキングに関する詳細情報
エンコードスキームは節で説明されています。我々はまた、作成するためのツールを実装して14.3、
BQAデータセットからの雲のマスク。また、あなたは、公式のツールを使用することができます
例としてUSGS.3によって提供されるL-LDOPEツールベルト、我々は雲を作成します
宗派からツールbqa2cloudとマスク。14.3:
bqa2 -i LC 82070232013160 LGN 00_ BQA。TIF -o LC 82070232013160クラウドLGN 00_ CLDを。TIF
bqa 2クラウド-i LC 82070232013192 LGN 00_ BQA。TIF -o LC 82070232013192 LGN 00_ CLD。TIF
我々は、次に、GeoTIFFの形式で単一のマルチバンドラスタデータセット内のすべてのバンドスタック
ユーティリティgdal_merge.pyを使用します。パンクロマチックバンドは少し余分な含まれている
スペクトル情報を、サーマルバンドが粗い空間で取得された
解像度。したがって、我々はバンド8、10と11が含まれていません。
gdal_ merge.py -o LC82070232013160LGN00_ merge.tif -of GTiff - PS →別個-CO INTERLEAVE = BAND -CO COMPRESS = - 30 30 LZW →LC82070232013160LGN00_B1.TIF LC82070232013160 LGN00_ B2.TIF →LC82070232013160LGN00_B3.TIF LC82070232013160 LGN00_ B4.TIF →LC82070232013160LGN00_B5 .TIF LC82070232013160 LGN00_ B6.TIF →LC82070232013160LGN00_B7.TIF LC82070232013160 LGN00_ B9.TIF →LC82070232013160LGN00_CLD.TIF
//
- O LC82070232013160LGN00_merge.tif 出力複数の名前 - バンドラスタデータセット。 - GTiffの 私たちは、明示的にGeoTIFFの画像フォーマット設定用(出力ではない、これは、ここで必要なのであるデフォルト)。 30 30 -ps 設定し、共通の空間分解能のための全ての入力が。ほとんどのバンドはこの解像度で既にあるが、パンクロマティックバンドを再サンプリングする必要があるから、 15〜30 メートル。 - 別 置き別積み上げ帯域に各入力ファイル。 -CO INTERLEAVE = BAND 作成オプションについてシーケンシャルエンコーディングをバンド(セクト。3.5を参照)。 -CO COMPRESS = LZW 作成オプションのための LZW圧縮された出力(セクト。5.5.4を参照されたいです)。
2つのマルチバンドに基づいて、図17.2cクラウドフリー複合作成する
図17.2a、Bのラスタデータセットを、我々はpktoolsからユーティリティpkcompositeを使用することができる
(セクト。12.2を参照します)。
-i LC82070232013160LGN00をpkcomposite。TIF - I →LC82070232013192LGN00.tif - O LC82070232013_コンポジット。TIF → -bndnodata 8 -srcnodata 1 -dstnodata 0 -CR maxndvi -cb 3 -cb 4
//
-i LC82070232013160LGN00.tif - I ... 入力ラスタデータセットのシーケンス。 - O LC82070232013_composite.tif つの入力データセットに基づいて、出力コンポジットの名前。myinline NODATA値を定義するクラウド・バンドの8ランキング-bndnodata に入力ラスタデータセット(インデックス開始として9バンドから0)。 -srcnodata 1 NODATAを示す値、すなわち、曇った値で、この場合。 - dstnodata 0 書き込まれた値で出力ラスタデータセットに(すなわち、NODATAの場合ならば、両方のシーンが濁っています)。 - CR maxndvi 場合に複数の入力があり、我々は最大NDVI値を持つピクセルを選択し、クラウド自由。 -cb 3 -cb 4 複合バンド:各指標のための赤色及び近赤外スペクトル帯域におけるオンザフライNDVIを計算するマルチバンド入力ラスタ(表17.1を参照のこと)