1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
#!/bin/bash
# $Id: GMT_registration.sh 11490 2013-05-16 06:26:21Z pwessel $
#
# Gridline registration
gmt psxy -R0/3/0/3 -JX2.5i/1.25i -B1g1 -Bwesn -P -K -Wthinner -L -Glightred << EOF > GMT_registration.ps
0.5 1.5
1.5 1.5
1.5 2.5
0.5 2.5
EOF
gmt grdmath -R0/3/0/3 -I1 0 = tt.nc
gmt grd2xyz tt.nc | gmt psxy -R0/3/0/3 -JX2.5i/1.25i -Sc0.12i -N -G0 -O -K >> GMT_registration.ps
# Pixel registration
gmt psxy -R0/3/0/3 -JX2.5i/1.25i -B1g1 -Bwesn -O -K -W0p -L -Glightred -X2.75i << EOF >> GMT_registration.ps
1 1
2 1
2 2
1 2
EOF
gmt grdmath -R0/3/0/3 -I1 -r 0 = tt.nc
gmt grd2xyz tt.nc | gmt psxy -R0/3/0/3 -JX2.5i/1.25i -Sc0.12i -Gblack -O >> GMT_registration.ps
|
