【satellite image】 人工衛星によって得られた観測データを画像化したもの。
今回は衛星から捉えた雲のお話についてです。 気象衛星ひまわり 可視画像 赤外画像 水蒸気画像 【まとめ】学習の要点 参考図書・参考URL 気象衛星ひまわり 日本ではさまざまな観測機器を用いて気象の観測を行っています。 例えば日本全国に設置されているアメダスだったり，気象レーダーだったり，温度計などをつけて上空の観測をする気球などなど。 その中でも宇宙空間から地球の観測を行う気象衛星は個人的にはもっともロマンがあって，もっともスケールの大きい気象観測手法の一つでしょう。 2023年2月現在運用されている日本の気象衛星はひまわり9号です。高度36000kmと地上のはるか遠くを飛行しており（地球の直…
2022年12月6日、宮城県議会（令和4年11月定例会）本議会一般質問にて、遠藤伸幸議員（公明党県議団）が、「水道事業の諸課題について」ということで、水道管の老朽化による漏水状況把握にAIや衛星画像解析を活用してはどうか？ と提案しました。 宮城県議会HPより miyagi-pref.stream.jfit.co.jp 遠藤伸幸議員（公明党県議団） 最後に、大綱６点目、水道事業の諸課題について伺います。 近年、全国的に水道管の老朽化による漏水や断水といった事故が多発しています。 昨年10月、和歌山市では、水管橋が崩落し、１週間にわたり市内の４割にあたる６万個が断水。今年5月には、静岡県菊川市で…
フィルタリングの効果を見てみたい ちょっと前に公開されたGoogle Earth Engineで使えるSentinel-1の前処理法を使ってフィルタリングする。以前の続きです。 前処理法はこちら https://github.com/adugnag/gee_s1_ard フィルタリングの方法は5種類 スペックルノイズを除去する前処理法は色々あって、今でも新しい方法が提唱されたりしているようですが、ここで紹介されているのは５種類。それぞれ１枚の画像だけをつかうのか、複数枚で時間方向にもフィルタをかけるのか？で5x10種類の処理が可能。 Boxcar フィルタ Gamma maximum a-po…
SAR勉強シリーズを始める ちょっと勉強したいと思ったら結構なボリュームを勉強しないと分からないことが片付かないSAR。ちょっとずつやれば、そのうち理解が深まるんじゃないかと期待してシリーズを始めることに。まず、前回ノイズとエラーの原因として挙げたことから見てみようかと。で、まあ簡単に見られそうな入射角から取り組んでみます。 問題のおさらい Google Earth Engineで提供されているSARデータは、衛星から斜め下方に照射されるマイクロ波のうち、衛星に帰ってきたものの強度（とタイミング的なやつが本当はあるけど入っていない）を測った結果をきれいにして画像にしたもの。衛星からはある幅をも…
そんな簡単にはいかない 簡単そう、お手軽そうと思って始めましたが、予想より３倍くらいは大変でした。なんで大変だったんかなと思いだしてみます。 分析に時間がかかる サクサク進まないと感じた理由が時間がかかること。GGEで8000点を使って石川県サイズのランダムフォレストをして、書き出すとおよそ30分かかる。で、結果を都度見ながらやりたいなとMap.addLayerを使うとタイムアウトでできない。だからお手本論文ではRでいろいろ実験する流れになってたんだなと思った。 Ｒでもランダムフォレストにそれなりに時間はかかるけど、随分と時間の節約にはなる。必要なサンプル数の計算とかループで回してしまえば待っ…
やっとGEEに到達 あーだこーだとやって、これで行きましょうと行くまでの体感が長かった。実際はそんなにかかってないんだけど。 というわけで、GEEで分類します。見本があるのでそれに倣って書けば大丈夫。 developers.google.com 実際はこう。 //ポイントデータからラスタ情報を取得して教師データを作成 var trainingData = trainingLayers.sampleRegions({ collection: pointsInt, // 教師データのポイントのFeatureCollection（シェープファイル） properties: [label], // 分…
教師データはいくつ？ とりあえず10000のランダムなポイントでとしてはじめたんでした。で、サンプリングしてRに持ってきてみたんですけど、データの偏りがすごい。 開放水域 耕作地 高木林 湿地 人工物 水田 草地 低木林 裸地 2723 225 5949 8 627 1083 236 231 82 いやー、どう考えても8点（湿地）だけではどうしようもない。なので100,000くらいに増やしてみたんです。が、そしたら後々のランダムフォレストの分類機を作る手順でデータ多すぎですエラーがでた。ちゃんと説明で、そういうことは止めてねと書いてあるので、するべきではありませんでした。 これを踏まえて、元の…
オーストラリア公共放送ABCニュース4/5 2017年、オックスフォード大学で中東歴史専門家の2人の考古学者は イスラエルとヨルダン川西岸の衛星画像で奇妙なことに気づきました イスラエルで撮影された写真にアクセスしようとしましたが、 それらがすべてぼやけていて、探していた地面の詳細を理解するのに 十分な品質ではないことに気づきました 彼らが見つけたのは、イスラエルの核施設でした 衛星画像の取り決めで、アメリカのキルビンガマン改正・Kyl-Bingaman Amendmentでアメリカと同盟国の重要な軍事施設は、衛星画像は鮮明に できなくなっていたのです 同じ衛星画像を使うGoogleEarth…
衛星画像のデータが大きいと画像処理に何かと負荷がかかったり、時間がかかる。そこで、観測地域が限られている時は画像の切り取りを行う。 衛星画像の切り取り方法は主に二つある。 ①緯度経度を使用する場合 ul_lon=135 #画像左上経度 ul_lat= 40 #画像左上緯度 lr_lon=145 #画像右下経度 lr_lat= 30 #画像右下緯度 img_extent = [ul_lon, lr_lon, lr_lat, ul_lat] plt.imshow(img, extent=img_extent) ②緯度経度を使用しない場合（画素の数で切り取り） 下の画像で、青部分は元画像、橙部分は切…
Sentinel2Aの処理レベル Sentinel2Aは「1Cレベル」という処理レベルになっている 簡単にまとめると「ジオメトリック、ラジオメトリック、幾何補正などを行なったのちに 大気上端反射率を計算している」状態。 そのため、ほとんど補正を自分の手で行わずに次の段階に進むことができる。 読み込むときにはバイナリデータ（0-255）になっている。 Sentinel2Aのデータダウンロード 私は普段earthexplorerで衛星データのダウンロードを行なっている。 ダウンロード方法は多くのネット記事があるので省略。 アドレスは下記の通りである。 https://earthexplorer.u…
note.com
www.47news.jp
www.cnn.co.jp
www.afpbb.com
mainichi.jp
www.cnn.co.jp
www.cnn.co.jp
mainichi.jp
www.itmedia.co.jp
