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="49.4/57/43.4/48.4"
gmt begin youtian_1 png E1000
    
    gmt basemap -JM8c -R%R% -Bxa2df60m -Bya60mf60m -BWSne
    gmt grdcut @earth_relief_15s.grd -R%R% -Gyt.grd
    gmt grdgradient yt.grd -A40 -Ne0.6 -Gyt.grad
    gmt makecpt -Cgray -T-1600/600
    gmt grdimage yt.grd -C -I+d

    REM 海岸线
    gmt coast -JM -R%R% -W0.2p -Clightblue
    
    REM 绘制自定义指北针
    echo 49.8 47.9 -90 | gmt plot -Skcompass/0.5 -W0.6p 
    echo 49.8 48.2 N | gmt text -F+f8p,4

    REM 绘制图中图
    gmt inset begin -DjBL+w0.72i+o0.04i/0.04i -F+gwhite+p0.2p+c0.1c
            rem gmt coast -JG53/46S/? -Rg -Bg -Wfaint -Gbrown -EKAZ+gbisque -A5000
            gmt coast -R46/87/40/56 -JM? -EKZ+g+p0.2p -A5000
            echo 49.4 48.4 57 43.75 | gmt plot -Sr+s -W0.2p,blue
    gmt inset end
        
        

    REM 绘制SAR影像范围
    gmt kml2gmt asc57.kml -V > asc57.txt
    gmt plot asc57.txt -W0.8p,black,-
    gmt kml2gmt asc159.kml -V > asc159.txt
    gmt plot asc159.txt -W0.8p,blue,-
    gmt kml2gmt des.kml -V > des.txt
    gmt plot des.txt -W0.8p,red,-


    REM 绘制油田位置
    echo 53.335211 46.073415 | gmt plot -Sa0.3c -W0.3p -Gyellow
    rem gmt kml2gmt youtian.kml -V > youtian.txt
    rem gmt plot youtian.txt -W0.4p,black -Gyellow

    REM 绘制LOS向箭头
    rem echo 99.5 30.2 103 1 | gmt plot -Sv0.25c+e -W0.6p -Gblack
    rem echo 99.5 30.2 13 1 | gmt plot -Sv0.25c+e -W0.6p -Gblack
    rem echo 99.5 30.2 LOS | gmt text -F+f7p+jBL+a13 -Dj0.2c/-0.4c
    rem echo 99.5 30.2 Azimuth | gmt text -F+f7p+jBL+a103 -Dj0c/0.2c

    REM 绘制比例尺
    
    gmt coast -A100000 -I0 -Lg53/43.8+c43.8+w100k+u --FONT_TITLE=7p,4

    REM 绘制colorbar
    rem gmt colorbar -DjBR+w1.6c/0.2c+o0.7c/0.2c -Cdem3 -Bx2000 -By+l(m) --FONT_ANNOT_PRIMARY=7p
    
    echo N 1 > legend.txt
    echo S 0.1i - 0.15i - 0.5p,black,- 0.3i Asc_path57_frame148 >> legend.txt
    echo G 0.1 >> legend.txt
    echo S 0.1i - 0.15i - 0.5p,blue,- 0.3i  Asc_path159_frame149 >> legend.txt
    echo G 0.1 >> legend.txt
    echo S 0.1i - 0.15i - 0.5p,red,- 0.3i  Des_path35_frame441 >> legend.txt

    gmt legend legend.txt -DjBR+w3c+o0c/0c -F+p0p+gwhite@30  --FONT_ANNOT_PRIMARY=6p,4
    
    rem gmt colorbar -DjBR+o0.6c/0.04c+w0.7c/0.15c -C -Bx2000 -By+l"m" --FONT_ANNOT_PRIMARY=6p,4 --MAP_FRAME_PEN=0.5p
    

gmt end

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

gmt begin yt1 png E1000

    set R="53.14/53.64/45.9/46.23"
    gmt basemap -JM8c -R%R% -Bxa10mf5m -Bya5mf5m -BWSne
    gmt grdimage guang.tif
    REM 绘制自定义指北针
    echo 53.18 46.19 -90 | gmt plot -Skcompass/0.5 -W0.6p 
    echo 53.18 46.21 N | gmt text -F+f8p,4

    REM 绘制比例尺
    gmt coast -A100000 -I0 -Lg53.2/45.93+c45.93+w5k+u --FONT_TITLE=7p,4

    REM 绘制研究区位置
    gmt kml2gmt fw.kml -V > fw.txt
    gmt plot fw.txt -W0.6p,black 
    
gmt end
del gmt.conf