带有公制网格叠加的墨卡托投影（从中心以米为增量） - GMT Generic Mapping Tools
Posted
技术标签:
【中文标题】带有公制网格叠加的墨卡托投影（从中心以米为增量） - GMT Generic Mapping Tools【英文标题】：Mercator projection with metric grid overlay (metre increments from centre) - GMT Generic Mapping Tools 【发布时间】：2016-05-16 13:34:36 【问题描述】：我有一张蒙特塞拉特图，以 -62.2°、16.75° 为中心，使用通用映射工具 (GMT) 创建了墨卡托投影，如下所示。
我想在此之上添加一个网格，以中心为原点。但是，我想以米为单位而不是度数来生成覆盖网格。 有没有一种方法可以将 -62.2°、16.75° 指定为原点，然后从原点（1000m、2000m、5000m 等）生成 1000m 倍数的网格线（可能通过 psxy？）？
用于生成此图的代码是
set loc_area = -R-63.2/-61.2/15.75/17.75
set loc_proj = -Jm3.25i
set outfile = bath_map.ps
psbasemap $loc_area $loc_proj -O -Bf1a1WSen:.'Montserrat Local': -X7.5i -Y0 -K >> $outfile
grdimage $loc_area $loc_proj -O -Clant.cpt lant_topo.nc -K >> $outfile
grdcontour lant_topo.nc -Ccont.cpt -A- -Q100 -W $loc_area $loc_proj -O -K >> $outfile
【问题讨论】：
【参考方案1】：只需将度数转换为米即可；你的网格会超级小
gmt pscoast -R-62.25/-62.12/16.65/16.84 -JM6i -P -B.01666666g.00898311174 -Ggray -Dh -Slightblue -Wred > GMT_tut_4.ps
【讨论】：
以上是关于带有公制网格叠加的墨卡托投影（从中心以米为增量） - GMT Generic Mapping Tools的主要内容，如果未能解决你的问题，请参考以下文章