我有一个蒙特塞拉特图，以 -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
