衛星データ解析の全体探索

何を計画し、何をやり、何が見えたか

2026年3月 / 松島 宏佑

Space Tech Accelerator 共有資料

背景と問題意識

衛星データ解析の全体像を理解したいと思いました。

衛星解析は、データの取得・加工・解析・可視化まで多くのステップがあり、使うセンサーやモデル、前処理の方法によって結果が大きく変わります。パームオイルなど特定テーマに絞る前に、まず全体像を身体感覚として掴みたい

特に知りたかったのは、AIがある時代に、これらがどの程度民主化されるのか、されないのかということです。

AIエージェントの能力が急速に上がっている。衛星データは無料で大量に手に入る。この2つが交差したとき、衛星解析はどこまで「誰でもできるもの」になるのか。それとも、専門家でなければ触れない領域が依然として残るのか。

座学ではなく、実際にデータを動かして一つずつ確認していく方法で、1週間探索しました。

学習計画

設計原則

原則内容
変数を一度に変えない各日で「固定する変数」と「動かす変数」を分け、差が何に起因するかを明確にする
2軸で広げる縦軸: データにどこまで遡るか / 横軸: 解析をどこまで深めるか

5日間の計画

Dayテーマ固定動かす時間
1光学画像分類(EuroSAT)データ解析手法 L1→L2→L35.5h
2SAR洪水検出(Sen1Floods11)解析の型センサー(光学→SAR)6.5h
3データパイプライン解析(L1固定)データ側(前処理・指標・センサー種)5-8h
4AI手法総当たりデータ(EuroSAT)モデル(古典ML→DL→基盤モデル)8.5h
5統合・バッファ未消化 + 時系列合成9h

実行体制

役割担当
問いを立てる / 結果を判断する人間(松島)
カリキュラム設計人間 + Claude
コード生成・実行・可視化・レポートAIエージェント(Codex)

全37タスク、合計33-38時間。 Google Colab(無料枠 T4 GPU)で実行。人間はコードを一行も書いていない。

到達目標チェックリスト

到達状態Day
L1→L2→L3の精度差を数値で説明できる1
光学とSARの違いを説明し、使い分けられる2
前処理の影響を「ありなし」の比較画像で説明できる2-3
STAC APIでデータを検索・取得できる3
20種の指数を用途別に使い分けられる3
各モデルの特徴と精度を一通り体験し、違いを説明できる4
教師データの形式と作成コストを理解している4
解像度・観測条件が精度にどう影響するか体感している4-5
P1〜P8の全パイプラインを1本通せる5

やったこと一覧

4

日間

41

タスク

155+

ビジュアル

278

ファイル

#タスク名主な学び
Day 1: 光学画像分類(5タスク)
1-1NDVI閾値分類(L1)精度64.4%。
1-2ResNet50 事前学習モデル推論(L2)精度95.7%(+31.3pt)。
1-3ResNet50 ファインチューニング(L3)精度98.5%(+2.8pt)。
1-4L1 / L2 / L3 精度比較 最大ジャンプはL1→L2(+31pt)。
1-5全モデル比較(8手法)ViT=ResNet50 FT=98.5%(同率首位)。
Day 2: SAR洪水検出・InSAR(19タスク)
2-1SAR閾値法による洪水検出(L1)F1=0.396。
2-2事前学習モデル洪水検出(L2)F1=0.167。
2-3U-Net学習による洪水検出(L3)F1=0.240。
2-4SAR洪水検出 L1/L2/L3 比較 Day1と真逆: L1(0.396) > L3(0.240) > L2(0.167)。
2-5NDVI 時系列解析農地・森林・都市の年間NDVI変動を季節分解。
2-6SAR 変化検出(3手法)差分法・比率法・CVAの3手法でPre/Post変化を検出。
2a-1Sentinel-1 データ取得と理解asf_searchによるSentinel-1データ検索。
2a-2HyP3 RTC前処理パイプラインHyP3 RTCによるSAR前処理(地形補正等)。
2a-3JRC恒常水域マスクJRC Global Surface Waterで恒常水域と洪水域を分離。
2a-4大規模SAR処理タイリング・並列処理・モザイク化による広域SAR解析パイプライン。
2b-1SAR可視化・スペックルフィルタ比較リニア vs dBスケール、スペックルフィルタ5種の効果比較。
2b-2Otsu閾値+水域マスクOtsu自動閾値・固定dB閾値最適化・恒常水域マスクの3手法を比較。
2b-3偏波活用+Pre/Post変化検出VV/VH偏波情報を活用した変化検出。
2b-4SAR 3ch U-Net洪水検出VV/VH/VV-VHの3チャンネル入力U-Net。
2b-5全手法比較+判断フローチャート 5手法を統一条件で比較。
2c-1InSAR基礎SLC・干渉縞・コヒーレンスの概念理解。
2c-2コヒーレンス変化検出Pre/Postコヒーレンス差分で地震被害域を推定。
2c-3DInSAR変位解析Unwrapped PhaseからLOS変位マップを生成。
2c-4SAR統合まとめ Day2全13タスクを「3つの物理量」(後方散乱/コヒーレンス/位相差)で統合整理。
Day 3: データパイプライン(8タスク)
3-0前処理比較実験 大気補正なしだとNDVIが254.7%過小評価。
3-1スペクトル指数20種 同じ画像からバンドの組み合わせを変えるだけで20通りの情報。
3-2バンド合成・雲マスク・大気補正10種バンド合成、3種雲マスク、L1C vs L2A大気補正を比較。
3-3DEM地形解析SRTM 30mで関東平野の地形解析。
3-4水域・海洋解析琵琶湖の水域抽出、東京湾のクロロフィル推定、Sentinel-3 SSTとの比較。
3-5大気解析(TROPOMI 7プロダクト)★ 地表ではなく大気の化学成分を測定する全く別種の衛星データ。
3-6熱赤外・ヒートアイランドLandsat TIRSで東京の地表面温度を実測。
3-7都市拡大検出+夜間光2018→2024年のNDBI差分で東京の都市拡大を検出。
Day 4: AI手法比較(8タスク)
4-1ファウンデーションモデル4種比較 ResNet50(96.0%) ≈ DINOv2(95.0%) >> CLIP(34.6%)。
4-2セグメンテーション3手法比較U-Net mIoU 28.1%。
4-3物体検出DeepForest 55本検出 vs YOLOv8 0検出。
4-4古典的機械学習6手法比較RF/SVM/k-NN等をPCA圧縮特徴で比較。
4-4bピクセルレベル古典MLピクセル単位の13バンドスペクトル分類。
4-5深層学習4手法比較ViT 99.1% > ViT+LoRA 98.8% > ResNet FT 98.4% > Scratch 91.…
4-6データ効率性実験 教師データ0枚(35%)→5枚(75%,+40pt!)→100枚(94%)→全量(98%)。
4-7解像度感度実験 分類は30mで95%維持、物体検出は0.5mで信頼度半減。

