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.pspsscale -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.psgrdimage phasefilt_mask_ll.grd -JM6i -P -Ba -Cphasefilt.cpt -V -K > phasefilt_mask_ll.pspsscale -DjTC+w3i/0.2i+o-7/7i -Rphasefilt_mask_ll.grd -J -Cphasefilt.cpt -I0.4 -By+lm -O >> phasefilt_mask_ll.ps