gmt set FORMAT_GEO_MAP ddd:mm:ssF
gmt set FONT_TITLE 5p 
gmt set FONT_ANNOT_PRIMARY 7p,Times-Roman
gmt set FONT_LABEL 7p 
gmt set MAP_FRAME_WIDTH 1.3p
gmt set MAP_FRAME_TYPE fancy 
gmt set MAP_FRAME_PEN 0.4p
gmt set MAP_TITLE_OFFSET 2p
gmt set MAP_ANNOT_OFFSET_PRIMARY 2p 
gmt set MAP_ANNOT_OFFSET_SECONDARY 2p 
gmt set MAP_TICK_LENGTH_PRIMARY 1p 

set GMT_SESSION_NAME=97401
set R="99/103.5/28.5/32.7"
gmt begin xsh png E1000
    gmt basemap -JM8c -R%R% -Bxa50mf25m -Bya50mf25m -BWSne
    rem gmt grdcut @earth_relief_15s.grd -R%R% -Gxsh.grd
    rem gmt grdgradient xsh.grd -A40 -Ne0.6 -Gxsh.grad
    gmt makecpt -Cdem3 -T419/7039/10 -D
    gmt grdimage xsh.grd -C -Ixsh.grad
    REM 绘制指北针
    gmt coast -Df -I0 -Tf103.1/32.36/0.3i/1 --FONT_TITLE=7p,4
    REM 绘制断层
    gmt kml2gmt fault-line.kml -Fl -V > fault-line.txt
    gmt plot fault-line.txt -W1p,black
    REM 绘制SAR影像范围
    gmt plot alo2_asc0.txt -W1p,yellow,-
    gmt plot alo2_asc1.txt -W1p,yellow,-
    gmt plot alo2_asc2.txt -W1p,yellow,-
    gmt plot alo2_asc3.txt -W1p,yellow,-
    gmt plot alo2_asc4.txt -W1p,yellow,-
    gmt plot alo2_asc5.txt -W1p,yellow,-

    gmt plot alo2_des1.txt -W1p,magenta,-
    gmt plot alo2_des2.txt -W1p,magenta,-
    gmt plot alo2_des3.txt -W1p,magenta,-
    gmt plot alo2_des4.txt -W1p,magenta,-
    gmt plot alo2_des5.txt -W1p,magenta,-

    gmt plot envisat_des.txt -W1p,orange,-

    gmt kml2gmt sen_asc_up.kml -V > sentinel_1_asc.txt
    gmt plot sentinel_1_asc.txt -W1p,blue,-
    gmt kml2gmt sen_asc_down.kml -V > sentinel_2_asc.txt
    gmt plot sentinel_2_asc.txt -W1p,blue,-

    gmt kml2gmt sen_des_1.kml -V > sen_des_1.txt
    gmt plot sen_des_1.txt -W1p,green,-
    gmt kml2gmt sen_des_2.kml -V > sen_des_2.txt
    gmt plot sen_des_2.txt -W1p,green,-
    gmt kml2gmt sen_des_3.kml -V > sen_des_3.txt
    gmt plot sen_des_3.txt -W1p,green,-

    gmt kml2gmt alo1_asc_1.kml -V > alo1_asc_1.txt
    gmt plot alo1_asc_1.txt -W1p,darkblue,-
    gmt kml2gmt alo1_asc_2.kml -V > alo1_asc_2.txt
    gmt plot alo1_asc_2.txt -W1p,darkblue,-
    gmt kml2gmt alo1_asc_3.kml -V > alo1_asc_3.txt
    gmt plot alo1_asc_3.txt -W1p,darkblue,-
    gmt kml2gmt alo1_asc_4.kml -V > alo1_asc_4.txt
    gmt plot alo1_asc_4.txt -W1p,darkblue,-
    gmt kml2gmt alo1_asc_5.kml -V > alo1_asc_5.txt
    gmt plot alo1_asc_5.txt -W1p,darkblue,-
    gmt kml2gmt alo1_asc_6.kml -V > alo1_asc_6.txt
    gmt plot alo1_asc_6.txt -W1p,darkblue,-

    gmt kml2gmt alo1_des_1.kml -V > alo1_des_1.txt
    gmt plot alo1_des_1.txt -W1p,red,-


    REM 绘制图例
    echo N 4 > legend.txt
    echo S 0.12i - 0.15i - 0.5p,blue,- 0.2i Sentinel A >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,darkblue,- 0.2i ALOS-1 A >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,yellow,- 0.2i ALOS-2 A >> legend.txt
    rem echo S 0.12i - 0.15i - 0.5p,red,- 0.2i Envisat A >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,black 0.2i  Fault >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,green,- 0.2i Sentinel D >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,red,- 0.2i ALOS-1 D >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,magenta,- 0.2i ALOS-2 D >> legend.txt
    echo S 0.12i - 0.15i - 0.5p,orange,- 0.2i Envisat D >> legend.txt
    gmt legend legend.txt -DjBL+w8c+o0c/0c -F+p0p+gwhite@30   --FONT_ANNOT_PRIMARY=7p,Times-Roman


gmt end show