Gridleme ve filtreleme
From GMT Türkiye Wiki
(Sürümler arası farklar)
| (Bir kullanıcı tarafından yapılan 17 ara revizyon gösterilmiyor) | |||
| 1. satır: | 1. satır: | ||
| − | + | * [[ batimetri verisini gridleme ve spike'ları maskeleme]] | |
| − | |||
| − | |||
| − | + | Batimetri verisi indirmek için [[Media:Bathy.txt|Bathy.txt]] | |
| + | <bash> | ||
| + | gmtset MEASURE_UNIT inch | ||
| + | set dosya = Bathy.txt | ||
########get minmax coordinates | ########get minmax coordinates | ||
set R = `minmax $dosya -C | awk '{print int($1/10)*10-20,int($2/10)*10+20,int($3/10)*10-20,int($4/10)*10+20}'` | set R = `minmax $dosya -C | awk '{print int($1/10)*10-20,int($2/10)*10+20,int($3/10)*10-20,int($4/10)*10+20}'` | ||
| 71. satır: | 72. satır: | ||
### 3D plotting see the spikes !!!!! | ### 3D plotting see the spikes !!!!! | ||
echo $R | echo $R | ||
| − | set nam = | + | set nam = 3dspykes.ps |
set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] | set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] | ||
set depth = ` grdinfo -C spiky.grd | awk '{print $6,$7}'` | set depth = ` grdinfo -C spiky.grd | awk '{print $6,$7}'` | ||
| 77. satır: | 78. satır: | ||
xv $nam | xv $nam | ||
| + | </bash> | ||
| + | |||
| + | [[Resim:3dspykes.jpg]] | ||
| + | <bash> | ||
############# remove the spikes now | ############# remove the spikes now | ||
| 99. satır: | 104. satır: | ||
### 3D plotting | ### 3D plotting | ||
set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] | set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] | ||
| + | set name = 3dnospykes.ps | ||
| + | grdview nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -Jz0.08 -Qi300 -V -Z0 -B250 > ! $name | ||
| + | xv $name & | ||
| + | </bash> | ||
| − | grdview nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -InospikeI.grd -Jz0.08 -Qi300 -V -Z0 -B250 > ! $ | + | |
| − | xv $ | + | [[Resim:3dnospykes.jpg]] |
| + | |||
| + | |||
| + | <bash> | ||
| + | ## with intensity | ||
| + | grdview nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -InospikeI.grd -Jz0.08 -Qi300 -V -Z0 -B250 > ! $name | ||
| + | xv $name & | ||
</bash> | </bash> | ||
15:49, 18 Şubat 2010 itibarı ile sayfanın şu anki hâli
Batimetri verisi indirmek için Bathy.txt
gmtset MEASURE_UNIT inch set dosya = Bathy.txt ########get minmax coordinates set R = `minmax $dosya -C | awk '{print int($1/10)*10-20,int($2/10)*10+20,int($3/10)*10-20,int($4/10)*10+20}'` #Math: aliased to set !:1 = `echo "!:3-$" | bc -l ` #Math x = $R[2]- $R[1] #Math y = $R[4]- $R[3] #echo $x $y ########### butun noktalari gridlemek icin # gridding with 0.5 metre pixel size xyz2grd -R$R[1]/$R[2]/$R[3]/$R[4] -I0.5 -F -V -Gtmp05.grd $dosya # 1 metre with 1 metre pixel size xyz2grd -R$R[1]/$R[2]/$R[3]/$R[4] -I1 -F -V -Gtmp1.grd $dosya # gridding with 2 metre pixel size xyz2grd -R$R[1]/$R[2]/$R[3]/$R[4] -I2 -F -V -Gtmp2.grd $dosya # gridding with 5 metre pixel size xyz2grd -R$R[1]/$R[2]/$R[3]/$R[4] -I5 -F -V -Gtmp5.grd $dosya ######### plot each to see the difference #grd2cpt tmp05.grd -Crainbow -V -Z > ! palet.cpt #grdimage tmp05.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.007 -P -B200 > ! x.ps #xv x.ps #grd2cpt tmp1.grd -Crainbow -V -Z > ! palet.cpt #grdimage tmp1.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.007 -P -B200 > ! x.ps #xv x.ps #grd2cpt tmp2.grd -Crainbow -V -Z > ! palet.cpt #grdimage tmp2.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.007 -P -B200 > ! x.ps #xv x.ps #grd2cpt tmp5.grd -Crainbow -V -Z > ! palet.cpt #grdimage tmp5.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.007 -P -B200 > ! x.ps #xv x.ps #### resample and make all with the same pixel size grdsample tmp1.grd -I0.5 -R$R[1]/$R[2]/$R[3]/$R[4] -Gtmp1S.grd grdsample tmp2.grd -I0.5 -R$R[1]/$R[2]/$R[3]/$R[4] -Gtmp2S.grd grdsample tmp5.grd -I0.5 -R$R[1]/$R[2]/$R[3]/$R[4] -Gtmp5S.grd grd2cpt tmp5S.grd -Crainbow -V -Z > ! palet.cpt grdimage tmp5S.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.08 -P -B200 > ! x.ps xv x.ps ##### stack them with highest resolution on top ## fill the gaps in the data with 0.5 m pixel size ## from the with 1 m, 2m and then 5 m grdmath tmp05.grd tmp1S.grd AND = tmpA.grd grdmath tmpA.grd tmp2S.grd AND = tmpB.grd grdmath tmpB.grd tmp5S.grd AND = spiky.grd ######### plot the stacked data with 0.5 m pixel size. ### create an intensity (shade) file with sun angle of 315 (NW) grdgradient spiky.grd -GspikyI.grd -Nt -A315 grd2cpt spiky.grd -Crainbow -V -Z > ! palet.cpt grdimage spiky.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.07 -P -B200 > ! x.ps xv x.ps & ### 3D plotting see the spikes !!!!! echo $R set nam = 3dspykes.ps set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] set depth = ` grdinfo -C spiky.grd | awk '{print $6,$7}'` grdview spiky.grd -E40/20 $ranges -Jx0.08 -Cpalet.cpt -IspikyI.grd -Jz0.08 -Qi300 -V > ! $nam xv $nam
############# remove the spikes now ### make a shade and a slope file grdgradient spiky.grd -Gintensity.grd -Dc -Sslope.grd -V ## convert pixels with a slope value higer than 1.2 to NaN (Not a Value) grdclip slope.grd -Gmask.grd -Sa1.2/NaN -V ## mask spikes now grdmath spiky.grd mask.grd OR = nospike.grd ## have a look at it grdgradient nospike.grd -GtmpI.grd -Nt -A315 grd2cpt nospike.grd -Crainbow -V -Z > ! palet.cpt grdimage nospike.grd -ItmpI.grd -Cpalet.cpt -R$R[1]/$R[2]/$R[3]/$R[4] -Jx0.08 -P -B200 > ! x.ps xv x.ps & # make intensity file grdgradient nospike.grd -GnospikeI.grd -Nt -A315 ### 3D plotting set ranges = -R$R[1]/$R[2]/$R[3]/$R[4] set name = 3dnospykes.ps grdview nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -Jz0.08 -Qi300 -V -Z0 -B250 > ! $name xv $name &
## with intensity grdview nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -InospikeI.grd -Jz0.08 -Qi300 -V -Z0 -B250 > ! $name xv $name &
