GMT绘制研究区与光学影像
2022-6-6 10:22
阅读:2961
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
转载本文请联系原作者获取授权,同时请注明本文来自吴远昆科学网博客。
链接地址:https://wap.sciencenet.cn/blog-3473451-1341767.html?mobile=1
收藏
当前推荐数:0
推荐到博客首页
网友评论0 条评论