★ = 特に重要な発見があったタスク。+ ドローン×衛星統合リサーチ(6ドキュメント、参考文献100件超)

プロセスマップ: 衛星解析の全体フロー

全工程を手を動かして通した結果、衛星解析の統一プロセスフローとして整理しました。

PRE-WORKFLOW(ワークフロー前段階) 0a. ユースケース定義(モニタリング / マッピング / 単発分析) 0b. アーカイブ存在確認(STAC照会) 0c. タスキング判断(アーカイブ利用 or 新規撮影) 0d. ライセンス・倫理確認 ① 目的定義(要件仕様の正式化) ↓ ② センサー選定最初の大きな分岐(データ融合計画を含む) ↓ ③ データ取得 3a. データ探索(STAC API) 3b. アーカイブDL or 商用タスキング 3c. データ品質初期確認 ↓ [ARDか?] → Yes → ④の大部分をスキップ → ⑤へ → No → ④前処理をフルに実施④ 前処理 4a. 幾何補正(オルソ化、DEM選択、リサンプリング) 4b. 放射補正(DN→放射輝度→TOA反射率) 4c. 大気補正(BOA反射率) 4d. 雲/影マスク 4e. BRDF正規化(多時期合成の場合) 4f. クロスセンサーハーモナイズ(Landsat+Sentinel混用の場合) 4g. SAR固有: 軌道補正→熱雑音→スペックルフィルタ→RTC 4h. 時系列補間・ギャップフィル 4i. マルチセンサー融合(SAR+光学等) 4j. 時系列合成(L3: 季節合成・年合成) ↓ ⑤ 解析最大の分岐点 5a. 特徴量エンジニアリング(スペクトル指数、時間メトリクス) 5b. ラベル収集・アノテーション 5c. チップ生成 + 空間的 train/val/test 分割 5d. 解析実行(→ 次スライドで詳細) 5e. タイルスティッチング(推論結果の統合) ↓ ⑥ 後処理・検証 6a. ノイズ除去(モルフォロジー、最小面積フィルタ) 6b. ラスタ→ベクタ変換 6c. ピクセル別品質フラグ・不確実性レイヤー 6d. 精度評価(混同行列、面積バイアス補正) 6e. フィードバックループ: 精度不足 → ⑤に戻る⑦ 地図化・重ね合わせ(既存地図/行政区画/インフラとオーバーレイ) ↓ ⑧ Web配信(タイルサーバー / ダッシュボード / API / アラート) ↓ ⑨ 運用モニタリング(定期ワークフローの場合) 9a. レイテンシ・SLA監視 9b. 自動変化アラート 9c. データ分布ドリフト検知 9d. モデル精度監視・再学習トリガー

重要な事実: 解析(⑤)は全体の一部に過ぎず、実際の作業時間の6-7割は①〜④のデータ側が占める。

センサーの全体像

