grdproject¶
- 官方文档
- 简介
对网格数据做地图变换和逆变换
该命令可以将地理坐标下的网格数据投影到一个矩形网格中,也可以将一个矩形坐标系下的网格数据反投影到地理坐标下。
必选选项¶
<in_grdfile>
要进行变换的2D网格数据
-G<out_grdfile>
输出的网格文件名
可选选项¶
-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
决定。
示例¶
将地理网格数据变换成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