gmt grdimage phasefilt.grd -JX6.5i -Bxaf+lRange -Byaf+lAzimuth -BWSen -Cphase.cpt -X1.3i -Y3i -P -K > phasefilt.ps
-JX线性投影
-B坐标轴设置
gmt psscale -Rphasefilt.grd -J -DJTC+w5i/0.2i+h -Cphase.cpt -Bxa1.57+l"Phase" -By+lrad -O >> phasefilt.ps
二、gacos使用方法
1、命令gacos_operation.csh master_ztd master_ztd.rsc slave_ztd slave_ztd.rsc reference_point incidence_angle
master_ztd、master_ztd.rsc、slave_ztd、slave_ztd.rsc 在GACOS网站下载
refence_point首先需要将corr.grd文件绘图，查找一个相干点高的点位（或者多个点位），然后创建一个文档，将距离向和方位向坐标写入，然后运行
最后将得到2个unwrap_GACOS_corrected.grd and unwrap_GACOS_corrected_trend.grd文件
2、将以上2个文件进行投影（proj_ra2ll.csh)
3、
gmt grdmath unwrap_mask_ll.grd 0.0554658 MUL -79.58 MUL = los_ll.grd
gmt grdedit -D//"mm"/1///"$PWD:t LOS displacement"/"equals negative range" los_ll.grd