センサー種別代表衛星バンド解像度回帰日コスト得意なこと
光学Sentinel-21310-60m5日無料植生、土地被覆、水域、都市、農業
光学+熱赤外Landsat 8/91130m(熱100m)16日無料温度測定、長期変化(1972年〜)
SAR (C-band)Sentinel-12偏波5-20m6日無料雲貫通、洪水、地盤変動
SAR (L-band)ALOS-24偏波10-25m14日有料*森林透過、バイオマス、水田
大気Sentinel-5P3.5-7km1日無料NO2、CO、メタン等
夜間光VIIRS375-750m1日無料経済活動、電化率、漁船
高解像度商用Planet/Maxar4-80.3-5m1-数日高額個別建物、車両、樹冠個体

見たいもの → 使うセンサー・指数 対応表

見たいもの使うバンド/指数最適センサー最低解像度
植生の健康度NDVI (NIR, Red)Sentinel-210m
植生ストレス・病害NDRE (Red Edge)Sentinel-220m
水域の境界NDWI/MNDWISentinel-210-20m
都市・不透水面NDBI (SWIR, NIR)Sentinel-2/Landsat10-30m
火災跡NBR/dNBRSentinel-2/Landsat20-30m
地表面温度TIR (B10/B11)Landsat 8/9100m
洪水域VV後方散乱低下Sentinel-110-20m
地盤沈下・変動InSAR位相差Sentinel-1 SLC5m
森林バイオマスL-band後方散乱ALOS-210-25m
大気汚染(NO2等)UV-SWIR吸収Sentinel-5P3.5-7km
個々の建物・樹冠RGB高解像度Maxar/Planet<1m
完全パイプライン
5段階の完全パイプライン(Day 2で作成)
SAR手法決定木
SAR: 目的別の手法選択ツリー

プロセスマップ: 解析の分岐判断フロー

⑤解析 は最大の分岐点。「何を使うか」の判断フローを整理しました。

Q1: その問題の指数(NDVI等)は存在するか? → Yes → L1で試す → 精度十分? → Yes → 完了 → No → Q2へ Q2: 学習済みモデルは存在するか? → Yes → L2 そのまま推論 → 精度十分? → Yes → 完了 → No → Q3へ → No → Q4へ Q3: 少量のGTで精度は回復するか? → 5-50枚のGTでFTしてみる → 精度十分? → Yes → 完了(L2.5) → No → Q4へ Q4: GT収集のコスト/期間は許容できるか? → Yes → GT収集 → L3 フルトレーニング → 学習・反復 → 完了 → No → 問題設定を見直す or 商用サービスを検討

解析手法の全体像

レベル手法条件精度目安自動化専門家
L1: 指数・閾値NDVI, NDWI, dB閾値等物理量に直結する指数が存在60-70%○ 完全可能不要
L2: 既存モデルpretrained ResNet, Prithvi等学習済みモデルが存在 AND 条件が近い80-96%○ 完全可能不要
L2.5: ドメイン適応Few-shot FT, 転移学習モデルはあるが地域・条件が異なる85-95%△ 半自動判断に必要
L3: フルトレーニングU-Net, ViT, 古典ML学習対応するモデルが存在しない90-99%△ 学習は自動全工程で必要
ゼロショットRemoteCLIP, SAM「まず試す」探索用30-50%○ 完全可能不要

★ L2.5が実務で最頻出。例: 欧州で学習した土地被覆モデルを東南アジアに適用→「水田」クラスがない、「パーム林」が「森林」に誤分類される。

教師データ(GT)の形式と作成コスト

GTの形式用途コスト
画像レベルラベル分類(この画像は「森林」)最も安い
ピクセルレベルラベルセグメンテーション(各ピクセルを塗り分け)高コスト
バウンディングボックス物体検出(この範囲に「建物」)中コスト
GPS+現地調査検証(実際に行って確認)最も高コスト

プロセスマップ: 自動化可能性マップ

8ステップの各工程で、どこまで自動化でき、どこに人間が必要かを整理しました。

ステップ自動化人間が必要な部分
① 目的定義人間が決める(ビジネス要件)
② センサー選定推薦は可能だが最終判断は人間
③ データ取得
④ 前処理
⑤ 解析(L1/L2)
⑤ 解析(L2.5)精度判断、モデル選定
⑤ 解析(L3)GT収集・設計、モデル選定
⑥ 後処理・検証
⑦ 地図化何と重ねるかは人間が決める
⑧ Web配信テンプレート化すれば完全自動

衛星解析でカバーできる領域の全体マップ

領域見えるもの主なセンサー代表的な手法
植生・農業活性度、作物種類、収穫予測、病害Sentinel-2, MODISNDVI, 時系列, 分類
森林森林減少、バイオマス、樹種、樹冠高Landsat, ALOS-2, GEDI変化検出, SAR, LiDAR
水域・海洋水体範囲、水質、海岸線変動Sentinel-2/3NDWI, 水色アルゴリズム
都市都市拡大、建物、不透水面、夜間光Sentinel-2, VIIRSNDBI, 物体検出
災害洪水、地震被害、土砂崩れ、干ばつSentinel-1/2SAR洪水, 変化検出
大気NO2, CO, メタン, エアロゾルSentinel-5Pカラム量マッピング
地形標高、傾斜、流域、浸水リスクSRTM DEMDEM解析
温度地表面温度、ヒートアイランドLandsat熱赤外LST物理モデル

