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

这里写图片描述

  • 3
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 6
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

IvanLJF

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值