13.26 grdproject¶
| 官方文档: | grdproject |
|---|---|
| 简介: | 对网格数据做地图变换和逆变换 |
该命令可以将地理坐标下的网格数据投影到一个矩形网格中，也可以将一个矩形坐标系下的网格数据反投影到地理坐标下。
13.26.1 必选选项¶
<in_grdfile>- 要进行变换的2D网格数据
-G<out_grdfile>- 输出的网格文件名
13.26.2 可选选项¶
-C[<dx>/<dy>]- 默认投影后的坐标是相对于区域的左下角，该选项使得投影后的坐标相对于投影的中心。
<dx>/<dy>是要加到投影后坐标的偏移量。 -D<xinc>[<unit>][+e|n][/<yinc>[<unit>][+e|n]]- 指定新网格的网格间隔。
-E<dpi>- 设置新网格的分辨率，即每英寸的点数。
-F[c|i|p|e|f|k|M|n|u]- 强制1:1比例，即输出数据的单位是真实的投影长度，默认单位为m。也可以指定为其他单位。
-I- 逆变换，将矩形区域变换成地理区域。
-Mc|i|p- 指定投影后的测量单位，默认值由参数
PROJ_LENGTH_UNIT决定。
13.26.3 示例¶
将地理网格数据变换成Mercator网格，分辨率为300dpi:
gmt grdproject dbdb5.nc -R20/50/12/25 -Jm0.25i -E300 -r -Gdbdb5_merc.nc
将网格数据逆变换为地理网格:
gmt grdproject topo_tm.nc -R-80/-70/20/40 -Jt-75/1:500000 -I -D5m -V -Gtopo.nc
将UTM（以米为单位）下的网格数据逆变换为地理网格:
gmt grdproject topo_utm.nc -R203/205/60/65 -Ju5/1:1 -I -Mm -Gtopo.nc -V