5日間の実験とプロセスマップの対応

Day体験した範囲プロセスフローの対応箇所
1光学分類 L1→L2→L3②光学 → ④大気補正なし → ⑤L1/L2/L3比較
2SAR洪水検出 L1→L2→L3②SAR → ④スペックルフィルタ → ⑤L1/L2/L3比較
3データパイプライン全体③STAC API → ④大気補正・雲マスク → ⑤L1(指数20種)
4モデル総当たり比較⑤L2(基盤モデル4種) + L3(古典ML6種/DL4種) + 解像度感度

主な発見 ①

転移学習の威力と「逆転現象」

光学(Day 1): 転移学習が圧勝

レベル手法精度
L1: 数式NDVI閾値64.4%
L2: 転移学習ResNet5095.7%(+31pt
L3: FTResNet50 FT98.5%(+2.8pt)

SAR(Day 2): 物理原理が勝つ

レベル手法F1
L1: 閾値後方散乱0.396(1位)
L2: 転移学習ResNet-UNet0.167(3位)
L3: FTU-Net学習0.240(2位)

なぜ逆転するか: ImageNetの知識(テクスチャ・形状)はSAR(後方散乱強度)に転用できない。「水面は鏡→マイクロ波が返ってこない→暗い=洪水」という物理原理のほうが強い。

教訓: AIの民主化は「既知のドメイン」には強力だが、ドメインが変わると汎用モデルは機能しない。

effort vs accuracy
労力 vs 精度のトレードオフ
SAR L1/L2/L3比較
SAR同一画像での検出結果

主な発見 ②

前処理がすべてを決める

大気補正なしだとNDVIが254.7%ズレる

L2A(補正済み)vs TOA(未補正)でNDVIの絶対値が信頼不可に。相関r=0.96でパターンは似るが、閾値分類は崩壊する。
→ どんな高精度モデルを使っても、前処理を間違えると結果が崩壊する。

1画像から20通りの情報が取れる

スペクトル指標20種(植生/水域/都市/災害/土壌/農業/合成)を同一画像から計算。実質5-6次元に圧縮可能。DL時代でも探索・説明性・業界標準として指標は残る。

5種のセンサーで見える世界が全く違う

20スペクトル指数
1画像 → 20の「見方」
TROPOMI 7プロダクト
大気7プロダクト — 地表が全く見えない別世界

主な発見 ③

データ効率とドメインギャップ

たった5枚で+40ポイント

教師データ量精度改善幅
0枚(zero-shot)35.4%
5枚/クラス75.2%+39.8pt
全量(2,700枚)98.0%+22.8pt

「良い表現(事前学習)+ 少量データ」 > 「悪い表現 + 大量データ」。

ドメインギャップの決定的実証

やりたいことで必要な解像度が2桁違う

分類は30mで95%維持(Landsatで十分)。物体検出は0.5mで信頼度半減(ドローンが必要)。

データ効率カーブ
教師データ量 vs 精度
解像度感度
分類精度: 30mでも95%維持

主な発見 ④

何が民主化され、何がされないか(3層構造)

内容具体例
第1層
完全自動化
AIエージェントに丸投げ データ取得(STAC API)、定型前処理(L2A)、スペクトル指標20種、事前学習モデル適用、DEM計算、レポート・可視化
第2層
AI+人間協働
AIが実装、人間が判断 閾値の最適化(-15dBは妥当か?)、モデル適合性判断、前処理パラメータ選択、結果の解釈(F1=0.4は許容か?)
第3層
専門家必須
AIでは代替困難 解析の目的定義、ドメインギャップの事前評価、物理原理に基づく妥当性検証、教師データの品質管理

実感としての結論: 80%の作業は自動化可能。残り20%が結果の質を決定的に左右する。

AIエージェント時代の衛星解析

消滅する専門性

  • コードを書く能力(rasterioでNDVI計算等)
  • データ取得・前処理のルーチンワーク
  • 既存モデルの適用スクリプト

残る / 高まる専門性

  • 問題定義 — 何を見たいか決める力
  • ドメイン適応 — 欧州モデルを東南アジアで使う調整
  • 妥当性保証 — この精度で実務に耐えるか
  • AIオペレーション — 指示→確認→FBの高速サイクル

「コードを書く専門家」→「問いを設計し、妥当性を保証する専門家」への不可逆なシフト

学習方法としての振り返り

「変数を一度に変えない」設計が効いた

Day固定動かす見えたもの
1データ手法手法の効果が明確に見える
2手法の型センサードメインの違いが浮き彫りに
3手法データ側前処理の重要性を実感
4データモデル手法間の差が数値で比較可能

AIエージェント駆動の学習の特性

特性内容
速い37タスクを4日間で実行。従来なら数ヶ月かかる量
網羅的全手法を総当たりで試せる。「とりあえず全部やる」が現実的に
再現可能全コードが残っている。2周目以降で即座に再利用
限界結果の妥当性判断は人間 / 「何が足りないか」のメタ認知はAI不可 / 専門家FBなしの盲点

議論したいこと

素人の初期仮説を1週間で形にしました。間違っているところも多いと思いますが、フィードバックをもらえたらすぐに修正して次にいけます。

#問い
1SARでドメイン知識がDLに勝つ — 実務でもそうですか?
25枚で+40ポイント — 教師データ戦略、実務ではどう考えますか?
33層構造(自動化/協働/専門家) — この整理は実感に合いますか?
4プロセスマップ — 抜けているステップや重要なのに触れていない領域は?
5この学習方法の有効性 — 素人がAIエージェントで学ぶアプローチの可能性と限界は?

STA が解析をするしないとはまったく別の次元として、AIエージェント時代にどう衛星解析がなっていくのかの仮説を持ちたくて探索しています。

Appendix — 全成果物

Day 1〜4の全タスク、全ビジュアルを収録

Day 1: 光学画像分類(EuroSAT / Sentinel-2)

Task 1-1: NDVI閾値分類(L1)

精度64.4%。NDVI/NDWI/NDBIの3指数だけで土地被覆を分類。都市の識別は完全に失敗(0.0%)。スペクトル指数の限界を3D散布図が示す。
L1 3d scatter
L1 3d scatter
L1 classification maps
L1 classification maps
L1 confusion matrix
L1 confusion matrix
eurosat 10classes rgb
eurosat 10classes rgb
index distributions
index distributions

Task 1-2: ResNet50 事前学習モデル推論(L2)

精度95.7%(+31.3pt)。ImageNetで犬猫を学んだモデルが、そのまま衛星画像の10クラス分類に有効。最もコスパの良いステップ。
L2 confusion matrix
L2 confusion matrix
L2 correct vs misclassified
L2 correct vs misclassified
L2 hard samples
L2 hard samples
L2 per class accuracy
L2 per class accuracy
L2 predictions
L2 predictions
L2 tsne
L2 tsne

Task 1-3: ResNet50 ファインチューニング(L3)

精度98.5%(+2.8pt)。衛星画像で再学習し全クラスF1≥0.96。データ50%でも98.0%、きれいなデータでは量より質。
L3 confusion diff
L3 confusion diff
L3 confusion matrix
L3 confusion matrix
L3 correct vs incorrect
L3 correct vs incorrect
L3 detailed learning curves
L3 detailed learning curves
L3 learning curves
L3 learning curves
L3 per class accuracy
L3 per class accuracy

Task 1-4: L1 / L2 / L3 精度比較 ★

最大ジャンプはL1→L2(+31pt)。学習済みモデルをそのまま使うだけで精度が劇的に向上する。effort vs accuracyが核心メッセージ。
3class comparison
3class comparison
L1 L2 L3 comparison bar
L1 L2 L3 comparison bar
L2 vs L3 confusion
L2 vs L3 confusion
L2 vs L3 radar
L2 vs L3 radar
comparison 6samples L1L2L3
comparison 6samples L1L2L3
comparison L3 error rate
comparison L3 error rate
comparison L3 failures
comparison L3 failures
comparison improvement map
comparison improvement map
effort vs accuracy
effort vs accuracy
same image comparison
same image comparison

Task 1-5: 全モデル比較(8手法)

ViT=ResNet50 FT=98.5%(同率首位)。SVM 93.8%で古典MLも健闘。アーキテクチャの差より事前学習の有無が重要。
all models class heatmap
all models class heatmap
all models comparison
all models comparison
all models cost vs accuracy
all models cost vs accuracy
all models ranking
all models ranking

Day 2: SAR洪水検出・時系列・InSAR(Sen1Floods11 / Sentinel-1)

Task 2-1: SAR閾値法による洪水検出(L1)

F1=0.396。水面はマイクロ波を鏡面反射→暗い=洪水。物理原理に基づくシンプルかつ高速な手法。
2-1 raw data explanation
2-1 raw data explanation
2-1 sar 6samples
2-1 sar 6samples
2-1 sar error map
2-1 sar error map
2-1 sar flood detection
2-1 sar flood detection
2-1 sar sample
2-1 sar sample
2-1 sar threshold sensitivity
2-1 sar threshold sensitivity
2-1 speckle filter comparison
2-1 speckle filter comparison

Task 2-2: 事前学習モデル洪水検出(L2)

F1=0.167。ImageNet/Prithviのzero-shotはSARで機能せず。光学と全く異なるドメインギャップ。
2-2 prithvi 6samples
2-2 prithvi 6samples
2-2 prithvi error map
2-2 prithvi error map
2-2 prithvi metrics comparison
2-2 prithvi metrics comparison
2-2 prithvi prob distribution
2-2 prithvi prob distribution
2-2 sar 6samples
2-2 sar 6samples
2-2 sar error map
2-2 sar error map
2-2 sar flood detection
2-2 sar flood detection
2-2 sar prob distribution
2-2 sar prob distribution
2-2 vs L1 comparison
2-2 vs L1 comparison

Task 2-3: U-Net学習による洪水検出(L3)

F1=0.240。CPU 20エポックでは学習不足。DLはSARで大量の学習データと計算リソースが必要。
2-3 confusion matrix
2-3 confusion matrix
2-3 unet learning curves
2-3 unet learning curves
2-3 unet predictions
2-3 unet predictions

Task 2-4: SAR洪水検出 L1/L2/L3 比較 ★

Day1と真逆: L1(0.396) > L3(0.240) > L2(0.167)。SARではドメイン知識がDLに勝つ。
2-4 L1 L2 L3 comparison same images
2-4 L1 L2 L3 comparison same images
2-4 flood comparison bar
2-4 flood comparison bar
2-4 flood comparison radar
2-4 flood comparison radar
2-4 flood comparison table
2-4 flood comparison table
2-4 flood f1 progression
2-4 flood f1 progression
2-4 flood l3 training with baselines
2-4 flood l3 training with baselines
2-4 flood precision recall
2-4 flood precision recall

Task 2-5: NDVI 時系列解析

農地・森林・都市の年間NDVI変動を季節分解。異常検知・フェノロジー分析で「いつ何が起きたか」を把握。
2-5 ndvi anomaly events
2-5 ndvi anomaly events
2-5 ndvi breakpoint detection
2-5 ndvi breakpoint detection
2-5 ndvi monthly heatmap
2-5 ndvi monthly heatmap
2-5 ndvi phenology
2-5 ndvi phenology
2-5 ndvi seasonal decomposition
2-5 ndvi seasonal decomposition
2-5 ndvi timeseries 3types
2-5 ndvi timeseries 3types
2-5 ndvi timeseries overlay
2-5 ndvi timeseries overlay
2-5 ndvi zscore anomaly
2-5 ndvi zscore anomaly

Task 2-6: SAR 変化検出(3手法)

差分法・比率法・CVAの3手法でPre/Post変化を検出。Day1(光学) vs Day2(SAR)の比較まとめ含む。
2-6 change detection 3methods
2-6 change detection 3methods
2-6 change detection comparison bar
2-6 change detection comparison bar
2-6 change detection error maps
2-6 change detection error maps
2-6 change detection inputs
2-6 change detection inputs
2-6 cva visualization
2-6 cva visualization
2-6 day1 vs day2 comparison
2-6 day1 vs day2 comparison

Task 2a-1: Sentinel-1 データ取得と理解

asf_searchによるSentinel-1データ検索。昇交/降交軌道、VV/VH偏波の特性理解。
2a-1 data structure
2a-1 data structure
2a-1 s1 orbit comparison
2a-1 s1 orbit comparison
2a-1 s1 polarization
2a-1 s1 polarization

Task 2a-2: HyP3 RTC前処理パイプライン

HyP3 RTCによるSAR前処理(地形補正等)。クラウド処理で月8,000クレジット無料。
2a-2 preprocessing pipeline
2a-2 preprocessing pipeline

Task 2a-3: JRC恒常水域マスク

JRC Global Surface Waterで恒常水域と洪水域を分離。誤検出の大幅削減。
2a-3 jrc water occurrence
2a-3 jrc water occurrence
2a-3 mask effect
2a-3 mask effect
2a-3 permanent vs flood
2a-3 permanent vs flood

Task 2a-4: 大規模SAR処理

タイリング・並列処理・モザイク化による広域SAR解析パイプライン。
2a-4 gee architecture
2a-4 gee architecture
2a-4 parallel speedup
2a-4 parallel speedup
2a-4 tiling grid
2a-4 tiling grid
2a-4 tiling mosaic
2a-4 tiling mosaic

Task 2b-1: SAR可視化・スペックルフィルタ比較

リニア vs dBスケール、スペックルフィルタ5種の効果比較。ENLによるフィルタ定量評価。
2b-1 filter histogram
2b-1 filter histogram
2b-1 filter params
2b-1 filter params
2b-1 filter params enl
2b-1 filter params enl
2b-1 linear vs db
2b-1 linear vs db
2b-1 speckle filters
2b-1 speckle filters

Task 2b-2: Otsu閾値+水域マスク

Otsu自動閾値・固定dB閾値最適化・恒常水域マスクの3手法を比較。
2b-2 lee filter effect
2b-2 lee filter effect
2b-2 method comparison
2b-2 method comparison
2b-2 otsu threshold
2b-2 otsu threshold
2b-2 threshold sweep
2b-2 threshold sweep
2b-2 water mask effect
2b-2 water mask effect

Task 2b-3: 偏波活用+Pre/Post変化検出

VV/VH偏波情報を活用した変化検出。偏波比(VV/VH)で植生の有無を判別。
2b-3 change comparison
2b-3 change comparison
2b-3 change detection 3methods
2b-3 change detection 3methods
2b-3 polarimetric features
2b-3 polarimetric features

Task 2b-4: SAR 3ch U-Net洪水検出

VV/VH/VV-VHの3チャンネル入力U-Net。431ペアで学習。閾値法との比較。
2b-4 unet predictions
2b-4 unet predictions
2b-4 unet training
2b-4 unet training
2b-4 unet vs threshold
2b-4 unet vs threshold
2b-4 unet vs threshold bar
2b-4 unet vs threshold bar

Task 2b-5: 全手法比較+判断フローチャート ★

5手法を統一条件で比較。災害対応の手法選択フローチャートを作成。精度だけでなくスピード・データ制約も考慮。
2b-5 all methods comparison
2b-5 all methods comparison
2b-5 decision flowchart
2b-5 decision flowchart
2b-5 prediction comparison
2b-5 prediction comparison
2b-5 tradeoff
2b-5 tradeoff
2b-5 validation
2b-5 validation

Task 2c-1: InSAR基礎

SLC・干渉縞・コヒーレンスの概念理解。2016年熊本地震を題材にした合成データ演習。
2c-1 coherence
2c-1 coherence
2c-1 interferogram
2c-1 interferogram
2c-1 slc amplitude phase
2c-1 slc amplitude phase

Task 2c-2: コヒーレンス変化検出

Pre/Postコヒーレンス差分で地震被害域を推定。都市部 vs 農村部の応答差。
2c-2 coherence pre post
2c-2 coherence pre post
2c-2 damage map
2c-2 damage map
2c-2 urban rural comparison
2c-2 urban rural comparison

Task 2c-3: DInSAR変位解析

Unwrapped PhaseからLOS変位マップを生成。mm精度の地殻変動計測。
2c-3 deformation map
2c-3 deformation map
2c-3 displacement histogram
2c-3 displacement histogram
2c-3 displacement profile
2c-3 displacement profile
2c-3 processing flow
2c-3 processing flow
2c-3 wrapped phase
2c-3 wrapped phase

Task 2c-4: SAR統合まとめ ★

Day2全13タスクを「3つの物理量」(後方散乱/コヒーレンス/位相差)で統合整理。決定木とパイプライン図。
2c-4 full pipeline
2c-4 full pipeline
2c-4 sar decision tree
2c-4 sar decision tree
2c-4 sar summary table
2c-4 sar summary table
2c-4 sensor guide
2c-4 sensor guide

Day 3: データパイプライン(Sentinel-2/3, Landsat, SRTM, TROPOMI)

Task 3-0: 前処理比較実験 ★

大気補正なしだとNDVIが254.7%過小評価。L2A vs TOAで絶対値が信頼不可に。前処理は必須。
3-0 cloud mask visualization
3-0 cloud mask visualization
3-0 l2a vs l1c scatter
3-0 l2a vs l1c scatter
3-0 preprocessing comparison ndvi
3-0 preprocessing comparison ndvi
3-0 preprocessing stats table
3-0 preprocessing stats table

Task 3-1: スペクトル指数20種 ★

同じ画像からバンドの組み合わせを変えるだけで20通りの情報。実質5-6次元に圧縮可能。NDVI↔SAVI r≥0.98。
3-1 indices correlation heatmap
3-1 indices correlation heatmap
3-1 spectral indices 20 grid
3-1 spectral indices 20 grid
3-1 spectral indices stats
3-1 spectral indices stats

Task 3-2: バンド合成・雲マスク・大気補正

10種バンド合成、3種雲マスク、L1C vs L2A大気補正を比較。SCL+膨張処理が実務的。
3-2 atmospheric correction comparison
3-2 atmospheric correction comparison
3-2 atmospheric correction stats
3-2 atmospheric correction stats
3-2 band composites 10
3-2 band composites 10
3-2 cloud mask comparison
3-2 cloud mask comparison
3-2 s2cloudless diff map
3-2 s2cloudless diff map
3-2 s2cloudless mask comparison
3-2 s2cloudless mask comparison
3-2 s2cloudless stats table
3-2 s2cloudless stats table

Task 3-3: DEM地形解析

SRTM 30mで関東平野の地形解析。傾斜・流域・浸水シミュレーション。水位5m→浸水10km2。
3-3 dem 3d hillshade
3-3 dem 3d hillshade
3-3 flood inundation graph
3-3 flood inundation graph
3-3 terrain analysis 5maps
3-3 terrain analysis 5maps

Task 3-4: 水域・海洋解析

琵琶湖の水域抽出、東京湾のクロロフィル推定、Sentinel-3 SSTとの比較。
3-4 sentinel2 vs sentinel3 comparison
3-4 sentinel2 vs sentinel3 comparison
3-4 sentinel3 ocean color
3-4 sentinel3 ocean color
3-4 sentinel3 sst
3-4 sentinel3 sst
3-4 water biwa extraction
3-4 water biwa extraction
3-4 water chiba coastline
3-4 water chiba coastline
3-4 water sst anomaly proxy
3-4 water sst anomaly proxy
3-4 water tokyo bay chlorophyll
3-4 water tokyo bay chlorophyll

Task 3-5: 大気解析(TROPOMI 7プロダクト)★

地表ではなく大気の化学成分を測定する全く別種の衛星データ。1ピクセル=面積50万倍の別世界。
3-5 atmosphere no2 comparison
3-5 atmosphere no2 comparison
3-5 atmosphere s5p overview
3-5 atmosphere s5p overview
3-5 atmosphere summary
3-5 atmosphere summary
3-5 tropomi aerosol
3-5 tropomi aerosol
3-5 tropomi all products
3-5 tropomi all products
3-5 tropomi ch4
3-5 tropomi ch4
3-5 tropomi co
3-5 tropomi co
3-5 tropomi hcho
3-5 tropomi hcho
3-5 tropomi no2
3-5 tropomi no2
3-5 tropomi o3
3-5 tropomi o3
3-5 tropomi so2
3-5 tropomi so2

Task 3-6: 熱赤外・ヒートアイランド

Landsat TIRSで東京の地表面温度を実測。LST 26-65度C。「温度」と「過去」はLandsat一択。
3-6 thermal lst 4panel
3-6 thermal lst 4panel
3-6 thermal lst highres
3-6 thermal lst highres
3-6 thermal uhi analysis
3-6 thermal uhi analysis

Task 3-7: 都市拡大検出+夜間光

2018→2024年のNDBI差分で東京の都市拡大を検出。6年で新規市街化108.5km2(+8.3%)。
3-7 urban analysis summary
3-7 urban analysis summary
3-7 urban expansion detection
3-7 urban expansion detection
3-7 urban greenery analysis
3-7 urban greenery analysis
3-7 urban nightlight proxy
3-7 urban nightlight proxy

Day 4: AI手法の深掘り比較(基盤モデル・セグメンテーション・DL)

Task 4-1: ファウンデーションモデル4種比較 ★

ResNet50(96.0%) ≈ DINOv2(95.0%) >> CLIP(34.6%)。ラベルなしDINOv2が教師あり学習と同精度。
4-1 foundation models comparison
4-1 foundation models comparison
4-1 foundation models samples
4-1 foundation models samples
4-1 foundation models summary
4-1 foundation models summary

Task 4-2: セグメンテーション3手法比較

U-Net mIoU 28.1%。分類(96%) vs セグメンテーション(28%)の難易度差。教師データ作成コストが桁違い。
4-2 segmentation comparison
4-2 segmentation comparison
4-2 segmentation metrics
4-2 segmentation metrics
4-2 segmentation summary
4-2 segmentation summary

Task 4-3: 物体検出

DeepForest 55本検出 vs YOLOv8 0検出。ドメインギャップの決定的実証。COCOの「建物」と上空の「屋根」は別物。
4-3 object detection comparison
4-3 object detection comparison
4-3 object detection deepforest
4-3 object detection deepforest
4-3 object detection summary
4-3 object detection summary

Task 4-4: 古典的機械学習6手法比較

RF/SVM/k-NN等をPCA圧縮特徴で比較。SVM 93.8%で手作り特徴でも健闘。
4-4 classical ml comparison
4-4 classical ml comparison
4-4 classical ml summary
4-4 classical ml summary

Task 4-4b: ピクセルレベル古典ML

ピクセル単位の13バンドスペクトル分類。実務に近いアプローチ。
task4 4b classification maps
task4 4b classification maps
task4 4b summary
task4 4b summary

Task 4-5: 深層学習4手法比較

ViT 99.1% > ViT+LoRA 98.8% > ResNet FT 98.4% > Scratch 91.3%。事前学習の効果が支配的。
4-5 deep learning comparison
4-5 deep learning comparison
4-5 deep learning confusion
4-5 deep learning confusion
4-5 deep learning summary
4-5 deep learning summary

Task 4-6: データ効率性実験 ★

教師データ0枚(35%)→5枚(75%,+40pt!)→100枚(94%)→全量(98%)。最初の数枚が最もコスパ高い。
4-6 data efficiency comparison
4-6 data efficiency comparison
4-6 data efficiency curve
4-6 data efficiency curve
4-6 data efficiency summary
4-6 data efficiency summary

Task 4-7: 解像度感度実験 ★

分類は30mで95%維持、物体検出は0.5mで信頼度半減。やりたいことで必要な解像度が2桁違う。
4-7 classification resolution
4-7 classification resolution
4-7 detection resolution
4-7 detection resolution
4-7 resolution summary
4-7 resolution summary

Generated by AGI Cockpit / Claude Code — 2026-03-15

Space Tech Accelerator / 松島 宏佑