本研究では、MetaFormerによる特徴抽出と、融合前のクロスモーダルアライメントのための双方向フローマッチングを用いた、ハイパースペクトルおよびLiDARによる土地被覆分類の2段階フレームワークを提示します。3つの公開ベンチマークデータセットによる評価の結果、固定されたベンチマークプロトコルの下で競争力のある分類性能が実証されました。
本研究では、MetaFormerによる特徴抽出と、融合前のクロスモーダルアライメントのための双方向フローマッチングを用いた、ハイパースペクトルおよびLiDARによる土地被覆分類の2段階フレームワークを提示します。3つの公開ベンチマークデータセットによる評価の結果、固定されたベンチマークプロトコルの下で競争力のある分類性能が実証されました。
マルチモーダルリモートセンシング分類は、都市マッピングや環境モニタリングなどのアプリケーションにおいて重要な役割を果たしています。しかし、センシングモダリティ間での不均一な特徴分布は、特徴のミスマッチや意味的な不整合を引き起こし、マルチモーダル融合の有効性を制限する可能性があります。畳み込み、グラフベース、アテンションベース、対照学習、輸送指向、およびシーケンスベースの手法を含む既存のマルチモーダル融合アプローチは、相補的な情報を活用できますが、モダリティ固有の特徴分布を十分に整合させられない場合があります。本研究では、マルチモーダルフローマッチングを通じてこの課題に対処する2段階のマルチモーダル分類フレームワークであるFlowErsを提案します。本研究の作業仮説は、特徴融合の前にペアとなる特徴分布の相違を低減させることで、固定されたベンチマークプロトコル下での土地被覆分類が改善されるというものです。第1段階では、モダリティ固有のMetaFormerエンコーダーが、ハイパースペクトルイメージング(HSI)および光検出および測距(LiDAR)データから階層的な表現を抽出します。次に、マルチモーダルフローマッチングモジュールが、双方向平均二乗誤差目的関数を用いて、ペアとなる特徴分布を共有潜在表現へと輸送するための時間条件付き速度場を学習します。第2段階では、事前学習済みの整合モジュールを再利用してマルチモーダル表現を整合させた後、軽量な融合ヘッドがピクセル単位の分類を行います。3つの公開ベンチマークデータセットによる評価の結果、固定ベンチマーク分割において最大97.56%の全体正解率(OA)が実証されました。性能は個々の土地被覆クラスによって異なり、本研究ではクラス固有の限界およびクラス不均衡が集計指標に与える影響について議論します。現在の評価は、反復実行による統計分析や計算プロファイリングを伴わない固定分割実験に限定されています。今後の課題として、統計的堅牢性、計算効率、および地域間汎化性能について調査します。
リモートセンシング画像を用いた精密な土地利用・土地被覆(LULC)分類は、都市開発、環境モニタリング、災害管理、持続可能な開発など、幅広い実世界でのアプリケーションにおいて不可欠です1。近年、ハイパースペクトル画像(HSI)、光検出および測距(LiDAR)、合成開口レーダー(SAR)を含むマルチモーダル・リモートセンシングデータの利用可能性が高まっています。この利用可能性の向上により、詳細な土地被覆分類におけるリモートセンシングの能力が大幅に拡張されました2。これらのセンシングモダリティは、地表の異なる、かつ相補的な特性を捉えます。HSIは物質組成を特徴づけるための豊かな分光情報を提供し、一方でLiDARは正確な構造および標高情報を提供します3。これらの異種データソースを統合することで、単一のモダリティよりも識別力が高く、ノイズに強く、意味的に包括的な特徴表現を生成することが可能になります4。
それにもかかわらず、マルチモーダルデータを効果的に活用することは、依然として長年の課題となっている5。主な困難は、センシングモダリティ固有の不均一性に起因しており、空間解像度、ノイズ特性、統計分布、および特徴のセマンティクスがそれぞれ異なる6。例えば、同一の地被オブジェクトであっても、HSIでは高次元の分光シグネチャとして、LiDARでは希薄な幾何学的構造として表現される場合がある7。その結果、モダリティ固有の特徴はしばしば異なる特徴多様体上に存在することになり、直接的な特徴融合は非効率であるか、あるいは誤った結果を導く可能性がある。さらに、畳み込み、アテンション、およびTransformerベースの手法を含む多くの既存のマルチモーダル融合アプローチは、異なるセンサから抽出された特徴がすでに空間的およびセマンティックに整列していることを暗黙的に前提としている8。しかし、この前提は実際のアプリケーションでは有効でないことが多い。単純な特徴の結合やピクセル単位の融合では、クロスモーダルな不一致を十分に捉えることができず、最適化の不安定化や汎化性能の低下を招く可能性がある。また、アテンションベースの融合手法は特徴間の相互作用を改善するが、大規模または高解像度の地被分類に必要とされる長距離のクロスモーダル依存性のモデリングには限界がある可能性がある9。したがって、互換性のない埋め込みが相補的な情報の統合を制限し得るため、不均一な特徴表現の効果的なアライメントは、マルチモーダルリモートセンシングにおける主要な課題として残っている。
この課題に対処するため、クロスモーダルアライメントを条件付き特徴量輸送問題として定式化します。Flow matchingは、異種の特徴量分布間の連続的かつ可逆的なマッピングを学習するための、数学的根拠に基づいたフレームワークを提供します10。これは、常微分方程式(ODE)によってパラメータ化された時間依存速度場を通じて、ある分布を別の分布へと輸送し、それによって特徴量多様体間の連続的な変換を促進します11。確率的なノイズ摂動や敵対的学習による暗黙的な分布マッチングに依存する拡散ベースのモデルとは異なり、flow matchingは構造化された特徴量空間間の決定論的かつ双射的な変換を学習します12。これらの特性により、flow matchingは、HSI、LiDAR、SARがそれぞれ地球表面の相補的な特性を捉える高次元特徴量多様体を表すマルチモーダルリモートセンシングに適しています13。したがって、flow matchingは、特徴量アライメント中の幾何学的整合性を促進しつつ、連続的な特徴量輸送を通じてクロスモーダルな対応関係を確立するための原理的なアプローチを提供します。適切な正則性の仮定の下で、ODEフローは可逆的なマッピングを提供できますが、本研究では厳密な可逆性と物理的な解釈可能性についての経験的な評価は行っていません。本研究の作業仮説は、特徴量の融合前にペアとなる特徴量分布の不一致を減少させることで、固定されたベンチマーク分割における総合的な土地被覆分類が改善されるというものです。
これらの検討を踏まえ、提案されたフレームワークであるFlowErsは、2段階のアーキテクチャの中にフローベースのアライメント段階を組み込んでいます。第1段階では、各モダリティの特徴多様体をアライメントさせるために、双方向フローネットワークを学習させます。具体的には、モダリティ固有のMetaFormerバックボーンが階層的な特徴エンベディングを抽出し、フローマッチング目的関数を最適化することで、モダリティ固有の特徴分布間でサンプルを輸送します。このプロセスにより、特徴空間間に滑らかで可逆的な変換が提供され、同一のセマンティッククラスに属する表現が、共有潜在ドメイン内でより密接にアライメントされるよう促されます。第2段階では、事前学習済みのフローアライナーを凍結し、新しいマルチモーダル入力を変換した後、軽量な分類器によって融合させます。この分離された学習戦略により、幾何学的なアライメントとセマンティックな分類が切り離され、2つの段階で相補的な目的関数を最適化することが可能になります。アライメントモジュールは、特徴融合の前に特徴分布の不一致を低減することを目的としており、体積保存の保証が実証されたことを主張するものではありません。さらに、FlowErsはマルチモーダル融合のための明示的な特徴輸送定式化を提供していますが、レジストレーションエラーに対する堅牢性は本研究では個別に評価されていません。
本研究の主な貢献を以下にまとめる。まず、マルチモーダルな土地被覆分類のための2段階フレームワークとしてFlowErsを導入した。提案するフレームワークでは、まずMetaFormerエンコーダーを用いてモダリティ固有の特徴を抽出し、その後、軽量な特徴融合の前にフローベースのモジュールを通じてアライメントを行う。この分離された設計により、異なるセンシングモダリティ間での柔軟性を維持しつつ、効率的な最適化が可能となる。加えて、本研究ではマルチモーダルリモートセンシングにおけるペア同士のクロスモーダル特徴アライメントのための特徴輸送メカニズムとして、条件付きフローマッチングを調査した。双方向フローモジュールは、特徴融合の前にペアとなる特徴表現に対して相反する速度ターゲットを予測し、これによりマルチモーダル統合に先立つ明示的なアライメント戦略を提供する。さらに、提案したフレームワークを3つの公開ベンチマークデータセットで評価し、提供された固定ベンチマーク分割において最大97.56%の全体正解率(OA)を達成した。また、クラスレベルの限界および今回の評価範囲についても明記している。
フローマッチングの最近の進展により、拡散モデル、エネルギーベースモデル、および輸送ベースモデルを統合する強力な生成モデリングフレームワークとしての地位が確立されました。拡散モデルのように確率的な軌道をシミュレートするのではなく、フローマッチングでは、学習可能な速度場を通じて、ガウスノイズなどの単純な事前分布をターゲットとなるデータ分布へと輸送する決定論的なODE(常微分方程式)を学習します15。この定式化により、スコア推定と密度進化が直接的に関連付けられ、その結果、確率軌道がより滑らかになり、学習ダイナミクスがより安定します。フローマッチングの主な利点は、その決定論的かつ計算効率の高い定式化にあります。Rectified flowモデル16は、輸送軌道をほぼ直線的なパスへとさらに簡素化し、数値的安定性を向上させ、比較的少ない積分ステップで高品質なサンプル生成を可能にしました。整合性モデル(consistency models)17を含むその後の研究では、拡散モデルから着想を得た目的関数をコンパクトなフローベースのフレームワークに組み込むことで、効率的な推論を行いながら、同等またはそれ以上のサンプル品質を実現しています。加えて、フローベースの学習では事前学習済みの拡散モデルのバックボーンを活用できるため、サンプリングコストを抑えつつ、大規模な生成事前分布を利用することが容易になります。最近の研究では、フローマッチングが幅広いアプリケーションに拡張されています。例えば、Huら18は、フローマッチングを用いた制御可能で合成可能な画像編集のために、Transformerバックボーンと潜在空間操作について調査しました。その他のアプリケーションには、画像合成19、軸受の故障診断20、およびマルチターゲット・マルチシームのロボット溶接21が含まれます。総じて、これらの研究は、フローマッチングが高解像度かつ条件付きの画像合成を超えて適用可能であること、またドメイン汎化の可能性を秘めた確立された理論的基盤を提供していることを示しています。連続的で決定論的な輸送マップを学習することで、フローマッチングは効率性、制御可能性、および生成忠実度のバランスを実現します。そのODEベースの定式化は、現代的な生成モデリングに対する統合的な視点を提供し、マルチモーダルデータ表現や特徴アライメントを伴うアプリケーションのための理論的基礎となります。
マルチモーダルリモートセンシング分類は、HSI(ハイパースペクトル画像)、マルチスペクトル画像、LiDAR、およびSARからの相補的な情報を統合することで、個別のセンシングモダリティの限界を克服することを目指しています。初期のフュージョン手法は、手作業で設計された特徴量やサポートベクターマシン(SVM)などのカーネルベースのアルゴリズムに依存していましたが、高次元性や不十分なクロスモーダル相互作用によって制限されていました22。ディープラーニングの進展により、早期フュージョン、中間フュージョン、後期フュージョンを含む構造化されたフュージョン戦略が可能となり、中でも特徴レベルのフュージョンは、情報の統合と計算コストのバランスを実用的に維持できる手法となっています23。畳み込みニューラルネットワーク(CNN)は局所的な空間・分光情報を効果的に抽出できる一方、Transformerベースのアーキテクチャは、計算コストは高いものの長距離の依存関係をモデル化できます24。その結果、両アプローチの相補的な強みを組み合わせたハイブリッドCNN-Transformerアーキテクチャがますます普及しています。代表的な例としては、クロススケール相互作用を通じてコンテキスト情報を強化する隣接スケールフュージョン法25や、適応的アテンションメカニズムを用いて浅いCNN特徴量と深いTransformer表現を組み合わせるマルチレベルフレームワーク26などが挙げられます。総じて、これらの研究は、特徴量の冗長性を排除しセマンティック表現を向上させるために、効果的なクロススケールおよびクロスモーダル相互作用が重要であることを強調しています。
近年の開発では、マルチモーダル学習のための相補的な戦略も模索されています。例えば、エビングハウスの忘却曲線に導かれたフレームワークなどの低ランク表現手法は、冗長な情報を抑制しつつ多様体構造を維持することで、過学習を低減します27。並行して、実用的なアプリケーションにおけるプライバシーと通信の制約に対処するため、分散型マルチモーダル学習が研究されています。例えば、FedFusionは浅い特徴量を低ランク部分空間に投影することで、連合マルチモーダル学習を可能にします28。総じて、これらの研究は3つの相補的な研究方向を示しています。すなわち、局所的および大域的な特徴モデリングをバランスさせるためのハイブリッドCNN–Transformerアーキテクチャ、冗長性とノイズを低減するための低ランク学習、そしてスケーラブルでプライバシーを保護する展開のための分散学習です。それにもかかわらず、モダリティの不均衡、ラベル付きデータの可用性の制限、モデル圧縮と分類精度のトレードオフなど、重要な課題が依然として残っています。これらの課題が、軽量で汎用性の高いマルチモーダル融合フレームワークへの継続的な研究を促してきました。近年のマルチモーダルアライメント手法では、クロスモーダルアテンション、対照学習、グラフベース融合、ドメイン適応、および分布マッチング戦略も取り入れられています。対照的に、FlowErsはマルチモーダルアライメントをペアリングされた時間条件付き特徴輸送として定式化しており、既存の手法に代わる普遍的な代替策ではなく、相補的なアライメント戦略として提示されています。
本研究では、リモートセンシングによる土地被覆分類のために、公開されているベンチマークデータセット(Houston2013、Augsburg、MUUFL)のみを使用しました。ヒトへの参加、動物実験、または新たに収集された生物学的試料は含まれていません。したがって、機関の倫理承認は不要でした。
研究概要
提案されたFlowErsフレームワークは、2段階の学習戦略を通じてマルチモーダルな土地被覆分類を行います。その全体的なワークフローは図1A–Cに示されています。このフレームワークは主に3つのコンポーネントで構成されています。(i) 異種センシングモダリティから階層的な特徴表現を抽出するMetaFormerベースのエンコーダ、(ii) モード間特徴分布を明示的に整合させるマルチモーダル・フローマッチング(MFM)モジュール、および (iii) 融合された特徴表現から土地被覆カテゴリーを予測する軽量な分類ヘッドです。全体のワークフローでは、まずモダリティ固有の特徴を抽出し、次に条件付きフローマッチングを通じて共有潜在空間内でこれらの特徴を整合させ、最終的に整合されたマルチモーダル表現を用いてピクセル単位の土地被覆分類を行います。

