2024年6月13日发(作者：乾致)
psbasemap GMT画底图的指令
一定要有的参数
-B 设定图边界尺标的间隔
a 文字批注的间隔
f 边框标记的间隔
g 隔线的间隔
如果要针对X轴或Y轴做不同的设定，就在参数间加斜线/
NEWS 每个轴加上标记
news 每个轴不加标记
:.title: 在图上面加title这个标题
:title: 在轴上加title这个标题（只有在-Jx的情况下有用）
例如-BNEWsf2a3g2:xx:/f3a3g2:yy::.title:
-Jx(-JX) 线性投影(Linear Projection)一般绘图用
-Jx后面接的数字代表每单位的长度（预设长度单位为英吋，加c的话就
变公分）
-JX后面接的数字代表轴的全长
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-Jm(-JM) 麦卡托投影(Mercator Projection)地图绘图用
-Jm后面接的数字代表每度的长度（预设长度单位为英吋，加c的话就变
公分）
-JM后面接的数字代表整张图的宽
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-R 地图的经纬度范围
西/东/南/北 依序填入经纬度
> 把画出来的图输出到这个档案
不一定要有的参数
-E 旋转地图的视野方向
-G 把图着色，用r/g/b调
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psbasemap -JM7 -R117/126/18/26 -BNEWsf2a3g2:xx:/f3a3g2:y::.test: -P
-V -K >
pscoast GMT画陆块海洋和海岸线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-E 旋转地图的视野方向
-G 把干地着色，用r/g/b调色
-S 把湿地着色，用r/g/b调色
-D 图的分辨率(f,h,i,l,c)左往右为分辨率高至低
-I 加河流上去(1-10)大河至小河(a,r,i,c)大小河都画
-W 画海岸线
-N 画边界(1.国界,2.州界,3.海洋,a.全部)
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
pscoast -JM7 -R117/126/18/26 -G50/50/50 -S150/150/150 -W1 -Di -V -O
-K >>
psxy GMT在地图上画点线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-A 画大圆线
-B 和先前psbasemap的一样
-C 对照颜色，后面要接文件名称
-G 把符号着色，用r/g/b调色
-K 在这一行下面还有后续的指令
-N 不要跳过超出范围的点
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-M 线分段
-Sa 在图上画星星，后面要接大小
-Sc 在图上画圆圈，后面要接大小
-Sd 在图上画钻石，后面要接大小
-Se 在图上画椭圆形，在资料的第3,4,5栏要有方向长轴短轴
-Sh 在图上画六角形，后面要接大小
-Ss 在图上画正方型，后面要接大小
-St 在图上画三角形，后面要接大小
-Si 在图上画反三角，后面要接大小
-Sv 在图上画向量，在数据的第3,4栏要有方向长度
-Sx 在图上画叉叉，后面要接大小
-V 画面上会显示本行指令的执行结果
-W 笔宽
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式一 把数据直接打在档案里
psxy -JM7 -R117/126/18/26 -W1 -O -K -V‹‹! >>
120 23
124 27
!
表示方式二 把数据另外打在数据文件里
psxy -JM7 -R117/126/18/26 -W1 -M -V -O -K >>
在的档案里是这样（这是画线条用的档案）
>
120 25
123 25
>
121 23
123 27
在的档案里是这样（这是点用的档案）
120 25
123 25
121 23
123 27
表示方式三 把更多的数据打在数据文件里
psxy -R117/126/18/26 -JM7 -Sc -Ccpt >>
在的档案里要有四栏数据，依序为X,Y,颜色代码,圆的大小
grdimage GMT画二维等高图的指令
一定要有的参数
grd grd档案的档名
-C 对照颜色，后面要接文件名称
-J 投影方式，和前面一样
不一定要有的参数
-B 和先前psbasemap的一样
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdimage -JOc122/21/22/69/4.8i -R117/126/18/26 -
-P -K -V >>
psscale GMT画比例尺的指令
一定要有的参数
-D 定义比例尺的位置和大小
如果比例尺是水平的 -D距中央/距顶端/长度/宽度h
如果比例尺是垂直的 -D距左端/距中央/长度/宽度
不一定要有的参数
-B 定义比例尺的间隔，和前面一样
-C 对照颜色，后面要接文件名称
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psscale - -D2.8/-1/4/0.22h -Bf1000a2000:”Elevation (m)”:
-O -K >>
grdcontour GMT画二维等高线图的指令
一定要有的参数
grd grd档案的档名
-C 等高线的间距
-J 投影方式，和前面一样
不一定要有的参数
-A 等高线数值标示的间隔
-B 和先前psbasemap的一样
-E 视角
-G 同一条等高线数值标示的间隔
-K 在这一行下面还有后续的指令
-L 等高线最低或最高的限制
-N 等高线的单位（如果有下-A的参数才有效）
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-T 高地（-Ｔ+）或低地（-Ｔ-）加上标记
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdcontour -JM7 -R117/126/18/26 -C0.025 -W1 -A1 -P -K -V >>
回到目录：回到最前面
grd数据的处理
转换成xyz值
grd2xyz -R117/126/18/26 >
块状平均（I是精确度每多少度一格）
blockmean -R117/126/18/26 -I0.1 >
转换成grd檔
surface -R117/126/18/26 -I0.1 -
查询档案信息
minmax
grdinfo
2024年6月13日发(作者：乾致)
psbasemap GMT画底图的指令
一定要有的参数
-B 设定图边界尺标的间隔
a 文字批注的间隔
f 边框标记的间隔
g 隔线的间隔
如果要针对X轴或Y轴做不同的设定，就在参数间加斜线/
NEWS 每个轴加上标记
news 每个轴不加标记
:.title: 在图上面加title这个标题
:title: 在轴上加title这个标题（只有在-Jx的情况下有用）
例如-BNEWsf2a3g2:xx:/f3a3g2:yy::.title:
-Jx(-JX) 线性投影(Linear Projection)一般绘图用
-Jx后面接的数字代表每单位的长度（预设长度单位为英吋，加c的话就
变公分）
-JX后面接的数字代表轴的全长
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-Jm(-JM) 麦卡托投影(Mercator Projection)地图绘图用
-Jm后面接的数字代表每度的长度（预设长度单位为英吋，加c的话就变
公分）
-JM后面接的数字代表整张图的宽
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-R 地图的经纬度范围
西/东/南/北 依序填入经纬度
> 把画出来的图输出到这个档案
不一定要有的参数
-E 旋转地图的视野方向
-G 把图着色，用r/g/b调
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psbasemap -JM7 -R117/126/18/26 -BNEWsf2a3g2:xx:/f3a3g2:y::.test: -P
-V -K >
pscoast GMT画陆块海洋和海岸线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-E 旋转地图的视野方向
-G 把干地着色，用r/g/b调色
-S 把湿地着色，用r/g/b调色
-D 图的分辨率(f,h,i,l,c)左往右为分辨率高至低
-I 加河流上去(1-10)大河至小河(a,r,i,c)大小河都画
-W 画海岸线
-N 画边界(1.国界,2.州界,3.海洋,a.全部)
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
pscoast -JM7 -R117/126/18/26 -G50/50/50 -S150/150/150 -W1 -Di -V -O
-K >>
psxy GMT在地图上画点线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-A 画大圆线
-B 和先前psbasemap的一样
-C 对照颜色，后面要接文件名称
-G 把符号着色，用r/g/b调色
-K 在这一行下面还有后续的指令
-N 不要跳过超出范围的点
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-M 线分段
-Sa 在图上画星星，后面要接大小
-Sc 在图上画圆圈，后面要接大小
-Sd 在图上画钻石，后面要接大小
-Se 在图上画椭圆形，在资料的第3,4,5栏要有方向长轴短轴
-Sh 在图上画六角形，后面要接大小
-Ss 在图上画正方型，后面要接大小
-St 在图上画三角形，后面要接大小
-Si 在图上画反三角，后面要接大小
-Sv 在图上画向量，在数据的第3,4栏要有方向长度
-Sx 在图上画叉叉，后面要接大小
-V 画面上会显示本行指令的执行结果
-W 笔宽
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式一 把数据直接打在档案里
psxy -JM7 -R117/126/18/26 -W1 -O -K -V‹‹! >>
120 23
124 27
!
表示方式二 把数据另外打在数据文件里
psxy -JM7 -R117/126/18/26 -W1 -M -V -O -K >>
在的档案里是这样（这是画线条用的档案）
>
120 25
123 25
>
121 23
123 27
在的档案里是这样（这是点用的档案）
120 25
123 25
121 23
123 27
表示方式三 把更多的数据打在数据文件里
psxy -R117/126/18/26 -JM7 -Sc -Ccpt >>
在的档案里要有四栏数据，依序为X,Y,颜色代码,圆的大小
grdimage GMT画二维等高图的指令
一定要有的参数
grd grd档案的档名
-C 对照颜色，后面要接文件名称
-J 投影方式，和前面一样
不一定要有的参数
-B 和先前psbasemap的一样
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdimage -JOc122/21/22/69/4.8i -R117/126/18/26 -
-P -K -V >>
psscale GMT画比例尺的指令
一定要有的参数
-D 定义比例尺的位置和大小
如果比例尺是水平的 -D距中央/距顶端/长度/宽度h
如果比例尺是垂直的 -D距左端/距中央/长度/宽度
不一定要有的参数
-B 定义比例尺的间隔，和前面一样
-C 对照颜色，后面要接文件名称
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psscale - -D2.8/-1/4/0.22h -Bf1000a2000:”Elevation (m)”:
-O -K >>
grdcontour GMT画二维等高线图的指令
一定要有的参数
grd grd档案的档名
-C 等高线的间距
-J 投影方式，和前面一样
不一定要有的参数
-A 等高线数值标示的间隔
-B 和先前psbasemap的一样
-E 视角
-G 同一条等高线数值标示的间隔
-K 在这一行下面还有后续的指令
-L 等高线最低或最高的限制
-N 等高线的单位（如果有下-A的参数才有效）
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-T 高地（-Ｔ+）或低地（-Ｔ-）加上标记
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdcontour -JM7 -R117/126/18/26 -C0.025 -W1 -A1 -P -K -V >>
回到目录：回到最前面
grd数据的处理
转换成xyz值
grd2xyz -R117/126/18/26 >
块状平均（I是精确度每多少度一格）
blockmean -R117/126/18/26 -I0.1 >
转换成grd檔
surface -R117/126/18/26 -I0.1 -
查询档案信息
minmax
grdinfo