第 490 回電波研連 F 分科会開催日 :2004 年 12 月 17 日 4 成分散乱電力分解による Pi-SAR 画像の解析石堂基山口芳雄山田寛喜新潟大学

Size: px

Start display at page:

Download "第 490 回電波研連 F 分科会開催日 :2004 年 12 月 17 日 4 成分散乱電力分解による Pi-SAR 画像の解析石堂基山口芳雄山田寛喜新潟大学"

ただきよひきぎ
5 years ago
Views:

1 第 9 回電波研連 F 分科会開催日 : 年月 7 日成分散乱電力分解による Pi-AR 画像の解析石堂基山口芳雄山田寛喜新潟大学

2 発表の流れ. 研究の背景目的 POLAR 画像解析を行う背景目的解析に使用した領域 POLAR 画像について. 三成分散乱モデル分解について Covrine mtrix について三成分散乱モデル分解三成分散乱モデルと Covrine mtrix の画像三成分散乱モデルの問題点と四成分散乱モデル分解の提案 3. 四成分散乱モデル分解について平均化 Covrine mtrix 基本ターゲットにおける平均化 Covrine mtrix 確率密度関数の変更. 解析結果市街地モデルの解析 ( 鳥屋野潟周辺画像 ) 植生モデルの解析 ( 苫小牧画像 ) 5. まとめ

3 . 研究の背景目的

4 背景地球観測地球環境観測大規模災害の迅速な把握地球の資源探査航空機に搭載された偏波合成開口レーダ (Pi-AR) により得られた POLAR 画像データを用いて地表面の解析を行う波情報を用いることにより地上ターゲットの識別がある程度可能偏波を用いた振幅情報位相情報の利用合成開口処理による高分解能の実現

5 目的市街地や植生の領域を検出するためには? 物理モデルを基本とした三成分散乱モデル分解を用いる定式化に問題がある均化 Covrine mtrix を用いた四成分散乱モデル分解の提案三成分散乱モデル分解と四成分散乱モデル分解を様々な方法で比較四成分散乱モデル分解の妥当性を評価 : Freemn, Durden らによる

6 解析領域 () Rnge diretion Azimuth diretion 画像データの諸元解析領域観測日時入射角偏波周波数波長画像サイズピクセルサイズ平均化サイズ新潟市鳥屋野潟周辺 3 年 8 月日 3.~7.(degree) VH.7[GHz] (L-nd) λ3.6(m) (pixel).5.5(m) 5 5 (pixel)

7 解析領域 () Rnge diretion Azimuth diretion 解析領域観測日時入射角偏波周波数波長画像サイズピクセルサイズ平均化サイズ画像データの諸元北海道苫小牧年月 8 日 36.~9.(degree) VH.7[GHz] (L-nd) λ3.6(m) (pixel).5.5(m) 5 5

8 POLAR 画像解析について集合平均 ttering mtrix [ ] Covrine mtrix Cohereny mtrix Kennugh mtrix [ ] C [ T ] [ K ] Physil Mthemtil Phenomenologil 三成分散乱モデル分解 Entropy, α Huynenn Prmeter

9 . 三成分散乱モデル分解について

10 Covrine mtrix() 散乱行列 [] を Covrine mtrix [C] に変換する集合平均をとれるため平均化処理が容易 [ ] [ C] VH N N T 独立要素によりベクトル化 : アンサンブル平均

11 Covrine mtrix() 特徴森林などの分布した自然ターゲットに対して偏波チャネルの電力によりほとんどの要素が決まる個々のターゲットに依存しない次統計量として用いることができる実験的にわかっているRefletion symmetry なターゲットに対して三成分散乱モデル分解を使うことができる C

12 三成分散乱モデル分解 () 測定散乱波を物理的な散乱過程に基づいた散乱モデルに分解 Covrine Mtrix の要素を用いて分解を行う表面散乱二回反射海域農地低植生域における散乱過程地表面に入射して樹幹や人工建造物に反射する散乱過程体積散乱ランダムに傾いた Wire が合成された散乱過程

13 三成分散乱モデル分解 () C C C 33 観測された Covrine Mtrix を以下のように 3 つの成分に展開する [ C] f [ C] f [ C] f [ C] s s β fs β f β f α f f f v s / 3 d d f v surfe f v d β f d doule α α Re( ) < β C 3 fsβ fdα fv / 3 Re( ) > C, α v α volume fv 3 3 P s f s P d f d 8 3 各成分による散乱電力 P v fv 3 ( ) β ( ) α 8

14 三成分散乱モデル分解 (3) Ps Pd Pv min mx

15 平均化 Covrine Mtrix の要素 () [ C]

16 平均化 Covrine Mtrix の要素 () Re [ C] m

