gmt绘制形变图-反缠绕

#! /bin/csh -f
rm gmt.history
#
gmt gmtset MAP_FRAME_TYPE=plain
#设置低图边框类型
gmt gmtset FONT_ANNOT_PRIMARY 9p
#FONT_ANNOT_PRIMARY设置字体大小为9p
#
gdal_translate -of GMT -b 2 filt_topophase.unw.geo filt_topophases.unw.geo.grd
gmt grdmath filt_topophases.unw.geo.grd -0.3 MUL = los.grd

Rg=80.5/82/41.5/43
#设置经纬度范围
J=m8c
#设置投影类型
outps=wrap.ps
#设置输出文件名

#gmt grdmath filt_topophase.unw.geo.grd -0.44 MUL = los.grd
#gmt makecpt -Crainbow -T-10/15/0.01 -Z >los.cpt
#makecpt生成cpt文件（cpt文件为调色板将灰度值转换为RGB），-C选择cpt主颜色表，rainbow为彩虹色（默认），-T为范围加min/max/间隔（cpt）采样间隔 -Z创建一个连续的 CPT [默认是不连续的，即， 每个间隔的恒定颜色] https://docs.generic-mapping-tools.org/5.4/makecpt.html
#gmt psbasemap -R$Rg -J$J -Ba0.5f0.25 -P -K > $outps
#绘制 PostScript 底图，-R选择感兴趣的区域，-J选择地图投影，前导 a 用于指定注释和主要刻度 间距 [默认]， f 表示小刻度间距， g 表示 网格线间距。-P选择“纵向”绘图方向。
#gmt grdimage los.grd -J$J -R$Rg -Clos.cpt  -K -O >> $outps
#grdimage投影网格或图像并将它们绘制在地图上。
#gmt pstext -J$J -R$Rg -F+f11p,Helvetica -K -O >>$outps << EOF
#82.23 44.8 (a)
#EOF
#pstext在地图上绘制或排版文本 -F调整字体，+f 设置 字体（大小、字体名称、颜色）
#gmt psxy -R$Rg -J$J -Sa0.5c -Gred  -W0.2p -K -O >>$outps << EOF
#100.008 25.761
#EOF
#psxy - 在地图上绘制线条、多边形和符号 -S绘制符号，-Sa是外接圆的直径。

gmt makecpt -Crainbow -T0/1/0.1 -Z >wrap.cpt
gmt grdmath los.grd 1 FMOD 1 ADD 1 FMOD = wrap.grd
gmt psbasemap -R$Rg -J$J -Ba0.25 -BWS -K >> $outps
gmt grdimage wrap.grd -J$J -R$Rg -Cwrap.cpt -Q  -B0 -Ben -K -O >> $outps
#gmt pstext -J$J -R$Rg -F+f11p,Helvetica -K -O >>$outps << EOF
#77.65 35.64 (b)
#EOF

#gmt psxy -R$Rg -J$J -Sa0.5c -Gred  -W0.2p -K -O >>$outps << EOF
#100.008 25.761
#EOF

gmt psscale -D8c/2.4c/5c/0.15c -Cwrap.cpt -Ba1f0.1 -O -X4.5c >> $outps
#psscale在地图上绘制灰色或彩色比例尺，-C cpt是要使用的 CPT。-B为颜色栏设置注释、刻度和网格线间隔。
gmt psconvert $outps -A -P -Tj
#psconvert使用 GhostScript 将 [E]PS 文件转换为其他格式 （BMP、EPS、JPEG、PDF、PNG、PPM、SVG、TIFF）-A 将 BoundingBox 和 HiResBoundingBox 调整到所需的最小值 通过图像内容。Tj 表示 JPEG