GMT5SAR--由*.grd文件生成*.ps

1.雷达坐标系下的DEM(topo_ra.grd),无地理坐标

#!/bin/csh
#生产颜色表
grd2cpt topo_ra.grd -Cgray -V -Z > topo_ra.cpt 
#
grdimage topo_ra.grd -JX8i -X0.2i -P -Ctopo_ra.cpt -V > topo_ra.ps

2.地理投影坐标系下的相干图(corr_ll.grd)

#!/bin/csh
#grdinfo los_ll.grd 查看x,y,z的取值范围
#颜色表去/usr/local/share/cpt文件夹下寻找
gmt makecpt -Cgray -T0/1/0.1 -Z > corr_ll.cpt

#可用-R指定画图范围
#grdimage phasefilt_mask_ll.grd -R103.2/104.3/32.5/33.3 -JM6i -P -Ba -Ctopo.cpt -V -K > GMT_tut_15.ps

#-J指定投影方式
grdimage corr_ll.grd -JM6i -P -Ba -Ccorr_ll.cpt -V -K > corr_ll.ps

#+w<length>/<width>颜色条的长度与宽度    +h画水平颜色条,默认垂直颜色条    +o<dx>/<dy>控制颜色条的水平垂直位置

psscale -DjTC+w2.5i/0.15i+o-6.5/7i -Rcorr_ll.grd -J -Ccorr_ll.cpt -I0.4 -By+lm -O >> corr_ll.ps

3.视线上形变速率图,带投影坐标(los_ll.grd)

#!/bin/csh

#grdinfo los_ll.grd 查看x,y,z的取值范围
#颜色表去/usr/local/share/cpt文件夹下寻找
#-I翻转颜色
gmt makecpt -Cjet -I -T-58/38/15 > los_ll.cpt

#可用-R指定画图范围
#grdimage phasefilt_mask_ll.grd -R103.2/104.3/32.5/33.3 -JM6i -P -Ba -Ctopo.cpt -V -K > GMT_tut_15.ps

#-J指定投影方式
grdimage los_ll.grd -JM6i -P -Ba -Clos_ll.cpt -V -K > los_ll.ps

#+w<length>/<width>颜色条的长度与宽度    +h画水平颜色条,默认垂直颜色条    +o<dx>/<dy>控制颜色条的水平垂直位置
#psscale -DjTC+w5i/0.25i+h+o0/-1i -Rphasefilt_mask_ll.grd -J -Ctopo.cpt -I0.4 -By+lm -O >> GMT_tut_15.ps

psscale -DjTC+w2.5i/0.15i+o5.5/6.3i -Rlos_ll.grd -J -Clos_ll.cpt -I0.4 -By+lm -O >> los_ll.ps

4.缠绕相位滤波掩膜图(phasefilt_mask_ll.grd)

#!/bin/csh
#颜色表去/usr/local/share/cpt文件夹下寻找
#gmt makecpt -Crainbow -T-3.14/3.14/0.8 > phasefilt.cpt
#选择使用jet颜色条
gmt makecpt -Cjet -T-3.14/3.14/0.8 > phasefilt.cpt

#可用-R指定画图范围
#grdimage phasefilt_mask_ll.grd -R103.2/104.3/32.5/33.3 -JM6i -P -Ba -Ctopo.cpt -V -K > GMT_tut_15.ps

grdimage phasefilt_mask_ll.grd -JM6i -P -Ba -Cphasefilt.cpt -V -K > phasefilt_mask_ll.ps

psscale -DjTC+w3i/0.2i+o-7/7i -Rphasefilt_mask_ll.grd -J -Cphasefilt.cpt -I0.4 -By+lm -O >> phasefilt_mask_ll.ps