17 平均化 Covrine Mtrix の要素 (3) Re ( ) ( m ) [ C]

18 平均化 Covrine Mtrix の要素 () Re ( ) ( ) m [ C]

19 平均化 Covrine Mtrix の要素 Re Re ( ) Re ( ) m m ( ) m ( )

20 四成分散乱モデル分解の提案三成分散乱モデル分解の欠点必ずしもが成立しないが必ず体積散乱成分にならない定式化において強制的ににして計算している部分があるこれらの問題を解決する分解法を提案すればより詳細に地表面を解析することが可能四成分散乱モデルを提案

21 3. 四成分散乱モデル分解について

22 平均化 Covrine mtrix () [ ] VH ( ) [ ] os sin sin os os sin sin os vv vh hv hh ( ) [ ] ( ) [ ] ( ) π d p C C 数式的な平均 Covrine mtrix ( ) [ ] vv hv vv hh vv vv hv hv hh hv vv hh hv hh hh C ( ) p は確率密度関数 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) } Re{ Re Re Re Re Re Re vv v hv h vv h hv vv hh ( ) ( ) ( ) ( ) ( ) π π π π π d p d p d p d p d p sin os os sin sin sin os 5 3 ( ) ( ) ( ) ( ) ( ) π π π π π d p d p d p d p d p os os os sin os sin os sin sin

23 平均化 Covrine mtrix () 確率密度関数を一様であると仮定し hh hh hv hh hh vv hv hv 8 p ( ) ( ) vv 8 vv hh Re 8 j hv vv m{ ( )} j hv vv m{ ( )} p π ( ) のように選べば実数実数実数虚数虚数散乱行列 [] を用いて基本ターゲットの平均化 Covrine mtrix を求める [ C ] VH [ ]?????????

24 基本ターゲットによる平均化 Covrine mtrix() ワイヤターゲット垂直ワイヤ水平ワイヤ [ ] vwire [ ] hwire [ C] wire 二回反射ターゲット ( 金属の面リフレクタ ) 垂直 diplne 水平 diplne [ ] vdiplne [ ] hdiplne [ C] diplne 二次統計量の利点

25 基本ターゲットによる平均化 Covrine mtrix() 奇数回反射ターゲット (Plte,phere,3 面リフレクタ ) [ ] plte プレート [ ] plte C 円偏波発生ターゲット (Helix) R-Helix L-Helix [ ] j j helix r [ ] j j helix l [ ] j j j j C helix r [ ] j j j j C helix l

26 基本ターゲットによる平均化 Covrine mtrix(3) hh hv j hv vv m{ ( )} の関係より m{ ( ) } m( ) Helix のみ存在 R-Helix j [ ] rhelix j j m L-Helix j [ ] lhelix j j m 円偏波発生成分 f P m f ( ) m{ ( )} P

27 確率密度関数の変更 () 垂直に立っている幹枝が多い植生を想定して確率密度関数を変更する sin p( ) ( < < π ) ( π < < π ) π p( ) d

28 確率密度関数の変更 () ( ) ( ) ( ) Re Re vv v hv h vv h hv vv hh sin ) ( のとき ( ) ( ) ( ) ( ) ( ) vv hv hv hh vv hh hv vv hh Re Re [ ] wire h 水平ワイヤ [ ] wire h C [ ] wire v 垂直ワイヤ [ ] wire v C

29 平均化 Covrine mtrix の展開 () 地面と木の幹により構成される回反射構造を想定誘電体 hh 成分と vv 成分の位相差を考慮した複素数 α を導入 e jγ α h gh th j γ v [ ] e R R tree R R e jγ v R gv R tv gh th α e j( γ γ ) h v R R gh gv R R th tv hh vv R t _ : 樹幹による反射係数 : 地面による反射係数 γ : 伝播定数 (_ はそれぞれ, 水平のhか垂直のv) R g この回反射構造は角度に対して植生のようにランダムとはならないので角度に対する積分は行わない [ C] doule α α α

30 平均化 Covrine mtrix の展開 () 表面散乱モデルの基本 Covrine mtrix 次 Brgg 反射モデルを用いて β [ ] surfe 円偏波発生電力 C surfe Covrine mtrixにおける, 3 の値を用いると f は f m m β β ( ) { } m { } 3 β

31 四成分散乱モデル分解 () Azimuth diretion Rnge diretion zimuth 方向に対して傾いた市街地偏波が多く発生全てが体積散乱成分の寄与となり得る Ps P d Pv P の四成分を用いた四成分散乱モデル分解を行う

32 四成分散乱モデル分解 () 観測された Covrine Mtrix を以下のようにつの成分に展開する [ C] f [ C] f [ C] f [ C] f [ C] f s β β s fs surfe d β f d doule α α v volume α 8 fv 5 α, β : 相対要素 ( 未定係数 ) irulr 3 f m j p( ) sin, fd, fv, f : 表面散乱, 二回反射, 体積散乱の円偏波発生成分の寄与 ± m j j ± j C C C 33 C 3 C, 8 fs β fd α fv 5 f fv 5 3 f fs fd fv 5 f fsβ fdα fv 5 f m{ f } Re( ) < β Re( ) > α 各成分による散乱電力 P s f s P d f d P v f v P f ( β ) ( α )

33 . 解析結果

34 解析結果 ( 鳥屋野潟周辺 ) Ps Pd 四成分散乱モデル分解

35 解析結果 ( 鳥屋野潟周辺 ) P v P 四成分散乱モデル分解

36 解析領域の説明 pth 3 pth pth

37 解析結果 (pth ) Ps P 三成分散乱モデル分解 d Pv Ground truth Ps Pd Pv P

38 pth Azimuth diretion Rnge diretion Ps Pd 三成分散乱モデル分解 Pv Ps Pd Pv P

39 解析結果 (pth ) Ps P 三成分散乱モデル分解 d P v Ground truth Ps Pd Pv P

40 pth Azimuth diretion Rnge diretion d Ground truth d : 平成大橋 : 水道橋 : 新潟バイパス d : 新幹線高架橋 Ps Pd 三成分散乱モデル分解 Pv d Ps Pd Pv 四成分散乱モデル分解 P

41 解析結果 (pth 3) Ps P 三成分散乱モデル分解 d Pv Ground truth Ps Pd Pv P

42 pth 3 Azimuth diretion Rnge diretion Ps Pd 三成分散乱モデル分解 Pv Ps Pd Pv 四成分散乱モデル分解 P

43 市街地モデルのまとめ四成分散乱モデル分解ではアジマス方向に平行な市街地であるほど二回反射成分が強くなるアジマス方向に対して平行な市街地で発生していた体積散乱成分を軽減し円偏波発生成分への寄与に分散三成分散乱モデル分解よりも物理現象に近いことを確認

44 解析領域の説明 pixel 偏波シグネチャの変化の様子解析領域 Diplne 7.5m おきに偏波シグネチャを作成 Horizontl wire

45 解析結果 ( 苫小牧 ) 三成分散乱モデル分解 Ps Pd 四成分散乱モデル分解

46 解析結果 ( 苫小牧 ) 三成分散乱モデル分解 P v 四成分散乱モデル分解

47 解析結果 ( 苫小牧 ) f β f α s d f v f v / 3 f f s d f v Ps Pd 三成分散乱モデル分解 (fv を除去しないもの ) Pv f β f α s d f v / 3 Ps Pd P v 四成分散乱モデル分解 P

48 5. まとめ

49 まとめ三成分散乱モデル分解の物理的な矛盾点を解消した四成分散乱モデル分解を提案矛盾点自然ターゲットで必ずしもが成立しない三成分散乱モデル分解と四成分散乱モデル分解を単純な領域で比較市街地領域と植生領域を選択三成分散乱モデル分解の時に生じた二回反射成分における水域の不要な電力が除去された人工ターゲットではアジマスに対するターゲットの傾きによって二回反射受信電力が変化する樹幹や沢の部分において二回反射成分と想定できる箇所が微量ではあるが検出された

50 補足資料 ( コンポジット画像 ) 三成分散乱モデル分解 Pd Pv Ps 四成分散乱モデル分解

51 補足資料 (α, β の分布 ) α β ( 白がα-の部分 ) ( 白がβの部分 )

52 補足資料 (α, β の分布 ) α β ( 緑がα-の部分 ) ( 緑がβの部分 )

53 補足資料 (X-nd, 三成分 ) Ps Pd Pv

54 補足資料 (X-nd, 四成分 ) Ps Pd P v P

画像処理工学

画像処理工学画像の空間周波数解析とテクスチャ特徴フーリエ変換の基本概念信号波形のフーリエ変換信号波形を周波数の異なる三角関数 ( 正弦波など ) に分解する逆に, 周波数の異なる三角関数を重ねあわせることにより, 任意の信号波形を合成できる正弦波の重ね合わせによる矩形波の表現フーリエ変換の基本概念フーリエ変換次元信号 f (t) のフーリエ変換変換 ( ω) ( ) ωt F f

第 490 回電波研連 F 分科会開催日 :2004 年 12 月 17 日 4 成分散乱電力分解による Pi-SAR 画像 の解析 石堂基 山口芳雄 山田寛喜 新潟大学

第 490 回電波研連 F 分科会開催日 :2004 年 12 月 17 日 4 成分散乱電力分解による Pi-SAR 画像の解析石堂基山口芳雄山田寛喜新潟大学