トレンド除去など

サンプル

重力ブーゲー異常

• grdtrendは双3次傾向面を決定する。
  • grdtrend data.grd -N10 -Ttrend.grd
    
  • grdtrend data.grd -N10 -Dres.grd
    
  • data.grd:入力ファイル，　N：fitting parameter 1 - 10 の整数を入れる
    
  • -T：トレンドファイルを出力，　-D：残差ファイルを出力
    

Script

オリジナルデータ

• grdimage data.grd -R135/142/25/33 -JM10 -Cbg267_grad.cpt -P -K -V > nama.eps
  
• grdcontour data.grd -R135/142/25/33 -JM10 -C10 -W1a5 -L-500/500 -K -O -V >> nama.eps
  
• psbasemap -R135/142/25/33 -JM10 -Bf1a5g5 -O -V >> nama.eps
  
• サンプル図
  

N=10でトレンド計算

• grdtrend data.grd -N10 -Ttrd.grd
  
• grdimage trd.grd -R135/142/25/33 -JM10 -Cbg267_grad.cpt -P -K -V > index-t.eps
  
• grdcontour trd.grd -R135/142/25/33 -JM10 -C10 -W1a5 -L-500/500 -K -O -V >> index-t.eps
  
• psbasemap -R135/142/25/33 -JM10 -Bf1a5g5 -O -V >> index-t.eps
  
• サンプル図
  

N=10の残差

• grdtrend data.grd -N10 -Dres.grd
  
• grdimage res.grd -R135/142/25/33 -JM10 -Cbg267_grad.cpt -P -K -V > index0.eps
  
• grdcontour res.grd -R135/142/25/33 -JM10 -C10 -W1a5 -L-500/500 -K -O -V >> index0.eps
  
• psbasemap -R135/142/25/33 -JM10 -Bf1a5g5 -O -V >> index0.eps
  
• サンプル図
  

N＝１の場合（トレンドは直流値）

• トレンド
  
• 残差