Challengeが今年も開催されました。COVID-19によるパンデミックを鑑み2年連続でオンライン開催となりました。Streamチーム(ソーシャルグッド賞受賞)として参加した筆者の体験を基に、衛星データを用いて海洋や河川のゴミが見つけられるのかという試みについて紹介いたします。とプログラミング教育を生業にしている者です)。枝葉末節において誤解に基づく記述があるかもしれません。1.宇宙から海洋デブリを見つけろ!海鳥や海獣に絡まる釣り糸、ウミガメの胃から出てきたプラスチック袋、人の目には見えないマイクロプラスチックなど心を痛める海のゴミ(海洋デブリ)の話は、悲しいことに数えきれないほどあります。その海洋デブリの発生源はどこなのでしょうか?海上や海岸での人為的な投棄は発生源のひとつです。人間は海岸だけに住んでいるわけではありません。内陸で行われたゴミの投棄はやがて河川によって運ばれ海洋デブリの一群に合流します。海洋デブリの約20%が河川を経由したものだという試算[3]がある一方、河床にまぎれ運ばれるものは正しく試算することが難しいという指摘[4]もあります。DEBRIS」(豊橋会場訳:海洋プラゴミに対するAI/機械学習の活用)に取り組みました。チームメンバーが2人とも内陸住まいでしたので、海だけだと難しく、河川も含めたデブリ検出を試みました。JAXA我々の目や衛星のカメラは、電磁波を捉えることで物を知覚します。X線、γ線に分けられます。我々は目で電磁波の波長の違いを色として認識しますが、それは可視光という非常に狭い波長帯域の中での事です。物理学者や衛星データ解析者達は便宜的に色という言葉を可視光を超えた範囲でも使用します。この様々な色の組み合わせによって物質や現象の特徴をより分かりやすく1つの値としたものが指標です[5]。今回は、先行研究[6]を参考に、Sentinel-2のデータを基に以下の3指標による海洋デブリ検出精度の比較を行いました。(NDVI)各指標の計算式は以下となります。NIRは近赤外線(SentinelのB8バンド)、RE2は赤色の長波長側(同B6バンド)、SWIR1は短波赤外線(同B11バンド)です。3指標とも赤色と赤外線の組み合わせから算出されます。FDIとPIはプラスチックゴミに感度があると提案されている指標、NDVIは植物の活性度を示す指標です。3指数とも可視光の赤から赤外線の複数の波長帯の観測量を組み合わせたものです[7]。プラスチックゴミと一見関係なさそうなNDVIを指標としたのは、藻などの自然な浮遊物と海洋デブリを見分けるためです。今回は赤をPI、青をFDI、緑をNDVIにそれぞれ割り当てることで、PIとFDIでは高い値を取り、かつNDVIでは低い値を取る、つまり赤と青が強く見える(=紫色になる場所)場所=海洋デブリがある場所、と仮説を立てて可視化に挑みます。この3指標をただ眺めるだけでは元々ゴミ以外の何らかの理由で元々高い値をとっていた高い場所をデブリか流れ藻がある場所と誤解してしまうかもしれません。NDVIと見比べる以外にも、調査対象期間と異なる年を基準として、その基準からの変化の割合を求めました。この3指数の変化をRGBの各チャンネルで表示することで、ゴミがありそうな場所を直感的に判断出来るようにしました。3.ゴミじゃないものは見えた手法に目途が立ったところで、実際に衛星画像でみてみましょう。最初はゴミではないと思われるものから始めます。JAXAのしきさい衛星による観測結果[8]にて流れ藻があると指摘される領域の近辺で、この3指標の解析結果を可視光と共に可視化してみました。図2:九州沖の流れ藻。https://code.earthengine.google.com/af6794541999ded1226aab6e07c1dc21図2でNDVIが局所的に高い緑色の領域は流れ藻だと推測されます。緑に混じり赤い部分が点在しますが、そこはPIのみが高くなっている場所です。この流れ藻の領域は所謂「潮目」にあり、可視光(図2上段左)で濃い青の領域と合致します。このような潮目は流れ藻だけでなく海洋ゴミも集積しやすいと言われていますが、この図からはFDIが局所的に高い領域は見当たりません。海上保安庁が公開している海洋状況表示システム「海しる」にて、海岸清掃時のゴミの量のデータ(環境省H25年)を確認した際に、たまたまカーソルが向かったゴミが比較的多かった地域です。先ほどと同様に可視光と3指標の変化で見てみました。図3:松江市島根町加賀の海岸。図2と同様に、上段左:2021年7月1日から8月31日、上段右:2021年7月1日から8月31日と2020年の同期間、下段:同。https://code.earthengine.google.com/e1826123020eb0166d4be84ffa1a5cee図3の黄色の破線は防波堤、矢印の先はPI(赤)が目立つ領域を指しています。PIが目立つ領域は可視光では白く映っています。湾の防波堤の内側なので、ゴミがたまりそうだと思ったのですが、可視光でこれほど白っぽいという事は砂地という可能性もあります。実際、衛星データとの大きさの比較は出来ていませんが、Googleストリートビューで小さな砂浜があることが確認できています。単純に新たな砂地に反応したのか、それともその砂浜が特にプラスチックゴミを多く含んでいるのかは実地調査が必要なところです。様々な所をこの手法で眺めたところ、新たな砂地に反応してそうな感想を持っています。海洋学研究室)に大雨などでデブリが多く流された後の衛星データを見ると海洋デブリ検出には良いとアドバイスを頂きました。最近の大雨で該当しそうな事案として思い出したのが2021年7月3日に発生した熱海の土石流でした。図4:熱海の土石流跡。https://code.earthengine.google.com/f174e2012d8d57730a28d5fae595b306土石流の流れた先に防波堤(黄破線)があり、その先に青い領域が広がっていることが認められます(矢印先)。土砂崩れ直後の空撮写真を確認したところ、土石流の濁った水と外海の水との切れ目あたりだと分かりました。問題はこれが本当にプラスチックゴミなのかということです。本当は衛星データと同時間に現地でサンプル調査をすべきです。メンバーは熱海近辺どころか海すら遠い内陸住まいでした。そこで海は諦め河川で同様の事が出来ないか挑戦してみました。筆者がアクセス可能な範囲でそれらしい変化の見えた場所として京都の嵐山を流れる桂川に注目しました。図6:京都嵐山の桂川。上段左:2021年10月3日正午頃に矢印の先あたりから見た桂川(筆者撮影)、以下図2と同様に、上段中:2021年7月1日から8月31日、上段右:2021年7月1日から8月31日と2020年の同期間、下段:同。https://code.earthengine.google.com/ced9745b8d20bf62e35b3f75ae4dee2e図6の3指標を確認すると、葛野大堰より上流(左側)で紫色が目立つ領域があります(図内の矢印の先)。Indexが共に前年よりも上がっているということです。先の加賀の例で砂地だろうと推測したシグナルと似たところもあります。Challenge期間にその場所に行きました(図6上段左)。このあたりは保津川の船着き場かつ手漕ぎの貸ボートが行き交う場所です。緊急事態宣言明けだったからか、観光地らしくゴミ1つ見当たらない水面で、特に新たな砂地が出現した様子もありませんでした。Indexの方が信頼できそうな手ごたえを得ました。これらを実証していくには、圧倒的に実地データが足りませんでした。日本の河川にゴミが溢れていればもう少し検証がしやすかったかもしれません。河川が綺麗な事はもちろん良い事なのですが、この様な研究では正解情報を得るのに少し難儀する環境です。熱海のケースの様に災害発生時の調査は難しいでしょうが、衛星が回ってくる晴れのタイミングで実際に海上のゴミのサンプル調査が出来れば、研究がより進められるだろうと思います。衛星データの拡充も大切ですが、同様に地上でのその場観測も重要だと実感したオンラインイベントでした。海洋学研究室)と牟田梓様(さくらインターネット)にこの場をお借りしまして御礼申し上げます。#マンガでわかる衛星データ」(https://sorabatake.jp/5192/)をご覧下さい。EngineでNDVIを計算するには、以下の様な関数を組み、`ee.ImageCollection.map(calcNDVI)`で適用します。選択した期間に複数の観測が存在する可能性がありますので、必ず`ImageCollection`の状態で指標を計算し、それから`median()`などで時間方向を縮退します。しなければ波長帯により異なる時間の観測を使用する可能性があります。Engineで画像を取得する事にかなり難儀しました。専門、通俗的には宇宙天気とも。研究の傍ら、宇宙をからめた様々な教育・アウトリーチ活動を行う。Challengeを通して地上観測も扱えるようになり、衛星データのオールラウンダーに。将棋のアマ二段目指し勉強中。2022/12/13衛星データを使って市区町村の緑地率を算出して平均寿命・自殺率との相関を調査してみた緑地の多さが認知症や脳卒中と関係があるとの論文を見つけました。Enginとオープンデータを使って、市区町村の緑地率と平均寿命・自殺率との相関関係を調査してみました。2022/7/28YAMAPのリアルタイム紅葉モニターの画像を参照して衛星データから紅葉予測チャレンジ!~データ確認編~春の桜と並ぶ、日本の四季の楽しみ方である秋の紅葉について、衛星データで見頃を予測できないかと考えた宙畑編集部。YAMAPのリアルタイム紅葉モニターの画像をお借りしながら、アルゴリズム制作に挑戦してみました。2022/11/18YAMAPのリアルタイム紅葉モニターの画像を参照して衛星データから紅葉予測チャレンジ!~データ解析編~春の桜と並ぶ、日本の四季の楽しみ方である秋の紅葉について、衛星データで見頃を予測できないかと考えた宙畑編集部。YAMAPのリアルタイム紅葉モニターの画像をお借りしながら、高尾山の紅葉時期予測を行いました。