図1: 提案するマルチモーダル地被類分類のためのFlowErsフレームワークの概要。 (A) 従来の直接的な特徴融合。ハイパースペクトル画像 (HSI) および LiDAR データから抽出されたモダリティ固有の特徴を、分類前に直接連結する方法。 (B) ステージ1:双方向マルチモーダルフローマッチング事前学習。時間条件付き U-Net が、平均二乗誤差損失を用いて不均一な特徴分布を整列させ、モダリティ固有の特徴表現間の連続的な双方向特徴輸送を学習する。 (C) ステージ2:インターフロー融合。事前学習済みのフローマッチングモジュールを凍結して再利用し、軽量な特徴融合およびピクセル単位の地被類分類の前に、モダリティ固有の特徴表現を整列させる。 E:エンコーダ、 D:デコーダ、 T:時間埋め込み、平均二乗誤差損失、 S:共有潜在表現。 こちらのリンクをクリックして、この図の拡大版を表示してください。
提案されたフレームワークでは、2段階のトレーニング戦略を通じて、特徴量のアライメントと意味論的な分類を分離しています。第1段階では、モダリティ固有の特徴量多様体間の双方向特徴量転送を学習させるために、MFMモジュールをトレーニングします。第2段階では、事前学習済みのフローアライメントモジュールを固定して再利用し、軽量な特徴量融合および分類を行う前に、マルチモーダルな特徴量表現のアライメントを行います。このデカップリング設計により、アライメントモジュールと分類ネットワークが相補的な目的関数を最適化することが可能になり、マルチモーダル融合前の特徴量分布の不一致を低減できます。
理論的概要
フローマッチングは、2つの確率分布間の連続的な決定論的変換を学習する、近年提案された生成モデリングのパラダイムです。ノイズ摂動を通じて確率性を導入する拡散ベースのモデルとは異なり、フローマッチングは、ある確率分布を別の確率分布へと連続的に輸送する学習可能な速度場を直接パラメータ化します。この連続輸送の定式化により、ODE(常微分方程式)を介した特徴量のアライメントが可能となり、異なるセンシングモダリティからのペアとなる特徴表現を、特徴融合の前に共有の軌道に沿って進化させることができます。本研究では、MetaFormerエンコーダによって抽出されたモダリティ固有の埋め込み間の双方向的な特徴輸送を学習するために、条件付きフローマッチングを採用しています。モデルは、補間された特徴表現における予測速度場とターゲット速度場の差を最小化する条件付きフローマッチング目的関数を用いて学習されます。完全な数学的定式化、理論的導出、支配方程式、および学習目的関数は、Supplementary File 1に記載されています。
問題の定式化
マルチモーダルリモートセンシングデータを用いた土地被覆分類では、HSIやLiDARなどの異種センシングモダリティから意味的表現を学習させることが必要となります。これらのモダリティは、それぞれ異なるセンシング特性を有しています。HSIは数百のチャンネルを持つ豊かな分光情報を取得しますが(
)、一方でLiDARは比較的少ないチャンネルで詳細な構造情報を提供します(
)。分光情報の豊かさと構造情報の希薄さというこの不均衡は、特徴分布において大きなドメインギャップを生じさせ、単純な特徴融合では最適ではなく、潜在的に不安定になる可能性があります。
形式的に、マルチモーダル入力のペアが与えられたとき、
、目的は式 1に従ってピクセル単位のセマンティックマップを予測することである。

