xiaoairen322的個人博客分享 http://www.ueservicedoffices.com/u/xiaoairen322

博文

GMT繪制中國區域的地形圖

已有 950 次閱讀 2019-5-7 21:41 |個人分類:GMT|系統分類:科研筆記

學習GMT也有幾天了 ,GMT作為地學領域的一門優秀繪圖語言,一直深深的吸引著我,聯系了幾天也繪制了一副中國區域的地形圖,分辨率為5m。在此感謝李建濤(濤)的文章http://www.ueservicedoffices.com/home.php?mod=space&uid=3391834,還有GMT中文社區https://gmt-china.org/#docs還有周老師的一些例子http://www.ueservicedoffices.com/home.php?mod=space&uid=350278,  非常感謝劉珠妹老師,九天學者郭志道博士的熱心講解,這次我是在linux平臺gmt5.4.5下繪制 代碼如下,在此不做注釋(雖說沒有注釋就是沒有靈魂的代碼)

#!/bin/bash

R=70/135/15/55

J=M110/16c

PS=china.ps

D=earth_relief_05m.grd


gmt set MAP_FRAME_TYPE=fancy


#gmt set MAP_GRID_PEN_PRIMARY 0.25p,black,2_2:1


#gmt set FROMAT_GEO_MAP ddd.xxF


gmt set MAP_FRAME_WIDTH 2.5p


gmt set FONT_ANNOT_PRIMARY 8p


gmt set MAP_FONT_OFFSET_PRIMARY 0.5c


gmt set FONT_TITLE 15P


gmt set FONT_TITLE_OFFSET_FRAME 2c


gmt set MAP_TICK_LENGTH_PRIMARY -2p


gmt grdcut $D -R65/140/12/65 -Gchina2.grd


gmt grdgradient china2.grd -A0 -Nt -Gint.grd


gmt psbasemap -R$R -J$J -B10g10 -BESwn+t"relief Map of China" -U"Made by Hongbin Li" -P -V -K > $PS


gmt grdimage -R$R -J$J china2.grd -Iint.grd  -Cetopo1.cpt -K -P -O >> $PS


gmt pscoast -R$R -J$J -B10g10   -Ia/0.15p,177/178/183 -I1/0.5p,61/99/172 -C81/174/254 -Lg102/12+c30+w1000+u+f -K -O >> $PS


gmt psxy CN-border-La.dat   -R$R -J$J -K -O >> $PS


gmt psxy CN-capitals.dat -J$J -R$R -Sc0.08c -Gblack -K -O >> $PS


gmt psxy CN-capitals.dat -J$J -R$R -Sc0.15c -W0.2p,black -K -O >> $PS


gmt pstext CN-capitals.dat -J$J -R$R -F+f7p,35+j -Dj0.15c/0.15c -K -O >> $PS

繪制的地形圖如下:

bbb.jpg



http://www.ueservicedoffices.com/blog-3387053-1177697.html

上一篇:關于ppp---online的一些認識與看法
下一篇:關于GNSS工程控制網投影變形的一些事

0

該博文允許注冊用戶評論 請點擊登錄 評論 (0 個評論)

數據加載中...
掃一掃,分享此博文

Archiver|手機版|科學網 ( 京ICP備14006957 )

GMT+8, 2019-8-19 11:45

Powered by ScienceNet.cn

Copyright © 2007- 中國科學報社

返回頂部
时时彩平台 衡东县  永吉县  光山县  湖口县  红安县  乐陵市   孝感市  塔河县  温州市  呼玛县   桦南县  房产  广东省  三门县  门源  若羌县   得荣县  泰州市  营山县  金秀