[Google Earth Engine入門]衛星画像のロードと地図へのオーバーレイ

概要

特定の場所や、特定期間の衛星画像をロードしてMap Viewへオーバーレイする方法です。

位置を指定する

後ほど衛星画像を絞り込むのに使う、位置情報を先に設定しておきます。

・ジオメトリツールを使う
緯度経度を調べてコード内に記述するのは手間ですが、Map Viewの右上にあるジオメトリーツールを使用すると簡単に地図上の緯度経度を取得することができます。

ジオメトリーツール

地図上にマーカーを設置

コードエディタにインポート

衛星データを参照するためのオブジェクトを生成する

LANDSAT 8号の衛星画像を参照するオブジェクト-イメージコレクションを生成します。
このイメージコレクションを通じて、Google が保持している2011年から現在までのLandsat 8号の衛星画像にアクセスすることができます。

print文を使用することでイメージコレクションの詳細情報をコンソールに出力できます。
(console.logみたいなもの)

ただし、データ量が多すぎる(5000エレメント以上)とコンソールに表示することができません。

sizeメソッドを使用することでイメージコレクションから参照できるエレメント-イメージの数を習得することができます。

約83万枚の衛星画像がこのイメージコレクションの中に収められているのがわかります。
(この数は、日々増加しているのであくまで現時点での数です)

このままでは衛星画像の数が多すぎるので、条件をつけてフィルタリングします。

フィルタリングする

まずは、衛星画像が撮影された期間で絞り込みを行ってみます。

・期間で絞り込む

フィルタリングを行った結果をコンソールに出力してみると、約83万枚あった衛星画像が約15万まで絞り込まれているのがわかります。

・場所で絞り込む

さらに冒頭で作成したジオメトリオブジェクトを使って、特定位置の衛星画像だけに絞り込みます

期間と場所を指定することで、一気に15枚までに絞り込むことができました。

5000枚以下になったのでコンソールに詳細情報を出力することもできます。

コレクションから画像を取り出し地図にオーバーレイする

15枚までに絞り込んだイメージコレクションの中から、とりあえず先頭の一枚を取り出して地図上にオーバーレイしてみます。

衛星画像を地図上にオーバーレイすることができましたが、ほぼ真っ暗でなにも見えません。
これは、衛星の観測バンドの値が適切なRGB値として適用されていないためです。
次の項では、可視域のバンド値を使ってトゥルーカラーと呼ばれる可視化を行います。

参考：
よくみかける画像へ可視化する

ビジュアライズする

landsat 8号の場合、バンド(2、3、4)がそれぞれ(青、緑、赤)の波長データとなっているので、bandsオプションで[“B4”, “B3”, “B2”]をRGBに適用します。

全サンプルコード

アドベントカレンダー

この記事は、Google Earth / Google Earth Engine Advent Calendar 2018への参加投稿です。