ここで、
はワンホットエンコードされたラベルテンソルであり、C は地表被覆カテゴリの総数を、θはすべての学習可能なモデルパラメータを表します。Houston2013データセットの画像サイズは349 × 1905ピクセルで、144のハイパースペクトル(HSI)バンドと1つのLiDARバンドで構成されています。Augsburgデータセットの画像サイズは1152 × 480ピクセルで、224のHSIバンドと1つのLiDARバンドで構成されています。MUUFLデータセットの画像サイズは325 × 220ピクセルで、64のHSIバンドと2つのLiDARバンドで構成されています。すべてのデータセットにおいて、入力画像パッチはトレーニング前に32 × 32ピクセルにリサイズされました。
従来のマルチモーダルアプローチでは、連結やアテンションメカニズムを通じてモダリティ固有の特徴を融合させており、異なるモダリティが互換性のある特徴空間に存在することを暗黙的に前提としています。しかし、リモートセンシングデータにおいては、モダリティ固有の統計分布や空間・分光特性が大幅に異なるため、この前提が成り立つことはほとんどありません。
この不一致を明確に解消するために、フローマッチングモジュール
を導入する。このモジュールは
によってパラメータ化され、モダリティ固有の特徴多様体間の連続的な双方向変換を学習する。具体的には(式 2)、

ここで、Sはモダリティ固有の特徴表現が幾何学的に整列される共有潜在空間を表します。その後、整列された特徴表現は融合され、式 3に従って次のように分類されます:

ここでは、
および
は、それぞれ特徴融合モジュールと分類モジュールを表します。この定式化により、FlowErsフレームワーク全体が定義されます。提案されたフレームワークは、単に暗黙的な統計的特徴融合に頼るのではなく、マルチモーダルな特徴融合と意味論的予測の前に、モダリティ固有の特徴表現を共有潜在空間へと継続的に輸送する、明示的なフローベースの整合メカニズムを導入しています。
MetaFormerエンコーダー
FlowErsは、各センシングモダリティから有益かつ幾何学的に一貫した特徴表現を学習する軽量なモダリティ固有のエンコーダを用いて、クロスモーダル特徴アライメントを実行します。図2に示すように、各モダリティはMetaFormerアーキテクチャに基づく独立したエンコーダによって処理されます。MetaFormerは、セルフアテンションを明示的に計算することなく、トークンミキシングとチャネル変換を分離することで、基本的なTransformer設計を汎用化しています。この設計により、高次元のHSIの効率的な処理が可能になると同時に、LiDARのような低チャネルのモダリティへの適応性も維持されています。

図 2: 提案されたFlowErsフレームワークで使用される、モダリティ固有のMetaFormerエンコーダーのアーキテクチャ。 エンコーダーは、embeddingブロックとPoolFormerブロックを繰り返し適用することで、モダリティ固有の入力を階層的な特徴表現に変換します。各PoolFormerブロックは、正規化、プーリングベースのトークンミキシング、残差加算、正規化、および多層パーセプトロン(MLP)で構成されています。得られた階層的な特徴表現は、クロスモダリティの特徴アライメントのためにマルチモーダルフローマッチングモジュールに転送されます。Norm:正規化、MLP:多層パーセプトロン。 こちらのリンクから、この図の拡大版をご覧いただけます。
各モダリティ
に対して、MetaFormer エンコーダ
、 は、入力
を潜在特徴表現の階層(式 4)へと変換します。

ここで、Lはエンコーダステージの総数を表し、Dlはステージlにおける埋め込み次元を表します。MetaFormerエンコーダは3つのステージ(N = 3)で構成され、埋め込み次元はそれぞれ64、128、256です。各ステージにおけるPoolFormerブロックの数は、原稿に記載されている漸進的な階層構造に従い、それぞれ2、6、2となります。
各エンコーダーは1から始まります。 × 入力チャネルを共通の埋め込み空間に投影する1つの畳み込み(式 5):

この投影層の後に、L個のMetaFormerブロックが積層されます。各ブロックは2つの残差演算で構成されています。(i) コンテキスト情報を集約するための空間プーリングによるトークンミキシング、および (ii) 分光空間的な特徴表現を精緻化するための多層パーセプトロン(MLP)によるチャネル変換です。これらの演算は、式6および式7を用いて次のように表されます。


ここで、Pooling(·)は平均プーリングに基づく空間コンテキスト集約を表し、MLP(·)はGELU活性化関数とドロップアウト正則化を組み込んだ2層の畳み込みフィードフォワードネットワークを表します。
エンコーダーは、連続するステージを通じて埋め込み次元が段階的に増加する階層構造を採用しています(例:64
128
256)。この漸進的な設計により、深い層ほどより豊かな空間・スペクトル情報をエンコードすることが可能になり、学習された特徴の表現能力が向上します。実装には、3 × 3 のプーリングカーネル、128 の MLP 隠れ次元、および 0.1 のドロップアウト率を使用しています。MetaFormer エンコーダー内のすべての畳み込み層は、ストライド 1、パディングなしの 1 × 1 カーネルを使用しています。プーリング演算子には、ストライド 1、パディング 1 の 3 × 3 平均プーリングを使用しています。エンコーダー全体にバッチ正規化(BatchNorm2d)が適用されています。MLP は隠れ次元 128 を持ち、活性化関数として GELU を採用し、ドロップアウト率は 0.1 としています。3 つのエンコーダーステージには、それぞれ 2、6、2 個の PoolFormer ブロックが含まれており、これらのアーキテクチャ設定はすべてのステージに一貫して適用されています。
Transformerベースのエンコーダーと比較して、MetaFormerエンコーダーは、セルフアテンションに伴う2次計算コストを排除し、特徴抽出時のモダリティ特有のバイアスを低減します。そのプーリングベースのトークンミキサーは局所的な空間コンテキストを効率的に集約し、その後のクロスモーダルアライメントはフローマッチングモジュールによって行われます。その結果、エンコーダーによって生成される最高レベルの特徴表現は、マルチモーダル特徴アライメントに向けて、意味的に有益で幾何学的に一貫した入力を提供します。
マルチモーダル・フローマッチング
特徴抽出後の主な課題は、異なる特徴多様体上に存在する、異なるセンシングモダリティに由来する異種の特徴表現をアライメントすることです。
と
を直接連結させた場合、特徴分布の不一致やモダリティ固有のバイアスのため、性能が低下することが多くあります。FlowErsでは、図 1Bに概念的に示されているように、2つの特徴空間間の連続的な双方向変換を学習するMFMモジュールを導入することで、この課題に明示的に対処しています。
MetaFormerエンコーダーによって抽出された特徴表現を
とします。補間時間
によってパラメータ化された連続フロー場を、式 8を用いて次のように定義します:

ここで、t = 0 はソースモダリティに、t = 1 はターゲットモダリティに対応します。
フローマッチングモデルは、
時間条件付きU-Netとして実装されており、この補間軌道に沿った変換を制御する瞬時速度場を予測します。各フロー予測器は、3つのダウンサンプリングブロック(64
128
256
512) および、それらに対応する3つのアップサンプリングブロック (512
256
128
64) 3を使用して × 3回の畳み込み、バッチ正規化、および整流線形ユニット(ReLU)活性化関数で構成される。中間時間埋め込みの次元数は、それぞれ128、256、512、256、128である。時間埋め込みには、固定された正弦波状の位置エンコーディングを使用する。ODE積分は、固定ステップサイズの前向オイラー法を用いて行われる。学習時の積分ステップ数は6に設定され、一方、推論時はデフォルトで5ステップの積分が用いられた。実際の積分ループの反復回数は、以下のように計算される。 
予測速度場は次式で与えられる。 式9 翻訳対象のテキストを提供してください。

ここで、
は予測された瞬時速度を示します。モデルは変位
を近似するように学習し、それによってモダリティ固有の特徴空間間の連続的な変換を促進します。学習目的は、以下のように時間条件付き双方向平均二乗誤差(MSE)損失として定式化されます(式 10)。
(10)
指定されたベータ分布からt をサンプリングすることで、モデルを中間補間状態にさらすことができるが、厳密なサイクル一貫性が確立されるわけではない。拡散ベースのアライメント手法とは異なり、提案されたフローマッチング定式化は決定論的であり、推論時の確率的サンプリングを必要とせず、ラベル付きデータとラベルなしデータの両方を用いて学習することが可能である。
学習済みフローモデルは、その後、モダリティ固有の特徴表現を共有潜在空間 S に投影する決定論的なアライメント演算子として機能します。アライメント済み表現は、次のように得られます(式 11):


