Gridleme ve filtreleme

From GMT Türkiye Wiki
(Sürümler arası farklar)
Jump to: navigation, search
 
 
(Bir kullanıcı tarafından yapılan 17 ara revizyon gösterilmiyor)
1. satır: 1. satır:
<bash>
+
* [[ batimetri verisini gridleme ve spike'ları maskeleme]]
  
### dosya ismi gir
 
set dosya = bathy.dat
 
  
gmtset MEASURE_UNIT inch
+
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 = 3d.ps
+
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 > ! $dosya:r3d.ps
+
 
xv $dosya:r3d.ps &
+
[[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
 

3dspykes.jpg

 
 
############# 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 &
 


3dnospykes.jpg


 
## with intensity
grdview   nospike.grd -E10/20 $ranges -Jx0.08 -Cpalet.cpt -InospikeI.grd -Jz0.08 -Qi300 -V -Z0 -B250 > ! $name
xv $name &
 
Personal tools