关于GMT绘图的几个投影

http://blog.163.com/jey_df/blog/static/182550161201381072156631/

NAME

pscoast ： To plot land-masses, water-masses, coastlines, borders, and rivers

SYNOPSIS

pscoast -Jparameters -Rwest/east/south/north[r] [ -Amin_area[/min_level/max_level][+r|l][ppercent] ] [ -B[p|s]parameters ] [-C[l|r/]fill ] [ -Dresolution[+] ] [ -Eazim/elev[+wlon/lat[/z]][+vx0/y0] ] [ -Gfill|c ] [ -Iriver[/pen] ] [ -Jz|Zparameters ] [ -K ] [ -L[f][x]lon0/lat0[/slon]/slat/length[m|n|k][+llabel][+jjust][+ppen][+ffill][+u] ] ] [ -O ] [ -Nborder[/pen] ] [ -O ] [ -P ] [ -Q ] [ -Sfill|c ] [ -T[f|m][x]lon0/lat0/size[/info][:w,e,s,n:][+gint[/mint]] ] [ -U[just/dx/dy/][c|label] ] [ -V ] [-W[level/]pen ] [ -X[a|c|r][x-shift[u]] ] [ -Y[a|c|r][y-shift[u]] ] [ -Zzlevel ] [ -ccopies ] [ -bo[s|S|d|D[ncol]|c[var1/...]] ] [ -m[flag] ]

1、绘制全球地图。

    pscoast -Rg -JN0/5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Robinson: >robinson.ps

 2、全球椭圆投影。

    pscoast -Rg -JH180/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Hammer: >Hammer.ps

 3、矩形投影（极区误差较大，一般选择±70°纬度范围）

    pscoast -R-180/180/-85/85 -JM4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Mercator: >Mercator.ps

     可见，极区误差较大，赤道附近变形较小。

4、全球柱面投影

pscoast -R-180/180/-90/90 -JQ80/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -B60/30:.Cylin-Equidistant: >cylin_equidistant.ps

5、立体投影 

pscoast -R-180/180/-90/90 -JG85/40/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Orthographic: >Orthographic.ps

 6、南极的投影。

    gmtset PLOT_DEGREE_FORMAT = ddd:mmF

    pscoast -R-180/180/-90/-60 -JA0/-90/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -B30/15:.LambertAzimutual: >LambertAzimutual.ps

 7、北极投影

    gmtset PLOT_DEGREE_FORMAT = ddd:mmF

    pscoast -R-180/180/45/90 -JA0/90/4.5i -Gwhite -Slightgray -Wthinnest -P -Dc -B30/15:.North_Polar_Stereographic: >North_Polar_Stereographic.ps

 

 

           这样看起来很不靠谱，不专业，换个投影方式。

     gmtset PLOT_DEGREE_FORMAT = ddd:mmF

     pscoast -R280/350/55/85 -Jl315/70/60/75/1:30000000  -Dl -N1/1p -N2/0.5p -W0.25p  -B10WESn:.Greenland: > GreenLand.ps  

     这样就靠谱好看得多。

    参考资料：

GMT地图投影研究

Fig_Related

makecpt -T-7.0/16.5/0.01 -Z >