Google Earth Engine

八郎潟のSAR画像が面白かった 前回見た八郎潟が面白かったので、詳しくみる。 何が面白かったって、あるタイミングで黒い部分が一気に増えたこと。これは田んぼに水を張ったからだと思う。 Sentinel-1 VV の八郎潟（再掲） 時系列グラフを見る 八郎潟の田んぼを３か所適当に選んでポリゴン→Featureにしてから、グラフを書く。 平らだからDESENDINGとAECSENDINGを一緒にしたら、ガタガタのグラフになった。DESとAECで値が違っているみたい。 左：VV、中：VH、右：VH/VV ならば、DESCENDINGだけ選択してグラフを書き直し。 Sentinel-1 VVの時系列グラ…

アニメーションにするといいらしい 前処理について理解が少し進んだので、データを見てみる。 が、何をみてるかよくわからない。教材動画ではアニメーションにするとわかりやすいといっているので、やってみる。 コードはスライドが公開されているので、そこのリンクからコピペ＋場所やら範囲やらを変更したらできる。 www.youtube.com スライド docs.google.com 先生のコードはこちら https://code.earthengine.google.com/3a7843212acec68df670188fdb29e814 アニメーション 岩木山 左：VV、右：VH 雪はあるかな？っていう…

Radiometric Terrain Flatteningがよくわかってない さて、次の補正ですが、何をどうする補正かよくわかっていません。そもそも日本語で何というのか、わかってません。超おおまかな理解としては、斜面だとマイクロ波の入射角が変わるから真の（水平面だった時の）値より大きくなったり小さくなったりする。それを衛星の角度のデータとDEMから補正しようということだと思います。 わかっていないから、まずはデータを見てみようと。そこで、以前からお世話になっているGoogle Earth EngineでSentinel-1の前処理をするスクリプトに登場いただき（下のリンク）、補正前後の様子を…

複数画像ImageCollectionのmean()とmedian()はどうだろう？ SARの勉強シリーズ１回目でも触れたけど、論文によっては複数画像の平均値だったり、一筆の畑の平均値だったりを使うことでスペックルノイズに対処しているのがある。どのくらいノイズ除去効果があって、エッジは残せるのか？ 中心極限定理だろ スペックル１回目で出てきた中心極限定理。分布の形を問わず、平均μ、標準偏差σの集団からN個サンプリングするとその平均はμ、標準偏差σ/√Nになる。 となると、一定期間でピクセルの散乱強度係数が変わらないと仮定したとき、その期間に得られたN枚の画像の平均から得られるピクセルの値は平均…

SAR勉強シリーズを始める ちょっと勉強したいと思ったら結構なボリュームを勉強しないと分からないことが片付かないSAR。ちょっとずつやれば、そのうち理解が深まるんじゃないかと期待してシリーズを始めることに。まず、前回ノイズとエラーの原因として挙げたことから見てみようかと。で、まあ簡単に見られそうな入射角から取り組んでみます。 問題のおさらい Google Earth Engineで提供されているSARデータは、衛星から斜め下方に照射されるマイクロ波のうち、衛星に帰ってきたものの強度（とタイミング的なやつが本当はあるけど入っていない）を測った結果をきれいにして画像にしたもの。衛星からはある幅をも…

教師データはいくつ？ とりあえず10000のランダムなポイントでとしてはじめたんでした。で、サンプリングしてRに持ってきてみたんですけど、データの偏りがすごい。 開放水域 耕作地 高木林 湿地 人工物 水田 草地 低木林 裸地 2723 225 5949 8 627 1083 236 231 82 いやー、どう考えても8点（湿地）だけではどうしようもない。なので100,000くらいに増やしてみたんです。が、そしたら後々のランダムフォレストの分類機を作る手順でデータ多すぎですエラーがでた。ちゃんと説明で、そういうことは止めてねと書いてあるので、するべきではありませんでした。 これを踏まえて、元の…

そう簡単にはいかない いや、わかってましたよ。そんなサクサクはいかないだろうなーと。でも思った以上に考えることやらケースが多くて、勉強・経験不足が露呈します。どこらへんがなのか、というのは追々として、これまでのあらすじを書いてみましょうか。 作業手順 QGISで植生図のカテゴリをリクラス QGISでランダムにポイント発生させて、ポイントに植生クラスを付与、シェープファイルに書き出し シェープファイルをGEEに読み込む GEEで特徴量として使う衛星データを整理して、読み込んだシェープファイルにそのデータをくっつける シェープファイルの各点のデータをcsvで書き出し。ダウンロードしてRでランダムフ…