SARデータから活用可能なマップへ：Pythonによるオープンソース洪水検知パイプラインの構築

洪水は壊滅的な被害をもたらす可能性があります。急速に発生し、地域社会の復旧を極めて困難にします。気候変動により極端な気象イベントの頻度が増加している現在、効果的な洪水検知が不可欠となっており、このテーマは特に重要です。

SARデータとは？

合成開口レーダ（SAR）は、地表の画像を生成するために使用されるレーダの一種です。快晴の空に依存する光学衛星とは異なり、SARは雲や雨を透過して観測することができます。そのため、洪水検知において特に有用です。ただし、SARデータにはノイズが多く複雑になる場合があり、分析が難しくなることがあります。

洪水検知パイプラインの構築

この記事では、生のSARデータを活用可能な洪水マップへと変換するパイプラインの作成手順を順を追って説明します。ここでは、対処すべき3つの主な課題を紹介します。

課題1：ゼロ値境界の処理

Sentinel Application Platform（SNAP）を使用してSARデータをエクスポートする際、画像内にゼロ値の境界が見られることがあります。これらの境界はデータ分析を歪める原因となります。

解決策：

1. ゼロ値をNaN（Not a Number）として処理する。
2. データをデシベル（dB）に変換する。
3. 有効な情報のみが含まれるようにデータをクロップする。

以下は、データをクロップするためのサンプルコードスニペットです：

rows = np.where(~np.all(np.isnan(dB_full), axis=1))[0]
cols = np.where(~np.all(np.isnan(dB_full), axis=0))[0]
dB = dB_full[rows.min():rows.max() + 1, cols.min():cols.max() + 1]

課題2：実際の浸水域の分離

-17 dBというシンプルな閾値で水域全体を検出できますが、河川や湖などの恒久的な水体も検出してしまいます。これにより、洪水検知において偽陽性が発生する可能性があります。

解決策：

1. OpenStreetMap（OSM）を使用して恒久的な水体データを取得する。
2. このデータをSARグリッドに合わせてラスタライズする。
3. 検出されたマスクから恒久的な水体を減算し、浸水域を分離する。

恒久的な水体データを取得する方法は以下のとおりです：

osm_water = ox.features_from_bbox(bbox=bbox, tags={"natural": "water"})
flood_only = mask_clean(~binary_dilation(permanent_water, iterations=2))

課題3：QGISの描画問題への対処

透過GeoTIFFをエクスポートする際、QGISが透過ピクセルを完全な黒として描画してしまうことがあります。これにより、データを正しく可視化するのが困難になる場合があります。

解決策：

1. RGBA（赤、緑、青、および透過用のアルファを含む）を使用する代わりに、シングルバンドラスタを作成する。
2. QGISが透過として扱う専用のNoData値（255）を使用する。

これを設定する例は以下のとおりです：

overlay = np.where(final_flood_raster == 1, 1, 255).astype(np.uint8)

結果

これらの手順に従うことで、パイプラインは以下を出力します：

• ベクターデータ： 正確な面積計算が含まれるGeoJSONまたはShapefile。
• インタラクティブマップ： 衛星写真とストリートビューを切り替えられるFolium HTMLマップ。
• QGIS対応ラスタ： インポート後すぐに機能するクリーンなオーバーレイ。

まとめ

SARデータを使用した洪水検知パイプラインの構築は価値あるスキルです。複雑なデータを災害対応や計画に役立つ有用な情報へと変換するのに役立ちます。

メリット

• 正確な洪水検知を提供する。
• オープンソースツールを活用しているため、アクセスしやすい。
• 災害管理の取り組みを強化する。

デメリット

• PythonおよびGISツールに対する習熟が必要。
• 処理に多くのリソースを消費する場合がある。
• 初心者にとっては初期設定が複雑な場合がある。

注意事項

この記事は教育目的で提供されています。コード内のプレースホルダー値は必ず実際のデータに置き換えてください。情報に依拠する前に、必ず元の情報源と照らし合わせて主張を検証してください。

よくある質問

• SARデータとは何ですか？ — SARデータとは地表の画像を生成するために使用されるレーダデータであり、曇天時にも有効です。
• 洪水時に光学衛星データが役立たないのはなぜですか？ — 光学衛星は快晴の空に依存しますが、洪水時には快晴でないことが多いためです。
• 洪水検知パイプラインの目的は何ですか？ — このパイプラインは、生のSARデータを活用可能な洪水マップに変換します。
• パイプラインにはどのプログラミング言語が使用されていますか？ — 洪水検知パイプラインの構築にはPythonが使用されています。
• このパイプラインにおけるOpenStreetMapの役割は何ですか？ — OSMは恒久的な水体に関するデータを提供し、浸水域の分離を支援します。
• QGISとは何ですか？ — QGISは、空間データのマッピングや分析に使用されるフリーでオープンソースの地理情報システムです。

タグ

#python #opensouce #flooddetection #SAR #GIS #QGIS #OpenStreetMap #climatechange #floodmapping #dataanalysis