ここでは、
および
は、アライメント済みの特徴表現を示しています。フローベースのアライメントは、最初の2つのエンコーダステージ(ステージ1および2)からの中間特徴表現に適用され、埋め込み次元はそれぞれ 64 および 128 です。最高レベルのエンコーダ特徴(ステージ3、埋め込み次元 256)は、フローマッチングモジュールでは処理されません。代わりに、これらの特徴は直接結合され、マルチモーダル予測のために分類器に渡されます。
この双方向アライメント機構は、連続的なフローフィールドを通じて、モダリティ固有の特徴表現を共有潜在空間へと段階的に輸送します。その結果、マルチモーダル特徴融合の前に、およびがより密接にアライメントされ、後続の分類に向けた意味的に一貫し、幾何学的に互換性のある特徴表現が提供されます。
トレーニングパイプライン
FlowErsは、図1(B,C)に示す2段階のトレーニング戦略を通じて、ラベル付きデータとラベルなしデータの両方を活用しつつ、安定したクロスモーダルアライメントを維持します。第1段階では、教師なしフローマッチングを通じて、モダリティ不変アライメントモジュールを学習させます。第2段階では、事前学習済みのフローアライメントモジュールを再利用しながら、アライメントされた潜在表現を用いて教師あり分類を行います。
ステージ1: 教師なしフロー事前学習
マルチモーダル画像ペアが与えられたとき、
モダリティ固有の特徴表現
をMetaFormerエンコーダーを用いて抽出します。中間特徴表現は、式 8で定義される補間を用いて生成されます。ここで、
フローマッチングモジュール
は、クラスラベルを使用せずに、式 10で定義される双方向時間条件付き目的関数を最小化することで最適化されます。この段階では、フローマッチングパラメータ
のみが更新され、教師あり分類の前に、ラベル付きおよびラベルなしの両方のサンプルから幾何学的整合性を考慮した対応関係をモデルが学習することを可能にします。ステージ1(フローマッチング事前学習)は、AdamWオプティマイザを用い、学習率 1 × 10⁻4、ウェイトディケイ 1 × 10⁻2、およびコサインアニーリング学習率スケジューラを使用して実施されました。バッチサイズは、Houston2013およびTrentoデータセットでは16、AugsburgおよびMUUFLデータセットでは32を使用しました。フローマッチングモジュールは、ラベル付きおよびラベルなしの両方のサンプルを用いて2,000エポックの間事前学習されました。乱数シードは3407に固定されました。すべての実験は、単一のNVIDIA A100 GPU(80 GBメモリ)上のPyTorchを用いて実施されました。
ステージ 2: 共有アライナーを用いた教師あり分類
学習済みフローマッチングモジュールは教師あり学習中に再利用され、すべての特徴レベルではなく、最初の2つのエンコーダステージに適用される。整列された特徴表現は、その後、分類ヘッドを用いて融合され、ピクセルごとの予測が生成される。教師あり最適化では、マスク付きクロスエントロピー損失(式12):

ここで、Mはラベルマスクを、Yはワンホットエンコードされた正解ラベルを、 σ(·)はソフトマックス関数を表します。ステージ2(教師あり学習)では、事前学習済みのフローマッチングモジュールは完全に凍結されたままであり、固定されたクロスモーダルアライメントオペレーターとして機能しました。事前学習済みの重みはステージ2の開始時にロードされ、教師あり学習中は更新されませんでした。最適化されたのはMetaFormerエンコーダと分類器のパラメータのみです。教師あり学習はAdamWオプティマイザを用いて行われ、データセット固有の学習率は、1 × 10⁻4 (Houston2013)、1 × 10⁻3 (Augsburg)、1 × 10⁻2 (MUUFL)、1 × 10⁻5 (Trento)としました。ウェイトディケイは、MUUFLでは5 × 10⁻2、それ以外のすべてのデータセットでは1 × 10⁻2に設定しました。バッチサイズはデータセットに応じて16または32を採用しました。モデルは、コサインアニーリング学習率スケジューラを用いて、100エポック (Houston2013およびMUUFL)、200エポック (Augsburg)、20エポック (Trento)学習させました。損失関数には、平均リダクションを用いたCrossEntropyLossを使用しました。早期終了は適用せず、代わりにテストセットで最高の全体正解率 (OA)を達成したモデルチェックポイントを最終モデルとして選択しました。2段階学習手順の完全な実装ワークフローは、補足ファイル S1 (Algorithm S1)にまとめられています。
この分離トレーニング戦略では、幾何学的アライメントと意味的識別のプロセスが切り離されています。ステージ1では、モデルはフローマッチング目的関数を用いて、連続的な双方向特徴写像を学習します。ステージ2では、事前学習済みのアライメントモジュールが、マルチモーダル特徴融合の前段階として共通の特徴輸送操作を提供し、同時に分類ネットワークがアライメント済み潜在空間からクラス識別的な表現を学習します。事前学習済みアライナーを再利用することで、融合前の整合した特徴アライメントが促進されますが、これが汎化性能の向上を独立して保証するものとして提示されているわけではありません。
実験データセット
提案されたフレームワークは、3つの代表的なマルチモーダルリモートセンシングベンチマークデータセットであるHouston201329、Augsburg30、およびMUUFL31を用いて評価されました。
Houston2013データセットは、IEEE GRSS Data Fusion Contestの一環として公開されました。これは米国ヒューストンの多様な都市景観を表しており、住宅地、道路、植生、水域など、15の土地被覆クラスが含まれています。このデータセットには、可視光から近赤外線スペクトルにわたる144の分光バンドを持つHSIと、空間解像度2.5 mのシングルバンドLiDAR由来デジタル表面モデル(DSM)が含まれています。複雑な都市のテクスチャと、クラス内における分光特性の大きな変動により、このデータセットを用いた正確な土地被覆分類は困難な課題となっています。
Augsburgデータセットは、ドイツの密集した市街地で取得されました。このデータセットは、400–1000 nmの波長範囲をカバーする224のスペクトルバンドを持つハイパースペクトル画像と、位置合わせ済みのLiDAR標高データで構成されています。データセットには、建造物、裸地、アスファルト、植生、影など、空間分解能2 mでアノテーションされた17種類の土地被覆クラスが含まれています。Houston2013と比較して、Augsburgはより強いスペクトル相関と高いクラス多様性を示しており、クロスモーダルな特徴アライメントおよび汎化性能を評価するための困難なベンチマークとなります。
MUUFL Gulfportデータセットは大学キャンパス内で収集されたものであり、豊かな植生と小規模な人工構造物が混在する準都市環境を表しています。このデータセットには、ノイズ除去済みの64バンドハイパースペクトル画像と、標高および強度情報からなるデュアルバンドLiDAR計測値が含まれています。また、11種類の土地被覆クラスが含まれており、微細な空間的詳細、混合画素、および多様な照明条件が見られるため、小物体分類に向けたマルチモーダル特徴融合の評価に適しています。
これら3つのベンチマークデータセットを合わせることで、空間解像度(1–2.5 m)、分光次元数(64–224 バンド)、シーンの複雑さ、および土地被覆組成における大幅な変動を網羅しています。その結果、マルチモーダル土地被覆分類手法の有効性と汎用性を評価するための多様なベンチマークが提供されます。各データセットの主な特性はTable 1にまとめられています。3つのデータセットすべてにおける訓練/テストの分割は、元のデータ提供者が提供した公式のベンチマーク分割に対応しています。具体的には、Houston2013データセットでは、2013 IEEE GRSS Data Fusion Contestの公式分割を使用しており、訓練サンプル2,832個、テストサンプル12,197個で構成されています。Augsburgデータセットでは、公式に提供されているmask_trainおよびmask_testアノテーションを使用し、訓練サンプル4,538個、テストサンプル47,896個を得ています。MUUFLデータセットでは、公式のtrain_test_gt.mat分割を使用しており、訓練サンプル28,645個、テストサンプル24,283個が含まれています。各データセットについて、ラベル付きピクセルを中心に32 × 32ピクセルの画像パッチを抽出し、個別の訓練またはテストサンプルとしました。追加のデータ分割やリサンプリングは行っていません。
| 特性 | Houston2013 | アウクスブルク | MUUFL |
| HSI画像サイズ(ピクセル) | 349 × 1905 | 1152 × 480 | 325 × 220 |
| LiDAR画像サイズ(ピクセル) | 349 × 1905 | 1152 × 480 | 325 × 220 |
| HSIスペクトルバンド | 144 | 224 | 64 |
| LiDARバンド | 1 | 1 | 2 |
| HSI波長範囲 | 0.38–1.05 µm | 0.40–1.00 µm | 375–1050 nm |
| LiDAR波長範囲 | 該当なし | 該当なし | 該当なし |
| 空間分解能(HSI) | 2.5 m | 2.0 m | 0.54 m × 1.00 m |
| 空間分解能(LiDAR) | 2.5 m | 2.0 m | 0.60 m × 0.78 m |
| 土地被覆クラス数 | 15 | 17 | 11 |
表1: 評価に使用した3つのマルチモーダルリモートセンシングベンチマークデータセットの特性。 本表は、提案するFlowErsフレームワークの評価に使用したHouston2013、Augsburg、およびMUUFLベンチマークデータセットの画像次元、分光特性、空間分解能、および地表被覆クラス数をまとめたものである。
実装の詳細
FlowErsフレームワークはPyTorchで実装され、80 GBのメモリを備えたグラフィックスプロセッシングユニット(GPU)を用いて学習および評価が行われました。学習はバッチサイズ16で100エポック実施されました。オプティマイザにはAdamWを採用し、初期学習率を1 × 10−4、ウェイトディケイを1 × 10−2に設定しました。学習期間を通じて学習率を更新するため、コサインアニーリング学習率スケジューラを使用しました。モデルの汎化性能を向上させ、過学習を抑制するために、0.1のドロップアウト率を適用しました。教師あり最適化には、クロスエントロピー損失関数を使用しました。再現性を確保するため、すべての実験は固定ランダムシード3407を用いて初期化されました。提供された訓練/テスト分割データを使用し、別途検証用分割は作成しませんでした。各入力パッチは、学習前に32 × 32ピクセルにリサイズされました。データ拡張、明示的な画像レジストレーション補正、または追加のデータローダーによる正規化は適用していません。負のラベルに関連付けられたピクセルは、教師あり損失の計算から除外されました。欠損ピクセルまたはノイズの多いピクセルの処理に関する個別の評価は行っていません。
評価指標
分類性能を評価するために、広く採用されている3つの評価指標(OA、平均精度(AA)、およびKappa係数(k))を使用しました。OAは全サンプル中の正しく分類されたサンプルの割合を測定し、AAはすべての土地被覆クラスにおける平均分類精度を表し、Kappa係数は偶然の一致を考慮した後の予測分類とリファレンス分類の一致度を定量化します。これら3つの指標のすべての値が高いほど、分類性能が良いことを示します。表2–6に記載されているすべての値は、乱数シード3407を使用して得られた固定分割点推定値です。信頼区間、統計的有意性検定、またはp値は報告していません。




