プロセシングツールボックスを使用して雪の分類を自動化する

大雪水力株式会社(SHL)は、ニューサウスウェールズ州、オーストラリアの大雪山に基盤を置くピーク時発電の会社です。大雪水力株式会社は、全国に多様な発電資産を所有しているにもかかわらず、オーストラリアのインフラの象徴的な作品、大雪山水力発電スキームの管理権が最も有名です。

大雪山水力発電スキームは、ニューサウスウェールズ州の山岳地域での季節的降雪によって可能になります。春の雪解けの水は国家電力市場のピーク需要に対処する電気を生成するために使用されますが、また同時に、南オーストラリア東部のより乾燥した内陸部に灌漑農業のために水を流用しています。

冬月にわたる積雪を分類することで、SHLでは積雪地域のライフサイクルを監視し、積雪地域の歴史的記録を保管し、水収率および貯水の予測のために積雪の統計情報を水路サービスチームへ提供できます。

Map 1.1 Snowy Water Catchment Map

地図1.1 雪に覆われた水の流域地図

背景

毎年冬の間(10月31日から6月1日まで)、SHLはMODISテラ(EOS AM)衛星からできるだけ多くの雲のない画像を得ます。これらの画像は、その後、雪が存在する領域を示すために使用される正規化差分降雪指数(NDSI)として知られているアルゴリズムに摂取されます。このアルゴリズムはいずれか1つの特定のセンサに固有のものではありませんが、MODISの装置がその十分な空間分解能とプロジェクトエリア(1.1雪に覆われた水の流域の地図)の日々のカバーのために選ばれました。

目的

このプロジェクトの目的は、NDSIのワークフローをQGIS処理ツールボックス内のモデルとして自動化することです。必要な出力は、雪水流域内の雪に覆われた範囲を示すベクターポリゴンで、雪に覆われた総面積を平方キロメートルで記述する属性と一緒にしたものです。

方法

NDSIアルゴリズム自体は非常に簡単です(以下式1.1に示す - 正規化差分降雪指数(NDSI))と積雪をベクター形式で抽出するために必要なその後の工程は、プロセッシングツールボックスモデルにうまく適しています。

Formula 1.1 - Normalised Difference Snow Index (NDSI)

式1.1 - 正規化差分降雪指数(NDSI)

このNDSIアルゴリズムは-1と1の間で結果を返し、雪のための閾値は典型的には0.4よりも大きな値です。2014年7月13日について、自然色の画像(地図1.2 -MODISトゥルーカラー)と一緒にNDSI出力の例(地図1.3-NDSI出力)を以下に示します。これらの画像からわかるように、NDSI出力は雪で覆われた領域を他の土地被覆タイプおよび雲から弁別するのは非常に得意です。しかしながら、スペクトルの近赤外部分に吸収特性を有する水によって引き起こされる大きな水域の一部を雪と誤分類することがありえます。NDSIのこの特性は、モデル内のいくつかの段階で水マスクを組み込む必要性を示しています。

Map 1.2 - MODIS true colour

地図1.2 - MODISトゥルーカラー

Map 1.3 - NDSI output

地図1.3 - NDSI出力

Map 1.4 - Recoded "Snow" vs "Non snow"

地図1.4 - 「非雪」対「雪」録画

NDSI出力はGRASSモジュール「r.mapcalculator」を用いて、いくつかの非常にシンプルなバンド・マップを使用して作成されます。NDSIグリッドが利用可能になったら0.4のNDSIしきい値が適用されることが必要で、この結果は「雪」と「非雪」の領域を描くブール値グリッドになるでしょう。その分類を取得するために、別のGRASSモジュール「r.recode」を使用し、地図1.4『「非雪」VS「雪」再コード化』に示す結果を作成しました。

再符号化グリッドは、次いで、ベクターポリゴンデータセットに変換され、雪で覆われた領域が抽出されました(以下に示す、地図1.5 - 抽出された雪)。雪に覆われた領域は、その後雪に覆われた水流域にクリップされた(下に示す地図1.6 - SWC内部の雪)、その後、大規模な水域を出てマスクされている(下に示す地図1.7 - 最終的な雪に覆われた領域)を、そして最後エリア雪に覆われた面積の統計が生成されます。

Map 1.5 - Extracted Snow

地図1.5 - 抽出された雪

Map 1.6 - Snow inside the SWC

地図1.6 - SWC内部の雪

Map 1.7 - Final snow covered area

地図1.7 - 最終雪に覆われた面積

モデルからの出力は、雪に覆われた水流域内の雪に覆われた領域を示すベクターポリゴンです。ポリゴンの属性は、雪に覆われた領域と覆われている流域の割合が含まれます。その全体モデルは下の「画像1.1 - 正規化後差分雪指数(NDSI)、QGIS処理ツールボックスモデル」に見ることができます。

Image 1.1 - Normalised Difference Snow Index (NDSI) QGIS processing toolbox model

画像1.1 - 正規化後差分雪指数(NDSI)QGIS処理ツールボックスモデル

結果/結論

QGIS処理ツールボックスはこの雪の分類のようなワークフローを自動化するための強力な方法です。上に示したモデルは、雪分類プロセスを上手く自動化し、さらに統計的アプリケーションに可視化し、入力に使用できるGISデータセットを提供します。再現性および一貫性の分類ツールを提供しながら、SHLで使用されていた以前の方法に関連づいた主観的なユーザエラーを上手く除去します。いくつかの季節の結果をここで見ることができます。

制作者

アンドリュー・ジェフリーは、オーストラリアのニューサウスウェールズ州クーマに拠点を置く大雪水力株式会社に勤務するGISアナリストです。氏は科学の学士(BSC)、空間情報システム専攻で、ワガワガのチャールズ・スタート大学を2004年に卒業しました。大雪水力株式会社にいる間、氏は人工降雨実験のようなエキサイティングなプロジェクト、および組織内の日常業務へに空間的なサポートを提供することができました。