Takuya Miyashita
This site
Web
Start:
* grdvector メモ [#jb30d7b6]
**公式doc [#bcec63ba]
[[Modules » grdvector>https://docs.generic-mapping-tools....
難しい…
**例 [#tb883a04]
2次元の平面分布の Z の色つきコンター図と u, v にもとづく...
ここでは z.grd, u.grd, v.grd の 3 つは全て同じサイズ・同...
#codeprettify{{
## output file
outps="test_grdvec.ps"
## makecpt
cpt="tmp.cpt"
gmt makecpt -Cpolar -T-1/1 -D > $cpt
## filenames
grd="z.grd"
grdu="u.grd"
grdv="v.grd"
## parameters
proj="X10/`gmt grdinfo $grd -Cn -o0,1,2,3 | awk '{print ...
inc=`gmt grdinfo $grd -Cn -o6 | awk '{print 50*$1}'` # 5...
## attributes
vecscale=0.5
vecatt="+a25+e+p1+h1+gblack+n1"
lw=1
## plot
gmt grdimage $grd -J$proj -Baf -BneSW -R$grd -C$cpt -K > ...
gmt psscale -C$cpt -Bxa0.5f0.1 -By+lm -DJMR+w5.0/0.3+o1.0...
gmt grdvector $grdu $grdv -W$lw -I$inc -Si$vecscale -Q$ve...
rm $cpt
}}
**オプション [#w5172c01]
*** -I [#u8af93ed]
x,y 方向の出力間隔.~
n 点おきということではなく,長さの単位で指定する.~
例えばデータが 100 m 間隔で入っていて,20 データごと(2000...
入力ファイルのデータ間隔で割り切れる値でないとエラーを返...
*** -W [#ub2a1014]
矢印の棒の部分の pen attributes.
*** -Q [#mef45938]
Vector arrtibutes の入力.一番重要で難しい.~
基本的にこのオプションは矢印の頭の部分を決めるオプション...
棒の長さは -Si オプションで決まる.
*** -S [#f829e123]
ベクトルの長さのスケール.-Si[length] で決める.~
どのような規模の図に対して適切な値が何なのかは全くわから...
u,v の値によらず長さを固定したい場合は -Sl[length] とする...
*** -T [#z3daea7f]
J オプションで指定している投影面が直交軸でないときもよし...
End:
* grdvector メモ [#jb30d7b6]
**公式doc [#bcec63ba]
[[Modules » grdvector>https://docs.generic-mapping-tools....
難しい…
**例 [#tb883a04]
2次元の平面分布の Z の色つきコンター図と u, v にもとづく...
ここでは z.grd, u.grd, v.grd の 3 つは全て同じサイズ・同...
#codeprettify{{
## output file
outps="test_grdvec.ps"
## makecpt
cpt="tmp.cpt"
gmt makecpt -Cpolar -T-1/1 -D > $cpt
## filenames
grd="z.grd"
grdu="u.grd"
grdv="v.grd"
## parameters
proj="X10/`gmt grdinfo $grd -Cn -o0,1,2,3 | awk '{print ...
inc=`gmt grdinfo $grd -Cn -o6 | awk '{print 50*$1}'` # 5...
## attributes
vecscale=0.5
vecatt="+a25+e+p1+h1+gblack+n1"
lw=1
## plot
gmt grdimage $grd -J$proj -Baf -BneSW -R$grd -C$cpt -K > ...
gmt psscale -C$cpt -Bxa0.5f0.1 -By+lm -DJMR+w5.0/0.3+o1.0...
gmt grdvector $grdu $grdv -W$lw -I$inc -Si$vecscale -Q$ve...
rm $cpt
}}
**オプション [#w5172c01]
*** -I [#u8af93ed]
x,y 方向の出力間隔.~
n 点おきということではなく,長さの単位で指定する.~
例えばデータが 100 m 間隔で入っていて,20 データごと(2000...
入力ファイルのデータ間隔で割り切れる値でないとエラーを返...
*** -W [#ub2a1014]
矢印の棒の部分の pen attributes.
*** -Q [#mef45938]
Vector arrtibutes の入力.一番重要で難しい.~
基本的にこのオプションは矢印の頭の部分を決めるオプション...
棒の長さは -Si オプションで決まる.
*** -S [#f829e123]
ベクトルの長さのスケール.-Si[length] で決める.~
どのような規模の図に対して適切な値が何なのかは全くわから...
u,v の値によらず長さを固定したい場合は -Sl[length] とする...
*** -T [#z3daea7f]
J オプションで指定している投影面が直交軸でないときもよし...
Page:
Edit with a page name which already exists