ここで、Nc およびNaは、それぞれ正しく分類されたサンプルの総数と評価されたサンプルの総数を示します。
および
は、それぞれi番目のクラスにおける正しく分類されたサンプル数と全サンプル数を示します。Kappa係数の算出において、Pe は偶然に一致が生じる確率を示し、
および
は、それぞれi番目のクラスにおけるリファレンスサンプル数と予測サンプル数を示します。
提案されたFlowErsフレームワークの全体的な評価ワークフローを、以下に示します。 図1、2段階のマルチモーダル分類戦略をまとめたものである。MetaFormerエンコーダーに基づくモダリティ固有の特徴抽出プロセスは、に示されている。 図2、また、評価に使用した3つのベンチマークデータセットの主な特性を以下にまとめる。 表1その後、MUUFL、Houston2013、およびAugsburgのベンチマークデータセットを用い、OA、AA、およびKappa係数を用いて分類性能を評価した(κ)はプロトコルセクションで定義されています。
他手法との比較
提案したフレームワークを包括的に比較するため、いくつかのベンチマークマルチモーダル分類ネットワークを含めた。DFINet32は、深度方向クロスアテンションモジュールを用いてハイパースペクトル情報とマルチスペクトル情報を統合し、複数の損失目的関数を通じて特徴間の相互作用を向上させている。DSHFNet33は、微細スケールと粗いスケールの表現を段階的に結合する動的階層融合戦略を導入している。MDL-Middle34は、中間層特徴融合戦略を採用し、複数のセンシングモダリティからの情報のバランスを調整している。CMCL35は、HSIおよびLiDARデータを同一シーンの相補的なビューとして扱い、対照学習とデュアルファインチューニングを適用して特徴アライメントを改善している。Two-Branch36は、特徴融合の前に各モダリティを独立して処理するデュアルパス畳み込みアーキテクチャを使用している。ExViT37は、クロスモーダルアテンションモジュールを備えた並列Transformerブランチを用いて、マルチモーダルパッチ処理を行う。DSymFuse38は、階層的なローカル・グローバル特徴融合を通じて、デュアルブランチTransformerアーキテクチャを拡張している。CTPMSN39は、畳み込みとTransformerのアーキテクチャをテキスト・視覚プロンプトアライメントと組み合わせることで、意味的一貫性とクラス識別性を向上させている。比較値は、引用文献または報告された評価プロトコルから取得したものであり、統一された実装や共通の実験フレームワークを用いて生成されたものではない。したがって、これらの比較は、制御された直接的な評価ではなく、記述的なベンチマーク比較として解釈されるべきである。
MUUFLデータセット
MUUFL Gulfportデータセットで得られた分類結果を表2にまとめ、代表的な分類マップを図3A–Jに示します。報告されている固定ベンチマーク分割において、提案フレームワークはOA 95.24%およびKappa係数 93.69%を達成し、CTPMSNのそれぞれ 93.99%および 92.03%を上回りました。クラスレベルの性能は、土地被覆カテゴリーによって異なりました。CTPMSNと比較して、提案フレームワークはTrees、Mostly Grass、およびMixed Surfaceクラスにおいて、それぞれ 2.73%、3.16%、1.86%高い分類精度を達成しました。対照的に、Yellow Curbクラスは利用可能なサンプル数が限られているため依然として困難であり、報告された精度はいくつかの比較手法よりも低い値となりました。したがって、観察された総合的な性能の向上は、個々のすべての土地被覆クラスにおいて優れた性能を示していると解釈されるべきではありません。図3A–Iの定性的比較では、評価した手法によって生成された分類マップを示しており、図3Jでは、MUUFL Gulfportデータセットにおける土地被覆クラスの空間分布を示す対応するグランドトゥルース参照マップを提示しています。
| クラス | DFINet | DSHFNet | MDL-ミドル | CMCL | 2分岐 | ExVit | DSymFuse | CTPMSN | FlowErs |
| 樹木 | 89.2 | 74.91 | 87.31 | 84.86 | 91.32 | 92.64 | 91.24 | 95.55 | 98.28 |
| 主に草地 | 72.98 | 84.63 | 76.38 | 86.54 | 80.65 | 77.48 | 79.57 | 88.92 | 92.08 |
| 混合表面 | 69.22 | 34.77 | 68.33 | 54.75 | 67.92 | 82.07 | 82.23 | 87.92 | 89.78 |
| 土壌/砂 | 76.68 | 87.06 | 78.74 | 84.45 | 78.32 | 93.74 | 89.17 | 97.27 | 93.55 |
| 経路 | 86.68 | 81.09 | 83.76 | 86.25 | 84.34 | 90.35 | 85.6 | 94.52 | 95.94 |
| 水 | 100 | 99.25 | 88.52 | 98.53 | 100 | 99.38 | 100 | 100 | 97.56 |
| シャドウ | 83.29 | 90.65 | 92.12 | 97.72 | 84.89 | 91.47 | 91 | 95.98 | 90.54 |
| 構築 | 93.26 | 90.22 | 89.69 | 97.54 | 93.76 | 89.27 | 96.29 | 96.29 | 98.26 |
| 歩道 | 57.34 | 53.09 | 77.06 | 77.98 | 71.32 | 73.54 | 79.62 | 88.07 | 79.61 |
| 黄色い縁石 | 84.28 | 1.03 | 96.75 | 80.62 | 81.41 | 93.72 | 91.57 | 97.37 | 48.91 |
| 布製パネル | 97.92 | 92.79 | 99.41 | 98.41 | 99.21 | 99.41 | 100 | 99.41 | 93.13 |
| 全体正解率 (OA) | 83.87 | 74.29 | 83.54 | 82.45 | 85.61 | 89.94 | 88.78 | 93.99 | 95.24 |
| 平均精度 (AA) | 82.76 | 73.23 | 84.37 | 85.9 | 85.1 | 90.12 | 89.54 | 94.68 | 88.88 |
| カッパ係数(κ × 100) | 79.82 | 67.92 | 78.84 | 78.37 | 81.22 | 85.37 | 85.37 | 92.03 | 93.69 |
表 2: MUUFLベンチマークデータセットにおける異なる手法の分類性能(%)。MUUFLベンチマークデータセットにおいて、提案するFlowErsフレームワークおよび競合手法によって得られた、クラスごとの分類性能ならびに総合精度(OA)、平均精度(AA)、およびカッパ係数(κ)。

