GMTver6で中国行政界を描く

バージョンの問題ではありません。Ver.5でも描画可能です。強いて言えばDCW(The Digital Chart of the World Data)の問題だったということです。東経105度北緯35度を投影中心とし、北緯20度と50度線を標準緯線に設定したランベルト正角図法で描画してみました。

f:id:S_E_Hyphen:20210821171804p:plain

 

#!/bin/bash
county="\
CN.34,CN.35,CN.62,CN.44,CN.52,CN.46,CN.13,CN.23,CN.41,CN.42,\
CN.43,CN.32,CN.36,CN.22,CN.21,CN.63,CN.61,CN.37,CN.14,CN.51,\
CN.53,CN.33,CN.45,CN.15,CN.64,CN.65,CN.54,CN.11,CN.50,CN.31,\
CN.12"

gmt begin MapofChina png
 gmt gmtset FONT_TITLE  12p,Song-Medium--UniGB-UTF8-H,black
 gmt gmtset FONT_LOGO  5p,STSong-Light--UniGB-UTF8-H
 gmt gmtset FONT_ANNOT_PRIMARY  3p,GothicBBB-Medium-UniJIS-UTF8-H,black
 gmt gmtset MAP_TITLE_OFFSET 2p
 gmt gmtset MAP_LABEL_OFFSET 1p
 gmt gmtset MAP_TICK_LENGTH 3p
 gmt gmtset MAP_ANNOT_OFFSET_PRIMARY 1p
 gmt gmtset MAP_LOGO true
 gmt gmtset MAP_LOGO_POS TR/+8.5c/-6.5p
 gmt gmtset FORMAT_TIME_STAMP %Y年%m月%d日

 gmt basemap \
    -X0.5c -Y0.5c \
    -JL105/35/20/50/8c \
    -R75/10/145/50r -BNEWS+t"中国省份地图"
 gmt coast -S0/0/0/20 -E${county}+p#000000
 gmt basemap -Ba10f10g10 -U
gmt end show