psmask

データの無いところにはマスクをかけてコンターを引く。


シェルスクリプトはこちら
#! /bin/tcsh -f
set file=test.data
set psfile=$file.ps
set gfile=$file.grd
set range=0/100/0/100
set scale=X16/16
set resol=150
set color=cpt.dat
## makecpt -T0/100/10 -Crainbow > $color
surface $file -G$gfile -I$resol -R$range -T1.0
psmask $file -J$scale -R$range -I$resol -K >! $psfile
grdimage $gfile -J$scale -R$range -C$color -O -K >> $psfile
grdcontour $gfile -J$scale -R$range -C10 -A10f13 -O -K >>$psfile
psmask -C -O >>$psfile
gs $psfile


surface
アスキーデータをならしてgridファイル(バイナリ)に落す
-T..ある関数で滑らかにデータを整理する。1だと弱く、0だと強く滑らかにする

psmask
アスキーデータを読んでマスクの準備?。

grdimage
グラデーションを描く。

grdcontour
コンターを描く。
最後のpsmask..マスクの実行?

注:psmaskコマンドでgrdimage,grdcontourを囲まないとうまく行かない。

出来上がりの図

モドル

Last modified: Sat May 11 18:49:55 JST 2002