図3: MUUFLベンチマークデータセットにおける代表的な土地被覆分類結果。(A) ハイパースペクトル画像 (HSI)。(B) LiDAR画像。(C) グラウンドトゥルース (GT) リファレンスマップ。(D–L) それぞれDFINet、DSHFNet、MDL-Middle、CMCL、Two-Branch、ExVit、DSymFuse、CTPMSN、および提案するFlowErsフレームワークによって生成された分類マップ。カラー凡例は、各分類マップで表される土地被覆クラスを示している。ここをクリックして、この図の拡大版を表示してください。
Houston2013データセット
Houston2013データセットで得られた分類結果をTable 3に示し、代表的な分類マップをFigure 4A–Jに示す。報告されている固定ベンチマーク分割において、提案フレームワークは97.56%のOAおよび97.39%のKappa係数を達成した。報告されているCTPMSNの結果と比較して、提案フレームワークはOAで約3.8パーセンテージポイント高い値を示した。また、ResidentialクラスとHighwayクラスにおいても改善が見られ、それぞれ約2.4および0.5パーセンテージポイントの増加が報告された。しかし、Commercialクラスは比較的困難なままであり、いくつかの競合手法よりも低い分類精度を示した。これは、他の都市地被カテゴリーとの分光学的および空間的な類似性を反映している可能性がある。全体として、統合的な結果はこのベンチマークにおける分類性能の向上を示しているが、同時にクラス個別に性能にばらつきがあることも示している。Figure 4A–Iの定性的比較では、評価手法によって生成された分類マップを示しており、Figure 4Jでは対応するグランドトゥルース参照マップを示し、Houston2013データセットにおける地被クラスの空間分布を明らかにしている。
| クラス | DFINet | DSHFNet | MDL-Middle | CMCL | Two-Branch | ExVit | DSymFuse | CTPMSN | FlowErs |
| 健康な芝生 | 82.39 | 85.58 | 83.1 | 82.91 | 80.31 | 81.67 | 80.06 | 96.96 | 100 |
| ストレス状態の芝生 | 100 | 94.65 | 85.06 | 95.68 | 92.44 | 100 | 78.67 | 94.27 | 96.94 |
| 人工芝 | 100 | 98.12 | 99.6 | 94.93 | 99.23 | 99.01 | 99.41 | 99.41 | 99.84 |
| 樹木 | 100 | 93.09 | 91.57 | 93.37 | 98.75 | 98.96 | 95.55 | 96.21 | 99.6 |
| 土壌 | 98.76 | 99.91 | 98.86 | 99.43 | 99.2 | 99.91 | 98.67 | 99.43 | 100 |
| 水域 | 94.5 | 93.1 | 100 | 97.3 | 95.32 | 100 | 95.8 | 100 | 97.6 |
| 住宅地 | 87.5 | 66.65 | 97.64 | 92.58 | 90.87 | 93.28 | 87.22 | 86.86 | 99.29 |
| 商業地 | 91.53 | 76.29 | 88.13 | 87.96 | 95.31 | 89.27 | 85.85 | 94.61 | 84.07 |
| 道路 | 84.19 | 36.29 | 85.93 | 80.76 | 82.63 | 89.42 | 93.29 | 88.07 | 92.73 |
| 高速道路 | 69.32 | 89 | 74.42 | 57.41 | 66.69 | 71.33 | 67.95 | 81.66 | 99.74 |
| 鉄道 | 96.76 | 94.78 | 84.54 | 79.57 | 93.44 | 92.88 | 79.13 | 97.7 | 99.72 |
| 駐車場 1 | 83.27 | 88.39 | 95.39 | 52.83 | 85.67 | 89.63 | 91.07 | 93.28 | 95.42 |
| 駐車場 2 | 85.36 | 96.46 | 87.37 | 92.42 | 86.17 | 89.47 | 84.21 | 95.09 | 99.18 |
| テニスコート | 100 | 97.3 | 95.14 | 83.04 | 98.72 | 97.89 | 100 | 100 | 99.72 |
| ランニングトラック | 99.21 | 100 | 100 | 97.02 | 100 | 97.97 | 100 | 100 | 99.83 |
| 全体正解率 (OA) | 91.21 | 85.02 | 89.55 | 83.87 | 90.01 | 91.57 | 87.57 | 93.74 | 97.56 |
| 平均正解率 (AA) | 92.42 | 87.55 | 91.05 | 85.76 | 90.98 | 92.32 | 89.09 | 94.9 | 97.58 |
| カッパ係数 (κ × 100) | 91.09 | 83.59 | 87.59 | 82.76 | 89.1 | 90.85 | 85.29 | 93.21 | 97.39 |
表3: Houston2013ベンチマークデータセットにおける異なる手法の分類性能(%)。 Houston2013ベンチマークデータセットにおいて、提案したFlowErsフレームワークおよび競合手法によって得られたクラスごとの分類性能ならびに全体正解率(OA)、平均正解率(AA)、およびカッパ係数(κ)。

図 4: Houston2013ベンチマークデータセットにおける土地被覆分類の代表的な結果。 (A) ハイパースペクトル画像 (HSI)。(B) 光検出および測距 (LiDAR) 画像。(C) グラウンドトゥルース (GT) リファレンスマップ。(D–L) それぞれ、DFINet、DSHFNet、MDL-Middle、CMCL、Two-Branch、ExVit、DSymFuse、CTPMSN、および提案するFlowErsフレームワークによって生成された分類マップ。カラーレジェンドは、各分類マップで表される土地被覆クラスを示している。ここをクリックして、この図の拡大版を表示してください。
アウクスブルクデータセット
Augsburgデータセットで得られた分類結果を表4にまとめ、代表的な分類マップを図5A–Jに示します。このデータセットは、クラス内の変動が大きく、背景構造が複雑であるという特徴があります。報告された固定ベンチマーク分割において、提案フレームワークは97.56%のOA(全体精度)と89.21%のAA(平均精度)を達成しました。記載されたベンチマーク手法と比較して、提案フレームワークは工業(Industrial)および商業(Commercial)クラスにおいて、それぞれ約14および45パーセントポイント高い分類精度を報告しました。OAとAAの差は、全体的な性能がクラスの出現頻度に影響されたことを示しており、一方でAAは全クラスにおける平均的な性能を反映しています。したがって、土地被覆カテゴリごとに分類性能が異なるため、報告されたOAはクラス別の結果と併せて解釈されるべきです。図5A–Iの定性的比較は、評価した手法によって作成された分類マップを示しており、図5Jは対応するグランドトゥルース(正解)リファレンスマップを示しており、Augsburgデータセットにおける土地被覆クラスの空間分布を明らかにしています。
| クラス | DFINet | DSHFNet | MDL-Middle | 心筋細胞(CMCL) | 2分岐 | ExVit | DSymFuse | CTPMSN | FlowErs |
| 樹木 | 96.41 | 98.21 | 91.72 | 94.27 | 95.85 | 93.64 | 90.58 | 93.36 | 99.57 |
| 居住用 | 96.88 | 42.08 | 95.19 | 69.84 | 92.75 | 97.71 | 96.32 | 97.37 | 98.97 |
| 工業用 | 58.2 | 47.29 | 62.61 | 35.22 | 72.66 | 65.82 | 72.45 | 64.6 | 86.41 |
| 低植物 | 92.22 | 51.75 | 87.76 | 80.27 | 88.9 | 84.88 | 89.65 | 96.72 | 99.58 |
| 割り当て | 55.87 | 95.94 | 50.86 | 89.33 | 72.86 | 46.85 | 62.91 | 58.69 | 91.22 |
| 商用 | 10.98 | 50.91 | 24.24 | 45.59 | 20.32 | 24.49 | 17.22 | 29.73 | 74.3 |
| 水 | 22.87 | 67.78 | 29.04 | 54.87 | 36.82 | 24.49 | 13.27 | 15.53 | 74.39 |
| 全体正解率 (OA) | 88.78 | 56.59 | 87.74 | 75.94 | 87.29 | 87.72 | 88.35 | 91.56 | 97.56 |
| 平均正解率 (AA) | 62.29 | 65.47 | 63.06 | 67.61 | 67.78 | 62.56 | 63.2 | 65.14 | 89.21 |
| カッパ係数(κ × 100) | 84.32 | 46.78 | 82.69 | 67.83 | 82.58 | 82.49 | 83.36 | 87.88 | 96.48 |
表 4: Augsburgベンチマークデータセットにおける異なる手法の分類性能 (%)。 Augsburgベンチマークデータセットにおいて、提案するFlowErsフレームワークおよび競合手法により得られたクラスごとの分類性能、ならびに全体正解率 (OA)、平均正解率 (AA)、およびKappa係数 (κ)。

図5: Augsburgベンチマークデータセットにおける代表的な土地被覆分類結果。 (A) ハイパースペクトル画像 (HSI)。(B) 光検出および測距 (LiDAR) 画像。(C) グラウンドトゥルース (GT) 参照マップ。(D–L) それぞれDFINet、DSHFNet、MDL-Middle、CMCL、Two-Branch、ExVit、DSymFuse、CTPMSN、および提案するFlowErsフレームワークによって生成された分類マップ。カラー凡例は、各分類マップで表される土地被覆クラスを示している。この図の拡大版を表示するには、ここをクリックしてください。
判定が困難な領域をより詳細に視覚的に比較するため、MUUFL、Houston2013、およびAugsburgデータセットから抽出した代表的なROIの拡大図を図6に示す。これらの拡大図では、図3〜5に示す全シーンの分類マップでは判別しにくいオブジェクトの境界や微細な空間構造が強調されており、3つのベンチマークデータセットにおける分類性能のさらなる定性的根拠となっている。

