0%

GMT极地投影

GMT的地图投影分成了4类,分别是圆锥投影、方位投影、圆柱投影和混合投影,共计有30多种投影方法。地球的极地区域多使用极地方位投影。

极地方位投影 Polar Stereographic Map

方位投影特点是保持角度无变形,面积和距离有变形,又称为等角投影。方位投影多用于单个半球以及极地区域地图,当用于极地时成为Polar Stereographic Map,一般区域则成为Stereographic Map

距离投影中心点越远,变形越大。

更多信息:

1
https://manifold.net/doc/mfd9/polar_stereographic_projection.htm

下面使用一行代码绘制北极区域地图。如果需要叠加其他点、线等要素,则使用对应的命令追加即可。

投影中心在北极点,展示极地区域

这种地图特点是边界为圆形,以北极点为中心,各纬圈方向变形均匀对称。

1
gmt pscoast -R0/360/60/90 -JS0/90/10c -Bag -Dl -A25000 -Groyalblue -Sseashell > GMT_stereographic_polar1.ps

投影中心点偏移,不在北极点

通过设置-JSlon/lat/size中的lon/lat调整投影中心点。地图仍旧是圆形边界,但是中心点已经不是北极点,投影的变形也不具备对称特点。

1
gmt pscoast -R120/260/10/75 -JS190/70/10c -Bag -Dl -A25000 -Groyalblue -Sseashell > GMT_stereographic_polar2.ps

投影中心点偏移,非全极效果

设置-JSlon/lat/size中的lat为90。但是-R的坐标范围是扇形区域,绘制出的地图也是扇形。

1
gmt pscoast -R100/260/10/75 -JS180/90/10c -Bag -Dl -A25000 -Groyalblue -Sseashell > GMT_stereographic_polar3.ps

矩形地图展示全极区

通过设置-R-45/50/135/50r,以左下角和右上角坐标为边界,并且加+r。这种地图也较为常用,但是GMT绘制技巧埋藏的比较深,同样是JS投影,改变的仅仅是坐标边界。

1
gmt pscoast -R-45/50/135/50r -JS0/90/11c -Ba30g30 -Dl -A25000 -Wthinnest -Groyalblue -Sseashell > GMT_stereographic_rect.ps

竖版地图

竖版地图对我们的海洋观念影响很大,在国内也一度热炒。GMT中有一种通用墨卡托全球投影(Transverse Mercator projection),可以将南北极同时绘制。

1
gmt pscoast -R0/360/-80/80 -JT180/-30/4i -A25000 -Groyalblue -Sseashell  -Bg30a30 >TMcontours.ps