Transcript StaMPS*** **** Step by Step

StaMPS講習会
解析処理 Step by Step
福島 洋・小澤 拓
2012/9/11-12
（最終更新日：2012/8/23）
StaMPS解析のコツ
• ひとつひとつのステップを確実に！間違えた処
理をした場合、その正しくステップをやり直しても
失敗することがあるため、最悪、最初からやり直
ししなければならなくなる。
• このスライド中、「正常な処理が行われればこう
なる」と書いてあるところは特に重要ポイントです。
そのステップの結果を注意して確認してください。
あと、Warning表示は無視してよい場合が多いで
すが、ERROR表示が出たら何かがおかしい可能
性が大きいです。
解析データ：九重火山
• Sensor: ALOS/PALSAR
• Path: 73, Frame:2950 (Descending, HH,
Offnadir:34.3deg.)
• Dates:
– 20070107, 20071125, 20080527, 20090530,
20100417, 20110420 の６枚
ソフトインストール状況の簡易チェック
•
•
•
•
•
•
roi
make_raw_alos.pl
doris
mt_prep
matlab
(in matlab) which stamps
• あとは、途中で問題が発生した時点で解決しま
しょう。
以下のようにパスを確認し、パスが
通っていなければ修正してください。
• echo $FFTW_LIB
/Users/yo/FUKU/programs/ROI_PAC_3_0_1/ROI_PAC/NetInst
/fftw-111104-1407/lib
• echo $MY_BIN
…./ALOS_preproc/bin/i386
• echo $STAMPS
/Users/yo/FUKU/programs/StaMPS_v3.2.1
• echo $MY_SCR
$STAMPS/ROI_PAC_SCR
• echo $DORIS_SCR
$STAMPS/DORIS_SCR
• echo $MATLABPATH
$STAMPS/matlab
データ・スクリプトの準備
•
•
データ（path-frame-yyyymmdd.tar.gz）を
~/kuju/RAW 内にまとめて置いてください。
suppl.tar.gz内のファイルを解凍＆展開してください（tar xzvf suppl.tar.gz）。ファイルを
次に指定する場所に移動させてください。ここに含まれているM関数ファイルを使え
ば、stampsはMatlab本体のみ（ツールボックス不要）で処理することができます。
– link_raw_pixel
$STAMPS/bin
– prep_alos_raw_pixel
$STAMPS/bin
– link_alos_frame_raw_yf
$STAMPS/bin
– make_cohs
$STAMPS/bin
– step_coh_sb
$STAMPS/bin
– ps_est_gamma_quick.m
$STAMPS/matlab
（このファイルは、同名のファイルの247行目あたりを修正したものです）
– mt_ml_select.m
$STAMPS/matlab
– repair_lonlat_5lraw.m
$STAMPS/matlab
– gausswin.m
$STAMPS/matlab
– 4931_ds.dehm
~/kuju/DEM
ロケールの変更
• 不具合が生じないように、ターミナルのlocale
を英語に変えておく。（日本語環境では、処理
中に用いられるdateコマンドの挙動が異なる
ため。）
• 解析を行うターミナル上で以下のいずれかを
実行。
export LANG=en_US.UTF-8（bashの場合）
setenv LANG en_US.UTF-8（csh, tcshの場合）
これを確実にやらないとハマります！
特に、複数のターミナルを開く場合に注意。
link_raw_pixel
データの解凍とディレクトリの作成
• （PIXELアーカイブのファイル名フォーマットを読み込
むように、オリジナルのlink_rawというシェルスクリプ
トを改変してある。）
• このコマンドにより、解析ディレクトリの下にフレーム
番号（2950）が作成され、 その下にyymmddという名
前のディレクトリが作成される。その中に、RAWデー
タ のリンクが置かれる。
使用法：
link_raw_pixel data_dir analysis_dir
使用例：
link_raw_pixel /Users/yo/kuju/RAW /Users/yo/kuju
（注意：絶対パスで与えること）
（参考）基線長
撮像日 垂直基線長（km）
2007/1/7
0
2007/11/25 -0.1407
2008/1/10 -0.3911
2008/2/25 -0.9280
2008/4/11 -1.2814
2008/5/27 -1.6176
2008/7/12
0.9758
2008/10/12 3.0673
2008/11/27 2.8823
2009/1/12
2.6156
2009/4/14
1.7866
2009/5/30
1.1259
2009/8/30
1.8113
2009/10/15 1.1608
2009/11/30 0.9332
2010/1/15
0.4644
2010/3/2
0.0161
2010/4/17 -0.4705
2010/12/3 -1.3912
2011/1/18 -1.6955
2011/3/5
-2.1100
2011/4/20 -2.6241
◯
◯
◯
◯
◯
◯
step_slc_alos
マスター画像のSLC作成
• 適当なマスターデータを決める。今回は、おおよそ
Bperpおよび観測日が中間値となるデータとして、
20090530とする。
• そこのディレクトリに入り、step_slc_alosを実行。
cd ~/kuju/2950/SLC/20090530
step_slc_alos
注）ALOS_preprocの`setenv MY_BIN (hoge)` (tcsh) 、`export MY_BIN=(hoge)`
(bash) を忘れずに！
以下のようなメッセージが最初に表示されるが気にしない。
ls: No match.
ls: IMG-HH-ALP*1.0__A: No such file or directory
SLCの強度画
像
（ルック数
5 x 15）
レーダ座標系なので、
左右逆であることに
注意。（北行軌道の
データだと、上下逆
になります。）
強度画像image.slc_5l.rasを開く
• GIMP (Linux/Mac) 、ToyViewer (Mac) など
• ここでやること：解析範囲（最終的な干渉画像
の範囲）を決定する
• そのために、まずras画像中で切り出したい領
域の画素番号をしらべる（今回は、X:230-630,
Y:1300-1600あたりを切り出すことにする）
• シングル・ルックの場合にどうなるかを計算。
X, Yそれぞれ5と15をかけた値をメモる。
（x=1150-3150, y=19550-24000)
../roi.proc の編集
切り出し領域を設定
• roi.procの中身
use1dopp=1
mean_pixel_rng = 2075 ←レンジ方向の解析領域中心
ymin = 18000 ←アジマス方向の解析領域下限
ymax = 25000 ←アジマス方向の解析領域上限
• アジマス方向のライン数は、1000以上余裕を与えて
設定する。つまり、最終的な解析領域を基準に、ymin
は1000以上小さい値、ymaxは1000以上大きい値に
設定する。
• ファイルの編集が終わったら、再度step_slc_alosを実
行。切り出された領域でSLCが作成される。
image.slc_5lk.ras
アジマス方向に切り出された画像ができる。
（合成開口長より短くはならなさそう）
master_crop.inの作成・編集
最終的な解析領域を設定
• cp $MY_SCR/master_crop.in .
• vi master_crop.in
（もちろんviでなくgedit等でも可）
first_l 1871
last_l 6000
first_p 1201
last_p 3100
• first_lとlast_lは、アジマス方向のライン数。
• 切り出された後のライン数を基準に設定。
step_master_setup
マスターの準備
cd ~/kuju/2950/SLC/20090530
step_master_setup
• パラメタファイルの作成などとともに、20090530_crop.ras
ができる。この画像が、最終的な解析領域。解析したい領
域が適切に切り出されているかチェックし、必要に応じ
master_crop.inを編集してやり直す。
• もし、20090530_crop.rasが真黒な場合、切り出しに失敗し
ている。master_crop.inで設定した領域がはじめに切り出
した領域外になっているとこうなるので、master_crop.inを
修正してやり直すこと。
make_slc_alos
すべてのスレーブ画像のSLCを作成
cd ~/kuju/2950/SLC
make_slcs_alos
• しばらく時間がかかる。（今回のデータセット（スレーブ
５枚）で20〜30分）
• 各データのディレクトリ内にras形式の強度画像が作
成されているので、 問題なくできているかを確認する。
• Toyviewerを使うなら、rasファイルはToyViewerで開くこ
とを定義しておき、open */*.rasとやると簡単。
Tips
• StaMPSのシェルスクリプトには、step… という
名前と make… という名前のものがある。
• step… のスクリプトは、あるひとつの画像に対
する処理をおこなう。
• make… は、すべての画像に対する処理をお
こなう。step… を順番に実行している。
DEM（DEHM）の準備
• 今回は、国土地理院基盤地図情報10m（0.4秒）
メッシュ標高とEGM96から作成したDEHM(Digital
Ellipsoildal Height Model)を解析用に整形したも
のを用いる。
• DEHMについて：pixelメンバー（共有データアクセ
ス可能な方）は、pixelサーバーからダウンロード
可能。
• 今回使うDEHM（DEM/4931_ds.dehm）は、データ
を半分に間引き（ピクセル間隔を倍にし）、範囲
を狭く切り取ったもの。
DEHM-10m について
http://pixel.eri.u-tokyo.ac.jp/data/DEHM-10m/
小澤さんが数年前に準備。
（[pixel: 986] のメールでアナウンス）
ここは、小澤さんに解説して
もらうといいと思います
用いるDEHM
timing.dorisin の編集
DEMパラメタの設定
• ~/kuju/2950/INSAR_20090530内にtiming.dorisinと
いうファイルがある。これはDORIS（M_SIMAMP,
M_TIMING）用のパラメータファイル。この中のDEM
に関するパラメータを書き換える。
SAM_IN_FORMAT real4
SAM_IN_DEM /Users/yo/kuju/DEM/4931_ds.dehm
SAM_IN_SIZE 2000 2500
// rows cols
SAM_IN_DELTA 0.0002222222 // in degrees
SAM_IN_UL
33.3332778 131.0000556 // lat and lon of upper left
SAM_IN_NODATA 0
SAM_OUT_FILE master_sam.raw // synthetic amplitude
SAM_OUT_DEM dem_sam.raw // cropped DEM for
step_master_timing
Timing（アジマス方向）Errorの推定
• DEMと位置合わせをすることにより、軌道
データから推定されるアジマス方向の位置の
ずれを推定。
cd ~/kuju/2950/INSAR_20090530
step_master_timing
step_master_timing中、下のようなWarningが出るとうまくいかな
い場合がある。もし、次のページのような結果が出なければ、
失敗していることになるので、修正してstep_master_timingをや
り直す。
こうなる原因は二つ（かそれ以上）ある。１：データの「時間」が
おかしい（ロケールの変更をしていないとこうなる）、２：DEMの
範囲が狭い。これは、DEMの範囲を拡げれば解決する。
•
•
•
•
•
•
•
•
•
•
•
WARNING : DEM does not cover entire interferogram.
WARNING : input DEM should be extended to the North.
WARNING : indexphiNDEM: 2085
WARNING : DEM does not cover entire interferogram.
WARNING : input DEM should be extended to the South.
WARNING : indexlambda0DEM: -126
WARNING : DEM does not cover entire interferogram.
WARNING : input DEM should be extended to the West.
WARNING : indexlambdaNDEM: 2459
WARNING : DEM does not cover entire interferogram.
WARNING : input DEM should be extended to the East.
step_master_timing最後に出るメッセージ
Number posL
posP
offsetL offsetP correlation
0 2914
1732
24 0 0.466933
1 2914
1991
22 0 0.449156
2 2914
2251
22 0 0.416776
.
.
.
26 4955
1789
22 0 0.271683
27 4955
2048
22 0 0.296915
28 4955
2308
24 0 0.325166
29 4955
2567
22 0 0.385394
Estimated total offset (l,p):
22, 0
Coherence NaN values are disregarded in the analysis.
***************************************************************
****
correlationが0.1以下だったり、Estimated total offset (l,p)が数百以上？になっている
場合、多分うまくいっていない。DEMのWarningが出ていないか確かめ、出ていたら
前頁に従い問題を解決する。
make_orbits
INSAR_20090530のディレクトリ内に各ス
レーブ用のディレクトリが作成される
make_coarse
大雑把な位置合わせ
cd ~/kuju/2950/INSAR_20090530
make_orbits
make_coarse
make_coreg
精密位置合わせ
• 通常、このステップは時間がかかる。
• 時間節約のためにパッチ数を調整してから実行してもよい。
• パッチ数の変更：$STAMPS/matlab/coreg_pos.mのファイルを
編集。nWinxとnWinyをそれぞれ20と50等に変更する。
• 注）MATLABPATHの設定`setenv MATLABPATH
/home/hoge/StaMPS_v3.2.1/matlab`を忘れずに！
cd ~/kuju/2950/INSAR_20090530
make_coreg
MacでR2011?以降のmatlabを使用している人は、以下の様なメッセージが出るかもしれ
ないが、無視してよい。
2012-08-27 00:03:55.953 MATLAB[16905:8103] This process is attempting to exclude an
item from Time Machine by path without administrator privileges. This is not supported.
make_dems
地形縞のシミュレート画像を作成
• マニュアルによると、地形の起伏が大きくなけ
れば、実行しなくてよいとのこと。
• refdem_11.raw などのファイルが各スレーブ
に対して作られる。
• この処理も結構時間がかかる。
cd ~/kuju/2950/INSAR_20090530
make_dems
make_resample
スレーブSLCのリサンプリング
• マスター画像と同じグリッドのスレーブ画像が
作成される。所要時間５分程度。
• このあと、リサンプルされた画像のファイルサ
イズがすべて同じになっているかを確認する。
cd ~/kuju/2950/INSAR_20090530
make_resample
ls –l */slave_res.slc
make_ifgs
干渉画像の作成
• 干渉画像 cint.minrefdem_5l.ras が各スレー
ブディレクトリ内にできる。すぐ終わる。
• できた干渉画像は、Matlabでplot_all_ifgsを使
うか、gimpやToyViewerなどで確認する。
• 全く縞が確認できないようなものは、失敗して
いる可能性が高い。本演習の例題の場合、
処理がうまくいけば干渉性は高くなるはず。
DEHM with the 20090530 as the master
plot_all_ifgs(0,3)
step_geo
ジオコーディング
• これは、後のPS処理のために必要。 スレーブ
画像のディレクトリのどれか一つだけに移動
し、コマンドを実行する。
• 実行後、InSAR_20090530にlat.raw、lon.raw
が作られる。
cd ~/kuju/2950/INSAR_20090530/20070107
step_geo
ここからPS処理
mt_prep
PS候補点の選定
cd ~/kuju/2950/INSAR_20090530/
mt_prep 0.4
（mt_prepのオプション）
mt_prep 0.4 2 2 50 200
（意味）
0.4: Amplitude dispersionのしきい値（これ以下を用いる）
2: レンジ方向のパッチ数（デフォルト1）
2: アジマス方向のパッチ数（デフォルト1）
50: レンジ方向のパッチ間のオーバーラップピクセル数（デフォルト50）
200: アジマス方向のパッチ間のオーバーラップピクセル数（デフォルト
200）
StaMPSの処理はパッチごとに行われ、あとでマージされる。解析エリアに
対してメモリに余裕があれば、パッチに分ける必要はない。
以下の処理は、matlab上で実行す
る。matlabを起動して、処理ディレク
トリに移動する。その後、処理を実
行するという流れが基本。
(getparm)
処理パラメタの確認
• ~/kuju/2950/INSAR_20090530/ で、
>> getparm
と打つと、処理パラメタが出力される。
• コマンド setparm でパラメタを変更することが
できる。通常、変更する必要はない。パラメタ
の説明はマニュアルにあるが、すべてではな
いかもしれない。
R2012aだとエラーになる部分あり。
その回避方法：
• matlabのコマンドラインで、”edit uw_interp”と
打ち、uw_interp.mをエディターで開く
• ５５行目を以下のように変更。
（変更前）
Z=dsearch(x,y,ele(:,2:4),X,Y); %index from grid to pixel node
（変更後）
Z=dsearchn([x,y],ele(:,2:4),[X(:),Y(:)]); %index from grid to pixel node
Z = reshape(Z,nrow,ncol);
stamps
PS処理の実行
• 何も指定しなければ、全ステップの処理を自
動的にしてくれる。
• stampsの処理は８つのステップに分かれてお
り、オプションを指定することにより各々のス
テップのみを実行することも可能。詳しくは、
help stamps 参照。
>> cd ~/kuju/2950/INSAR_20090530
>> stamps
新しいMATLABのバージョンだと、以下のエラーが出る（少し古いバージョンだと
Warning）。
Error using griddata (line 62)
GRIDDATA no longer requires Qhull-specific options.
Please remove these options when calling GRIDDATA.
Error in ps_load_initial (line 258)
la0_ps=griddata(gridX,gridY,la0,ij(:,3),ij(:,2),'linear',{'QJ'});
Error in stamps (line 116)
sb_load_initial
このエラーを回避するために、ps_load_initial.m でgriddataを呼んでいる行で、オプ
ションの ,{‘QJ’} を消す。
（修正前）griddata(gridX,gridY,la0,ij(:,3),ij(:,2),'linear',{'QJ'});
（修正後）griddata(gridX,gridY,la0,ij(:,3),ij(:,2),'linear');
ps_plot
結果のプロット・エクスポート
• このmatlab関数には、非常の多くの機能があ
る。詳しくは、 help ps_plot や、StaMPSマニュ
アルの Chapter 9 参照。
ps_plot(‘u-dm’,0,[-4 4],1) の結果
ps_plot 使用例：kmlファイル出力
>> ps_plot('v-d',-1,[-20 20])
>> load ps_plot_v-d
>> ps_gescatter('kuju-psinsar.kml',ph_disp,1,1)
Tips:
• ps_plotの２つ目の変数に-1を入れると、プロットされ
るはずのデータがmatファイルにエクスポートされる。
matファイルの名前は、ps_plot_処理方法指定.mat
となる。
• この例では、一度エクスポートしたものをロードし、
ロードされたph_dispという変数をさらにkmlファイル
にエクスポートしている。
Google Earthで見た結果（平均速度）
[-20 20] mm/year
Small Baseline Processing
短基線長（SB）解析
• PS解析とは別の時系列解析を行う手法
• 「StaMPS」は、正確にはPS解析の部分を指す。
• StaMPS/MTI（Multi-temporal InSAR）は、StaMPSの拡張版で、SB解
析や、PS解析とSB解析を合わせた解析を行うことができる。
• StaMPS/MTIのSB解析は、single-look干渉画像で処理を行うので、
基本的にPS解析と類似した結果になる。（multi-look干渉画像で行
うプログラムも公開されていないがあるらしく、問合せ中。）
• SB解析のメリット：アンラッピングの信頼度評価ができる。PS解析と
違うアルゴリズムで解くので、結果の比較が可能。PS+SB解析する
ことにより、信頼性が向上、ピクセル数増加。
（PS解析をせずSB解析をする場合）
• MATLABで基線長に関する情報をロードする
ために以下を実行。
mt_extract_info
matlab
>>ps_load_info
sb_find
短基線長ペアの選定
• コヒーレンス、時間基線長、垂直基線長を基
準に短基線長ペアを選定する。
(in matlab)
>>sb_find(0.5, 5000, 4000)
（0.5: コヒーレンスのしきい値、5000: 臨界時間基線長、
4000: 臨界垂直基線長）
コヒーレンス ＝ (1 – Bt/Btc)*(1-Bp/Bpc) というモデルから計算
したコヒーレンスがしきい値を超えるものを選出する。
表示される図を見てみて、線でつながっていない点があった
り複数のクラスタに分かれたりする場合は、パラメタを調整し
てつながるようにする。
sb_find（続き）
• sb_findを実行すると、small_baselines.listとい
うファイルができる。その中にペアが書かれ
ている。
• 自分でこのファイルにペアを書き込んでもよ
い。
• plot_sb_baselinesを実行すると、
small_baselines.listに含まれるペアが図示化
される。
make_small_baselines
干渉処理
• PS干渉画像との違いは、common bandpath
filterをかけるというところ。range and/or azimuth
フィルタをかけないようにするというオプションも
ある（マニュアル参照）。
• この処理により、INSAR_20090530の下に
SMALL_BASELINESというディレクトリができ、そこ
に干渉画像が格納される。
cd ~/kuju/2950/INSAR_20090530
make_small_baselines
in SMALL_BASELINES
plot_all_ifgs(0,4,10)
mt_prep
候補点の選定
cd ~/kuju/2950/INSAR_20090530/SMALL_BASELINES
mt_prep 0.6
mt_prepのオプション:SB解析編
(PSの場合とほとんど意味は同じ）
mt_prep 0.6 2 2 50 200
（意味）
0.6: Amplitude difference dispersionのしきい値（これ以下を用いる）
2: レンジ方向のパッチ数（デフォルト1）
2: アジマス方向のパッチ数（デフォルト1）
50: レンジ方向のパッチ間のオーバーラップピクセル数（デフォルト50）
200: アジマス方向のパッチ間のオーバーラップピクセル数（デフォルト200）
StaMPSの処理はパッチごとに行われ、あとでマージされる。解析エリアに
対してメモリに余裕があれば、パッチに分ける必要はない。
stamps
SB処理の実行
• 何も指定しなければ、全ステップの処理を自動的に
してくれる。
• stampsの処理は８つのステップに分かれており、オ
プションを指定することにより各々のステップのみを
実行することも可能。詳しくは、help stamps 参照。
• マニュアルによると、PS処理よりメモリをだいぶ食う
とのこと。パッチに分けて解析する必要性が生じる
かもしれない。
>> cd ~/kuju/2950/INSAR_20090530/SMALL_BASELINES
>> stamps
新しいMATLABのバージョンだと、以下のエラーが出る（少し古いバージョンだと
Warning）。
Error using griddata (line 62)
GRIDDATA no longer requires Qhull-specific options.
Please remove these options when calling GRIDDATA.
Error in sb_load_initial (line 227)
la0_ps=griddata(gridX,gridY,la0,ij(:,3),ij(:,2),'linear',{'QJ'});
Error in stamps (line 116)
sb_load_initial
このエラーを回避するために、sb_load_initial 116行目、117行目、254行目
の ,{‘QJ’} を消す。
la0_ps=griddata(gridX,gridY,la0,ij(:,3),ij(:,2),'linear');
la1000_ps=griddata(gridX,gridY,la1000,ij(:,3),ij(:,2),'linear');
bp0_ps=griddata(gridX,gridY,bp0,ij(:,3),ij(:,2),'linear');
ps_plot(‘u-dm’,0,[-4 4],1) の結果
(default parameters)
ps_plot(‘usb’,0,[-4 4],1)
（干渉画像）
ps_plot(‘rsb’,0,[-4 4],1)
（アンラップ誤差）
ピクセル数を増やすため、やり直し。
• SBの場合、’density_rand’がデフォルトで2に設定され
ている（PSは20）。
• setparm(‘density_rand’,20) とすると、PSと同じパラメタ
値でやることになる。
• このパラメタは、km^2あたりに許されるピクセル数。
• このパラメタが関係するのは、Step 3なので、その先
からやり直せばよい（stamps(3,8))。
• なお、先に”setparm”をしてから”stamps”をしても、自
動的にdensity_randが2に書き換えられてしまった。
stamps(1,2)→setparm(‘density_rand’,20)→stamps(3,8)
とやれば多分OK。
ps_plot(‘u-dm’,0,[-4 4],1) の結果
(density_rand = 20)
ps_plot(‘u-dm’,0,[-4 4],1) の結果
(default parameters)
PSとSB解析結果の融合
cd ~/kuju/2950/INSAR_20090530
ps_sb_merge
これにより、INSAR_20090530の下に MERGED というディレクトリができ、
そこにマージしたデータが置かれる。これは、stamps処理のステップ５ま
で（アンラッピングの手前まで）の処理結果を融合したものである。よって、
最終結果を得るためには残りのステップ６〜８を実行する。
cd MERGED
stamps(6,8)
これらの処理は、すぐ終わります。
ps_plot(0,[-4 4],1) in MERGED
マルチルック処理をしたSB解析
• この機能は、公開パッケージには含まれてい
ませんが、開発者（Andy Hooper）の許可をも
らいました。（←打診中）
• 通常のSB解析と同様に、マルチルック処理を
した干渉画像を作成し、それを用いて最小二
乗法で解きます。
ルック処理版SB解析
• cd ~/kuju/2950/INSAR_20090530/SMALL_BASELINES
• ここで make_ifgs.listを作成
/Users/yo/kuju/2950/INSAR_20090530/20070107_20071125
/Users/yo/kuju/2950/INSAR_20090530/20070107_20080527
…（省略）…
/Users/yo/kuju/2950/INSAR_20090530/20080527_20110420
/Users/yo/kuju/2950/INSAR_20090530/20090530_20100417
• make_cohs （！今講習会で配布するバージョン）
• mt_extract_info
(in matlab)
• >> cd ~/kuju/2950/INSAR_20090530
• >> repair_lonlat_5lraw （今講習会で配布。バグ？に対応）
• >> cd SMALL_BASELINES
• >> mt_ml_select （今講習会で配布）
• >> stamps(6,8)
ps_plot(0,[-4 4],1)
SMALL_BASELINES (マルチルック版）
ps_plot(0,[-4 4],1) in MERGED
ps_plot(‘u-dm’,0,[-4 4],1) の結果
(SB only, density_rand = 20)
ps_plot(‘u-dm’,0,[-4 4],1) の結果
(PS only, density_rand = 20)
LOSベクトルの計算
• heading角と入射角が必要。
• heading角は、stamps処理の最中に、
heading.1.in という名前のファイルに書き出され
る。ROI_PACのパラメタファイル（*.rsc）にも書いて
ある。
• 入射角は、干渉画像が置かれているスレーブの
ディレクトリ（INSAR_20090530/yyyymmdd、どれ
でもよい）にあるdorisのファイルcoreg.out,
dem.out, geocode.outに書かれている。
（grep incidence coreg.out などとやると簡単）
その他諸々
• ps_plotには、時系列をプロットするオプション
がある。少し試してみたところ、プロット領域を
大きくしたときは異常終了してしまった。
• データのエクスポートは、ps_plotだけでなく、
ps_outputでもできる。
• ref_lon, ref_lat（またはref_centre_lonなど）で
不動域を指定することができます（ちゃんとで
きるかは未確認です）
参考文献など
•
•
•
•
•
•
•
StaMPSマニュアル（StaMPSパッケージに含まれる。下のウェブサイトからも入手可能）
StaMPSウェブサイト
http://radar.tudelft.nl/~ahooper/stamps/index.html
ユーザーグループ（過去ログの検索可）
http://groups.google.com/group/mainsar
StaMPSアルゴリズム
Hooper, A., P. Segall, and H. Zebker (2007): Persistent scatterer InSAR for crustal
deformation analysis, with application to Volcan Alcedo, Galpagos, J. Geophys. Res., 112,
B07407, doi:10.1029/2006JB004763.
StaMPSを用いた解析処理の流れの説明
福島 洋(2011): StaMPS パッケージを用いた PS 干渉 SAR 解析, 測地学会誌, 57, 41-48.
PSとSB解析の融合アルゴリズム
Hooper, A. (2008): A multi-temporal InSAR method incorporating both persistent
scatterer and small baseline approaches. Geophys. Res. Lett., 35:L16302.
アンラッピングアルゴリズム
Hooper, A. (2009): A statistical-cost approach to unwrapping the phase of InSAR time
series, Proceedings FRINGE Workshop, Frascati.
http://earth.eo.esa.int/workshops/fringe09/proceedings/papers/p1_26hoop.pdf