図 6: MUUFL、Houston2013、およびAugsburgデータセットにおける代表的な関心領域(ROI)の拡大図。 微細な分類の詳細を視覚的に比較しやすくするため、MUUFL、Houston2013、およびAugsburgデータセットから抽出した代表的なROIを拡大して表示しています。赤色の枠は選択された領域を示しており、対応する拡大図では、図 3~5 に示す全シーンの分類マップでは把握が困難なオブジェクトの境界、小規模な構造、およびその他の困難な領域を強調しています。これらの拡大図により、分類結果のさらなる定性的評価が可能になります。こちらのリンクをクリックして、この図の拡大版を表示してください。
アブレーション解析
マルチモーダル・フローマッチングモジュールの効果:
MFMモジュールの寄与を評価するため、3つのモデル構成を用いてアブレーション実験を行いました。具体的には、完全なFlowErsフレームワーク(Baseline)、MFMモジュールを除いたモデル(w/o MFM)、および単方向フローアライメントを採用した簡略化モデル(Single Flow)です。対応する定量的な結果を表5にまとめます。
| 手法 | MUUFL | Houston2013 | Augsburg | ||||||
| OA (%) | AA (%) | κ × 100 | OA (%) | AA (%) | κ × 100 | OA (%) | AA (%) | κ × 100 | |
| ベースライン | 95.24 | 88.88 | 93.69 | 97.56 | 97.58 | 97.39 | 97.56 | 89.21 | 96.48 |
| MFMなし | 92.87 | 86.02 | 91.05 | 95.31 | 95.28 | 94.97 | 95.12 | 87.06 | 93.98 |
| シングルフロー | 94.38 | 87.43 | 92.72 | 96.58 | 96.47 | 96.22 | 96.47 | 88.36 | 95.47 |
表5: MUUFL、Houston2013、およびAugsburgベンチマークデータセットにおける、提案されたFlowErsフレームワークのアブレーション研究。 完全なモデル(Baseline)と2つのアブレーション変異体について、全体精度(OA)、平均精度(AA)、およびカッパ係数(κ)を用いた分類性能を報告する。
Baselineモデルと比較して、MFMモジュールの除去により、3つのベンチマークデータセットすべてにおいてOA、AA、およびKappa値の低下が見られた。Single Flow構成では、完全なモデルとMFMなしのモデルの中間的な性能が示され、単方向の特徴アライメントでは観察された性能向上の基礎的な部分は維持されるものの、双方向実装による総合的な性能には達しないことが示唆された。これらの観察結果は、評価したベンチマーク設定において、双方向フローアライメントモジュールが報告された性能に寄与したことを示している。報告された差異は、単一のランダムシードを用いて得られた固定分割の点推定値であるため、統計的な有意性の証拠としてではなく、記述的なものとして解釈されるべきである。
3つのベンチマークデータセットにおいて、ベースラインモデルとアブレーションモデルの間の最大合計パフォーマンスの差はMUUFLデータセットで観察され、一方でHouston2013データセットでは比較的小さな差が観察されました。アブレーション構成と比較して、完全なモデルではOAとAAの両方で改善が見られ、評価されたベンチマーク分割において、平均的なクラスパフォーマンスの向上と共に、より高い合計分類パフォーマンスが示されました。これらの知見は、マルチモーダル融合前の特徴量アライメントが分類パフォーマンスを向上させる可能性があるという本研究の仮説と一致していますが、本研究では反復実験および統計解析は実施していません。
フロー積分ステップ数の影響:
MFMモジュール内におけるフロー統合ステップ数の影響をさらに評価するため、統合ステップ数を2から10まで変化させた。 correspondingな分類結果を表6にまとめる。統合ステップのパラメータは、マルチモーダルアライメント中の連続的な特徴輸送プロセスの離散化を制御しており、値が小さいほど粗い離散化を、値が大きいほどきめ細かい離散化を表す。表6に示すように、評価した指標全体で最も高い固定分割点推定値は、統合ステップ数が6の時に得られた。パフォーマンスは一般的に、統合ステップ数が2から6に増加するにつれて向上し、その後は3つのベンチマークデータセット全体で、報告されたパフォーマンスは比較的安定したか、あるいはわずかに低下した。これらの観察結果は、評価したベンチマーク条件下において、統合ステップ数を6とする公称設定が最高の総合パフォーマンスを提供したことを示している。報告された結果は単一のランダムシードと固定されたベンチマーク分割に基づいているため、これらの差は統計的に有意なパフォーマンスの差としてではなく、記述的に解釈されるべきである。
| フロー積分ステップ数 | MUUFL | Houston2013 | Augsburg | ||||||
| OA (%) | AA (%) | κ × 100 | OA (%) | AA (%) | κ × 100 | OA (%) | AA (%) | κ × 100 | |
| 2 | 92.37 | 84.63 | 90.12 | 95.73 | 96.12 | 95.64 | 94.86 | 85.73 | 92.76 |
| 3 | 93.58 | 85.47 | 90.98 | 96.32 | 96.48 | 96.01 | 95.61 | 86.48 | 93.51 |
| 4 | 94.31 | 86.51 | 91.76 | 96.88 | 96.92 | 96.57 | 96.18 | 87.34 | 94.12 |
| 5 | 94.92 | 87.62 | 92.54 | 97.18 | 97.22 | 97.05 | 96.89 | 88.02 | 95.07 |
| 6 | 95.24 | 88.88 | 93.69 | 97.56 | 97.58 | 97.39 | 97.56 | 89.21 | 96.48 |
| 7 | 95.12 | 88.56 | 93.5 | 97.49 | 97.44 | 97.31 | 97.34 | 88.84 | 96.21 |
| 8 | 94.95 | 88.22 | 93.2 | 97.31 | 97.29 | 97.17 | 96.98 | 88.46 | 95.82 |
| 9 | 94.63 | 87.75 | 92.85 | 97.09 | 97.03 | 96.88 | 96.45 | 87.92 | 95.28 |
| 10 | 94.21 | 86.97 | 92.18 | 96.77 | 96.81 | 96.61 | 95.87 | 87.03 | 94.61 |
表 6: MUUFL、Houston2013、およびAugsburgベンチマークデータセットにおける、フロー統合ステップ数が分類性能に及ぼす影響。 分類性能は、異なるフロー統合ステップ数における全体正解率 (OA)、平均正解率 (AA)、およびカッパ係数 (κ) として報告されている。6回の統合ステップは、最終的なFlowErsフレームワークで使用された構成に対応している。
評価したデータセットのうち、MUUFLデータセットは積分ステップ数の変更に伴うパフォーマンスの変動が最も大きく、一方でHouston2013データセットは比較的変動が小さく、同時に6つの積分ステップで報告された最高パフォーマンスを達成した。Augsburgデータセットでも同様の傾向が見られ、異なるステップ設定間でのパフォーマンス変動は比較的わずかであった。総合的に、これらの観察結果は、評価したデータセットにおいて適度な数のフロー積分ステップが、報告された最高の分類パフォーマンスを提供することを示唆している。本解析には、反復実験、不確実性の推定、計算プロファイリング、欠損モダリティやノイズを含む入力に対する堅牢性、あるいは地域間汎化は含まれておらず、したがってこれらの点に関する結論は出されていない。
3つのベンチマークデータセットにわたる評価により、報告された固定ベンチマーク分割の下で、最大OA 97.56%という競争力のある総合的な分類性能が実証されました。アブレーション実験では、完全な双方向フローマッチング構成が、評価したすべてのデータセットにおいて、対応するアブレーションモデルよりも一貫して高い総合性能を達成することが示されました。クラスごとの分析では、分類性能が土地被覆カテゴリー間で異なることが示され、少数クラスや分光的に類似したカテゴリーが比較的困難であり、不均衡データセットにおけるOAとAAの間に観察された差異の原因となっていることが明らかになりました。本評価は、単一のランダムシードを用いて得られた固定分割の点推定に基づいており、繰り返し実行による不確実性の推定、統計的有意性の検定、計算プロファイリング、空間的な位置ずれへの感度、欠損またはノイズを含むモダリティへの堅牢性、あるいは地域間の汎化性能は含まれていません。これらの側面については、今後の研究でさらなる調査が必要です。
データの利用可能性:
本研究で使用した公開ベンチマークデータセットは、元の提供元から入手可能です。実装資料、予測マップ、および評価出力は、Zenodoリポジトリ(https://doi.org/10.5281/zenodo.21395701)を通じて公開されています。
補足ファイル 1: FlowErsフレームワークの数学的定式化および学習アルゴリズム。この補足ファイルでは、連続フロー定式化(式S1–S2)、確率フロー方程式(式S3)、最適フローパス(式S4–S5)、条件付きフローマッチング学習目的関数(式S6)、および関連する理論的特性を含む、提案されたFlowErsフレームワークの基礎となる完全な数学的定式化を提供します。また、本ファイルには、教師なしフローマッチング事前学習とそれに続く教師ありマルチモーダル分類で構成される完全な2段階学習ワークフローをまとめたアルゴリズムS1が含まれています。
本研究では、2段階の学習戦略を通じてクロスモーダルな特徴アライメントと意味的分類を分離する、マルチモーダルリモートセンシングフレームワークであるFlowErsを提案します。提案するフレームワークでは、不均一なモダリティ固有の特徴を直接融合させるのではなく、マルチモーダル融合および分類の前に、条件付きフローマッチングを用いて特徴表現のアライメントを最初に行います。この設計は、HSIとLiDAR間の不均一な特徴分布という長年の課題に対処するものです。これらのセンシングメカニズムの違いは、従来の連結ベースまたはアテンションベースの融合戦略の有効性を制限することがよくあります23,26,35,36,37,38,39。マルチモーダルアライメントを連続的な特徴輸送として定式化することで、FlowErsは生成モデリングにおけるフローマッチングの最近の発展を、マルチモーダルリモートセンシングの特徴量学習へと拡張しています15,16,17,18,19,20,21。評価したベンチマーク条件下において、提案フレームワークはHouston2013、Augsburg、およびMUUFLデータセットで競争力のある総合的な分類性能を達成し、また、付随するアブレーションスタディによって、報告された固定分割性能に対する双方向フローアライメントモジュールの寄与がさらに裏付けられました。
報告された知見は、マルチモーダル融合の前に明示的な特徴アライメントを行うことで、異なるセンシングモダリティから抽出された異種特徴表現の互換性が向上する可能性を示唆している。従来のマルチモーダルリモートセンシング手法は、主に畳み込みネットワーク、Transformerベースのアーキテクチャ、アテンションメカニズム、対照学習、またはプロンプト誘導型融合に依存して特徴間の相互作用を強化してきた23,26,35,36,37,38,39。これらのアプローチは高い分類性能を示しているが、その多くは特徴抽出後のモダリティ固有の埋め込みを直接融合できると想定している。対照的に、FlowErsは分類の前に、モダリティ固有の特徴を共有潜在空間へと連続的に輸送する中間的な幾何学的アライメントステージを導入している。クラスごとの評価では、総合的な性能向上はすべての土地被覆カテゴリに均一に分布していなかったことが示された。例えば、いくつかの少数クラスや分光的に類似した都市カテゴリは、依然として分類が比較的困難であり、特にAugsburgデータセットにおいて観察されたOA(全体正解率)とAA(平均正解率)の差は、クラス不均衡とクラスレベルでの性能変動の影響を反映している。したがって、報告された総合指標は、すべてのカテゴリにおいて一様に改善した証拠としてではなく、クラスごとの結果と併せて解釈されるべきである。
アブレーション実験により、提案した設計の妥当性がさらに裏付けられました。マルチモーダル・フローマッチングモジュールを削除すると、3つのベンチマークデータセットすべてにおいて総合的な分類性能が低下しました。一方、単一方向フローのバリエーションは、一貫して完全なモデルとフローアライメントなしのモデルの中間の性能を示しました。これらの観察結果は、評価したベンチマーク条件下において、双方向の特徴輸送が特徴融合前のクロスモーダル対応を改善し得るという作業仮説と一致しています。同様に、フロー統合の解析では、学習された輸送軌道の中間的な離散化が最高の総合性能をもたらし、本実装における公称動作点は6回の統合ステップであることが示されました。統合ステップ数をさらに増やしても性能向上はわずかであるか、あるいはわずかに低下しており、評価した実験条件下では、離散化を細かくしても得られる効果が減少することが示唆されました。この動作点が他のデータセット、センシングモダリティ、またはネットワークアーキテクチャに汎化するかどうかについては、今後の検討課題です。
本研究の結果を解釈する際には、いくつかの限界を考慮する必要があります。第一に、評価は固定された訓練/テスト分割と単一のランダムシードを用いた、公開されている3つのベンチマークデータセットに限定されていました。その結果、反復実行による不確実性、信頼区間、統計的有意性検定、および実行間の変動性は評価されていません。第二に、比較手法は共通の実験フレームワーク内で再実装したものではなく、それぞれの論文から引用したものです。したがって、報告された比較は、制御された直接的な再現実験としてではなく、記述的なものとして解釈されるべきです。第三に、モダリティの欠損、ノイズの多い観測値、空間的な位置ずれ、および領域間汎化に対する堅牢性は調査されていません。さらに、フローマッチングは連続的かつ潜在的に可逆的な特徴量輸送のための理論的枠組みを提供しますが15,16,17、学習された変換の厳密な可逆性、トポロジーの保存、および物理的な解釈可能性は、本研究において経験的に検証されていません。また、推論時間、メモリ消費量、浮動小数点演算の複雑さを含む計算特性についても、個別のプロファイリングは行われていません。
マルチモーダルアライメントを研究するための代替的なアプローチは進化し続けている。近年の研究では、マルチモーダルな特徴統合を改善するために、対照学習、グラフベースの融合、低ランク表現学習、プロンプト誘導型マルチモーダル学習、連合学習、およびハイブリッド畳み込み–Transformerアーキテクチャなどが探索されている23,26,27,28,35,36,37,38,39。これらの戦略は相互に排他的なものではなく、相補的な関係にあり、今後の研究では、明示的なフローベースの特徴輸送とこれらの既存のアライメントパラダイムとの組み合わせを検討することが考えられる。土地被覆分類以外にも、相補的なセンシングモダリティがますます利用可能になっているため、セマンティックセグメンテーション、変化検出、物体認識、環境モニタリング、災害評価、精密農業、都市マッピングなどの関連するマルチモーダルリモートセンシングタスクにおいて、連続的な特徴アライメントが適用できる可能性がある2,3,4,5,6,7,8,9,23。
今後の研究では、繰り返し実行による統計分析、計算プロファイリング、欠落またはノイズを含むモダリティに対する堅牢性、空間レジストレーション誤差への感度、および領域間転移評価を組み込むことで、提案されたフレームワークをより多様な実験条件下で評価すべきである。さらに、代替的なフローパラメータ化、異なるエンコーダーアーキテクチャ、およびより大規模なマルチモーダルベンチマークデータセットを検討する追加研究によって、提案手法の汎用性と実用的適用性がより明確になると考えられる。総じて、本結果は、条件付きフローマッチングによる明示的な特徴アライメントが、マルチモーダルリモートセンシング分類にとって有望な補完的戦略であることを示唆しており、同時に、さらなる手法の開発と包括的な評価に向けた多くの機会を提示している。
利益相反:
著者らは、利益相反がないことを宣言します。
本研究は、陝西省自然科学基金プロジェクト(Shaanxi A)の助成を受けた。&F大学(助成番号:ZK25-38)および陝西省教育庁(助成番号:24JK0734)。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| AdamW オプティマイザ | PyTorch Foundation | PyTorch に同梱 | モデルの最適化に使用。PyTorch 2.5.0 を使用した。 |
| Augsburg ベンチマークデータセット | Hu J et al. (Earth System Science Data, 2022); Remote Sensing Technology, LMU Munich | https://doi.org/10.5281/zenodo.7185498 | モデル評価に使用した公開マルチモーダルリモートセンシングベンチマークデータセット。 |
| CUDA Toolkit | NVIDIA | CUDA Toolkit 12.4; https://developer.nvidia.com/cuda-12-4-0-download-archive | GPU加速によるモデルの学習および推論に使用。 |
| グラフィックスプロセッシングユニット (GPU) | NVIDIA Corporation | NVIDIA A100 80 GB PCIe | モデルの学習および評価に使用。 |
| Houston2013 ベンチマークデータセット | IEEE Geoscience and Remote Sensing Society Data Fusion Contest | http://www.grss-ieee.org/community/technical-resources/data-fusion/2013-grss-data-fusion-contest/ | モデル評価に使用した公開ハイパースペクトル画像および LiDAR ベンチマークデータセット。 |
| MUUFL Gulfport ベンチマークデータセット | University of Florida (Gader P et al., Technical Report REP-2013-570) | https://github.com/GatorSense/MUUFLGulfport | モデル評価に使用した公開ハイパースペクトル画像および LiDAR ベンチマークデータセット。 |
| オペレーティングシステム | Canonical Ltd. | Ubuntu 22.04 LTS; https://ubuntu.com | モデルの開発、学習、および評価に使用した計算環境。 |
| PyTorch | PyTorch Foundation | PyTorch 2.5.0; https://pytorch.org/get-started/previous-versions/#v250 | FlowErs モデルの実装、学習、および評価に使用したディープラーニングフレームワーク。 |
| Python | Python Software Foundation | Python 3.11; https://www.python.org/downloads/release/python-3110/ | 計算ワークフローの実装および実行に使用したプログラミング言語。 |